Helix motion in protein C12A‐p8MTCP1: Comparison of molecular dynamics simulations and multifield NMR relaxation data

The human p8MTCP1 protein is constituted by an original disulfide bridged α‐hairpin motif, and a third hydrophilic helix that appeared mobile and independent in NMR analysis. To get atomic scale description of the possible motions involved, a 12‐ns molecular dynamics simulation of the C12A‐p8MTCP1 mutated protein and 10 400‐ps simulations were performed in explicit water at 298 K. Analyses of the essential dynamics subspace and of the time‐evolution of secondary structures indicate large displacements and internal motions of the third helix. Simulated order parameters, S,2 are consistent to those obtained from 15N‐NMR relaxation data recorded at five different magnetic fields (proton resonance frequencies: 400, 500, 600, 700, and 800 MHz). The simulations show a fairly rigid α‐hairpin motif with only very small local conformational changes, although the loop between helices undergoes significant correlated motions. The displacements and internal motions of the third helix can be essentially accounted for by α‐helical to 310‐helical transitions near its N‐terminus and by a hinge‐bending motion near glycine 52. The loss of the α‐helical conformation for a large part of the third helix during about 25% of the simulation time supports the possible occurrence of a slow motion as suggested by previous analyses of the 3JNH‐Hα coupling constant values and of the proton‐deuterium NH exchange kinetics. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1577–1586, 2002

[1]  M. Karplus,et al.  Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin , 1998, Protein science : a publication of the Protein Society.

[2]  D. Case,et al.  Molecular Dynamics Analysis of NMR Relaxation in a Zinc-Finger Peptide , 1992 .

[3]  Alexey K. Mazur,et al.  The Hierarchy of Fast Motions in Protein Dynamics , 1998, physics/9801022.

[4]  L. Li,et al.  Role of gamma-carboxyglutamic acid in the calcium-induced structural transition of conantokin G, a conotoxin from the marine snail Conus geographus. , 1997, Biochemistry.

[5]  A. Pohorille,et al.  The development/application of a ‘minimalist’ organic/biochemical molecular mechanic force field using a combination of ab initio calculations and experimental data , 1997 .

[6]  G. Lipari Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules , 1982 .

[7]  M. Sundaralingam,et al.  Water-inserted alpha-helical segments implicate reverse turns as folding intermediates. , 1989, Science.

[8]  K. Wüthrich,et al.  Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long‐time molecular dynamics calculations in explicit water , 1996, Proteins.

[9]  J B Findlay,et al.  Protein dynamics derived from clusters of crystal structures. , 1997, Biophysical journal.

[10]  Eric Westhof,et al.  MULTIPLE MOLECULAR DYNAMICS SIMULATIONS OF THE ANTICODON LOOP OF TRNAASP IN AQUEOUS SOLUTION WITH COUNTERIONS , 1995 .

[11]  H. Kaback,et al.  Engineering conformational flexibility in the lactose permease of Escherichia coli: use of glycine-scanning mutagenesis to rescue mutant Glu325-->Asp. , 2001, Biochemistry.

[12]  A. Szabó,et al.  Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity , 1982 .

[13]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[14]  A M Lesk,et al.  Domain closure in mitochondrial aspartate aminotransferase. , 1992, Journal of molecular biology.

[15]  C. Brooks,et al.  310 Helixes in Peptides and Proteins As Studied by Modified Zimm-Bragg Theory , 1995 .

[16]  D. Cafiso,et al.  The role of proline and glycine in determining the backbone flexibility of a channel-forming peptide. , 1999, Biophysical journal.

[17]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[18]  J. Lefèvre,et al.  Refined solution structure and backbone dynamics of 15N-labeled C12A-p8MTCP studied by NMR relaxation , 1999, Journal of biomolecular NMR.

[19]  M Gerstein,et al.  Analysis of protein loop closure. Two types of hinges produce one motion in lactate dehydrogenase. , 1991, Journal of molecular biology.

[20]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[21]  F. Sigaux,et al.  MTCP-1: a novel gene on the human chromosome Xq28 translocated to the T cell receptor alpha/delta locus in mature T cell proliferations. , 1993, Oncogene.

[22]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[23]  A. Amadei,et al.  The essential dynamics of Cu, Zn superoxide dismutase: suggestion of intersubunit communication. , 1997, Biophysical journal.

[24]  E. Westhof,et al.  RNA hydration: three nanoseconds of multiple molecular dynamics simulations of the solvated tRNA(Asp) anticodon hairpin. , 1997, Journal of molecular biology.

[25]  P. Barthe,et al.  Synthesis and NMR solution structure of an α‐helical hairpin stapled with two disulfide bridges , 2000, Protein science : a publication of the Protein Society.

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

[27]  P. Kollman,et al.  The application of different solvation and electrostatic models in molecular dynamics simulations of ubiquitin: How well is the x‐ray structure “maintained”? , 1996, Proteins.

[28]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[29]  L. Nilsson,et al.  Structural fluctuations between two conformational states of a transmembrane helical peptide are related to its channel-forming properties in planar lipid membranes. , 1993, European journal of biochemistry.

[30]  C. Toniolo,et al.  The polypeptide 310-helix. , 1991, Trends in biochemical sciences.

[31]  M. Dyer,et al.  The chromosomal translocation t(X;14)(q28;q11) in T-cell pro-lymphocytic leukaemia breaks within one gene and activates another. , 1993, Oncogene.

[32]  J. Soulier,et al.  Solution structure of human p8MTCP1, a cysteine-rich protein encoded by the MTCP1 oncogene, reveals a new alpha-helical assembly motif. , 1997, Journal of molecular biology.

[33]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[34]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[35]  J. Thornton,et al.  Helix geometry in proteins. , 1988, Journal of molecular biology.