Position of helical kinks in membrane protein crystal structures and the accuracy of computational prediction.

The structural features of helical transmembrane (TM) proteins, such as helical kinks, tilts, and rotational orientations are important in modulation of their function and these structural features give rise to functional diversity in membrane proteins with similar topology. In particular, the helical kinks caused by breaking of the backbone hydrogen bonds lead to hinge bending flexibility in these helices. Therefore it is important to understand the nature of the helical kinks and to be able to reproduce these kinks in structural models of membrane proteins. We have analyzed the position and extent of helical kinks in the transmembrane helices of all the crystal structures of membrane proteins taken from the MPtopo database, which are about 405 individual helices of length between 19 and 35 residues. 44% of the crystal structures of TM helices showed a significant helical kink, and 35% of these kinks are caused by prolines. Many of the non-proline helical kinks are caused by other residues like Ser and Gly that are located at the center of helical kinks. The side chain of Ser makes a hydrogen bond with the main chain carbonyl of the i - 4th or i + 4th residue thus making a kink. We have also studied how well molecular dynamics (MD) simulations on isolated helices can reproduce the position of the helical kinks in TM helices. Such a method is useful for structure prediction of membrane proteins. We performed MD simulations, starting from a canonical helix for the 405 TM helices. 1 ns of MD simulation results show that we can reproduce about 79% of the proline kinks, only 59% of the vestigial proline kinks and 18% of the non-proline helical kinks. We found that similar results can be obtained from choosing the lowest potential energy structure from the MD simulation. 4-14% more of the vestigial prolines were reproduced by replacing them with prolines before performing MD simulations, and changing the amino acid back to proline after the MD simulations. From these results we conclude that the position of the helical kinks is dependent on the TM sequence. However the extent of helical kinking may depend on the packing of the rest of the protein and the lipid bilayer.

[1]  R. Friesner,et al.  Generalized Born Model Based on a Surface Integral Formulation , 1998 .

[2]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[3]  Peter L. Freddolino,et al.  Prediction of structure and function of G protein-coupled receptors , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[4]  D. Engelman,et al.  The GxxxG motif: a framework for transmembrane helix-helix association. , 2000, Journal of molecular biology.

[5]  S. L. Mayo,et al.  DREIDING: A generic force field for molecular simulations , 1990 .

[6]  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.

[7]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[8]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[9]  Thomas L. Madden,et al.  PowerBLAST: a new network BLAST application for interactive or automated sequence analysis and annotation. , 1997, Genome research.

[10]  S H White,et al.  MPtopo: A database of membrane protein topology , 2001, Protein science : a publication of the Protein Society.

[11]  Rodrigo Lopez,et al.  Multiple sequence alignment with the Clustal series of programs , 2003, Nucleic Acids Res..

[12]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

[13]  H Weinstein,et al.  Prokink: a protocol for numerical evaluation of helix distortions by proline. , 2000, Protein engineering.

[14]  Adrian A Canutescu,et al.  Access the most recent version at doi: 10.1110/ps.03154503 References , 2003 .

[15]  P. Y. Chou,et al.  Prediction of protein conformation. , 1974, Biochemistry.

[16]  Sharon Brunett,et al.  Molecular dynamics for very large systems on massively parallel computers: The MPSim program , 1997 .

[17]  Nagarajan Vaidehi,et al.  First principles predictions of the structure and function of g-protein-coupled receptors: validation for bovine rhodopsin. , 2004, Biophysical journal.

[18]  L. Pardo,et al.  Ser and Thr residues modulate the conformation of pro-kinked transmembrane alpha-helices. , 2004, Biophysical journal.

[19]  Duan Yang,et al.  The evolution of transmembrane helix kinks and the structural diversity of G protein-coupled receptors. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[20]  J. Valjakka,et al.  Unreliability of the Chou-Fasman parameters in predicting protein secondary structure. , 1998, Protein engineering.

[21]  Nagarajan Vaidehi,et al.  Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors. , 2008, Journal of molecular biology.

[22]  Giovanni Colonna,et al.  Amino acid propensities for secondary structures are influenced by the protein structural class. , 2006, Biochemical and biophysical research communications.

[23]  W. DeGrado,et al.  Helix-packing motifs in membrane proteins , 2006, Proceedings of the National Academy of Sciences.

[24]  Leonardo Pardo,et al.  Ser and Thr Residues Modulate the Conformation of Pro-Kinked Transmembrane α-Helices , 2004 .

[25]  R. Horuk,et al.  I want a new drug: G-protein-coupled receptors in drug development. , 2006, Drug discovery today.