1
|
Kunzmann P, Krumbach JH, Saponaro A, Moroni A, Thiel G, Hamacher K. Anisotropic Network Analysis of Open/Closed HCN4 Channel Advocates Asymmetric Subunit Cooperativity in cAMP Modulation of Gating. J Chem Inf Model 2024; 64:4727-4738. [PMID: 38830626 PMCID: PMC11203669 DOI: 10.1021/acs.jcim.4c00360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 05/07/2024] [Accepted: 05/10/2024] [Indexed: 06/05/2024]
Abstract
Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are opened in an allosteric manner by membrane hyperpolarization and cyclic nucleotides such as cAMP. Because of conflicting reports from experimental studies on whether cAMP binding to the four available binding sites in the channel tetramer operates cooperatively in gating, we employ here a computational approach as a promising route to examine ligand-induced conformational changes after binding to individual sites. By combining an elastic network model (ENM) with linear response theory (LRT) for modeling the apo-holo transition of the cyclic nucleotide-binding domain (CNBD) in HCN channels, we observe a distinct pattern of cooperativity matching the "positive-negative-positive" cooperativity reported from functional studies. This cooperativity pattern is highly conserved among HCN subtypes (HCN4, HCN1), but only to a lesser extent visible in structurally related channels, which are only gated by voltage (KAT1) or cyclic nucleotides (TAX4). This suggests an inherent cooperativity between subunits in HCN channels as part of a ligand-triggered gating mechanism in these channels.
Collapse
Affiliation(s)
- Patrick Kunzmann
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Jan H. Krumbach
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Andrea Saponaro
- Department
of Pharmacology and Biomolecular Sciences, University of Milano, via Balzaretti 9, 20133 Milano, Italy
| | - Anna Moroni
- Department
of Biosciences, Ion Channel Biophysics, University of Milan, via Celoria 26, 20133 Milan, Italy
| | - Gerhard Thiel
- Department
of Biology, Membrane Biophysics, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| | - Kay Hamacher
- Department
of Biology, Computational Biology & Simulation, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
- Centre
for Synthetic Biology, TU Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt, Germany
| |
Collapse
|
2
|
Khan MW, Murali A. Normal mode analysis and comparative study of intrinsic dynamics of alcohol oxidase enzymes from GMC protein family. J Biomol Struct Dyn 2023:1-16. [PMID: 37676256 DOI: 10.1080/07391102.2023.2255275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 08/30/2023] [Indexed: 09/08/2023]
Abstract
Glucose-Methanol-Choline (GMC) family enzymes are very important in catalyzing the oxidation of a wide range of structurally diverse substrates. Enzymes that constitute the GMC family, share a common tertiary fold but < 25% sequence identity. Cofactor FAD, FAD binding signature motif, and similar structural scaffold of the active site are common features of oxidoreductase enzymes of the GMC family. Protein functionality mainly depends on protein three-dimensional structures and dynamics. In this study, we used the normal mode analysis method to search the intrinsic dynamics of GMC family enzymes. We have explored the dynamical behavior of enzymes with unique substrate catabolism and active site characteristics from different classes of the GMC family. Analysis of individual enzymes and comparative ensemble analysis of enzymes from different classes has shown conserved dynamic motion at FAD binding sites. The present study revealed that GMC enzymes share a strong dynamic similarity (Bhattacharyya coefficient >90% and root mean squared inner product >52%) despite low sequence identity across the GMC family enzymes. The study predicts that local deformation energy between atoms of the enzyme may be responsible for the catalysis of different substrates. This study may help that intrinsic dynamics can be used to make meaningful classifications of proteins or enzymes from different organisms.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Mohammad Wahab Khan
- Department of Bioinformatics, School of Life Science, Pondicherry University, Puducherry, India
| | - Ayaluru Murali
- Department of Bioinformatics, School of Life Science, Pondicherry University, Puducherry, India
| |
Collapse
|
3
|
Sabei A, Caldas Baia TG, Saffar R, Martin J, Frezza E. Internal Normal Mode Analysis Applied to RNA Flexibility and Conformational Changes. J Chem Inf Model 2023; 63:2554-2572. [PMID: 36972178 DOI: 10.1021/acs.jcim.2c01509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
Abstract
We investigated the capability of internal normal modes to reproduce RNA flexibility and predict observed RNA conformational changes and, notably, those induced by the formation of RNA-protein and RNA-ligand complexes. Here, we extended our iNMA approach developed for proteins to study RNA molecules using a simplified representation of the RNA structure and its potential energy. Three data sets were also created to investigate different aspects. Despite all the approximations, our study shows that iNMA is a suitable method to take into account RNA flexibility and describe its conformational changes opening the route to its applicability in any integrative approach where these properties are crucial.
Collapse
|
4
|
Costa MGS, Batista PR, Gomes A, Bastos LS, Louet M, Floquet N, Bisch PM, Perahia D. MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments. J Chem Theory Comput 2023; 19:412-425. [PMID: 36622950 DOI: 10.1021/acs.jctc.2c00599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.
Collapse
Affiliation(s)
- Mauricio G S Costa
- Programa de Computação Científica, Vice-Presidência de Educação Informação e Comunicação, Fundação Oswaldo Cruz, Av. Brasil 4365, Residência Oficial, Manguinhos, 21040-900Rio de Janeiro, Brasil
- Laboratoire de Biologie et de Pharmacologie Appliquée (LBPA), UMR 8113, CNRS, École Normale Supérieure Paris-Saclay, 4 Avenue des Sciences, 91190Gif-sur-Yvette, France
| | - Paulo R Batista
- Programa de Computação Científica, Vice-Presidência de Educação Informação e Comunicação, Fundação Oswaldo Cruz, Av. Brasil 4365, Residência Oficial, Manguinhos, 21040-900Rio de Janeiro, Brasil
| | - Antoniel Gomes
- Laboratório de Física Biológica, Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro21941-902, Brasil
| | - Leonardo S Bastos
- Programa de Computação Científica, Vice-Presidência de Educação Informação e Comunicação, Fundação Oswaldo Cruz, Av. Brasil 4365, Residência Oficial, Manguinhos, 21040-900Rio de Janeiro, Brasil
| | - Maxime Louet
- Institut des Biomolecules Max Mousseron, UMR5247, CNRS, Université De Montpellier, ENSCM, 1919 Route de Mende, Montpellier, Cedex 0534095, France
| | - Nicolas Floquet
- Institut des Biomolecules Max Mousseron, UMR5247, CNRS, Université De Montpellier, ENSCM, 1919 Route de Mende, Montpellier, Cedex 0534095, France
| | - Paulo M Bisch
- Laboratório de Física Biológica, Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro21941-902, Brasil
| | - David Perahia
- Laboratoire de Biologie et de Pharmacologie Appliquée (LBPA), UMR 8113, CNRS, École Normale Supérieure Paris-Saclay, 4 Avenue des Sciences, 91190Gif-sur-Yvette, France
| |
Collapse
|
5
|
Dasgupta B, Tiwari SP. Explicit versus implicit consideration of binding partners in protein-protein complex to elucidate intrinsic dynamics. Biophys Rev 2022; 14:1379-1392. [PMID: 36659985 PMCID: PMC9842844 DOI: 10.1007/s12551-022-01026-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Accepted: 11/27/2022] [Indexed: 12/23/2022] Open
Abstract
The binding of many proteins to their protein partners is tightly regulated via control of their relative intrinsic dynamics during the binding process, a phenomenon which can in turn be modulated. Therefore, investigating the intrinsic dynamics of proteins is necessary to understand function in a comprehensive way. By intrinsic dynamics herein, we principally refer to the vibrational signature of a protein molecule popularly obtained from normal modes or essential modes. For normal modes, one often considers that the molecule under investigation is a collection of springs in a solvent-free or implicit-solvent medium. In the context of a protein-binding partner, the analysis of vibration of the target protein is often complicated due to molecular interaction within the complex. Generally, it is assumed that the isolated bound conformation of the target protein captures the implicit effect of the binding partner on the intrinsic dynamics, therefore suggesting that any influence of the partner molecule is also already integrated. Such an assumption allows large-scale studies of the conservation of protein flexibility. However, in cases where a partner protein directly influences the vibration of the target via critical contacts at the protein-protein interface, the above assumption falls short of providing a detailed view. In this review article, we discuss the implications of considering the dynamics of a protein in a protein-protein complex, as modelled implicitly and explicitly with methods dependent on elastic network models. We further propose how such an explicit consideration can be applied to understand critical protein-protein contacts that can be targeted in future studies.
Collapse
Affiliation(s)
- Bhaskar Dasgupta
- Research Center for Advanced Science and Technology, University of Tokyo, 4-6-1 Komaba, Meguro-Ku, Tokyo, 153-8904 Japan
| | - Sandhya P. Tiwari
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashi-Hiroshima City, 1-3-1 Kagamiyama, Hiroshima, 739-8526 Japan
- Present Address: Institute of Protein Research, Osaka University, 3-2 Yamadaoka, Suita-Shi, Osaka, 565-0871 Japan
| |
Collapse
|
6
|
Singh O, Venugopal PP, Mathur A, Chakraborty D. Exploring the multiple conformational states of RNA genome through interhelical dynamics and network analysis. J Mol Graph Model 2022; 116:108264. [PMID: 35820344 DOI: 10.1016/j.jmgm.2022.108264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2022] [Revised: 06/24/2022] [Accepted: 06/27/2022] [Indexed: 12/14/2022]
Abstract
The structural variation of RNA is often very transient and can be easily missed in experiments. Molecular dynamics simulation studies along with network analysis can be an effective tool to identify prominent conformations of such dynamic biomolecular systems. Here we describe a method to effectively sample different RNA conformations at six different temperatures based on the changes in the interhelical orientations. This method gives the information about prominent states of the RNA as well as the probability of the existence of different conformations and their interconnections during the process of evolution. In the case of the SARS-CoV-2 genome, the change of prominent structures was found to be faster at 333 K as compared to higher temperatures due to the formation of the non-native base pairs. ΔΔG calculated between 288 K and 363 K are found to be 10.31 kcal/mol (88 nt) considering the contribution from the multiple states of the RNA which agrees well with the experimentally reported denaturation energy for E. coli α mRNA pseudoknot (∼16 kcal/mol, 112 nt) determined by calorimetry/UV hyperchromicity and human telomerase RNA telomerase (4.5-6.6 kcal/mol, 54 nt) determined by FRET analysis.
Collapse
Affiliation(s)
- Omkar Singh
- Biophysical and Computational Chemistry Laboratory, Department of Chemistry, National Institute of Technology, Karnataka, 575025, India
| | - Pushyaraga P Venugopal
- Biophysical and Computational Chemistry Laboratory, Department of Chemistry, National Institute of Technology, Karnataka, 575025, India
| | - Apoorva Mathur
- Biophysical and Computational Chemistry Laboratory, Department of Chemistry, National Institute of Technology, Karnataka, 575025, India
| | - Debashree Chakraborty
- Biophysical and Computational Chemistry Laboratory, Department of Chemistry, National Institute of Technology, Karnataka, 575025, India.
| |
Collapse
|
7
|
Yu CC, Raj N, Chu JW. Edge weights in a protein elastic network reorganize collective motions and render long-range sensitivity responses. J Chem Phys 2022; 156:245105. [DOI: 10.1063/5.0095107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effects of inter-residue interactions on protein collective motions are analyzed by comparing two elastic network models (ENM)—structural contact ENM (SC-ENM) and molecular dynamics (MD)-ENM—with the edge weights computed from an all-atom MD trajectory by structure-mechanics statistical learning. A theoretical framework is devised to decompose the eigenvalues of ENM Hessian into contributions from individual springs and to compute the sensitivities of positional fluctuations and covariances to spring constant variation. Our linear perturbation approach quantifies the response mechanisms as softness modulation and orientation shift. All contacts of C α positions in SC-ENM have an identical spring constant by fitting the profile of root-of-mean-squared-fluctuation calculated from an all-atom MD simulation, and the same trajectory data are also used to compute the specific spring constant of each contact as an MD-ENM edge weight. We illustrate that the soft-mode reorganization can be understood in terms of gaining weights along the structural contacts of low elastic strengths and loosing magnitude along those of high rigidities. With the diverse mechanical strengths encoded in protein dynamics, MD-ENM is found to have more pronounced long-range couplings and sensitivity responses with orientation shift identified as a key player in driving the specific residues to have high sensitivities. Furthermore, the responses of perturbing the springs of different residues are found to have asymmetry in the action–reaction relationship. In understanding the mutation effects on protein functional properties, such as long-range communications, our results point in the directions of collective motions as a major effector.
Collapse
Affiliation(s)
- Chieh Cheng Yu
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, 75 Bo-Ai Street, Hsinchu 30010, Taiwan
| | - Nixon Raj
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, 75 Bo-Ai Street, Hsinchu 30010, Taiwan
| | - Jhih-Wei Chu
- Institute of Bioinformatics and Systems Biology, Department of Biological Science and Technology, Institute of Molecular Medicine and Bioengineering, Center for Intelligent Drug Systems and Smart Bio-devices (IDS2B), National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
| |
Collapse
|
8
|
Kaynak BT, Krieger JM, Dudas B, Dahmani ZL, Costa MGS, Balog E, Scott AL, Doruker P, Perahia D, Bahar I. Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods. Front Mol Biosci 2022; 9:832847. [PMID: 35187088 PMCID: PMC8855042 DOI: 10.3389/fmolb.2022.832847] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 01/12/2022] [Indexed: 12/17/2022] Open
Abstract
Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.
Collapse
Affiliation(s)
- Burak T. Kaynak
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - James M. Krieger
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - Balint Dudas
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- Laboratoire de Biologie et Pharmacologie Appliquée, Ecole Normale Supérieure Paris-Saclay, Gif-sur-Yvette, France
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Zakaria L. Dahmani
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
| | - Mauricio G. S. Costa
- Programa de Computação Científica, Vice-Presiden̂cia de Educação, Informação e Comunicação, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Erika Balog
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Ana Ligia Scott
- Laboratory of Bioinformatics and Computational Biology, Center of Mathematics, Computation and Cognition, Federal University of ABC-UFABC, Santo André, Brazil
| | - Pemra Doruker
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| | - David Perahia
- Laboratoire de Biologie et Pharmacologie Appliquée, Ecole Normale Supérieure Paris-Saclay, Gif-sur-Yvette, France
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- *Correspondence: Ivet Bahar, ; David Perahia, ; Pemra Doruker,
| |
Collapse
|
9
|
Koehl P, Orland H, Delarue M. Parameterizing elastic network models to capture the dynamics of proteins. J Comput Chem 2021; 42:1643-1661. [PMID: 34117647 DOI: 10.1002/jcc.26701] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Revised: 12/14/2020] [Accepted: 05/23/2021] [Indexed: 11/09/2022]
Abstract
Coarse-grained normal mode analyses of protein dynamics rely on the idea that the geometry of a protein structure contains enough information for computing its fluctuations around its equilibrium conformation. This geometry is captured in the form of an elastic network (EN), namely a network of edges between its residues. The normal modes of a protein are then identified with the normal modes of its EN. Different approaches have been proposed to construct ENs, focusing on the choice of the edges that they are comprised of, and on their parameterizations by the force constants associated with those edges. Here we propose new tools to guide choices on these two facets of EN. We study first different geometric models for ENs. We compare cutoff-based ENs, whose edges have lengths that are smaller than a cutoff distance, with Delaunay-based ENs and find that the latter provide better representations of the geometry of protein structures. We then derive an analytical method for the parameterization of the EN such that its dynamics leads to atomic fluctuations that agree with experimental B-factors. To limit overfitting, we attach a parameter referred to as flexibility constant to each atom instead of to each edge in the EN. The parameterization is expressed as a non-linear optimization problem whose parameters describe both rigid-body and internal motions. We show that this parameterization leads to improved ENs, whose dynamics mimic MD simulations better than ENs with uniform force constants, and reduces the number of normal modes needed to reproduce functional conformational changes.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Sciences and Genome Center, University of California, Davis, California, USA
| | - Henri Orland
- Institut de Physique Théorique, Université Paris-Saclay, Gif sur Yvette, France
| | - Marc Delarue
- Unité de Dynamique Structurale des Macromolécules, Institut Pasteur, UMR 3528 du CNRS, Paris, France
| |
Collapse
|
10
|
Na H, Hinsen K, Song G. The amounts of thermal vibrations and static disorder in protein X-ray crystallographic B-factors. Proteins 2021; 89:1442-1457. [PMID: 34174110 DOI: 10.1002/prot.26165] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Revised: 05/31/2021] [Accepted: 06/06/2021] [Indexed: 12/20/2022]
Abstract
Crystallographic B-factors provide direct dynamical information on the internal mobility of proteins that is closely linked to function, and are also widely used as a benchmark in assessing elastic network models. A significant question in the field is: what is the exact amount of thermal vibrations in protein crystallographic B-factors? This work sets out to answer this question. First, we carry out a thorough, statistically sound analysis of crystallographic B-factors of over 10 000 structures. Second, by employing a highly accurate all-atom model based on the well-known CHARMM force field, we obtain computationally the magnitudes of thermal vibrations of nearly 1000 structures. Our key findings are: (i) the magnitude of thermal vibrations, surprisingly, is nearly protein-independent, as a corollary to the universality for the vibrational spectra of globular proteins established earlier; (ii) the magnitude of thermal vibrations is small, less than 0.1 Å2 at 100 K; (iii) the percentage of thermal vibrations in B-factors is the lowest at low resolution and low temperature (<10%) but increases to as high as 60% for structures determined at high resolution and at room temperature. The significance of this work is that it provides for the first time, using an extremely large dataset, a thorough analysis of B-factors and their thermal and static disorder components. The results clearly demonstrate that structures determined at high resolution and at room temperature have the richest dynamics information. Since such structures are relatively rare in the PDB database, the work naturally calls for more such structures to be determined experimentally.
Collapse
Affiliation(s)
- Hyuntae Na
- Department of Computer Science, Penn State Harrisburg, Middletown, Pennsylvania, USA
| | - Konrad Hinsen
- Centre de Biophysique Moleculaire, CNRS, Orleans, France.,Synchrotron SOLEIL, Division Expériences, Gif sur Yvette, France
| | - Guang Song
- Department of Computer Science, Program of Bioinformatics and Computational Biology, Iowa State University, Ames, Iowa, USA
| |
Collapse
|
11
|
Song G. Bridging between material properties of proteins and the underlying molecular interactions. PLoS One 2021; 16:e0247147. [PMID: 33951045 PMCID: PMC8099097 DOI: 10.1371/journal.pone.0247147] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2021] [Accepted: 04/19/2021] [Indexed: 11/18/2022] Open
Abstract
In this work, we develop a novel method that bridges between material properties of proteins, particularly the modulus of elasticity, and the underlying molecular interactions. To this end, we employ both an all-atom normal mode analysis (NMA) model with the CHARMM force field and an elastic solid model for proteins and protein interfaces. And the "bridge" between the two models is a common physical property predictable by both models: the magnitude of thermal vibrations. This connection allows one to calibrate the Young's moduli of proteins and protein interface regions. We find that the Young's moduli of proteins are in the range of a few Gpa to 10 Gpa, while the Young's moduli of the interface regions are several times smaller. The work is significant as it represents the first attempt to systematically compute the elastic moduli of proteins from molecular interactions.
Collapse
Affiliation(s)
- Guang Song
- Department of Computer Science, Iowa State University, Ames, IA, United States of America
- Program of Bioinformatics and Computational Biology, Iowa State University, Ames, IA, United States of America
| |
Collapse
|
12
|
Procyk J, Poppleton E, Šulc P. Coarse-grained nucleic acid-protein model for hybrid nanotechnology. SOFT MATTER 2021; 17:3586-3593. [PMID: 33398312 DOI: 10.1039/d0sm01639j] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The emerging field of hybrid DNA-protein nanotechnology brings with it the potential for many novel materials which combine the addressability of DNA nanotechnology with the versatility of protein interactions. However, the design and computational study of these hybrid structures is difficult due to the system sizes involved. To aid in the design and in silico analysis process, we introduce here a coarse-grained DNA/RNA-protein model that extends the oxDNA/oxRNA models of DNA/RNA with a coarse-grained model of proteins based on an anisotropic network model representation. Fully equipped with analysis scripts and visualization, our model aims to facilitate hybrid nanomaterial design towards eventual experimental realization, as well as enabling study of biological complexes. We further demonstrate its usage by simulating DNA-protein nanocage, DNA wrapped around histones, and a nascent RNA in polymerase.
Collapse
Affiliation(s)
- Jonah Procyk
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, Arizona 85281, USA.
| | | | | |
Collapse
|
13
|
Laine E, Grudinin S. HOPMA: Boosting Protein Functional Dynamics with Colored Contact Maps. J Phys Chem B 2021; 125:2577-2588. [PMID: 33687221 DOI: 10.1021/acs.jpcb.0c11633] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
In light of the recent very rapid progress in protein structure prediction, accessing the multitude of functional protein states is becoming more central than ever before. Indeed, proteins are flexible macromolecules, and they often perform their function by switching between different conformations. However, high-resolution experimental techniques such as X-ray crystallography and cryogenic electron microscopy can catch relatively few protein functional states. Many others are only accessible under physiological conditions in solution. Therefore, there is a pressing need to fill this gap with computational approaches. We present HOPMA, a novel method to predict protein functional states and transitions by using a modified elastic network model. The method exploits patterns in a protein contact map, taking its 3D structure as input, and excludes some disconnected patches from the elastic network. Combined with nonlinear normal mode analysis, this strategy boosts the protein conformational space exploration, especially when the input structure is highly constrained, as we demonstrate on a set of more than 400 transitions. Our results let us envision the discovery of new functional conformations, which were unreachable previously, starting from the experimentally known protein structures. The method is computationally efficient and available at https://github.com/elolaine/HOPMA and https://team.inria.fr/nano-d/software/nolb-normal-modes.
Collapse
Affiliation(s)
- Elodie Laine
- CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), Sorbonne Université, 75005 Paris, France
| | - Sergei Grudinin
- CNRS, Inria, Grenoble INP, LJK, Univ. Grenoble Alpes, 38000 Grenoble, France
| |
Collapse
|
14
|
Gong W, Liu Y, Zhao Y, Wang S, Han Z, Li C. Equally Weighted Multiscale Elastic Network Model and Its Comparison with Traditional and Parameter-Free Models. J Chem Inf Model 2021; 61:921-937. [PMID: 33496590 DOI: 10.1021/acs.jcim.0c01178] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dynamical properties of proteins play an essential role in their function exertion. The elastic network model (ENM) is an effective and efficient tool in characterizing the intrinsic dynamical properties encoded in biomacromolecule structures. The Gaussian network model (GNM) and anisotropic network model (ANM) are the two often-used ENM models. Here, we introduce an equally weighted multiscale ENM (equally weighted mENM) based on the original mENM (denoted as mENM), in which fitting weights of Kirchhoff/Hessian matrixes in mENM are removed since they neglect the details of pairwise interactions. Then, we perform its comparison with the mENM, traditional ENM, and parameter-free ENM (pfENM) in reproducing dynamical properties for the six representative proteins whose molecular dynamics (MD) trajectories are available in http://mmb.pcb.ub.es/MoDEL/. In the results, for B-factor prediction, mENM performs best, while the equally weighted mENM performs also well, better than the traditional ENM and pfENM models. As to the dynamical cross-correlation map calculation, mENM performs worst, while the results produced from the equally weighted mENM and pfENM models are close to those from MD trajectories with the latter a little better than the former. Furthermore, encouragingly, the equally weighted mANM displays the best performance in capturing the functional motional modes, followed by pfANM and traditional ANM models, while the mANM fails in all the cases. This work is helpful for strengthening the understanding of the elastic network model and provides a valuable guide for researchers to utilize the model to explore protein dynamics.
Collapse
Affiliation(s)
- Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yang Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yanpeng Zhao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Shihao Wang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
15
|
Rosário-Ferreira N, Marques-Pereira C, Gouveia RP, Mourão J, Moreira IS. Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View. Methods Mol Biol 2021; 2315:3-28. [PMID: 34302667 DOI: 10.1007/978-1-0716-1468-6_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.
Collapse
Affiliation(s)
- Nícia Rosário-Ferreira
- Coimbra Chemistry Center, Department of Chemistry, University of Coimbra, Coimbra, Portugal.,Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal
| | - Catarina Marques-Pereira
- Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal.,PhD Programme in Experimental Biology and Biomedicine, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, Coimbra, Portugal
| | - Raquel P Gouveia
- Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal
| | - Joana Mourão
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal
| | - Irina S Moreira
- Department of Life Sciences, University of Coimbra, Coimbra, Portugal.
| |
Collapse
|
16
|
Scaramozzino D, Khade PM, Jernigan RL, Lacidogna G, Carpinteri A. Structural compliance: A new metric for protein flexibility. Proteins 2020; 88:1482-1492. [PMID: 32548853 DOI: 10.1002/prot.25968] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Revised: 04/29/2020] [Accepted: 06/06/2020] [Indexed: 12/17/2022]
Abstract
Proteins are the active players in performing essential molecular activities throughout biology, and their dynamics has been broadly demonstrated to relate to their mechanisms. The intrinsic fluctuations have often been used to represent their dynamics and then compared to the experimental B-factors. However, proteins do not move in a vacuum and their motions are modulated by solvent that can impose forces on the structure. In this paper, we introduce a new structural concept, which has been called the structural compliance, for the evaluation of the global and local deformability of the protein structure in response to intramolecular and solvent forces. Based on the application of pairwise pulling forces to a protein elastic network, this structural quantity has been computed and sometimes is even found to yield an improved correlation with the experimental B-factors, meaning that it may serve as a better metric for protein flexibility. The inverse of structural compliance, namely the structural stiffness, has also been defined, which shows a clear anticorrelation with the experimental data. Although the present applications are made to proteins, this approach can also be applied to other biomolecular structures such as RNA. This present study considers only elastic network models, but the approach could be applied further to conventional atomic molecular dynamics. Compliance is found to have a slightly better agreement with the experimental B-factors, perhaps reflecting its bias toward the effects of local perturbations, in contrast to mean square fluctuations. The code for calculating protein compliance and stiffness is freely accessible at https://jerniganlab.github.io/Software/PACKMAN/Tutorials/compliance.
Collapse
Affiliation(s)
- Domenico Scaramozzino
- Department of Structural, Geotechnical and Building Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, Torino, Italy
| | - Pranav M Khade
- Bioinformatics and Computational Biology Program, Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa, USA
| | - Robert L Jernigan
- Bioinformatics and Computational Biology Program, Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, Iowa, USA
| | - Giuseppe Lacidogna
- Department of Structural, Geotechnical and Building Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, Torino, Italy
| | - Alberto Carpinteri
- Department of Structural, Geotechnical and Building Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, Torino, Italy
| |
Collapse
|
17
|
Shape-preserving elastic solid models of macromolecules. PLoS Comput Biol 2020; 16:e1007855. [PMID: 32407309 PMCID: PMC7297265 DOI: 10.1371/journal.pcbi.1007855] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Revised: 06/09/2020] [Accepted: 04/07/2020] [Indexed: 02/07/2023] Open
Abstract
Mass-spring models have been a standard approach in molecular modeling for the last few decades, such as elastic network models (ENMs) that are widely used for normal mode analysis. In this work, we present a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as ENMs in producing the equilibrium dynamics and moreover, offers some significant new features that may greatly benefit the research community. ESM is different from ENM in that it treats macromolecules as elastic solids. Our particular version of ESM presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape and thus makes normal mode computations and visualization of extremely large complexes more manageable. Secondly, as a solid model, ESM’s close link to finite element analysis renders it ideally suited for studying mechanical responses of macromolecules under external force. Lastly, we show that ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy. The complete MATLAB code of αESM is provided. Mass-spring models have been a standard approach in classical molecular modeling where atoms are modeled as spheres with a mass and their interactions modeled as springs. The models have been extremely successful. Thinking ahead, however, as molecular systems of our interest grow more quickly in size or dimension than what our computation resources can keep up with, some adjustments in methodology are timely. This work presents a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as mass-spring models in producing the equilibrium dynamics and moreover, offers some unique features that make it suitable for much larger systems. ESM is different from ENMs in that it treats macromolecules as elastic solids. Our particular version of ESM model presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape. ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy.
Collapse
|
18
|
Tang QY, Kaneko K. Long-range correlation in protein dynamics: Confirmation by structural data and normal mode analysis. PLoS Comput Biol 2020; 16:e1007670. [PMID: 32053592 PMCID: PMC7043781 DOI: 10.1371/journal.pcbi.1007670] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 02/26/2020] [Accepted: 01/21/2020] [Indexed: 11/18/2022] Open
Abstract
Proteins in cellular environments are highly susceptible. Local perturbations to any residue can be sensed by other spatially distal residues in the protein molecule, showing long-range correlations in the native dynamics of proteins. The long-range correlations of proteins contribute to many biological processes such as allostery, catalysis, and transportation. Revealing the structural origin of such long-range correlations is of great significance in understanding the design principle of biologically functional proteins. In this work, based on a large set of globular proteins determined by X-ray crystallography, by conducting normal mode analysis with the elastic network models, we demonstrate that such long-range correlations are encoded in the native topology of the proteins. To understand how native topology defines the structure and the dynamics of the proteins, we conduct scaling analysis on the size dependence of the slowest vibration mode, average path length, and modularity. Our results quantitatively describe how native proteins balance between order and disorder, showing both dense packing and fractal topology. It is suggested that the balance between stability and flexibility acts as an evolutionary constraint for proteins at different sizes. Overall, our result not only gives a new perspective bridging the protein structure and its dynamics but also reveals a universal principle in the evolution of proteins at all different sizes. The long-range correlated fluctuations are closely related to many biological processes of the proteins, such as catalysis, ligand binding, biomolecular recognition, and transportation. In this paper, we elucidate the structural origin of the long-range correlation and describe how native contact topology defines the slow-mode dynamics of the native proteins. Our result suggests an evolutionary constraint for proteins at different sizes, which may shed light on solving many biophysical problems such as structure prediction, multi-scale molecular simulations, and the design of molecular machines. Moreover, in statistical physics, as the long-range correlations are notable signs of the critical point, unveiling the origin of such criticality can extend our understanding of the organizing principle of a large variety of complex systems.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
- * E-mail:
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
| |
Collapse
|
19
|
Analyzing Fluctuation Properties in Protein Elastic Networks with Sequence-Specific and Distance-Dependent Interactions. Biomolecules 2019; 9:biom9100549. [PMID: 31575003 PMCID: PMC6843209 DOI: 10.3390/biom9100549] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2019] [Revised: 09/20/2019] [Accepted: 09/24/2019] [Indexed: 01/26/2023] Open
Abstract
Simple protein elastic networks which neglect amino-acid information often yield reasonable predictions of conformational dynamics and are broadly used. Recently, model variants which incorporate sequence-specific and distance-dependent interactions of residue pairs have been constructed and demonstrated to improve agreement with experimental data. We have applied the new variants in a systematic study of protein fluctuation properties and compared their predictions with those of conventional anisotropic network models. We find that the quality of predictions is frequently linked to poor estimations in highly flexible protein regions. An analysis of a large set of protein structures shows that fluctuations of very weakly connected network residues are intrinsically prone to be significantly overestimated by all models. This problem persists in the new models and is not resolved by taking into account sequence information. The effect becomes even enhanced in the model variant which takes into account very soft long-ranged residue interactions. Beyond these shortcomings, we find that model predictions are largely insensitive to the integration of chemical information, at least regarding the fluctuation properties of individual residues. One can furthermore conclude that the inherent drawbacks may present a serious hindrance when improvement of elastic network models are attempted.
Collapse
|
20
|
Zhang PF, Su JG. Identification of key sites controlling protein functional motions by using elastic network model combined with internal coordinates. J Chem Phys 2019; 151:045101. [PMID: 31370540 DOI: 10.1063/1.5098542] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
The elastic network model (ENM) is an effective method to extract the intrinsic dynamical properties encoded in protein tertiary structures. We have proposed a new ENM-based analysis method to reveal the motion modes directly responsible for a specific protein function, in which an internal coordinate related to the specific function was introduced to construct the internal/Cartesian hybrid coordinate space. In the present work, the function-related internal coordinates combined with a linear perturbation method were applied to identify the key sites controlling specific protein functional motions. The change in the fluctuations of the internal coordinate in response to residue perturbation was calculated in the hybrid coordinate space by using the linear response theory. The residues with the large fluctuation changes were identified to be the key sites that allosterically control the specific protein function. Two proteins, i.e., human DNA polymerase β and the chaperonin from Methanococcus maripaludis, were investigated as case studies, in which several collective and local internal coordinates were applied to identify the functionally key residues of these two studied proteins. The calculation results are consistent with the experimental observations. It is found that different collective internal coordinates lead to similar results, where the predicted functionally key sites are located at similar positions in the protein structure. While for the local internal coordinates, the predicted key sites tend to be situated at the region near to the coordinate-involving residues. Our studies provide a starting point for further exploring other function-related internal coordinates for other interesting proteins.
Collapse
Affiliation(s)
- Peng Fei Zhang
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao 066004, China
| | - Ji Guo Su
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao 066004, China
| |
Collapse
|
21
|
Poma AB, Li MS, Theodorakis PE. Generalization of the elastic network model for the study of large conformational changes in biomolecules. Phys Chem Chem Phys 2019; 20:17020-17028. [PMID: 29904772 DOI: 10.1039/c8cp03086c] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Gō-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Gō-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for α-helix and β-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland.
| | | | | |
Collapse
|
22
|
Echave J. Beyond Stability Constraints: A Biophysical Model of Enzyme Evolution with Selection on Stability and Activity. Mol Biol Evol 2018; 36:613-620. [DOI: 10.1093/molbev/msy244] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Julian Echave
- Escuela de Ciencia y Tecnología, Universidad Nacional de San Martín (UNSAM), Buenos Aires, Argentina
| |
Collapse
|
23
|
Wang WB, Liang Y, Zhang J, Wu YD, Du JJ, Li QM, Zhu JZ, Su JG. Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model. Sci Rep 2018; 8:9487. [PMID: 29934573 PMCID: PMC6015066 DOI: 10.1038/s41598-018-27745-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2017] [Accepted: 06/08/2018] [Indexed: 11/28/2022] Open
Abstract
Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.
Collapse
Affiliation(s)
- Wei Bu Wang
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao, 066004, China
| | - Yu Liang
- Beijing Institute of Biological Products Co., Ltd, Beijing, 101111, China
| | - Jing Zhang
- Beijing Institute of Biological Products Co., Ltd, Beijing, 101111, China
| | - Yi Dong Wu
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao, 066004, China
| | - Jian Jun Du
- Beijing Key Lab of Digital Plant, Beijing Research Center for Information Technology in Agriculture, Beijing, 100097, China
| | - Qi Ming Li
- Beijing Institute of Biological Products Co., Ltd, Beijing, 101111, China
| | - Jian Zhuo Zhu
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao, 066004, China.
| | - Ji Guo Su
- Key Laboratory for Microstructural Material Physics of Hebei Province, College of Science, Yanshan University, Qinhuangdao, 066004, China.
| |
Collapse
|
24
|
Tiwari SP, Reuter N. Conservation of intrinsic dynamics in proteins — what have computational models taught us? Curr Opin Struct Biol 2018; 50:75-81. [DOI: 10.1016/j.sbi.2017.12.001] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2017] [Revised: 11/24/2017] [Accepted: 12/08/2017] [Indexed: 12/12/2022]
|
25
|
Dehouck Y, Bastolla U. The maximum penalty criterion for ridge regression: application to the calibration of the force constant in elastic network models. Integr Biol (Camb) 2018; 9:627-641. [PMID: 28555214 DOI: 10.1039/c7ib00079k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Tikhonov regularization, or ridge regression, is a popular technique to deal with collinearity in multivariate regression. We unveil a formal analogy between ridge regression and statistical mechanics, where the objective function is comparable to a free energy, and the ridge parameter plays the role of temperature. This analogy suggests two novel criteria for selecting a suitable ridge parameter: specific-heat (Cv) and maximum penalty (MP). We apply these fits to evaluate the relative contributions of rigid-body and internal fluctuations, which are typically highly collinear, to crystallographic B-factors. This issue is particularly important for computational models of protein dynamics, such as the elastic network model (ENM), since the amplitude of the predicted internal motion is commonly calibrated using B-factor data. After validation on simulated datasets, our results indicate that rigid-body motions account on average for more than 80% of the amplitude of B-factors. Furthermore, we evaluate the ability of different fits to reproduce the amplitudes of internal fluctuations in X-ray ensembles from the B-factors in the corresponding single X-ray structures. The new ridge criteria are shown to be markedly superior to the commonly used two-parameter fit that neglects rigid-body rotations and to the full fits regularized under generalized cross-validation. In conclusion, the proposed fits ensure a more robust calibration of the ENM force constant and should prove valuable in other applications.
Collapse
Affiliation(s)
- Yves Dehouck
- Machine Learning Group, Université Libre de Bruxelles (ULB), Boulevard du Triomphe CP 212, 1050 Brussels, Belgium.
| | | |
Collapse
|
26
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 327] [Impact Index Per Article: 54.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
27
|
Abstract
Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.
Collapse
Affiliation(s)
- Jeremy C Smith
- UT/ORNL Center for Molecular Biophysics, Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6309, USA; .,Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Pan Tan
- School of Physics and Astronomy and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Loukas Petridis
- UT/ORNL Center for Molecular Biophysics, Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6309, USA; .,Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Liang Hong
- School of Physics and Astronomy and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
28
|
Sankar K, Mishra SK, Jernigan RL. Comparisons of Protein Dynamics from Experimental Structure Ensembles, Molecular Dynamics Ensembles, and Coarse-Grained Elastic Network Models. J Phys Chem B 2018; 122:5409-5417. [DOI: 10.1021/acs.jpcb.7b11668] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Kannan Sankar
- Bioinformatics and Computational Biology Interdepartmental Graduate Program, Iowa State University, Ames, Iowa 50011-1178, United States
- Roy J. Carver Department of Biochemistry, Biophysics, and Molecular Biology, Iowa State University, Ames, Iowa 50011-1178, United States
| | - Sambit K. Mishra
- Bioinformatics and Computational Biology Interdepartmental Graduate Program, Iowa State University, Ames, Iowa 50011-1178, United States
- Roy J. Carver Department of Biochemistry, Biophysics, and Molecular Biology, Iowa State University, Ames, Iowa 50011-1178, United States
| | - Robert L. Jernigan
- Bioinformatics and Computational Biology Interdepartmental Graduate Program, Iowa State University, Ames, Iowa 50011-1178, United States
- Roy J. Carver Department of Biochemistry, Biophysics, and Molecular Biology, Iowa State University, Ames, Iowa 50011-1178, United States
| |
Collapse
|
29
|
Tirion MM, ben-Avraham D. PDB-NMA of a protein homodimer reproduces distinct experimental motility asymmetry. Phys Biol 2018; 15:026004. [DOI: 10.1088/1478-3975/aaa277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
30
|
Kaynak BT, Findik D, Doruker P. RESPEC Incorporates Residue Specificity and the Ligand Effect into the Elastic Network Model. J Phys Chem B 2017; 122:5347-5355. [DOI: 10.1021/acs.jpcb.7b10325] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Burak T. Kaynak
- Department of Physics, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Doga Findik
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, 34342, Bebek, Istanbul, Turkey
| | - Pemra Doruker
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, 34342, Bebek, Istanbul, Turkey
| |
Collapse
|
31
|
Patterns of coevolving amino acids unveil structural and dynamical domains. Proc Natl Acad Sci U S A 2017; 114:E10612-E10621. [PMID: 29183970 DOI: 10.1073/pnas.1712021114] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Patterns of interacting amino acids are so preserved within protein families that the sole analysis of evolutionary comutations can identify pairs of contacting residues. It is also known that evolution conserves functional dynamics, i.e., the concerted motion or displacement of large protein regions or domains. Is it, therefore, possible to use a pure sequence-based analysis to identify these dynamical domains? To address this question, we introduce here a general coevolutionary coupling analysis strategy and apply it to a curated sequence database of hundreds of protein families. For most families, the sequence-based method partitions amino acids into a few clusters. When viewed in the context of the native structure, these clusters have the signature characteristics of viable protein domains: They are spatially separated but individually compact. They have a direct functional bearing too, as shown for various reference cases. We conclude that even large-scale structural and functionally related properties can be recovered from inference methods applied to evolutionary-related sequences. The method introduced here is available as a software package and web server (spectrus.sissa.it/spectrus-evo_webserver).
Collapse
|
32
|
Putz I, Brock O. Elastic network model of learned maintained contacts to predict protein motion. PLoS One 2017; 12:e0183889. [PMID: 28854238 PMCID: PMC5576689 DOI: 10.1371/journal.pone.0183889] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2016] [Accepted: 08/14/2017] [Indexed: 12/21/2022] Open
Abstract
We present a novel elastic network model, lmcENM, to determine protein motion even for localized functional motions that involve substantial changes in the protein's contact topology. Existing elastic network models assume that the contact topology remains unchanged throughout the motion and are thus most appropriate to simulate highly collective function-related movements. lmcENM uses machine learning to differentiate breaking from maintained contacts. We show that lmcENM accurately captures functional transitions unexplained by the classical ENM and three reference ENM variants, while preserving the simplicity of classical ENM. We demonstrate the effectiveness of our approach on a large set of proteins covering different motion types. Our results suggest that accurately predicting a "deformation-invariant" contact topology offers a promising route to increase the general applicability of ENMs. We also find that to correctly predict this contact topology a combination of several features seems to be relevant which may vary slightly depending on the protein. Additionally, we present case studies of two biologically interesting systems, Ferric Citrate membrane transporter FecA and Arachidonate 15-Lipoxygenase.
Collapse
Affiliation(s)
- Ines Putz
- Robotics and Biology Laboratory, Department of Computer Science and Electrical Engineering, Technische Universität Berlin, Berlin, Berlin, Germany
| | - Oliver Brock
- Robotics and Biology Laboratory, Department of Computer Science and Electrical Engineering, Technische Universität Berlin, Berlin, Berlin, Germany
| |
Collapse
|
33
|
Gaussian network model can be enhanced by combining solvent accessibility in proteins. Sci Rep 2017; 7:7486. [PMID: 28790346 PMCID: PMC5548781 DOI: 10.1038/s41598-017-07677-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2016] [Accepted: 06/29/2017] [Indexed: 01/03/2023] Open
Abstract
Gaussian network model (GNM), regarded as the simplest and most representative coarse-grained model, has been widely adopted to analyze and reveal protein dynamics and functions. Designing a variation of the classical GNM, by defining a new Kirchhoff matrix, is the way to improve the residue flexibility modeling. We combined information arising from local relative solvent accessibility (RSA) between two residues into the Kirchhoff matrix of the parameter-free GNM. The undetermined parameters in the new Kirchhoff matrix were estimated by using particle swarm optimization. The usage of RSA was motivated by the fact that our previous work using RSA based linear regression model resulted out higher prediction quality of the residue flexibility when compared with the classical GNM and the parameter free GNM. Computational experiments, conducted based on one training dataset, two independent datasets and one additional small set derived by molecular dynamics simulations, demonstrated that the average correlation coefficients of the proposed RSA based parameter-free GNM, called RpfGNM, were significantly increased when compared with the parameter-free GNM. Our empirical results indicated that a variation of the classical GNMs by combining other protein structural properties is an attractive way to improve the quality of flexibility modeling.
Collapse
|
34
|
Mahajan S, Sanejouand YH. Jumping between protein conformers using normal modes. J Comput Chem 2017; 38:1622-1630. [PMID: 28470912 DOI: 10.1002/jcc.24803] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 02/03/2017] [Accepted: 03/19/2017] [Indexed: 12/27/2022]
Abstract
The relationship between the normal modes of a protein and its functional conformational change has been studied for decades. However, using this relationship in a predictive context remains a challenge. In this work, we demonstrate that, starting from a given protein conformer, it is possible to generate in a single step model conformers that are less than 1 Å (Cα -RMSD) from the conformer which is the known endpoint of the conformational change, particularly when the conformational change is collective in nature. Such accurate model conformers can be generated by following either the so-called robust or the 50 lowest-frequency modes obtained with various Elastic Network Models (ENMs). Interestingly, the quality of many of these models compares well with actual crystal structures, as assessed by the ROSETTA scoring function and PROCHECK. The most accurate and best quality conformers obtained in the present study were generated by using the 50 lowest-frequency modes of an all-atom ENM. However, with less than ten robust modes, which are identified without any prior knowledge of the nature of the conformational change, nearly 90% of the motion described by the 50 lowest-frequency modes of a protein can be captured. Such results strongly suggest that exploring the robust modes of ENMs may prove efficient for sampling the functionally relevant conformational repertoire of many proteins. © 2017 Wiley Periodicals, Inc.
Collapse
|
35
|
Abstract
We present a new conceptually simple and computationally efficient method for nonlinear normal-mode analysis called NOLB. It relies on the rotations-translations of blocks (RTB) theoretical basis developed by Y.-H. Sanejouand and colleagues [ Durand et al. Biopolymers 1994 , 34 , 759 - 771 . Tama et al. Proteins: Struct., Funct., Bioinf . 2000 , 41 , 1 - 7 ]. We demonstrate how to physically interpret the eigenvalues computed in the RTB basis in terms of angular and linear velocities applied to the rigid blocks and how to construct a nonlinear extrapolation of motion out of these velocities. The key observation of our method is that the angular velocity of a rigid block can be interpreted as the result of an implicit force, such that the motion of the rigid block can be considered as a pure rotation about a certain center. We demonstrate the motions produced with the NOLB method on three different molecular systems and show that some of the lowest frequency normal modes correspond to the biologically relevant motions. For example, NOLB detects the spiral sliding motion of the TALE protein, which is capable of rapid diffusion along its target DNA. Overall, our method produces better structures compared to the standard approach, especially at large deformation amplitudes, as we demonstrate by visual inspection, energy, and topology analyses and also by the MolProbity service validation. Finally, our method is scalable and can be applied to very large molecular systems, such as ribosomes. Standalone executables of the NOLB normal-mode analysis method are available at https://team.inria.fr/nano-d/software/nolb-normal-modes/ . A graphical user interface created for the SAMSON software platform will be made available at https://www.samson-connect.net .
Collapse
|
36
|
Koehl P, Poitevin F, Navaza R, Delarue M. The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems. J Chem Theory Comput 2017; 13:1424-1438. [PMID: 28170254 DOI: 10.1021/acs.jctc.6b01136] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Sciences and Genome Center, University of California, Davis , Davis, California 95616, United States
| | - Frédéric Poitevin
- Department of Structural Biology, Stanford University , Stanford, California 94305, United States.,Stanford PULSE Institute, SLAC National Accelerator Laboratory, Standford University , Menlo Park, California 94025, United States
| | - Rafael Navaza
- Platform of Crystallogenesis and Crystallography, CiTech, Institut Pasteur , 75015 Paris, France
| | - Marc Delarue
- Unité de Dynamique Structurale des Macromolécules, UMR 3528 du CNRS, Institut Pasteur , 75015 Paris, France
| |
Collapse
|
37
|
Hsieh YC, Poitevin F, Delarue M, Koehl P. Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids. Front Mol Biosci 2016; 3:85. [PMID: 28083537 PMCID: PMC5187361 DOI: 10.3389/fmolb.2016.00085] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2016] [Accepted: 12/12/2016] [Indexed: 11/13/2022] Open
Abstract
Key steps in the life cycle of a virus, such as the fusion event as the virus infects a host cell and its maturation process, relate to an intricate interplay between the structure and the dynamics of its constituent proteins, especially those that define its capsid, much akin to an envelope that protects its genomic material. We present a comprehensive, comparative analysis of such interplay for the capsids of two viruses from the flaviviridae family, Dengue (DENV) and Zika (ZIKV). We use for that purpose our own software suite, DD-NMA, which is based on normal mode analysis. We describe the elements of DD-NMA that are relevant to the analysis of large systems, such as virus capsids. In particular, we introduce our implementation of simplified elastic networks and justify their parametrization. Using DD-NMA, we illustrate the importance of packing interactions within the virus capsids on the dynamics of the E proteins of DENV and ZIKV. We identify differences between the computed atomic fluctuations of the E proteins in DENV and ZIKV and relate those differences to changes observed in their high resolution structures. We conclude with a discussion on additional analyses that are needed to fully characterize the dynamics of the two viruses.
Collapse
Affiliation(s)
- Yin-Chen Hsieh
- Department of Computer Science and Genome Center, University of California, Davis Davis, CA, USA
| | - Frédéric Poitevin
- Department of Structural Biology, Stanford UniversityStanford, CA, USA; SLAC National Accelerator Laboratory, Stanford PULSE InstituteMenlo Park, CA, USA
| | - Marc Delarue
- Unit of Structural Dynamics of Macromolecules, UMR 3528 du Centre National de la Recherche Scientifique, Institut Pasteur Paris, France
| | - Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis Davis, CA, USA
| |
Collapse
|
38
|
Calligari P, Gerolin M, Abergel D, Polimeno A. Decomposition of Proteins into Dynamic Units from Atomic Cross-Correlation Functions. J Chem Theory Comput 2016; 13:309-319. [DOI: 10.1021/acs.jctc.6b00702] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Paolo Calligari
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
| | - Marco Gerolin
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
- Département
de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules (LBM), 24 rue Lhomond, 75005 Paris, France
| | - Daniel Abergel
- Département
de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules (LBM), 24 rue Lhomond, 75005 Paris, France
| | - Antonino Polimeno
- Dipartimento
di Scienze Chimiche, Università di Padova, via Marzolo, 1, I-35131 Padova, Italy
| |
Collapse
|
39
|
Yao XQ, Skjærven L, Grant BJ. Rapid Characterization of Allosteric Networks with Ensemble Normal Mode Analysis. J Phys Chem B 2016; 120:8276-88. [PMID: 27056373 DOI: 10.1021/acs.jpcb.6b01991] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Allosteric regulation is a primary means of controlling protein function. By definition, allostery involves the propagation of structural dynamic changes between distal protein sites that yields a functional change. Gaining improved knowledge of these fundamental mechanisms is important for understanding many biomolecular processes and for guiding protein engineering and drug design efforts. In this work we compare and contrast a range of normal mode analysis (NMA) approaches together with network analysis for the prediction of structural dynamics and allosteric sites. Application to heterotrimeric G proteins, hemoglobin, and caspase 7 indicates that atomistic elastic network models provide improved predictions of experimental allosteric mutation sites. Results for G proteins also display an improved consistency with those derived from more computationally demanding MD simulations. Application of this approach across available experimental structures for a given protein family in a unified manner, that we refer to as ensemble NMA, yields the best overall predictive performance. We propose that this atomistic ensemble NMA approach represents an efficient and powerful tool for guiding the exploration of coupled motions and allosteric mechanisms in cases where multiple structures are available and where MD may prove prohibitively expensive.
Collapse
Affiliation(s)
- Xin-Qiu Yao
- Department of Computational Medicine and Bioinformatics, University of Michigan , 100 Washtenaw Avenue, 2017 Palmer Commons Building, Ann Arbor, Michigan 48109-2218, United States
| | - Lars Skjærven
- Department of Biomedicine, University of Bergen , Bergen N-5020, Norway
| | - Barry J Grant
- Department of Computational Medicine and Bioinformatics, University of Michigan , 100 Washtenaw Avenue, 2017 Palmer Commons Building, Ann Arbor, Michigan 48109-2218, United States
| |
Collapse
|
40
|
Hati S, Bhattacharyya S. Incorporating modeling and simulations in undergraduate biophysical chemistry course to promote understanding of structure-dynamics-function relationships in proteins. BIOCHEMISTRY AND MOLECULAR BIOLOGY EDUCATION : A BIMONTHLY PUBLICATION OF THE INTERNATIONAL UNION OF BIOCHEMISTRY AND MOLECULAR BIOLOGY 2016; 44:140-159. [PMID: 26801683 DOI: 10.1002/bmb.20942] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2015] [Accepted: 11/15/2015] [Indexed: 06/05/2023]
Abstract
A project-based biophysical chemistry laboratory course, which is offered to the biochemistry and molecular biology majors in their senior year, is described. In this course, the classroom study of the structure-function of biomolecules is integrated with the discovery-guided laboratory study of these molecules using computer modeling and simulations. In particular, modern computational tools are employed to elucidate the relationship between structure, dynamics, and function in proteins. Computer-based laboratory protocols that we introduced in three modules allow students to visualize the secondary, super-secondary, and tertiary structures of proteins, analyze non-covalent interactions in protein-ligand complexes, develop three-dimensional structural models (homology model) for new protein sequences and evaluate their structural qualities, and study proteins' intrinsic dynamics to understand their functions. In the fourth module, students are assigned to an authentic research problem, where they apply their laboratory skills (acquired in modules 1-3) to answer conceptual biophysical questions. Through this process, students gain in-depth understanding of protein dynamics-the missing link between structure and function. Additionally, the requirement of term papers sharpens students' writing and communication skills. Finally, these projects result in new findings that are communicated in peer-reviewed journals.
Collapse
Affiliation(s)
- Sanchita Hati
- Department of Chemistry, University Wisconsin, Eau Claire, Wisconsin
| | | |
Collapse
|
41
|
López-Blanco JR, Chacón P. New generation of elastic network models. Curr Opin Struct Biol 2015; 37:46-53. [PMID: 26716577 DOI: 10.1016/j.sbi.2015.11.013] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Revised: 11/23/2015] [Accepted: 11/26/2015] [Indexed: 12/16/2022]
Abstract
The intrinsic flexibility of proteins and nucleic acids can be grasped from remarkably simple mechanical models of particles connected by springs. In recent decades, Elastic Network Models (ENMs) combined with Normal Model Analysis widely confirmed their ability to predict biologically relevant motions of biomolecules and soon became a popular methodology to reveal large-scale dynamics in multiple structural biology scenarios. The simplicity, robustness, low computational cost, and relatively high accuracy are the reasons behind the success of ENMs. This review focuses on recent advances in the development and application of ENMs, paying particular attention to combinations with experimental data. Successful application scenarios include large macromolecular machines, structural refinement, docking, and evolutionary conservation.
Collapse
Affiliation(s)
- José Ramón López-Blanco
- Department of Biological Chemical Physics, Rocasolano Physical Chemistry Institute C.S.I.C., Serrano 119, 28006 Madrid, Spain
| | - Pablo Chacón
- Department of Biological Chemical Physics, Rocasolano Physical Chemistry Institute C.S.I.C., Serrano 119, 28006 Madrid, Spain.
| |
Collapse
|
42
|
Carvalho HF, Roque ACA, Iranzo O, Branco RJF. Comparison of the Internal Dynamics of Metalloproteases Provides New Insights on Their Function and Evolution. PLoS One 2015; 10:e0138118. [PMID: 26397984 PMCID: PMC4580569 DOI: 10.1371/journal.pone.0138118] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2015] [Accepted: 08/25/2015] [Indexed: 11/20/2022] Open
Abstract
Metalloproteases have evolved in a vast number of biological systems, being one of the most diverse types of proteases and presenting a wide range of folds and catalytic metal ions. Given the increasing understanding of protein internal dynamics and its role in enzyme function, we are interested in assessing how the structural heterogeneity of metalloproteases translates into their dynamics. Therefore, the dynamical profile of the clan MA type protein thermolysin, derived from an Elastic Network Model of protein structure, was evaluated against those obtained from a set of experimental structures and molecular dynamics simulation trajectories. A close correspondence was obtained between modes derived from the coarse-grained model and the subspace of functionally-relevant motions observed experimentally, the later being shown to be encoded in the internal dynamics of the protein. This prompted the use of dynamics-based comparison methods that employ such coarse-grained models in a representative set of clan members, allowing for its quantitative description in terms of structural and dynamical variability. Although members show structural similarity, they nonetheless present distinct dynamical profiles, with no apparent correlation between structural and dynamical relatedness. However, previously unnoticed dynamical similarity was found between the relevant members Carboxypeptidase Pfu, Leishmanolysin, and Botulinum Neurotoxin Type A, despite sharing no structural similarity. Inspection of the respective alignments shows that dynamical similarity has a functional basis, namely the need for maintaining proper intermolecular interactions with the respective substrates. These results suggest that distinct selective pressure mechanisms act on metalloproteases at structural and dynamical levels through the course of their evolution. This work shows how new insights on metalloprotease function and evolution can be assessed with comparison schemes that incorporate information on protein dynamics. The integration of these newly developed tools, if applied to other protein families, can lead to more accurate and descriptive protein classification systems.
Collapse
Affiliation(s)
- Henrique F. Carvalho
- UCIBIO-REQUIMTE, Department of Chemistry, Faculty of Science and Technology, Universidade NOVA de Lisboa, 2829-516 Caparica, Portugal
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, 2780–157 Oeiras, Portugal
| | - Ana C. A. Roque
- UCIBIO-REQUIMTE, Department of Chemistry, Faculty of Science and Technology, Universidade NOVA de Lisboa, 2829-516 Caparica, Portugal
| | - Olga Iranzo
- Aix Marseille Université, Centrale Marseille, CNRS, iSm2 UMR 7313, 13397, Marseille, France
| | - Ricardo J. F. Branco
- UCIBIO-REQUIMTE, Department of Chemistry, Faculty of Science and Technology, Universidade NOVA de Lisboa, 2829-516 Caparica, Portugal
| |
Collapse
|
43
|
Pinamonti G, Bottaro S, Micheletti C, Bussi G. Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments. Nucleic Acids Res 2015; 43:7260-9. [PMID: 26187990 PMCID: PMC4551938 DOI: 10.1093/nar/gkv708] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2015] [Accepted: 06/30/2015] [Indexed: 11/23/2022] Open
Abstract
Elastic network models (ENMs) are valuable and efficient tools for characterizing the collective internal dynamics of proteins based on the knowledge of their native structures. The increasing evidence that the biological functionality of RNAs is often linked to their innate internal motions poses the question of whether ENM approaches can be successfully extended to this class of biomolecules. This issue is tackled here by considering various families of elastic networks of increasing complexity applied to a representative set of RNAs. The fluctuations predicted by the alternative ENMs are stringently validated by comparison against extensive molecular dynamics simulations and SHAPE experiments. We find that simulations and experimental data are systematically best reproduced by either an all-atom or a three-beads-per-nucleotide representation (sugar-base-phosphate), with the latter arguably providing the best balance of accuracy and computational complexity.
Collapse
Affiliation(s)
- Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Cristian Micheletti
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
44
|
Tiberti M, Invernizzi G, Papaleo E. (Dis)similarity Index To Compare Correlated Motions in Molecular Simulations. J Chem Theory Comput 2015; 11:4404-14. [PMID: 26575932 DOI: 10.1021/acs.jctc.5b00512] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Molecular dynamics (MD) simulations are widely used to complement or guide experimental studies in the characterization of protein dynamics, thanks to improvements in force-field accuracy, along with in the software and hardware to sample the conformational landscape of proteins. Among the different applications of MD simulations, the study of correlated motions is largely employed for different purposes. Several metrics have been developed to describe correlated motions in the MD ensemble, such as methods based on Pearson Correlation or Mutual Information. Cross-correlation analysis of MD trajectories is indeed appealing not only to identify residues characterized by coupled fluctuations in protein structures but also since it can be used to extrapolate motions along directions in which major conformational changes should occur, for example on longer time scales than the ones that are actually simulated. Nevertheless, most of the MD studies employ average correlation maps and mostly in a qualitative way, even when different systems or different replicates of the same system are compared. The broad application of correlation metrics in the analysis of MD simulations, especially for comparative purposes, requires a step forward toward more quantitative and accurate comparisons. We thus here employed a simple but effective index, which is based on a normalized Frobenius norm of the differences between protein correlation maps, to compare correlated motions. We applied this index for a quantitative comparison of correlated motions from MD simulations of seven proteins of different size and fold. We also employed the index to assess the robustness of correlation description when multi-replicate MD simulations of a same system are used, and we compared our index to metrics for comparison of structural ensembles such as Root Mean Square Inner Product and the Bhattacharyya Coefficient.
Collapse
Affiliation(s)
- Matteo Tiberti
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Gaetano Invernizzi
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| | - Elena Papaleo
- Department of Biotechnology and Biosciences, University of Milano-Bicocca , Piazza della Scienza 2, 20126 Milan, Italy
| |
Collapse
|
45
|
Dorner ME, McMunn RD, Bartholow TG, Calhoon BE, Conlon MR, Dulli JM, Fehling SC, Fisher CR, Hodgson SW, Keenan SW, Kruger AN, Mabin JW, Mazula DL, Monte CA, Olthafer A, Sexton AE, Soderholm BR, Strom AM, Hati S. Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis. Protein Sci 2015; 24:1495-507. [PMID: 26130403 DOI: 10.1002/pro.2737] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2015] [Revised: 04/18/2015] [Accepted: 06/14/2015] [Indexed: 12/24/2022]
Abstract
Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.
Collapse
Affiliation(s)
- Mariah E Dorner
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Ryan D McMunn
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Thomas G Bartholow
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Brecken E Calhoon
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Michelle R Conlon
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Jessica M Dulli
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Samuel C Fehling
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Cody R Fisher
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Shane W Hodgson
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Shawn W Keenan
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Alyssa N Kruger
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Justin W Mabin
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Daniel L Mazula
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Christopher A Monte
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Augustus Olthafer
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Ashley E Sexton
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Beatrice R Soderholm
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Alexander M Strom
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| | - Sanchita Hati
- Department of Chemistry, University of Wisconsin-Eau Claire, Eau Claire, Wisconsin, 54702
| |
Collapse
|
46
|
Ponzoni L, Polles G, Carnevale V, Micheletti C. SPECTRUS: A Dimensionality Reduction Approach for Identifying Dynamical Domains in Protein Complexes from Limited Structural Datasets. Structure 2015; 23:1516-1525. [PMID: 26165596 DOI: 10.1016/j.str.2015.05.022] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2015] [Revised: 05/23/2015] [Accepted: 05/29/2015] [Indexed: 02/06/2023]
Abstract
Identifying dynamical, quasi-rigid domains in proteins provides a powerful means for characterizing functionally oriented structural changes via a parsimonious set of degrees of freedom. In fact, the relative displacements of few dynamical domains usually suffice to rationalize the mechanics underpinning biological functionality in proteins and can even be exploited for structure determination or refinement purposes. Here we present SPECTRUS, a general scheme that, by solely using amino acid distance fluctuations, can pinpoint the innate quasi-rigid domains of single proteins or large complexes in a robust way. Consistent domains are usually obtained by using either a pair of representative structures or thousands of conformers. The functional insights offered by the approach are illustrated for biomolecular systems of very different size and complexity such as kinases, ion channels, and viral capsids. The decomposition tool is available as a software package and web server at spectrus.sissa.it.
Collapse
Affiliation(s)
- Luca Ponzoni
- Molecular and Statistical Biophysics, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy.
| | - Guido Polles
- Molecular and Statistical Biophysics, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, College of Science and Technology, Temple University, SERC, 1925 North 12th Street, Philadelphia, PA 19122, USA
| | - Cristian Micheletti
- Molecular and Statistical Biophysics, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
47
|
Mahajan S, Sanejouand YH. On the relationship between low-frequency normal modes and the large-scale conformational changes of proteins. Arch Biochem Biophys 2015; 567:59-65. [PMID: 25562404 DOI: 10.1016/j.abb.2014.12.020] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2014] [Revised: 12/16/2014] [Accepted: 12/20/2014] [Indexed: 11/15/2022]
Abstract
Normal mode analysis is a computational technique that allows to study the dynamics of biological macromolecules. It was first applied to small protein cases, more than thirty years ago. The interest in this technique then raised when it was realized that it can provide insights about the large-scale conformational changes a protein can experience, for instance upon ligand binding. As it was also realized that studying highly simplified protein models can provide similar insights, meaning that this kind of analysis can be both quick and simple to handle, several applications were proposed, in the context of various structural biology techniques. This review focuses on these applications, as well as on how the functional relevance of the lowest-frequency modes of proteins was established.
Collapse
|
48
|
Tiwari SP, Fuglebakk E, Hollup SM, Skjærven L, Cragnolini T, Grindhaug SH, Tekle KM, Reuter N. WEBnm@ v2.0: Web server and services for comparing protein flexibility. BMC Bioinformatics 2014; 15:427. [PMID: 25547242 PMCID: PMC4339738 DOI: 10.1186/s12859-014-0427-6] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2014] [Accepted: 12/11/2014] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins. RESULTS We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma . CONCLUSION WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.
Collapse
Affiliation(s)
- Sandhya P Tiwari
- Department of Molecular Biology, University of Bergen, Bergen, Norway.
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Edvin Fuglebakk
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Siv M Hollup
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Lars Skjærven
- Department of Biomedicine, University of Bergen, Bergen, Norway.
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Tristan Cragnolini
- Department of Molecular Biology, University of Bergen, Bergen, Norway.
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
- Present address: University Chemical Laboratories, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| | - Svenn H Grindhaug
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Kidane M Tekle
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen, Bergen, Norway.
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| |
Collapse
|
49
|
Perica T, Kondo Y, Tiwari SP, McLaughlin SH, Kemplen KR, Zhang X, Steward A, Reuter N, Clarke J, Teichmann SA. Evolution of oligomeric state through allosteric pathways that mimic ligand binding. Science 2014; 346:1254346. [PMID: 25525255 PMCID: PMC4337988 DOI: 10.1126/science.1254346] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Evolution and design of protein complexes are almost always viewed through the lens of amino acid mutations at protein interfaces. We showed previously that residues not involved in the physical interaction between proteins make important contributions to oligomerization by acting indirectly or allosterically. In this work, we sought to investigate the mechanism by which allosteric mutations act, using the example of the PyrR family of pyrimidine operon attenuators. In this family, a perfectly sequence-conserved helix that forms a tetrameric interface is exposed as solvent-accessible surface in dimeric orthologs. This means that mutations must be acting from a distance to destabilize the interface. We identified 11 key mutations controlling oligomeric state, all distant from the interfaces and outside ligand-binding pockets. Finally, we show that the key mutations introduce conformational changes equivalent to the conformational shift between the free versus nucleotide-bound conformations of the proteins.
Collapse
Affiliation(s)
- Tina Perica
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK. Medical Research Council (MRC) Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge CB2 0QH, UK
| | - Yasushi Kondo
- Medical Research Council (MRC) Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge CB2 0QH, UK
| | - Sandhya P Tiwari
- Department of Molecular Biology, University of Bergen University of Bergen, P.O. Box 7803, N-5020 Bergen, Norway. Computational Biology Unit, Department of Informatics, University of Bergen, P.O. Box 7803, N-5020 Bergen, Norway
| | - Stephen H McLaughlin
- Medical Research Council (MRC) Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge CB2 0QH, UK
| | - Katherine R Kemplen
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Xiuwei Zhang
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Annette Steward
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen University of Bergen, P.O. Box 7803, N-5020 Bergen, Norway. Computational Biology Unit, Department of Informatics, University of Bergen, P.O. Box 7803, N-5020 Bergen, Norway
| | - Jane Clarke
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Sarah A Teichmann
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK.
| |
Collapse
|
50
|
Skjærven L, Yao XQ, Scarabelli G, Grant BJ. Integrating protein structural dynamics and evolutionary analysis with Bio3D. BMC Bioinformatics 2014; 15:399. [PMID: 25491031 PMCID: PMC4279791 DOI: 10.1186/s12859-014-0399-6] [Citation(s) in RCA: 249] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2014] [Accepted: 11/26/2014] [Indexed: 12/01/2022] Open
Abstract
Background Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution. Results Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case. Conclusions The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/. Electronic supplementary material The online version of this article (doi:10.1186/s12859-014-0399-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Lars Skjærven
- Department of Biomedicine, University of Bergen, Bergen, Norway. .,Structural and Computational Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany.
| | - Xin-Qiu Yao
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA.
| | - Guido Scarabelli
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA.
| | - Barry J Grant
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA.
| |
Collapse
|