Computational and Mathematical Methods in Medicine

Volume 2012 (2012), Article ID 242903, 3 pages

doi:10.1155/2012/242903

Editorial

Functional Dynamics of Proteins

Lee-Wei Yang,¹ Silvina Matysiak,² Shang-Te Danny Hsu,³ Gabriela Mustata Wilson,⁴ and Yasumasa Joti⁵

¹Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu 300, Taiwan

²Fischell Department of Bioengineering, University of Maryland, College Park, MD 20742, USA

³Institute of Biological Chemistry, Academia Sinica, Taipei 115, Taiwan

⁴Health Services and Health Administration, University of Southern Indiana, Evansville, IN 47712, USA

⁵SPring-8/JASRI, Hyogo 679-5198, Japan

Received 17 December 2012; Accepted 17 December 2012

Copyright © 2012 Lee-Wei Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

X-ray crystallography captures protein conformational states that are stable enough to be captured. These observable states are simply part of the functional event in the context of protein dynamics. Life at the molecular level has a vibrant nature at physiological temperature, governed by physics laws.

Molecular dynamics (MDs) simulations have been used to study biomolecular systems since 1970s [1]. Advancing to this day, it has been a common practice for researchers who access supercomputers to obtain trajectories for hundreds of nanoseconds to two microseconds for small proteins [2, 3] and a few tens of nanoseconds for supramolecular assemblies such as the ribosome of a size larger than 200 Å in explicit solvent with full atomic details [4, 5]. For the few privileged who are able to access machines with build-in hard-wired forcefields, a millisecond trajectory for proteins of ordinary sizes can be made possible [3]. The advancement in parallel computing and hardware speeds enables MD to access timescales of molecular motions that are verifiable by spectroscopy such as IR [6], UVRR [7], FRET [8], and NMR [9]. On the other hand, experimental techniques, especially in X-ray crystallography that has advanced to resolve atomic details for molecules of ever increasing sizes [10], have provided theoreticians a reliable ground to initiate tracking of molecular motions over lengthened time and extended space by leveraging the modern computing environment.

In this issue, we report dynamics stories for biomolecules of different sizes, from crystal waters near the enzyme active sites, “Implication of crystal water molecules in inhibitor binding at ALR2 active site,” HIV-inhibiting peptides “Molecular dynamics simulation of HIV fusion inhibitor T-1249: insights on peptide-lipid interaction,” V2 vasopressin receptor “Membrane protein stability analyses by means of protein energy profiles in case of nephrogenic diabetes insipidus,” ribosomes “Revealing-1 programmed ribosomal frameshifting mechanisms by single-molecule techniques and computational methods,” to “walking” myosins “Coarse-grained simulation of myosin-V movement.” None of these studies can be made complete with static structures alone. In Hymavati’s work, theoretical evidence is presented that proteins bind different ligands (or drug-like molecules) by engaging a variety of water-bridging patterns “Implication of crystal water molecules in inhibitor binding at ALR2 active site.” In the MD study led by Luis Loura, potent HIV-inhibitory peptides are found to concentrate near the cell membrane by engaging hydrogen bond contacts with lipids and cholesterol so as to prevent the fusion of HIV envelope with the membrane and therefore increase the antivirial efficacy “Molecular dynamics simulation of HIV fusion inhibitor T-1249: insights on peptide-lipid interaction.” Heinke and Labudde demonstrate here how MD can help produce reliable receptor models, whereby stability-related energy profiles can be derived and show how genetic mutations in V2R and aquaporin-2 impair their stability and functions, which leads to the disordered phenotype, nephrogenic diabetes insipidus (NDI) “Membrane protein stability analyses by means of protein energy profiles in case of nephrogenic diabetes insipidus”. In similar vein, Zhang et al. offer a comprehensive review on how missense mutations could impact protein stability and dynamics and therefore function “Analyzing effects of naturally occurring missense mutations.”

Despite of the great potential of MD, protein size can easily limit its applicability. It is certainly challenging to simulate macromolecules for timescales of immediate biological interest. The problem can be worsened with inherited sampling problems [11]. The directions of initial momenta of atoms are not specifically assigned (although the magnitude of velocity is deterministic at given temperatures), and therefore a few trajectories have to be obtained to ensure an unbiased sampling.

Coarse-grained (CG) modeling and simulation techniques emerged in recent years as a promising alternative to obtain biologically relevant dynamics for large macromolecules [12–19]. Dynamic nature and functional details of ribosomes, myosins or even virus capsids of a size as large as 675 Å “Coarse-grained simulation of myosin-V movement,” “Analyzing effects of naturally occurring missense mutations” [20] can now be theoretically studied thanks to the development of CG techniques. The computational speed of elastic network models (ENMs), a normal mode analysis for proteins in their CG-ed presentation, can be made five orders of magnitude faster than MD simulations [12–14, 21, 22]. The elastic description of residue-residue interaction reflects the physical nature of folded proteins that sample their surrounding energy landscape [14, 22]. Due to the elimination of fine interaction details, the potential energy surface underlying the conformational transition is smoothened, and the sampling for global dynamics at the expense of atomistic details is enhanced [22, 23]. ENMs have demonstrated that the observed “bound” structures have been readily accessible to the “unbound” structures by ENM-inferred low-frequency normal modes [14, 22, 24], while the “induced fit” can be made apparent to dominate the conformational transitions in strong ligand-protein interactions [25] (Figure 1).

Figure 1: Induced fit versus conformational selection. Protein (in cyan) exists preexisting conformers and “selects” a given conformer by stabilizing it with the incoming ligand (in blue). ENM reveals such transition via its low-frequency normal modes. KNF, Koshland-Néméthy-Filmer; MWC, Monod-Wyman-Changeux.

Aforementioned physical nature can be well explored for supramolecular assemblies with a size of the ribosome by atomistic MD simulations and CG-techniques. Chang gives a comprehensive review in this issue on the possible helicase activity in ribosomes and a phenomenon that is called −1 programmed ribosomal frameshifting (−1 PRF), the functional dynamics of which is a realization of “conformational selection” and “induced fit.” The review covers recent progresses on single molecular and modeling techniques, especially the optical tweezers, fluorescence (Förster) resonance energy transfer (FRET), and ENM, followed by the active uses of these techniques to unravel the mystery of −1 PRF, an elaborate and efficient use of one messenger RNA to produce two or more gene products of totally different amino acid sequences and structures so as to maintain life.

Katsimitsoulia and Taylor herein report a study using coarse-grained hierarchical models, at three different structural resolutions, to simulate Myosin’s “walking” on the actin filaments. The simulations reproduce observed mean length of a processive run for myosin V, estimated at 2.4 m, or approximately 66 steps of 36 nm each. The simulations can potentially model more than a hundred of such macromolecules using only 8–12 cpu cores, as the article suggested. The progress in simulation techniques along with fast and parallel computing facilities is gradually bridging the gaps between molecular and cellular simulations. Before another decade, we may expect to see married molecular dynamics and systems biology simulations developed to bring unprecedented mechanistic insights of life at the molecular and cellular levels.

Lee-Wei Yang

Silvina Matysiak

Shang-Te Danny Hsu

Gabriela Mustata Wilson

Yasumasa Joti

References

J. A. McCammon, B. R. Gelin, and M. Karplus, “Dynamics of folded proteins,” Nature, vol. 267, no. 5612, pp. 585–590, 1977.

P. L. Freddolino, F. Liu, M. Gruebele, and K. Schulten, “Ten-microsecond molecular dynamics simulation of a fast-folding WW domain,” Biophysical Journal, vol. 94, no. 10, pp. L75–L77, 2008.

R. O. Dror, R. M. Dirks, J. P. Grossman, H. Xu, and D. E. Shaw, “Biomolecular simulation: a computational microscope for molecular biology,” Annual Review of Biophysics, vol. 41, pp. 429–452, 2012.

K. Y. Sanbonmatsu, S. Joseph, and C. S. Tung, “Simulating movement of tRNA into the ribosome during decoding,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 44, pp. 15854–15859, 2005.

P. C. Whitford, P. Geggier, R. B. Altman, S. C. Blanchard, J. N. Onuchic, and K. Y. Sanbonmatsu, “Accommodation of aminoacyl-tRNA into the ribosome involves reversible excursions along multiple pathways,” RNA, vol. 16, no. 6, pp. 1196–1204, 2010.

R. Nakamura, K. Machida, M. Oobatake, and S. Hayashi, “Molecular dynamics simulation of infrared spectra and average structure of benzoic acid crystal,” Molecular Physics, vol. 64, no. 2, pp. 215–227, 1988.

A. Sato, Y. Gao, T. Kitagawa, and Y. Mizutani, “Primary protein response after ligand photodissociation in carbonmonoxy myoglobin,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 23, pp. 9627–9632, 2007.

J. B. Munro, A. Vaiana, K. Y. Sanbonmatsu, and S. C. Blanchard, “A new view of protein synthesis: mapping the free energy landscape of the ribosome using single-molecule FRET,” Biopolymers, vol. 89, no. 7, pp. 565–577, 2008.

P. Maragakis, K. Lindorff-Larsen, M. P. Eastwood et al., “Microsecond molecular dynamics simulation shows effect of slow loop dynamics on backbone amide order parameters of proteins,” Journal of Physical Chemistry B, vol. 112, no. 19, pp. 6155–6158, 2008.

F. Schlünzen, R. Zarivach, J. Harms et al., “Structural basis for the interaction of antibiotics with the peptidyl transferase centre in eubacteria,” Nature, vol. 413, no. 6858, pp. 814–821, 2001.

J. B. Clarage, T. Romo, B. K. Andrews, B. M. Pettitt, and G. N. Phillips Jr., “A sampling problem in molecular dynamics simulations of macromolecules,” Proceedings of the National Academy of Sciences of the United States of America, vol. 92, no. 8, pp. 3288–3292, 1995.

M. M. Tirion, “Large amplitude elastic motions in proteins from a single-parameter, atomic analysis,” Physical Review Letters, vol. 77, no. 9, pp. 1905–1908, 1996.

I. Bahar, A. R. Atilgan, and B. Erman, “Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential,” Folding and Design, vol. 2, no. 3, pp. 173–181, 1997.

I. Bahar, T. R. Lezon, L. W. Yang, and E. Eyal, “Global dynamics of proteins: bridging between structure and function,” Annual Review of Biophysics, vol. 39, no. 1, pp. 23–42, 2010.

C. P. Chng and L. W. Yang, “Coarse-grained models reveal functional dynamics—II. Molecular dynamics simulation at the coarse-grained level—theories and biological applications,” Bioinformatics and Biology Insights, vol. 2, pp. 171–185, 2008.

H. Kenzaki, N. Koga, N. Hori et al., “CafeMol: a coarse-grained biomolecular simulator for simulating proteins at work,” Journal of Chemical Theory and Computation, vol. 7, no. 6, pp. 1979–1989, 2011.

P. Das and S. Matysiak, “Direct characterization of hydrophobic hydration during cold and pressure denaturation,” The Journal of Physical Chemistry B, vol. 116, no. 18, pp. 5342–5348, 2012.

S. Matysiak and C. Clementi, “Minimalist protein model as a diagnostic tool for misfolding and aggregation,” Journal of Molecular Biology, vol. 363, no. 1, pp. 297–308, 2006.

S. Matysiak and C. Clementi, “Optimal combination of theory and experiment for the characterization of the protein folding landscape of S6: how far can a minimalist model go?” Journal of Molecular Biology, vol. 343, no. 1, pp. 235–248, 2004.

M. K. Kim, R. L. Jernigan, and G. S. Chirikjian, “An elastic network model of HK97 capsid maturation,” Journal of Structural Biology, vol. 143, no. 2, pp. 107–117, 2003.

A. J. Rader, D. H. Vlad, and I. Bahar, “Maturation dynamics of bacteriophage HK97 capsid,” Structure, vol. 13, no. 3, pp. 413–421, 2005.

L. W. Yang and C. P. Chng, “Coarse-grained models reveal functional dynamics—I. Elastic network models–theories, comparisons and perspectives,” Bioinformatics and Biology Insights, vol. 2, pp. 25–45, 2008.

J. Kim, J. E. Straub, and T. Keyes, “Statistical-temperature Monte Carlo and molecular dynamics algorithms,” Physical Review Letters, vol. 97, Article ID 050601, 2006.

D. Tobi and I. Bahar, “Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 52, pp. 18908–18913, 2005.

K. I. Okazaki and S. Takada, “Dynamic energy landscape view of coupled binding and protein conformational change: induced-fit versus population-shift mechanisms,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 32, pp. 11182–11187, 2008.