1
|
Kim MH, Seo S, Jeong JI, Kim BJ, Liu WK, Lim BS, Choi JB, Kim MK. A mass weighted chemical elastic network model elucidates closed form domain motions in proteins. Protein Sci 2013; 22:605-13. [PMID: 23456820 DOI: 10.1002/pro.2244] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2012] [Revised: 02/15/2013] [Accepted: 02/24/2013] [Indexed: 12/18/2022]
Abstract
An elastic network model (ENM), usually Cα coarse-grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass-weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well-known proteins of which both closed and open conformations are available as well as three α-helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.
Collapse
Affiliation(s)
- Min Hyeok Kim
- SKKU Advanced Institute of Nanotechnology-SAINT, Sungkyunkwan University, Suwon 440-746, Korea
| | | | | | | | | | | | | | | |
Collapse
|
2
|
Qian P, Seo S, Kim J, Kim S, Lim BS, Liu WK, Kim BJ, LaBean TH, Park SH, Kim MK. DNA nanotube formation based on normal mode analysis. NANOTECHNOLOGY 2012; 23:105704. [PMID: 22361575 DOI: 10.1088/0957-4484/23/10/105704] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Ever since its inception, a popular DNA motif called the cross tile has been recognized to self-assemble into addressable 2D templates consisting of periodic square cavities. Although this may be conceptually correct, in reality certain types of cross tiles can only form planar lattices if adjacent tiles are designed to bind in a corrugated manner, in the absence of which they roll up to form 3D nanotube structures. Here we present a theoretical study on why uncorrugated cross tiles self-assemble into counterintuitive 3D nanotube structures and not planar 2D lattices. Coarse-grained normal mode analysis of single and multiple cross tiles within the elastic network model was carried out to expound the vibration modes of the systems. While both single and multiple cross tile simulations produce results conducive to tube formations, the dominant modes of a unit of four cross tiles (one square cavity), termed a quadruplet, fully reflect the symmetries of the actual nanotubes found in experiments and firmly endorse circularization of an array of cross tiles.
Collapse
Affiliation(s)
- PengFei Qian
- School of Mechanical Engineering, Sungkyunkwan University, Suwon, Republic of Korea
| | | | | | | | | | | | | | | | | | | |
Collapse
|
4
|
Schuyler AD, Jernigan RL, Qasba PK, Ramakrishnan B, Chirikjian GS. Iterative cluster-NMA: A tool for generating conformational transitions in proteins. Proteins 2009; 74:760-76. [PMID: 18712827 DOI: 10.1002/prot.22200] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.
Collapse
Affiliation(s)
- Adam D Schuyler
- Department of Neurology, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | | | | | | |
Collapse
|
5
|
Peeters K, Taormina A. Group theory of icosahedral virus capsid vibrations: a top-down approach. J Theor Biol 2008; 256:607-24. [PMID: 19014954 DOI: 10.1016/j.jtbi.2008.10.019] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2008] [Revised: 09/11/2008] [Accepted: 10/12/2008] [Indexed: 11/18/2022]
Abstract
We explore the use of a top-down approach to analyse the dynamics of icosahedral virus capsids and complement the information obtained from bottom-up studies of viral vibrations available in the literature. A normal mode analysis based on protein association energies is used to study the frequency spectrum, in which we reveal a universal plateau of low-frequency modes shared by a large class of Caspar-Klug capsids. These modes break icosahedral symmetry and are potentially relevant to the genome release mechanism. We comment on the role of viral tiling theory in such dynamical considerations.
Collapse
Affiliation(s)
- Kasper Peeters
- Institute for Theoretical Physics, Utrecht University, P.O. Box 80.195, 3508 TD Utrecht, The Netherlands.
| | | |
Collapse
|
6
|
Englert F, Peeters K, Taormina A. Twenty-four near-instabilities of Caspar-Klug viruses. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:031908. [PMID: 18851066 DOI: 10.1103/physreve.78.031908] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2008] [Indexed: 05/26/2023]
Abstract
Group theoretical arguments combined with normal mode analysis techniques are applied to a coarse-grained approximation of icosahedral viral capsids which incorporates areas of variable flexibility. This highlights a remarkable structure of the low-frequency spectrum in this approximation, namely, the existence of a plateau of 24 near zero modes with universal group theory content.
Collapse
Affiliation(s)
- François Englert
- Service de Physique Théorique and The International Solvay Institutes, Université Libre de Bruxelles, Belgium
| | | | | |
Collapse
|
7
|
Maguid S, Fernandez-Alberti S, Echave J. Evolutionary conservation of protein vibrational dynamics. Gene 2008; 422:7-13. [PMID: 18577430 DOI: 10.1016/j.gene.2008.06.002] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
The aim of the present work is to study the evolutionary divergence of vibrational protein dynamics. To this end, we used the Gaussian Network Model to perform a systematic analysis of normal mode conservation on a large dataset of proteins classified into homologous sets of family pairs and superfamily pairs. We found that the lowest most collective normal modes are the most conserved ones. More precisely, there is, on average, a linear correlation between normal mode conservation and mode collectivity. These results imply that the previously observed conservation of backbone flexibility (B-factor) profiles is due to the conservation of the most collective modes, which contribute the most to such profiles. We discuss the possible roles of normal mode robustness and natural selection in the determination of the observed behavior. Finally, we draw some practical implications for dynamics-based protein alignment and classification and discuss possible caveats of the present approach.
Collapse
Affiliation(s)
- Sandra Maguid
- Centro de Estudios e Investigaciones, Universidad Nacional de Quilmes, Bernal, Argentina
| | | | | |
Collapse
|
9
|
Rueda M, Chacón P, Orozco M. Thorough validation of protein normal mode analysis: a comparative study with essential dynamics. Structure 2007; 15:565-75. [PMID: 17502102 DOI: 10.1016/j.str.2007.03.013] [Citation(s) in RCA: 133] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2007] [Revised: 03/27/2007] [Accepted: 03/29/2007] [Indexed: 11/19/2022]
Abstract
The deformation patterns of a large set of representative proteins determined by essential dynamics extracted from atomistic simulations and coarse-grained normal mode analysis are compared. Our analysis shows that the deformational space obtained with both approaches is quite similar when taking into account a representative number of modes. The results provide not only a comprehensive validation of the use of a low-frequency modal spectrum to describe protein flexibility, but also a complete picture of normal mode limitations.
Collapse
Affiliation(s)
- Manuel Rueda
- Molecular Modeling and Bioinformatics Unit, Institut de Recerca Biomèdica, Parc Cientific de Barcelona, 08028 Barcelona, Spain
| | | | | |
Collapse
|
10
|
Franklin J, Koehl P, Doniach S, Delarue M. MinActionPath: maximum likelihood trajectory for large-scale structural transitions in a coarse-grained locally harmonic energy landscape. Nucleic Acids Res 2007; 35:W477-82. [PMID: 17545201 PMCID: PMC1933200 DOI: 10.1093/nar/gkm342] [Citation(s) in RCA: 68] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The non-linear problem of simulating the structural transition between two known forms of a macromolecule still remains a challenge in structural biology. The problem is usually addressed in an approximate way using 'morphing' techniques, which are linear interpolations of either the Cartesian or the internal coordinates between the initial and end states, followed by energy minimization. Here we describe a web tool that implements a new method to calculate the most probable trajectory that is exact for harmonic potentials; as an illustration of the method, the classical Calpha-based Elastic Network Model (ENM) is used both for the initial and the final states but other variants of the ENM are also possible. The Langevin equation under this potential is solved analytically using the Onsager and Machlup action minimization formalism on each side of the transition, thus replacing the original non-linear problem by a pair of linear differential equations joined by a non-linear boundary matching condition. The crossover between the two multidimensional energy curves around each state is found numerically using an iterative approach, producing the most probable trajectory and fully characterizing the transition state and its energy. Jobs calculating such trajectories can be submitted on-line at: http://lorentz.dynstr.pasteur.fr/joel/index.php.
Collapse
Affiliation(s)
- Joel Franklin
- Department of Physics, Reed College, Portland, OR 97202, USA, Department of Computer Science and Genome Center, UC Davis, Davis, CA 95616, USA, Departments of Physics and Applied Physics, Stanford University, Stanford, CA 94305-4045, USA and Department of Structural Biology and Chemistry and URA 2185 du C.N.R.S., Institut Pasteur, Paris, France
| | - Patrice Koehl
- Department of Physics, Reed College, Portland, OR 97202, USA, Department of Computer Science and Genome Center, UC Davis, Davis, CA 95616, USA, Departments of Physics and Applied Physics, Stanford University, Stanford, CA 94305-4045, USA and Department of Structural Biology and Chemistry and URA 2185 du C.N.R.S., Institut Pasteur, Paris, France
| | - Sebastian Doniach
- Department of Physics, Reed College, Portland, OR 97202, USA, Department of Computer Science and Genome Center, UC Davis, Davis, CA 95616, USA, Departments of Physics and Applied Physics, Stanford University, Stanford, CA 94305-4045, USA and Department of Structural Biology and Chemistry and URA 2185 du C.N.R.S., Institut Pasteur, Paris, France
| | - Marc Delarue
- Department of Physics, Reed College, Portland, OR 97202, USA, Department of Computer Science and Genome Center, UC Davis, Davis, CA 95616, USA, Departments of Physics and Applied Physics, Stanford University, Stanford, CA 94305-4045, USA and Department of Structural Biology and Chemistry and URA 2185 du C.N.R.S., Institut Pasteur, Paris, France
- *To whom correspondence should be addressed. +33-1-45-688605+33-1-40-613793
| |
Collapse
|