DAI, JIAN | 24 Jun 2012 06:27
Picon
Favicon

Factors that determine the tilt angle of a peptide

Hello everyone:
We are calculating the structure of a membrane peptide with two helices linked with a loop with PISEMA data
(chemical shift, dipolar coupling) with dihedral angle, hydrogen bond restraints for the helical
part, we also have some distance restraints for atoms from both helices.
>From PISEMA experiment we expect to get models with the right tilt angle with respect to the membrane
normal (i.e., the z-axis), but we got models with a variety of tilt angles. For example, we expect to have a
tilt angle of 12 to 17 degrees, but some of the models have tilt angle that is 26, which is way too large. 
I'm wondering what are the factors that determine the tilt angle of the peptide in XPLOR calculation.
Thanks.
Jian
Charles Schwieters | 25 Jun 2012 17:56

Re: Factors that determine the tilt angle of a peptide


Hello Jian--

> We are calculating the structure of a membrane peptide with two
> helices linked with a loop with PISEMA data (chemical shift, dipolar
> coupling) with dihedral angle, hydrogen bond restraints for the
> helical part, we also have some distance restraints for atoms from
> both helices.
>
> From PISEMA experiment we expect to get models with the right tilt
> angle with respect to the membrane normal (i.e., the z-axis), but we
> got models with a variety of tilt angles. For example, we expect to
> have a tilt angle of 12 to 17 degrees, but some of the models have
> tilt angle that is 26, which is way too large. 
> I'm wondering what are the factors that determine the tilt angle of
> the peptide in XPLOR calculation.

Are your restraints well-satisfied? There may be some problem such
that there are large violations causing convergence problems. If your
PISEMA-based restraints are satisfied, then the tilt angle should be
correct. No?

best regards--
Charles

--
Charles Schwieters     email:   Charles <at> Schwieters.org
                       www:     http://schwieters.org/cds
(Continue reading)

DAI, JIAN | 26 Jun 2012 19:56
Picon
Favicon

Re: Factors that determine the tilt angle of a peptide

Charles:
You are right, there are large violations. I'm confused about the values that should be used in the CS and DC tables.
Our experimental collaborators have given us a table of the CS values (let's call it Set 1), but another
computational collaborator uses the value that is subtracted by the isotropic value, 119.3 ppm for
non-Glycine residues (let's call this Set 2. In this table, each value is the corresponding value from Set
1 minus 119.3).
It seems strange for me that when Set 1 was used, the structures look OK, but with huge violations for CS;  when
Set 2 was used, the structure look strange but violations are small. What do you think might be the problem?
Do the simulation parameters in refine.py matter much? Is there a script specifically designed for
membrane proteins?
Thanks a lot.
Jian

________________________________________
From: Charles Schwieters [charles <at> schwieters.org]
Sent: Monday, June 25, 2012 11:56 AM
To: DAI, JIAN
Cc: xplor-nih <at> nmr.cit.nih.gov
Subject: Re: [Xplor-nih] Factors that determine the tilt angle of a peptide

Hello Jian--

> We are calculating the structure of a membrane peptide with two
> helices linked with a loop with PISEMA data (chemical shift, dipolar
> coupling) with dihedral angle, hydrogen bond restraints for the
> helical part, we also have some distance restraints for atoms from
> both helices.
>
> From PISEMA experiment we expect to get models with the right tilt
> angle with respect to the membrane normal (i.e., the z-axis), but we
(Continue reading)

Charles Schwieters | 26 Jun 2012 23:26

Re: Factors that determine the tilt angle of a peptide


Hello Jian--

> You are right, there are large violations. I'm confused about the
> values that should be used in the CS and DC tables.  Our
> experimental collaborators have given us a table of the CS values
> (let's call it Set 1), but another computational collaborator uses
> the value that is subtracted by the isotropic value, 119.3 ppm for
> non-Glycine residues (let's call this Set 2. In this table, each
> value is the corresponding value from Set 1 minus 119.3).  It seems
> strange for me that when Set 1 was used, the structures look OK, but
> with huge violations for CS; when Set 2 was used, the structure look
> strange but violations are small. What do you think might be the
> problem?

You might check out

http://nmr.cit.nih.gov/xplor-nih/doc/current/python/ref/csa-methods-marvin.pdf

(referenced from)
http://nmr.cit.nih.gov/xplor-nih/doc/current/python/ref/csaPotTools.html

> Do the simulation parameters in refine.py matter much?  Is
> there a script specifically designed for membrane proteins?

please see eginput/assignFit/02refine.py in the Xplor-NIH
distribution for an example.

I hope this helps--
Charles
(Continue reading)

Jakob Toudahl Nielsen | 13 Jul 2012 10:45
Picon
Picon

psf generation for at cyclic peptide

Dear all,

I am trying to generate a psf file for a small (13 aa) cyclic peptide  
based on the template, "eginput/PSF_generation/genCircPep.py".
However, it seems that it is not working: all HT1/HT2/HT3 and OT1/OT2  
are still present in the generated psf file and from careful reading  
of the numbers of the bond definitions in the psf-file it seems that  
there is defined a bond between resi 1 C and resi 13 N (it should be  
the oposite).
I appeared to be that the references in the patch definition should be  
swapped:

before: patch PEPT reference=-=(resid 1) reference=+=(resid 13) end
after:  patch PEPT reference=+=(resid 1) reference=-=(resid 13) end

but after this change xplor complaints about missing bond parameters:

  %CODBON-ERR: missing bond parameters %%%%%%%%%%%%%%%%%%%%%%%%%%
   bond energy constant missing.
   target bond length missing.
   ATOM1: SEGId="    ",  RESId="13  ",  NAME="C   ",  CHEMical="C   "
   ATOM2: SEGId="    ",  RESId="1   ",  NAME="N   ",  CHEMical="NH3 "

I hope someone can help me figure out this problem!

best regards,

Jakob

--

-- 
(Continue reading)

Charles Schwieters | 13 Jul 2012 16:25

Re: psf generation for at cyclic peptide


Hello Jakob--

>
> I am trying to generate a psf file for a small (13 aa) cyclic peptide
> based on the template, "eginput/PSF_generation/genCircPep.py".
> However, it seems that it is not working: all HT1/HT2/HT3 and OT1/OT2
> are still present in the generated psf file and from careful reading
> of the numbers of the bond definitions in the psf-file it seems that
> there is defined a bond between resi 1 C and resi 13 N (it should be
> the oposite).

Well, that's embarrassing. Clearly, that script was never
tested. Please find attached a working version of the script

my apologies--
Charles

--
Charles Schwieters     email:   Charles <at> Schwieters.org
                       www:     http://schwieters.org/cds
phone: (301) 402-4914  PGP key: http://schwieters.org/cds/pgp.txt

#
# generate PSF and initial pdb w/ correct covalent geometry for
#
(Continue reading)


Gmane