1
|
Krieger JM, Sorzano COS, Carazo JM. Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing. Int J Mol Sci 2023; 24:14245. [PMID: 37762547 PMCID: PMC10532346 DOI: 10.3390/ijms241814245] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 09/10/2023] [Accepted: 09/15/2023] [Indexed: 09/29/2023] Open
Abstract
Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.
Collapse
Affiliation(s)
- James M. Krieger
- Biocomputing Unit, National Centre for Biotechnology (CNB CSIC), Campus Universidad Autónoma de Madrid, Darwin 3, Cantoblanco, 28049 Madrid, Spain
| | | | - Jose Maria Carazo
- Biocomputing Unit, National Centre for Biotechnology (CNB CSIC), Campus Universidad Autónoma de Madrid, Darwin 3, Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
2
|
Yu CC, Raj N, Chu JW. Statistical Learning of Protein Elastic Network from Positional Covariance Matrix. Comput Struct Biotechnol J 2023; 21:2524-2535. [PMID: 37095762 PMCID: PMC10121796 DOI: 10.1016/j.csbj.2023.03.033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 03/20/2023] [Accepted: 03/20/2023] [Indexed: 03/30/2023] Open
Abstract
Positional fluctuation and covariance during protein dynamics are key observables for understanding the molecular origin of biological functions. A frequently employed potential energy function for describing protein structural variation at the coarse-gained level is elastic network model (ENM). A long-standing issue in biomolecular simulation is thus the parametrization of ENM spring constants from the components of positional covariance matrix (PCM). Based on sensitivity analysis of PCM, the direct-coupling statistics of each spring, which is a specific combination of position fluctuation and covariance, is found to exhibit prominent signal of parameter dependence. This finding provides the basis for devising the objective function and the scheme of running through the effective one-dimensional optimization of every spring by self-consistent iteration. Formal derivation of the positional covariance statistical learning (PCSL) method also motivates the necessary data regularization for stable calculations. Robust convergence of PCSL is achieved in taking an all-atom molecular dynamics trajectory or an ensemble of homologous structures as input data. The PCSL framework can also be generalized with mixed objective functions to capture specific property such as the residue flexibility profile. Such physical chemistry-based statistical learning thus provides a useful platform for integrating the mechanical information encoded in various experimental or computational data.
Collapse
|
3
|
Mhashal AR, Yoluk O, Orellana L. Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations. Front Mol Biosci 2022; 9:890851. [PMID: 35836931 PMCID: PMC9275627 DOI: 10.3389/fmolb.2022.890851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Accepted: 05/05/2022] [Indexed: 11/13/2022] Open
Abstract
Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.
Collapse
|
4
|
Castillo F, Corbi-Verge C, Murciano-Calles J, Candel AM, Han Z, Iglesias-Bexiga M, Ruiz-Sanz J, Kim PM, Harty RN, Martinez JC, Luque I. Phage display identification of nanomolar ligands for human NEDD4-WW3: Energetic and dynamic implications for the development of broad-spectrum antivirals. Int J Biol Macromol 2022; 207:308-323. [PMID: 35257734 DOI: 10.1016/j.ijbiomac.2022.03.010] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Revised: 03/01/2022] [Accepted: 03/02/2022] [Indexed: 11/29/2022]
Abstract
The recognition of PPxY viral Late domains by the third WW domain of the human HECT-E3 ubiquitin ligase NEDD4 (NEDD4-WW3) is essential for the budding of many viruses. Blocking these interactions is a promising strategy to develop broad-spectrum antivirals. As all WW domains, NEDD4-WW3 is a challenging therapeutic target due to the low binding affinity of its natural interactions, its high conformational plasticity, and its complex thermodynamic behavior. In this work, we set out to investigate whether high affinity can be achieved for monovalent ligands binding to the isolated NEDD4-WW3 domain. We show that a competitive phage-display set-up allows for the identification of high-affinity peptides showing inhibitory activity of viral budding. A detailed biophysical study combining calorimetry, nuclear magnetic resonance, and molecular dynamic simulations reveals that the improvement in binding affinity does not arise from the establishment of new interactions with the domain, but is associated to conformational restrictions imposed by a novel C-terminal -LFP motif in the ligand, unprecedented in the PPxY interactome. These results, which highlight the complexity of WW domain interactions, provide valuable insight into the key elements for high binding affinity, of interest to guide virtual screening campaigns for the identification of novel therapeutics targeting NEDD4-WW3 interactions.
Collapse
Affiliation(s)
- Francisco Castillo
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Carles Corbi-Verge
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain; Donnelly Centre for Cellular and Biomolecular Research, Department of Molecular Genetics & Department of Computer Science, University of Toronto, Toronto, ON M5S 3E1, Canada
| | - Javier Murciano-Calles
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Adela M Candel
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Ziying Han
- Department of Pathobiology, School of Veterinary Medicine, University of Pennsylvania, 3800 Spruce St., Philadelphia, PA 19104, USA
| | - Manuel Iglesias-Bexiga
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Javier Ruiz-Sanz
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Philip M Kim
- Donnelly Centre for Cellular and Biomolecular Research, Department of Molecular Genetics & Department of Computer Science, University of Toronto, Toronto, ON M5S 3E1, Canada
| | - Ronald N Harty
- Department of Pathobiology, School of Veterinary Medicine, University of Pennsylvania, 3800 Spruce St., Philadelphia, PA 19104, USA
| | - Jose C Martinez
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain
| | - Irene Luque
- Department of Physical Chemistry, Institute of Biotechnology and Excelence Unit in Chemistry Applied to Biomedicine and Environment, School of Sciences, University of Granada, Campus Fuentenueva s/n 18071, Granada, Spain.
| |
Collapse
|
5
|
Protein Fluctuations in Response to Random External Forces. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12052344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Elastic network models (ENMs) have been widely used in the last decades to investigate protein motions and dynamics. There the intrinsic fluctuations based on the isolated structures are obtained from the normal modes of these elastic networks, and they generally show good agreement with the B-factors extracted from X-ray crystallographic experiments, which are commonly considered to be indicators of protein flexibility. In this paper, we propose a new approach to analyze protein fluctuations and flexibility, which has a more appropriate physical basis. It is based on the application of random forces to the protein ENM to simulate the effects of collisions of solvent on a protein structure. For this purpose, we consider both the Cα-atom coarse-grained anisotropic network model (ANM) and an elastic network augmented with points included for the crystallized waters. We apply random forces to these protein networks everywhere, as well as only on the protein surface alone. Despite the randomness of the directions of the applied perturbations, the computed average displacements of the protein network show a remarkably good agreement with the experimental B-factors. In particular, for our set of 919 protein structures, we find that the highest correlation with the B-factors is obtained when applying forces to the external surface of the water-augmented ANM (an overall gain of 3% in the Pearson’s coefficient for the entire dataset, with improvements up to 30% for individual proteins), rather than when evaluating the fluctuations obtained from the normal modes of a standard Cα-atom coarse-grained ANM. It follows that protein fluctuations should be considered not just as the intrinsic fluctuations of the internal dynamics, but also equally well as responses to external solvent forces, or as a combination of both.
Collapse
|
6
|
Matsuoka R, Fudim R, Jung S, Zhang C, Bazzone A, Chatzikyriakidou Y, Robinson CV, Nomura N, Iwata S, Landreh M, Orellana L, Beckstein O, Drew D. Structure, mechanism and lipid-mediated remodeling of the mammalian Na +/H + exchanger NHA2. Nat Struct Mol Biol 2022; 29:108-120. [PMID: 35173351 PMCID: PMC8850199 DOI: 10.1038/s41594-022-00738-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 12/14/2021] [Indexed: 11/09/2022]
Abstract
The Na+/H+ exchanger SLC9B2, also known as NHA2, correlates with the long-sought-after Na+/Li+ exchanger linked to the pathogenesis of diabetes mellitus and essential hypertension in humans. Despite the functional importance of NHA2, structural information and the molecular basis for its ion-exchange mechanism have been lacking. Here we report the cryo-EM structures of bison NHA2 in detergent and in nanodiscs, at 3.0 and 3.5 Å resolution, respectively. The bison NHA2 structure, together with solid-state membrane-based electrophysiology, establishes the molecular basis for electroneutral ion exchange. NHA2 consists of 14 transmembrane (TM) segments, rather than the 13 TMs previously observed in mammalian Na+/H+ exchangers (NHEs) and related bacterial antiporters. The additional N-terminal helix in NHA2 forms a unique homodimer interface with a large intracellular gap between the protomers, which closes in the presence of phosphoinositol lipids. We propose that the additional N-terminal helix has evolved as a lipid-mediated remodeling switch for the regulation of NHA2 activity.
Collapse
Affiliation(s)
- Rei Matsuoka
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | - Roman Fudim
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden.
| | - Sukkyeong Jung
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | - Chenou Zhang
- Center for Biological Physics, Department of Physics, Arizona State University, Tempe, AZ, USA
| | | | | | | | - Norimichi Nomura
- Graduate School of Medicine, Kyoto University, Konoe-cho, Yoshida, Sakyo-ku, Kyoto, Japan
| | - So Iwata
- Graduate School of Medicine, Kyoto University, Konoe-cho, Yoshida, Sakyo-ku, Kyoto, Japan
| | - Michael Landreh
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Laura Orellana
- Department of Oncology-Pathology, Karolinska Institutet, Stockholm, Sweden
| | - Oliver Beckstein
- Center for Biological Physics, Department of Physics, Arizona State University, Tempe, AZ, USA.
| | - David Drew
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden.
| |
Collapse
|
7
|
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
|
8
|
Almeida BC, Kaczmarek JA, Figueiredo PR, Prather KLJ, Carvalho ATP. Transcription factor allosteric regulation through substrate coordination to zinc. NAR Genom Bioinform 2021; 3:lqab033. [PMID: 33987533 PMCID: PMC8092373 DOI: 10.1093/nargab/lqab033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Revised: 03/30/2021] [Accepted: 04/08/2021] [Indexed: 11/14/2022] Open
Abstract
The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.
Collapse
Affiliation(s)
- Beatriz C Almeida
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| | - Jennifer A Kaczmarek
- MIT-Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Pedro R Figueiredo
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| | - Kristala L J Prather
- MIT-Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Alexandra T P Carvalho
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| |
Collapse
|
9
|
Structure and Dynamics of an Archeal Monoglyceride Lipase from Palaeococcus ferrophilus as Revealed by Crystallography and In Silico Analysis. Biomolecules 2021; 11:biom11040533. [PMID: 33916727 PMCID: PMC8065475 DOI: 10.3390/biom11040533] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 03/31/2021] [Accepted: 04/01/2021] [Indexed: 01/06/2023] Open
Abstract
The crystallographic analysis of a lipase from Palaeococcus ferrophilus (PFL) previously annotated as a lysophospholipase revealed high structural conservation with other monoglyceride lipases, in particular in the lid domain and substrate binding pockets. In agreement with this observation, PFL was shown to be active on various monoacylglycerols. Molecular Dynamics (MD) studies performed in the absence and in the presence of ligands further allowed characterization of the dynamics of this system and led to a systematic closure of the lid compared to the crystal structure. However, the presence of ligands in the acyl-binding pocket stabilizes intermediate conformations compared to the crystal and totally closed structures. Several lid-stabilizing or closure elements were highlighted, i.e., hydrogen bonds between Ser117 and Ile204 or Asn142 and its facing amino acid lid residues, as well as Phe123. Thus, based on this complementary crystallographic and MD approach, we suggest that the crystal structure reported herein represents an open conformation, at least partially, of the PFL, which is likely stabilized by the ligand, and it brings to light several key structural features prone to participate in the closure of the lid.
Collapse
|
10
|
Kosek D, Hickman AB, Ghirlando R, He S, Dyda F. Structures of ISCth4 transpososomes reveal the role of asymmetry in copy-out/paste-in DNA transposition. EMBO J 2021; 40:e105666. [PMID: 33006208 PMCID: PMC7780238 DOI: 10.15252/embj.2020105666] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Revised: 08/07/2020] [Accepted: 09/10/2020] [Indexed: 01/23/2023] Open
Abstract
Copy-out/paste-in transposition is a major bacterial DNA mobility pathway. It contributes significantly to the emergence of antibiotic resistance, often by upregulating expression of downstream genes upon integration. Unlike other transposition pathways, it requires both asymmetric and symmetric strand transfer steps. Here, we report the first structural study of a copy-out/paste-in transposase and demonstrate its ability to catalyze all pathway steps in vitro. X-ray structures of ISCth4 transposase, a member of the IS256 family of insertion sequences, bound to DNA substrates corresponding to three sequential steps in the reaction reveal an unusual asymmetric dimeric transpososome. During transposition, an array of N-terminal domains binds a single transposon end while the catalytic domain moves to accommodate the varying substrates. These conformational changes control the path of DNA flanking the transposon end and the generation of DNA-binding sites. Our results explain the asymmetric outcome of the initial strand transfer and show how DNA binding is modulated by the asymmetric transposase to allow the capture of a second transposon end and to integrate a circular intermediate.
Collapse
Affiliation(s)
- Dalibor Kosek
- Laboratory of Molecular BiologyNational Institute of Diabetes and Digestive and Kidney DiseasesNational Institutes of HealthBethesdaMDUSA
| | - Alison B Hickman
- Laboratory of Molecular BiologyNational Institute of Diabetes and Digestive and Kidney DiseasesNational Institutes of HealthBethesdaMDUSA
| | - Rodolfo Ghirlando
- Laboratory of Molecular BiologyNational Institute of Diabetes and Digestive and Kidney DiseasesNational Institutes of HealthBethesdaMDUSA
| | - Susu He
- Laboratory of Molecular BiologyNational Institute of Diabetes and Digestive and Kidney DiseasesNational Institutes of HealthBethesdaMDUSA
- Present address:
State Key Laboratory of Pharmaceutical BiotechnologyMedical School of Nanjing UniversityNanjingJiangsuChina
| | - Fred Dyda
- Laboratory of Molecular BiologyNational Institute of Diabetes and Digestive and Kidney DiseasesNational Institutes of HealthBethesdaMDUSA
| |
Collapse
|
11
|
Winklemann I, Matsuoka R, Meier PF, Shutin D, Zhang C, Orellana L, Sexton R, Landreh M, Robinson CV, Beckstein O, Drew D. Structure and elevator mechanism of the mammalian sodium/proton exchanger NHE9. EMBO J 2020; 39:e105908. [PMID: 33118634 PMCID: PMC7737618 DOI: 10.15252/embj.2020105908] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Revised: 09/09/2020] [Accepted: 09/10/2020] [Indexed: 12/21/2022] Open
Abstract
Na+ /H+ exchangers (NHEs) are ancient membrane-bound nanomachines that work to regulate intracellular pH, sodium levels and cell volume. NHE activities contribute to the control of the cell cycle, cell proliferation, cell migration and vesicle trafficking. NHE dysfunction has been linked to many diseases, and they are targets of pharmaceutical drugs. Despite their fundamental importance to cell homeostasis and human physiology, structural information for the mammalian NHE was lacking. Here, we report the cryogenic electron microscopy structure of NHE isoform 9 (SLC9A9) from Equus caballus at 3.2 Å resolution, an endosomal isoform highly expressed in the brain and associated with autism spectrum (ASD) and attention deficit hyperactivity (ADHD) disorders. Despite low sequence identity, the NHE9 architecture and ion-binding site are remarkably similar to distantly related bacterial Na+ /H+ antiporters with 13 transmembrane segments. Collectively, we reveal the conserved architecture of the NHE ion-binding site, their elevator-like structural transitions, the functional implications of autism disease mutations and the role of phosphoinositide lipids to promote homodimerization that, together, have important physiological ramifications.
Collapse
Affiliation(s)
- Iven Winklemann
- Department of Biochemistry and BiophysicsStockholm UniversityStockholmSweden
| | - Rei Matsuoka
- Department of Biochemistry and BiophysicsStockholm UniversityStockholmSweden
| | - Pascal F Meier
- Department of Biochemistry and BiophysicsStockholm UniversityStockholmSweden
| | - Denis Shutin
- Department of ChemistryUniversity of OxfordOxfordUK
| | - Chenou Zhang
- Department of PhysicsCenter for Biological PhysicsArizona State UniversityTempeAZUSA
| | - Laura Orellana
- Department of Biochemistry and BiophysicsStockholm UniversityStockholmSweden
| | - Ricky Sexton
- Department of PhysicsCenter for Biological PhysicsArizona State UniversityTempeAZUSA
| | - Michael Landreh
- Department of Microbiology, Tumor and Cell BiologyKarolinska InstituteStockholmSweden
| | | | - Oliver Beckstein
- Department of PhysicsCenter for Biological PhysicsArizona State UniversityTempeAZUSA
| | - David Drew
- Department of Biochemistry and BiophysicsStockholm UniversityStockholmSweden
| |
Collapse
|
12
|
Orellana L, Gustavsson J, Bergh C, Yoluk O, Lindahl E. eBDIMS server: protein transition pathways with ensemble analysis in 2D-motion spaces. Bioinformatics 2020; 35:3505-3507. [PMID: 30838394 PMCID: PMC6748756 DOI: 10.1093/bioinformatics/btz104] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2018] [Revised: 12/30/2018] [Accepted: 02/15/2019] [Indexed: 11/15/2022] Open
Abstract
SUMMARY Understanding how proteins transition between different conformers, and how conformers relate to each other in terms of structure and function, is not trivial. Here, we present an online tool for transition pathway generation between two protein conformations using Elastic Network Driven Brownian Dynamics Importance Sampling, a coarse-grained simulation algorithm, which spontaneously predicts transition intermediates trapped experimentally. In addition to path-generation, the server provides an interactive 2D-motion landscape graphical representation of the transitions or any additional conformers to explore their structural relationships. AVAILABILITY AND IMPLEMENTATION eBDIMS is available online: http://ebdims.biophysics.se/ or as standalone software: https://github.com/laura-orellana/eBDIMS, https://github.com/cabergh/eBDIMS. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Laura Orellana
- Department of Biochemistry and Biophysics, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Johan Gustavsson
- Department of Computational Science and Technology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Cathrine Bergh
- Department of Applied Physics, Swedish e-Science Research Center, KTH Royal Institute of Technology, Solna, Sweden
| | - Ozge Yoluk
- Department of Applied Physics, Swedish e-Science Research Center, KTH Royal Institute of Technology, Solna, Sweden
| | - Erik Lindahl
- Department of Biochemistry and Biophysics, Science for Life Laboratory, Stockholm University, Stockholm, Sweden.,Department of Applied Physics, Swedish e-Science Research Center, KTH Royal Institute of Technology, Solna, Sweden
| |
Collapse
|
13
|
Leherte L, Haufroid M, Mirgaux M, Wouters J. Investigation of bound and unbound phosphoserine phosphatase conformations through elastic network models and molecular dynamics simulations. J Biomol Struct Dyn 2020; 39:3958-3974. [PMID: 32448044 DOI: 10.1080/07391102.2020.1772883] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
The human phosphoserine phosphatase (hPSP) catalyses the last step in the biosynthesis of L-serine. It involves conformational changes of the enzyme lid once the substrate, phosphoserine (PSer), is bound in the active site. Here, Elastic Network Model (ENM) is applied to the crystal structure of hPSP to probe the transition between open and closed conformations of hPSP. Molecular Dynamics (MD) simulations are carried out on several PSer-hPSP systems to characterise the intermolecular interactions and their effect on the dynamics of the enzyme lid. Systems involving either Ca++ or Mg++ are considered. The first ENM normal mode shows that an open-closed transition can be explained from a simple description of the enzyme in terms of harmonic potentials. Principal Component Analyses applied to the MD trajectories also highlight a trend for a closing/opening motion. Different PSer orientations inside the enzyme cavity are identified, i.e. either the carboxylate, the phosphate group of PSer, or both, are oriented towards the cation. The interaction patterns are analysed in terms of hydrogen bonds, electrostatics, and bond critical points of the electron density distributions. The latter approach yields a global description of the bonding intermolecular interactions. The PSer orientation determines the content of the cation coordination shell and the mobility of the substrate, while Lys158 and Thr182, involved in the reaction mechanism, are always in interaction with the substrate. Closed enzyme conformations involve Met52-Gln204, Arg49-Glu29, and Arg50-Glu29 interactions. Met52, as well as Arg49 and Arg50, also stabilize PSer inside the cavity. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Laurence Leherte
- Laboratoire de Chimie Biologique Structurale, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), NAmur MEdicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Namur, Belgium
| | - Marie Haufroid
- Laboratoire de Chimie Biologique Structurale, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), NAmur MEdicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Namur, Belgium
| | - Manon Mirgaux
- Laboratoire de Chimie Biologique Structurale, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), NAmur MEdicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Namur, Belgium
| | - Johan Wouters
- Laboratoire de Chimie Biologique Structurale, Unité de Chimie Physique Théorique et Structurale, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), NAmur MEdicine & Drug Innovation Center (NAMEDIC), Namur Institute of Structured Matter (NISM), University of Namur, Namur, Belgium
| |
Collapse
|
14
|
Orellana L. Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier. Front Mol Biosci 2019; 6:117. [PMID: 31750315 PMCID: PMC6848229 DOI: 10.3389/fmolb.2019.00117] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2019] [Accepted: 10/14/2019] [Indexed: 12/16/2022] Open
Abstract
Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.
Collapse
Affiliation(s)
- Laura Orellana
- Institutionen för Biokemi och Biofysik, Stockholms Universitet, Stockholm, Sweden.,Science for Life Laboratory, Solna, Sweden
| |
Collapse
|
15
|
Aydınkal RM, Serçinoğlu O, Ozbek P. ProSNEx: a web-based application for exploration and analysis of protein structures using network formalism. Nucleic Acids Res 2019; 47:W471-W476. [PMID: 31114881 PMCID: PMC6602423 DOI: 10.1093/nar/gkz390] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 04/17/2019] [Accepted: 05/09/2019] [Indexed: 01/14/2023] Open
Abstract
ProSNEx (Protein Structure Network Explorer) is a web service for construction and analysis of Protein Structure Networks (PSNs) alongside amino acid flexibility, sequence conservation and annotation features. ProSNEx constructs a PSN by adding nodes to represent residues and edges between these nodes using user-specified interaction distance cutoffs for either carbon-alpha, carbon-beta or atom-pair contact networks. Different types of weighted networks can also be constructed by using either (i) the residue-residue interaction energies in the format returned by gRINN, resulting in a Protein Energy Network (PEN); (ii) the dynamical cross correlations from a coarse-grained Normal Mode Analysis (NMA) of the protein structure; (iii) interaction strength. Upon construction of the network, common network metrics (such as node centralities) as well as shortest paths between nodes and k-cliques are calculated. Moreover, additional features of each residue in the form of conservation scores and mutation/natural variant information are included in the analysis. By this way, tool offers an enhanced and direct comparison of network-based residue metrics with other types of biological information. ProSNEx is free and open to all users without login requirement at http://prosnex-tool.com.
Collapse
Affiliation(s)
- Rasim Murat Aydınkal
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Ali Nihat Gokyigit Foundation, Etiler, Istanbul 34340, Turkey
| | - Onur Serçinoğlu
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Department of Bioengineering, Faculty of Engineering, Recep Tayyip Erdoğan University, Rize 53100, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| |
Collapse
|
16
|
Lawless MJ, Pettersson JR, Rule GS, Lanni F, Saxena S. ESR Resolves the C Terminus Structure of the Ligand-free Human Glutathione S-Transferase A1-1. Biophys J 2019; 114:592-601. [PMID: 29414705 DOI: 10.1016/j.bpj.2017.12.016] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2017] [Revised: 12/11/2017] [Accepted: 12/18/2017] [Indexed: 01/05/2023] Open
Abstract
Nitroxide- and Cu2+-based electron spin resonance (ESR) are combined to provide insight into the conformational states of the functionally important α-helix of the human glutathione S-transferase A1. Distance measurements on various spin-labeled dimeric human glutathione S-transferase A1-1 all result in bimodal distance distributions, indicating that the C-terminus exists in two distinct conformations in solution, one of which closely matches that found in the crystal structure of the ligand-bound enzyme. These measurements permit the generation of a model of the unliganded conformation. Room temperature ESR indicates that the second conformation has high mobility, potentially enabling the enzyme's high degree of substrate promiscuity. This model is then validated using computational modeling and further Cu2+-based ESR distance measurements. Cu2+-based ESR also provides evidence that the secondary structure of the second conformation is of helical nature. Addition of S-hexyl glutathione results in a shift in relative populations, favoring the state that is similar to the previously known structure of the ligand-bound enzyme.
Collapse
Affiliation(s)
- Matthew J Lawless
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania
| | - John R Pettersson
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Gordon S Rule
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Frederick Lanni
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Sunil Saxena
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania.
| |
Collapse
|
17
|
Advances in coarse-grained modeling of macromolecular complexes. Curr Opin Struct Biol 2018; 52:119-126. [PMID: 30508766 DOI: 10.1016/j.sbi.2018.11.005] [Citation(s) in RCA: 83] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2018] [Revised: 11/05/2018] [Accepted: 11/17/2018] [Indexed: 01/12/2023]
Abstract
Recent progress in coarse-grained (CG) molecular modeling and simulation has facilitated an influx of computational studies on biological macromolecules and their complexes. Given the large separation of length-scales and time-scales that dictate macromolecular biophysics, CG modeling and simulation are well-suited to bridge the microscopic and mesoscopic or macroscopic details observed from all-atom molecular simulations and experiments, respectively. In this review, we first summarize recent innovations in the development of CG models, which broadly include structure-based, knowledge-based, and dynamics-based approaches. We then discuss recent applications of different classes of CG models to explore various macromolecular complexes. Finally, we conclude with an outlook for the future in this ever-growing field of biomolecular modeling.
Collapse
|
18
|
Sladek V, Tokiwa H, Shimano H, Shigeta Y. Protein Residue Networks from Energetic and Geometric Data: Are They Identical? J Chem Theory Comput 2018; 14:6623-6631. [PMID: 30500196 DOI: 10.1021/acs.jctc.8b00733] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein residue networks (PRN) from energetic and geometric data are probably not identical. PRNs constructed from ab initio pair interaction energies are analyzed for the first time and compared to PRN based on center of mass separation. We use modern, previously unused algorithms such as global and local efficiencies to quantitatively confirm that both types of PRNs do exhibit small-world character. The main novelty finding is that interaction energy-based PRNs preserve small-world character even when clustered. A node hierarchy independent of the cutoff energy used for the edge creation is characteristic for them. Efficiency centrality identifies hubs responsible for such behavior. The interaction energy-based PRNs seem to comply with the scale-free network model with respect to efficiency centrality distribution as opposed to distance based PRNs. Community detection is introduced into protein network research as an extension beyond cluster analysis to study tertiary and quaternary structures.
Collapse
Affiliation(s)
- Vladimir Sladek
- Institute of Chemistry - Centre for Glycomics , Dubravska cesta 9 , 84538 Bratislava , Slovakia.,Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan
| | - Hiroaki Tokiwa
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Chemistry , Rikkyo University , Nishi-Ikebukuro , Toshima, Tokyo 171-8501 , Japan
| | - Hitoshi Shimano
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Internal Medicine, Faculty of Medicine , University of Tsukuba , 1-1-1 Tennodai , Tsukuba, Ibaraki 305-8575 , Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences , University of Tsukuba , Tennodai 1-1-1 , Tsukuba, Ibaraki 305-8577 , Japan
| |
Collapse
|
19
|
Koehl P. Large Eigenvalue Problems in Coarse-Grained Dynamic Analyses of Supramolecular Systems. J Chem Theory Comput 2018; 14:3903-3919. [DOI: 10.1021/acs.jctc.8b00338] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Patrice Koehl
- Department of Computer Sciences and Genome Center, University of California, Davis, California 95616, United States
| |
Collapse
|
20
|
Cossins BP, Lawson ADG, Shi J. Computational Exploration of Conformational Transitions in Protein Drug Targets. Methods Mol Biol 2018; 1762:339-365. [PMID: 29594780 DOI: 10.1007/978-1-4939-7756-7_17] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2022]
Abstract
Protein drug targets vary from highly structured to completely disordered; either way dynamics governs function. Hence, understanding the dynamical aspects of how protein targets function can enable improved interventions with drug molecules. Computational approaches offer highly detailed structural models of protein dynamics which are becoming more predictive as model quality and sampling power improve. However, the most advanced and popular models still have errors owing to imperfect parameter sets and often cannot access longer timescales of many crucial biological processes. Experimental approaches offer more certainty but can struggle to detect and measure lightly populated conformations of target proteins and subtle allostery. An emerging solution is to integrate available experimental data into advanced molecular simulations. In the future, molecular simulation in combination with experimental data may be able to offer detailed models of important drug targets such that improved functional mechanisms or selectivity can be accessed.
Collapse
Affiliation(s)
- Benjamin P Cossins
- Computer-Aided Drug Design and Structural Biology, UCB Pharma, Slough, UK.
| | | | - Jiye Shi
- Computer-Aided Drug Design and Structural Biology, UCB Pharma, Slough, UK
| |
Collapse
|
21
|
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
|
22
|
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
|
23
|
Fogarty AC, Potestio R, Kremer K. A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site. Proteins 2016; 84:1902-1913. [PMID: 27699855 DOI: 10.1002/prot.25173] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2016] [Revised: 09/07/2016] [Accepted: 09/16/2016] [Indexed: 11/11/2022]
Abstract
In multi-resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902-1913. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Aoife C Fogarty
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| |
Collapse
|
24
|
Prediction and validation of protein intermediate states from structurally rich ensembles and coarse-grained simulations. Nat Commun 2016; 7:12575. [PMID: 27578633 PMCID: PMC5013691 DOI: 10.1038/ncomms12575] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Accepted: 07/13/2016] [Indexed: 12/28/2022] Open
Abstract
Protein conformational changes are at the heart of cell functions, from signalling to ion transport. However, the transient nature of the intermediates along transition pathways hampers their experimental detection, making the underlying mechanisms elusive. Here we retrieve dynamic information on the actual transition routes from principal component analysis (PCA) of structurally-rich ensembles and, in combination with coarse-grained simulations, explore the conformational landscapes of five well-studied proteins. Modelling them as elastic networks in a hybrid elastic-network Brownian dynamics simulation (eBDIMS), we generate trajectories connecting stable end-states that spontaneously sample the crystallographic motions, predicting the structures of known intermediates along the paths. We also show that the explored non-linear routes can delimit the lowest energy passages between end-states sampled by atomistic molecular dynamics. The integrative methodology presented here provides a powerful framework to extract and expand dynamic pathway information from the Protein Data Bank, as well as to validate sampling methods in general. Protein conformational changes are key to a wide range of cellular functions but remain difficult to access experimentally. Here the authors describe eBDIMS, a novel approach to predict intermediates observed in structural transition pathways from experimental ensembles.
Collapse
|
25
|
González ÀL, Teixidó J, Borrell JI, Estrada-Tejedor R. On the Applicability of Elastic Network Models for the Study of RNA CUG Trinucleotide Repeat Overexpansion. PLoS One 2016; 11:e0152049. [PMID: 27010216 PMCID: PMC4806922 DOI: 10.1371/journal.pone.0152049] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2015] [Accepted: 03/08/2016] [Indexed: 11/18/2022] Open
Abstract
Non-coding RNAs play a pivotal role in a number of diseases promoting an aberrant sequestration of nuclear RNA-binding proteins. In the particular case of myotonic dystrophy type 1 (DM1), a multisystemic autosomal dominant disease, the formation of large non-coding CUG repeats set up long-tract hairpins able to bind muscleblind-like proteins (MBNL), which trigger the deregulation of several splicing events such as cardiac troponin T (cTNT) and insulin receptor’s, among others. Evidence suggests that conformational changes in RNA are determinant for the recognition and binding of splicing proteins, molecular modeling simulations can attempt to shed light on the structural diversity of CUG repeats and to understand their pathogenic mechanisms. Molecular dynamics (MD) are widely used to obtain accurate results at atomistic level, despite being very time consuming, and they contrast with fast but simplified coarse-grained methods such as Elastic Network Model (ENM). In this paper, we assess the application of ENM (traditionally applied on proteins) for studying the conformational space of CUG repeats and compare it to conventional and accelerated MD conformational sampling. Overall, the results provided here reveal that ANM can provide useful insights into dynamic rCUG structures at a global level, and that their dynamics depend on both backbone and nucleobase fluctuations. On the other hand, ANM fail to describe local U-U dynamics of the rCUG system, which require more computationally expensive methods such as MD. Given that several limitations are inherent to both methods, we discuss here the usefulness of the current theoretical approaches for studying highly dynamic RNA systems such as CUG trinucleotide repeat overexpansions.
Collapse
Affiliation(s)
- Àlex L. González
- Grup d’Enginyeria Molecular (GEM), Institut Químic de Sarrià (IQS) – Universitat Ramon Llull (URL), Barcelona, Catalonia, 08017, Spain
| | - Jordi Teixidó
- Grup d’Enginyeria Molecular (GEM), Institut Químic de Sarrià (IQS) – Universitat Ramon Llull (URL), Barcelona, Catalonia, 08017, Spain
| | - José I. Borrell
- Grup d’Enginyeria Molecular (GEM), Institut Químic de Sarrià (IQS) – Universitat Ramon Llull (URL), Barcelona, Catalonia, 08017, Spain
| | - Roger Estrada-Tejedor
- Grup d’Enginyeria Molecular (GEM), Institut Químic de Sarrià (IQS) – Universitat Ramon Llull (URL), Barcelona, Catalonia, 08017, Spain
- * E-mail:
| |
Collapse
|
26
|
Sfriso P, Duran-Frigola M, Mosca R, Emperador A, Aloy P, Orozco M. Residues Coevolution Guides the Systematic Identification of Alternative Functional Conformations in Proteins. Structure 2016; 24:116-126. [DOI: 10.1016/j.str.2015.10.025] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2015] [Revised: 10/13/2015] [Accepted: 10/17/2015] [Indexed: 12/12/2022]
|
27
|
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
|
28
|
Hospital A, Goñi JR, Orozco M, Gelpí JL. Molecular dynamics simulations: advances and applications. Adv Appl Bioinform Chem 2015; 8:37-47. [PMID: 26604800 PMCID: PMC4655909 DOI: 10.2147/aabc.s70333] [Citation(s) in RCA: 233] [Impact Index Per Article: 25.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.
Collapse
Affiliation(s)
- Adam Hospital
- Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, University of Barcelona, Barcelona, Spain
| | - Josep Ramon Goñi
- Joint BSC-IRB Research Program in Computational Biology, University of Barcelona, Barcelona, Spain ; Barcelona Supercomputing Center, University of Barcelona, Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine, The Barcelona Institute of Science and Technology, University of Barcelona, Barcelona, Spain ; Joint BSC-IRB Research Program in Computational Biology, University of Barcelona, Barcelona, Spain ; Barcelona Supercomputing Center, University of Barcelona, Barcelona, Spain ; Department of Biochemistry and Molecular Biology, University of Barcelona, Barcelona, Spain
| | - Josep L Gelpí
- Joint BSC-IRB Research Program in Computational Biology, University of Barcelona, Barcelona, Spain ; Barcelona Supercomputing Center, University of Barcelona, Barcelona, Spain ; Department of Biochemistry and Molecular Biology, University of Barcelona, Barcelona, Spain
| |
Collapse
|
29
|
Ahalawat N, Murarka RK. Conformational changes and allosteric communications in human serum albumin due to ligand binding. J Biomol Struct Dyn 2015; 33:2192-204. [DOI: 10.1080/07391102.2014.996609] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Navjeet Ahalawat
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal , Bhauri, Bhopal, Madhya Pradesh 462066, India
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal , Bhauri, Bhopal, Madhya Pradesh 462066, India
| |
Collapse
|
30
|
Heterogeneous elastic network model improves description of slow motions of proteins in solution. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2014.11.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
31
|
Fenwick RB, Orellana L, Esteban-Martín S, Orozco M, Salvatella X. Correlated motions are a fundamental property of β-sheets. Nat Commun 2014; 5:4070. [PMID: 24915882 DOI: 10.1038/ncomms5070] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2014] [Accepted: 05/08/2014] [Indexed: 01/19/2023] Open
Abstract
Correlated motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. The mechanisms that underlie these processes remain largely unknown due mainly to limitations in their direct detection. Here, based on a detailed analysis of protein structures deposited in the protein data bank, as well as on state-of-the art molecular simulations, we provide general evidence for the transfer of structural information by correlated backbone motions, mediated by hydrogen bonds, across β-sheets. We also show that the observed local and long-range correlated motions are mediated by the collective motions of β-sheets and investigate their role in large-scale conformational changes. Correlated motions represent a fundamental property of β-sheets that contributes to protein function.
Collapse
Affiliation(s)
- R Bryn Fenwick
- 1] Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), Baldiri Reixac 10, 08028 Barcelona, Spain [2]
| | - Laura Orellana
- 1] Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), Baldiri Reixac 10, 08028 Barcelona, Spain [2]
| | - Santi Esteban-Martín
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), Baldiri Reixac 10, 08028 Barcelona, Spain
| | - Modesto Orozco
- 1] Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), Baldiri Reixac 10, 08028 Barcelona, Spain [2] Departament de Bioquímica i Biologia Molecular, Facultat de Biologia, Universitat de Barcelona, Avinguda Diagonal 645, 08028 Barcelona, Spain
| | - Xavier Salvatella
- 1] Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine (IRB Barcelona), Baldiri Reixac 10, 08028 Barcelona, Spain [2] Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
| |
Collapse
|
32
|
López-Blanco JR, Aliaga JI, Quintana-Ortí ES, Chacón P. iMODS: internal coordinates normal mode analysis server. Nucleic Acids Res 2014; 42:W271-6. [PMID: 24771341 PMCID: PMC4086069 DOI: 10.1093/nar/gku339] [Citation(s) in RCA: 369] [Impact Index Per Article: 36.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Normal mode analysis (NMA) in internal (dihedral) coordinates naturally reproduces the collective functional motions of biological macromolecules. iMODS facilitates the exploration of such modes and generates feasible transition pathways between two homologous structures, even with large macromolecules. The distinctive internal coordinate formulation improves the efficiency of NMA and extends its applicability while implicitly maintaining stereochemistry. Vibrational analysis, motion animations and morphing trajectories can be easily carried out at different resolution scales almost interactively. The server is versatile; non-specialists can rapidly characterize potential conformational changes, whereas advanced users can customize the model resolution with multiple coarse-grained atomic representations and elastic network potentials. iMODS supports advanced visualization capabilities for illustrating collective motions, including an improved affine-model-based arrow representation of domain dynamics. The generated all-heavy-atoms conformations can be used to introduce flexibility for more advanced modeling or sampling strategies. The server is free and open to all users with no login requirement at http://imods.chaconlab.org.
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
| | - José I Aliaga
- Department of Computer Science and Engineering, University Jaume I, 12071 Castellón, Spain
| | | | - Pablo Chacón
- Department of Biological Chemical Physics, Rocasolano Physical Chemistry Institute C.S.I.C., Serrano 119, 28006 Madrid, Spain
| |
Collapse
|
33
|
Bastolla U. Computing protein dynamics from protein structure with elastic network models. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1186] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Ugo Bastolla
- Centro de Biologa Molecular Severo Ochoa (CSIC‐UAM)Universidad Autónoma de MadridMadridSpain
| |
Collapse
|
34
|
Abstract
Proteins are fascinating supramolecular structures, which are able to recognize ligands transforming binding information into chemical signals. They can transfer information across the cell, can catalyse complex chemical reactions, and are able to transform energy into work with much more efficiency than any human engine. The unique abilities of proteins are tightly coupled with their dynamic properties, which are coded in a complex way in the sequence and carefully refined by evolution. Despite its importance, our experimental knowledge of protein dynamics is still rather limited, and mostly derived from theoretical calculations. I will review here, in a systematic way, the current state-of-the-art theoretical approaches to the study of protein dynamics, emphasizing the most recent advances, examples of use and the expected lines of development in the near future.
Collapse
Affiliation(s)
- Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), Baldiri i Reixac 8, Barcelona 08028, Spain.
| |
Collapse
|
35
|
Identifying essential pairwise interactions in elastic network model using the alpha shape theory. J Comput Chem 2014; 35:1111-21. [DOI: 10.1002/jcc.23587] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2013] [Revised: 01/24/2014] [Accepted: 02/26/2014] [Indexed: 11/07/2022]
|
36
|
Setny P, Zacharias M. Elastic Network Models of Nucleic Acids Flexibility. J Chem Theory Comput 2013; 9:5460-70. [PMID: 26592282 DOI: 10.1021/ct400814n] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Elastic network models (ENMs) are a useful tool for describing large scale motions in protein systems. While they are well validated in the context of proteins, relatively little is known about their applicability to nucleic acids, whose different architecture does not necessarily warrant comparable performance. In this study we thoroughly evaluate and optimize the efficiency of popular ENMs for capturing RNA and DNA flexibility. We also introduce two alternative models in which the strength of elastic connections at a coarse-grained level is governed by distance distribution at atomic resolution. For each of the considered ENMs we report the optimal length of spring connections as well as the scaling of elastic force constants that provides the best agreement of vibrational frequencies with normal modes based on atomic force field. In order to determine the absolute values of force constants we introduce a novel method based on the overlap of pseudoinverse of Hessian matrices.
Collapse
Affiliation(s)
- Piotr Setny
- Centre for New Technologies, University of Warsaw , 00-927 Warsaw, Poland
| | - Martin Zacharias
- Physics Department T38, Technical University Munich , 85748 Garching, Germany
| |
Collapse
|
37
|
Xia F, Tong D, Lu L. Robust Heterogeneous Anisotropic Elastic Network Model Precisely Reproduces the Experimental B-factors of Biomolecules. J Chem Theory Comput 2013; 9:3704-14. [PMID: 26584122 DOI: 10.1021/ct4002575] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A computational method called the progressive fluctuation matching (PFM) is developed for constructing robust heterogeneous anisotropic network models (HANMs) for biomolecular systems. An HANM derived through the PFM approach consists of harmonic springs with realistic positive force constants, and yields the calculated B-factors that are basically identical to the experimental ones. For the four tested protein systems including crambin, trypsin inhibitor, HIV-1 protease, and lysozyme, the root-mean-square deviations between the experimental and the computed B-factors are only 0.060, 0.095, 0.247, and 0.049 Å(2), respectively, and the correlation coefficients are 0.99 for all. By comparing the HANM/ANM normal modes to their counterparts derived from both an atomistic force field and an NMR structure ensemble, it is found that HANM may provide more accurate results on protein dynamics.
Collapse
Affiliation(s)
- Fei Xia
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| | - Dudu Tong
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| | - Lanyuan Lu
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| |
Collapse
|
38
|
Sfriso P, Hospital A, Emperador A, Orozco M. Exploration of conformational transition pathways from coarse-grained simulations. ACTA ACUST UNITED AC 2013; 29:1980-6. [PMID: 23740746 DOI: 10.1093/bioinformatics/btt324] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
MOTIVATION A new algorithm to trace conformational transitions in proteins is presented. The method uses discrete molecular dynamics as engine to sample protein conformational space. A multiple minima Go-like potential energy function is used in combination with several enhancing sampling strategies, such as metadynamics, Maxwell Demon molecular dynamics and essential dynamics. The method, which shows an unprecedented computational efficiency, is able to trace a wide range of known experimental transitions. Contrary to simpler methods our strategy does not introduce distortions in the chemical structure of the protein and is able to reproduce well complex non-linear conformational transitions. The method, called GOdMD, can easily introduce additional restraints to the transition (presence of ligand, known intermediate, known maintained contacts, …) and is freely distributed to the community through the Spanish National Bioinformatics Institute (http://mmb.irbbarcelona.org/GOdMD). AVAILABILITY Freely available on the web at http://mmb.irbbarcelona.org/GOdMD.
Collapse
Affiliation(s)
- Pedro Sfriso
- Institute for Research in Biomedicine (IRB Barcelona), Joint IRB-BSC Program in Computational Biology, Baldiri Reixac 10, Barcelona, Spain
| | | | | | | |
Collapse
|
39
|
|
40
|
PIM: phase integrated method for normal mode analysis of biomolecules in a crystalline environment. J Mol Biol 2013; 425:1082-98. [PMID: 23333742 DOI: 10.1016/j.jmb.2012.12.026] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2012] [Accepted: 12/31/2012] [Indexed: 11/21/2022]
Abstract
In this study, a normal mode analysis, named phase integrated method (PIM), is developed for computing modes of biomolecules in a crystalline environment. PIM can calculate low-frequency modes on one or a few asymmetric units (AUs) and generate exact modes of a whole unit cell according to space group symmetry, while the translational symmetry between unit cells is maintained via the periodic boundary condition. Therefore, the method can dramatically reduce computational cost in mode calculation in the presence of crystal symmetry. PIM also has an option to map modes onto a single AU to form an orthonormalized mode set, which can be directly applied to normal-mode-based thermal parameter refinement in X-ray crystallography. The performance of PIM was tested on all 65 space groups available in protein crystals (one protein for each space group) and on another set of 83 ultra-high-resolution X-ray structures. The results showed that considering space group symmetry in mode calculation is crucial for accurately describing vibrational motion in a crystalline environment. Moreover, the optimal inter-AU packing stiffness was found to be about 60% of that of intra-AU interactions (non-bonded interaction only).
Collapse
|
41
|
Jeschke G. Optimization of Algorithms for Modeling Protein Structural Transitions from Sparse Long-Range Spin-Label Distance Constraints. ACTA ACUST UNITED AC 2012. [DOI: 10.1524/zpch.2012.0289] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Abstract
Function-related structural transitions of proteins are often large-scale conformational changes that are related to essential dynamics of a protein, which in turn proceeds along a small number of slow normal modes. Hence, it should be possible to characterize such transitions by an equally small number of distance constraints. These constraints should contain the information how the backbone coordinates move along the active normal modes. Such an approach based on residue-level elastic network models for the protein backbone is optimized with respect to the fit algorithm, constraint selection, and parametrization of the elastic network model. A stable fitting algorithm can be based on energy equipartitioning among the modes in the active space. This stabilization allows for extending active space dimension beyond the number of available constraints, which improves fit quality for transitions with a normal mode spectrum of only slowly increasing frequencies. In constraint selection, discrimination between the active normal modes appears to be more important than achieving large distance changes between initial and final structure. Parametrization of the network model has only a small influence on fit quality, as long as scaling of force constants with the inverse sixth power of the distance between network nodes is maintained. Elastic network models with a uniform force constant below a cutoff distance perform significantly worse. With 50 distance constraints, the optimized approach covers more than 50% of the structural change for 44% of all test cases, between 25 and 50% for 22% of the cases, and it fails for 33%.
Collapse
|
42
|
Micheletti C. Comparing proteins by their internal dynamics: exploring structure-function relationships beyond static structural alignments. Phys Life Rev 2012. [PMID: 23199577 DOI: 10.1016/j.plrev.2012.10.009] [Citation(s) in RCA: 70] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
The growing interest for comparing protein internal dynamics owes much to the realisation that protein function can be accompanied or assisted by structural fluctuations and conformational changes. Analogously to the case of functional structural elements, those aspects of protein flexibility and dynamics that are functionally oriented should be subject to evolutionary conservation. Accordingly, dynamics-based protein comparisons or alignments could be used to detect protein relationships that are more elusive to sequence and structural alignments. Here we provide an account of the progress that has been made in recent years towards developing and applying general methods for comparing proteins in terms of their internal dynamics and advance the understanding of the structure-function relationship.
Collapse
Affiliation(s)
- Cristian Micheletti
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste, Italy.
| |
Collapse
|
43
|
Gniewek P, Kolinski A, Jernigan RL, Kloczkowski A. Elastic network normal modes provide a basis for protein structure refinement. J Chem Phys 2012; 136:195101. [PMID: 22612113 DOI: 10.1063/1.4710986] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.
Collapse
Affiliation(s)
- Pawel Gniewek
- Laboratory of Theory of Biopolymers, Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | | | | | | |
Collapse
|
44
|
Sfriso P, Emperador A, Orellana L, Hospital A, Gelpí JL, Orozco M. Finding Conformational Transition Pathways from Discrete Molecular Dynamics Simulations. J Chem Theory Comput 2012; 8:4707-18. [DOI: 10.1021/ct300494q] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Pedro Sfriso
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
| | - Agusti Emperador
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
| | - Laura Orellana
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
| | - Adam Hospital
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
- Structural Bioinformatics Node,
Instituto Nacional De Bioinformática, Institute of Research
in Biomedicine, Josep Samitier 1-5, Barcelona, 08028, Spain
| | - Josep Lluis Gelpí
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
- Computational Bioinformatics Node,
Instituto Nacional De Bioinformática, Barcelona Supercomputing
Center, Jordi Girona 29, Barcelona, 08034, Spain
- Departament de Bioquímica,
Facultat de Biologia, Universtitat de Barcelona, Avgda Diagonal 647,
Barcelona, 08028, Spain
| | - Modesto Orozco
- Joint IRB-BSC Program in Computational
Biology, Institute of Research in Biomedicine, Josep Samitier 1-5,
Barcelona, 08028, Spain
- Structural Bioinformatics Node,
Instituto Nacional De Bioinformática, Institute of Research
in Biomedicine, Josep Samitier 1-5, Barcelona, 08028, Spain
- Departament de Bioquímica,
Facultat de Biologia, Universtitat de Barcelona, Avgda Diagonal 647,
Barcelona, 08028, Spain
| |
Collapse
|
45
|
Leioatts N, Romo TD, Grossfield A. Elastic Network Models are Robust to Variations in Formalism. J Chem Theory Comput 2012; 8:2424-2434. [PMID: 22924033 DOI: 10.1021/ct3000316] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Understanding the functions of biomolecules requires insight not only from structures, but from dynamics as well. Often, the most interesting processes occur on time scales too slow for exploration by conventional molecular dynamics (MD) simulations. For this reason, alternative computational methods such as elastic network models (ENMs) have become increasingly popular. These simple, coarse-grained models represent molecules as beads connected by harmonic springs; the system's motions are solved analytically by normal mode analysis. In the past few years, many different formalisms for performing ENM calculations have emerged, and several have been optimized using all-atom MD simulations. In contrast to other studies, we have compared the various formalisms in a systematic, quantitative way. In this study, we optimize many ENM functional forms using a uniform dataset containing only long (> 1 μs) all-atom MD simulations. Our results show that all models once optimized produce spring constants for immediate neighboring residues that are orders of magnitude stiffer than more distal contacts. In addition, the statistical significance of ENM performance varied with model resolution. We also show that fitting long trajectories does not improve ENM performance due to a problem inherent in all network models tested: they underestimate the relative importance of the most concerted motions. Finally, we characterize ENMs' resilience by tessellating the parameter space to show that broad ranges of parameters produce similar quality predictions. Taken together our data reveals that choice of spring function and parameters are not vital to performance of a network model and that simple parameters can by derived "by hand" when no data is available for fitting, thus illustrating the robustness of these models.
Collapse
Affiliation(s)
- Nicholas Leioatts
- Department of Biochemistry & Biophysics, University of Rochester Medical Center, Rochester, NY 14642, USA
| | | | | |
Collapse
|
46
|
Jeschke G. Characterization of Protein Conformational Changes with Sparse Spin-Label Distance Constraints. J Chem Theory Comput 2012; 8:3854-63. [PMID: 26593026 DOI: 10.1021/ct300113z] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The combination of site-directed spin labeling with pulse EPR distance measurements can provide a moderate number of distance constraints on the nanometer length scale for proteins in different states. By adapting an existing algorithm (Zheng, W.; Brooks, B. R. Biophys. J. 2006, 90, 4327) to the problem, we address the question to what extent conformational change can be characterized when the protein structure is known for one of the states, whereas only a sparse set of distance constraints between spin labels is available for the other state. We find that the type and general direction of the conformational change can be recognized, while the amplitude may be uncertain.
Collapse
Affiliation(s)
- G Jeschke
- Lab. Phys. Chem., ETH Zürich, Wolfgang-Pauli-Strasse 10, CH-8093 Zürich, Switzerland
| |
Collapse
|
47
|
Perez A, Yang Z, Bahar I, Dill KA, MacCallum JL. FlexE: Using elastic network models to compare models of protein structure. J Chem Theory Comput 2012; 8:3985-3991. [PMID: 25530735 DOI: 10.1021/ct300148f] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
It is often valuable to compare protein structures to determine how similar they are. Structure comparison methods such as RMSD and GDT-TS are based solely on fixed geometry and do not take into account the intrinsic flexibility or energy landscape of the protein. We propose a method, which we call FlexE, that is based on a simple elastic network model and uses the deformation energy as measure of the similarity between two structures. FlexE can distinguish biologically relevant conformational changes from random changes, while existing geometry-based methods cannot. Additionally, FlexE incorporates the concept of thermal energy, which provides a rational way to determine when two models are "the same". FlexE provides a unique measure of the similarity between protein structures that is complementary to existing methods.
Collapse
Affiliation(s)
- Alberto Perez
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794-5252
| | - Zheng Yang
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, 3064 BST3, 3501 Fifth Ave, Pittsburgh, PA 15213
| | - Ivet Bahar
- Department of Computational and Systems Biology, and Clinical & Translational Science Institute, School of Medicine, University of Pittsburgh, 3064 BST3, 3501 Fifth Ave, Pittsburgh, PA 15213
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794-5252
| | - Justin L MacCallum
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794-5252
| |
Collapse
|
48
|
Lopéz-Blanco JR, Garzón JI, Chacón P. iMod: multipurpose normal mode analysis in internal coordinates. ACTA ACUST UNITED AC 2011; 27:2843-50. [PMID: 21873636 DOI: 10.1093/bioinformatics/btr497] [Citation(s) in RCA: 147] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
MOTIVATION Dynamic simulations of systems with biologically relevant sizes and time scales are critical for understanding macromolecular functioning. Coarse-grained representations combined with normal mode analysis (NMA) have been established as an alternative to atomistic simulations. The versatility and efficiency of current approaches normally based on Cartesian coordinates can be greatly enhanced with internal coordinates (IC). RESULTS Here, we present a new versatile tool chest to explore conformational flexibility of both protein and nucleic acid structures using NMA in IC. Consideration of dihedral angles as variables reduces the computational cost and non-physical distortions of classical Cartesian NMA methods. Our proposed framework operates at different coarse-grained levels and offers an efficient framework to conduct NMA-based conformational studies, including standard vibrational analysis, Monte-Carlo simulations or pathway exploration. Examples of these approaches are shown to demonstrate its applicability, robustness and efficiency. CONTACT pablo@chaconlab.org SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- José Ramón Lopéz-Blanco
- Department of Biological Chemical Physics, Rocasolano Physical Chemistry Institute, CSIC, Serrano 119, Madrid 28006, Spain
| | | | | |
Collapse
|
49
|
Stein A, Rueda M, Panjkovich A, Orozco M, Aloy P. A Systematic Study of the Energetics Involved in Structural Changes upon Association and Connectivity in Protein Interaction Networks. Structure 2011; 19:881-9. [DOI: 10.1016/j.str.2011.03.009] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2010] [Revised: 02/22/2011] [Accepted: 03/13/2011] [Indexed: 01/26/2023]
|
50
|
Orozco M, Orellana L, Hospital A, Naganathan AN, Emperador A, Carrillo O, Gelpí JL. Coarse-grained representation of protein flexibility. Foundations, successes, and shortcomings. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2011; 85:183-215. [PMID: 21920324 DOI: 10.1016/b978-0-12-386485-7.00005-3] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Flexibility is the key magnitude to understand the variety of functions of proteins. Unfortunately, its experimental study is quite difficult, and in fact, most experimental procedures are designed to reduce flexibility and allow a better definition of the structure. Theoretical approaches have become then the alternative but face serious timescale problems, since many biologically relevant deformation movements happen in a timescale that is far beyond the possibility of current atomistic models. In this complex scenario, coarse-grained simulation methods have emerged as a powerful and inexpensive alternative. Along this chapter, we will review these coarse-grained methods, and explain their physical foundations and their range of applicability.
Collapse
Affiliation(s)
- Modesto Orozco
- Joint IRB-BSC Program on Computational Biology, Barcelona Supercomputing Center and Institute of Research in Biomedicine, Parc Científic de Barcelona, Barcelona, Spain
| | | | | | | | | | | | | |
Collapse
|