Hi all,<div>I have a 25 amino acid peptide that contains 4 phosphoserines. The forcefield that I have been using is the ff43ap which was modified to include these groups. So during pdb2gmx I have the rtp and itp files for the forcefield in the folder in order for the software to find the appropriate charges. Anyway, simulation is all honkey dory but the problem is in DSSP analysis which separates the SEP from the protein.</div>
<div><br></div><div>So the entire protein is 252 atoms but when I run the do_dssp program it is separated into two:</div><div>Protein 208 atoms &</div><div>SEP is 44 atoms </div><div><br></div><div>I can run the dssp separately on the two groups, but since dssp uses information about adjacent sequence of amino acids, I would say separating the analysis in the SEP region 16-19 would be wrong as secondary structure would not be additive like lets say dihedral angles or rmsf. </div>
<div><br></div><div>For results, in the SEP region I am getting the random coil (white as coloured by xpm2ps) which is equivalent to no secondary structure. Although I don't have a hard time believing this, just to be a devils advocate of my own results, I would like to process them together instead of separately. So my question is how to I join atoms into one group instead of having them separated?</div>
<div><br></div><div>Cordially,</div><div>Ali</div>