Biophys J, June 2002, p. 3063-3071, Vol. 82, No. 6

Contribution of Energy Values to the Analysis of Global Searching Molecular Dynamics Simulations of Transmembrane Helical Bundles

Jaume Torres,^* John A. G. Briggs, and Isaiah T. Arkin

 ^*Cambridge Centre for Molecular Recognition, Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom;  Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom; and  The Alexander Silberman Institute of Life Sciences, Department of Biological Chemistry, The Hebrew University of Jerusalem, Givat-Ram, Jerusalem 91904, Israel

Molecular interactions between transmembrane -helices can be explored using global searching molecular dynamics simulations (GSMDS), a method that produces a group of probable low energy structures. We have shown previously that the correct model in various homooligomers is always located at the bottom of one of various possible energy basins. Unfortunately, the correct model is not necessarily the one with the lowest energy according to the computational protocol, which has resulted in overlooking of this parameter in favor of experimental data. In an attempt to use energetic considerations in the aforementioned analysis, we used global searching molecular dynamics simulations on three homooligomers of different sizes, the structures of which are known. As expected, our results show that even when the conformational space searched includes the correct structure, taking together simulations using both left and right handednesses, the correct model does not necessarily have the lowest energy. However, for the models derived from the simulation that uses the correct handedness, the lowest energy model is always at, or very close to, the correct orientation. We hypothesize that this should also be true when simulations are performed using homologous sequences, and consequently lowest energy models with the right handedness should produce a cluster around a certain orientation. In contrast, using the wrong handedness the lowest energy structures for each sequence should appear at many different orientations. The rationale behind this is that, although more than one energy basin may exist, basins that do not contain the correct model will shift or disappear because they will be destabilized by at least one conservative (i.e. silent) mutation, whereas the basin containing the correct model will remain. This not only allows one to point to the possible handedness of the bundle, but can be used to overcome ambiguities arising from the use of homologous sequences in the analysis of global searching molecular dynamics simulations. In addition, because clustering of lowest energy models arising from homologous sequences only happens when the estimation of the helix tilt is correct, it may provide a validation for the helix tilt estimate.

Biophys J, June 2002, p. 3063-3071, Vol. 82, No. 6
© 2002 by the Biophysical Society   0006-3495/02/06/3063/09  $2.00

This article has been cited by other articles:

				S. W. Gan, L. Ng, X. Lin, X. Gong, and J. Torres Structure and ion channel activity of the human respiratory syncytial virus (hRSV) small hydrophobic protein transmembrane domain Protein Sci., May 1, 2008; 17(5): 813 - 820. [Abstract] [Full Text] [PDF]

				L. Bu, W. Im, and C. L. Brooks III Membrane Assembly of Simple Helix Homo-Oligomers Studied via Molecular Dynamics Simulations Biophys. J., February 1, 2007; 92(3): 854 - 863. [Abstract] [Full Text] [PDF]

				J. Torres, J. Wang, K. Parthasarathy, and D. X. Liu The Transmembrane Oligomers of Coronavirus Protein E Biophys. J., February 1, 2005; 88(2): 1283 - 1290. [Abstract] [Full Text] [PDF]

				W. Im, M. Feig, and C. L. Brooks III An Implicit Membrane Generalized Born Theory for the Study of Structure, Stability, and Interactions of Membrane Proteins Biophys. J., November 1, 2003; 85(5): 2900 - 2918. [Abstract] [Full Text] [PDF]