Implementation and application of helix–helix distance and crossing angle restraint potentials

Based on the definition of helix–helix distance and crossing angle introduced by Chothia et al. (J Mol Biol 1981, 145, 215), we have developed the restraint potentials by which the distance and crossing angle of two selected helices can be maintained around target values during molecular dynamics simulations. A series of assessments show that calculated restraint forces are numerically accurate. Since the restraint forces are only exerted on atoms which define the helical principal axes, each helix can rotate along its helical axis, depending on the helix–helix intermolecular interactions. Such a restraint potential enables us to characterize the helix–helix interactions at atomic details by sampling their conformational space around specific distance and crossing angle with (restraint) force‐dependent fluctuations. Its efficacy is illustrated by calculating the potential of mean force as a function of helix–helix distance between two transmembrane helical peptides in an implicit membrane model. © 2006 Wiley Periodicals, Inc. J Comput Chem 28: 669–680, 2007

[1]  David R. Westhead,et al.  Calculation of Helix Packing Angles in Protein Structures , 2003, Bioinform..

[2]  B. Roux,et al.  Lipid-mediated interactions between intrinsic membrane proteins: dependence on protein size and lipid composition. , 2001, Biophysical journal.

[3]  H. Bosshard,et al.  Fast folding of the HIV-1 and SIV gp41 six-helix bundles. , 2004, Journal of molecular biology.

[4]  W. DeGrado,et al.  How do helix–helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo‐oligomeric helical bundles , 2003, Protein science : a publication of the Protein Society.

[5]  M. Bansal,et al.  HELANAL: A Program to Characterize Helix Geometry in Proteins , 2000, Journal of biomolecular structure & dynamics.

[6]  S. O. Smith,et al.  Structure of the transmembrane dimer interface of glycophorin A in membrane bilayers. , 2001, Biochemistry.

[7]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[8]  P. Adams,et al.  Experimentally based orientational refinement of membrane protein models: A structure for the Influenza A M2 H+ channel. , 1999, Journal of molecular biology.

[9]  Benny G. Johnson,et al.  A rotationally invariant procedure for density functional calculations , 1994 .

[10]  D. Engelman,et al.  Helical membrane protein folding, stability, and evolution. , 2000, Annual review of biochemistry.

[11]  James H. Prestegard,et al.  A Transmembrane Helix Dimer: Structure and Implications , 1997, Science.

[12]  D. Raleigh,et al.  De novo design of helical bundles as models for understanding protein folding and function. , 2000, Accounts of chemical research.

[13]  D. Engelman,et al.  Membrane protein folding and oligomerization: the two-stage model. , 1990, Biochemistry.

[14]  R. J. Cushley,et al.  NMR studies of lipoprotein structure. , 2002, Annual review of biophysics and biomolecular structure.

[15]  D. Engelman,et al.  Motifs of two small residues can assist but are not sufficient to mediate transmembrane helix interactions. , 2004, Journal of molecular biology.

[16]  D. Engelman,et al.  Interhelical hydrogen bonding drives strong interactions in membrane proteins , 2000, Nature Structural Biology.

[17]  Walter L Ash,et al.  Direct simulation of transmembrane helix association: role of asparagines. , 2004, Biophysical journal.

[18]  G. Heijne,et al.  Recognition of transmembrane helices by the endoplasmic reticulum translocon , 2005, Nature.

[19]  Yang Zhang,et al.  Application of sparse NMR restraints to large-scale protein structure prediction. , 2004, Biophysical journal.

[20]  Ming-Ming Zhou.  Phosphothreonine recognition comes into focus , 2000, Nature Structural Biology.

[21]  M. Gerstein,et al.  Statistical analysis of amino acid patterns in transmembrane helices: the GxxxG motif occurs frequently and in association with beta-branched residues at neighboring positions. , 2000, Journal of molecular biology.

[22]  J. Thornton,et al.  PROMOTIF—A program to identify and analyze structural motifs in proteins , 1996, Protein science : a publication of the Protein Society.

[23]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[24]  C. Brooks,et al.  An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins. , 2003, Biophysical journal.

[25]  F E Cohen,et al.  Helix‐helix packing angle preferences for finite helix axes , 1998, Proteins.

[26]  Lukas K. Tamm,et al.  Membrane structure and fusion-triggering conformational change of the fusion domain from influenza hemagglutinin , 2001, Nature Structural Biology.

[27]  W. DeGrado,et al.  Polar side chains drive the association of model transmembrane peptides. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[28]  D. Engelman,et al.  Sequence motifs, polar interactions and conformational changes in helical membrane proteins. , 2003, Current opinion in structural biology.

[29]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[30]  P. Strevens.  Iii , 1985 .

[31]  Johan Åqvist,et al.  A simple way to calculate the axis of an -helix , 1986, Comput. Chem..

[32]  Alessandro Senes,et al.  Folding of helical membrane proteins: the role of polar, GxxxG-like and proline motifs. , 2004, Current opinion in structural biology.

[33]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[34]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[35]  T. Cross,et al.  Transmembrane four-helix bundle of influenza A M2 protein channel: structural implications from helix tilt and orientation. , 1997, Biophysical journal.

[36]  B. Roux,et al.  Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[37]  D. Engelman,et al.  Specificity in transmembrane helix–helix interactions can define a hierarchy of stability for sequence variants , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[38]  William F. DeGrado,et al.  Asparagine-mediated self-association of a model transmembrane helix , 2000, Nature Structural Biology.

[39]  G M Anantharamaiah,et al.  Correction for Genetic inactivation of AKT1, AKT2, and PDPK1 in human colorectal cancer cells clarifies their roles in tumor growth regulation , 2010 .

[40]  J Wang,et al.  Structure of the transmembrane region of the M2 protein H+ channel , 2001, Protein science : a publication of the Protein Society.

[41]  C. Chothia,et al.  Helix to helix packing in proteins. , 1981, Journal of molecular biology.

[42]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[43]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[44]  Alessandro Senes,et al.  The Cα—H⋅⋅⋅O hydrogen bond: A determinant of stability and specificity in transmembrane helix interactions , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[45]  D. Engelman,et al.  Polar residues drive association of polyleucine transmembrane helices , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Charles L. Brooks,et al.  Application of torsion angle molecular dynamics for efficient sampling of protein conformations , 2005, J. Comput. Chem..