51
|
Wu H, Dalal Y, Papoian GA. Binding Dynamics of Disordered Linker Histone H1 with a Nucleosomal Particle. J Mol Biol 2021; 433:166881. [PMID: 33617899 DOI: 10.1016/j.jmb.2021.166881] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 02/11/2021] [Accepted: 02/12/2021] [Indexed: 01/30/2023]
Abstract
Linker histone H1 is an essential regulatory protein for many critical biological processes, such as eukaryotic chromatin packaging and gene expression. Mis-regulation of H1s is commonly observed in tumor cells, where the balance between different H1 subtypes has been shown to alter the cancer phenotype. Consisting of a rigid globular domain and two highly charged terminal domains, H1 can bind to multiple sites on a nucleosomal particle to alter chromatin hierarchical condensation levels. In particular, the disordered H1 amino- and carboxyl-terminal domains (NTD/CTD) are believed to enhance this binding affinity, but their detailed dynamics and functions remain unclear. In this work, we used a coarse-grained computational model, AWSEM-DNA, to simulate the H1.0b-nucleosome complex, namely chromatosome. Our results demonstrate that H1 disordered domains restrict the dynamics and conformation of both globular H1 and linker DNA arms, resulting in a more compact and rigid chromatosome particle. Furthermore, we identified regions of H1 disordered domains that are tightly tethered to DNA near the entry-exit site. Overall, our study elucidates at near-atomic resolution the way the disordered linker histone H1 modulates nucleosome's structural preferences and conformational dynamics.
Collapse
Affiliation(s)
- Hao Wu
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, United States
| | - Yamini Dalal
- Laboratory of Receptor Biology and Gene Expression, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States.
| | - Garegin A Papoian
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, United States; Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, United States.
| |
Collapse
|
52
|
Bari KJ, Prakashchand DD. Fundamental Challenges and Outlook in Simulating Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins. J Phys Chem Lett 2021; 12:1644-1656. [PMID: 33555894 DOI: 10.1021/acs.jpclett.0c03404] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Intrinsically disordered proteins (IDPs) populate an ensemble of dynamic conformations, making their structural characterization by experiments challenging. Many IDPs undergo liquid-liquid phase separation into dense membraneless organelles with myriad cellular functions. Multivalent interactions in low-complexity IDPs promote the formation of these subcellular coacervates. While solution NMR, Förster resonance energy transfer (FRET), and small-angle X-ray scattering (SAXS) studies on IDPs have their own challenges, recent computational methods draw a rational trade-off to characterize the driving forces underlying phase separation. In this Perspective, we critically evaluate the scope of approximation-free field theoretic simulations, well-tempered ensemble methods, enhanced sampling techniques, coarse-grained force fields, and slab simulation approaches to offer an improved understanding of phase separation. A synergy between simulation length scale and model resolution would reduce the existing caveats and enable theories of polymer physics to elucidate finer details of liquid-liquid phase separation (LLPS). These computational advances offer promise for rigorous characterization of the IDP proteome and designing peptides with tunable material and self-assembly properties.
Collapse
Affiliation(s)
- Khandekar Jishan Bari
- Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Gopanpally, Hyderabad 500107, India
- Department of Chemical Sciences, Indian Institute of Science Education and Research, Berhampur, Odisha 760010, India
| | - Dube Dheeraj Prakashchand
- Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Gopanpally, Hyderabad 500107, India
| |
Collapse
|
53
|
Mu J, Liu H, Zhang J, Luo R, Chen HF. Recent Force Field Strategies for Intrinsically Disordered Proteins. J Chem Inf Model 2021; 61:1037-1047. [PMID: 33591749 PMCID: PMC8256680 DOI: 10.1021/acs.jcim.0c01175] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Intrinsically disordered proteins (IDPs) are widely distributed across eukaryotic cells, playing important roles in molecular recognition, molecular assembly, post-translational modification, and other biological processes. IDPs are also associated with many diseases such as cancers, cardiovascular diseases, and neurodegenerative diseases. Due to their structural flexibility, conventional experimental methods cannot reliably capture their heterogeneous structures. Molecular dynamics simulation becomes an important complementary tool to quantify IDP structures. This review covers recent force field strategies proposed for more accurate molecular dynamics simulations of IDPs. The strategies include adjusting dihedral parameters, adding grid-based energy correction map (CMAP) parameters, refining protein-water interactions, and others. Different force fields were found to perform well on specific observables of specific IDPs but also are limited in reproducing all available experimental observables consistently for all tested IDPs. We conclude the review with perspective areas for improvements for future force fields for IDPs.
Collapse
Affiliation(s)
- Junxi Mu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 20025, China
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Molecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, California 92697-3900, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| |
Collapse
|
54
|
Benayad Z, von Bülow S, Stelzl LS, Hummer G. Simulation of FUS Protein Condensates with an Adapted Coarse-Grained Model. J Chem Theory Comput 2021; 17:525-537. [PMID: 33307683 PMCID: PMC7872324 DOI: 10.1021/acs.jctc.0c01064] [Citation(s) in RCA: 100] [Impact Index Per Article: 33.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Indexed: 01/02/2023]
Abstract
Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physicochemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein-protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low-complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS-LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model. For a scaling factor α = 0.65, we recover experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS-LCD condenses. In the region of phase separation, we simulate FUS-LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS-LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01-0.4 mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02 Pa s for the FUS-LCD droplets with scaling factors α in the range of 0.625-0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
Collapse
Affiliation(s)
- Zakarya Benayad
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
- Département
de Chimie, École Normale Supérieure, PSL University, 75005 Paris, France
| | - Sören von Bülow
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Lukas S. Stelzl
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
| | - Gerhard Hummer
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Straße 3, 60438 Frankfurt am Main, Germany
- Institute
for Biophysics, Goethe University Frankfurt, 60438 Frankfurt
am Main, Germany
| |
Collapse
|
55
|
Klein F, Barrera EE, Pantano S. Assessing SIRAH's Capability to Simulate Intrinsically Disordered Proteins and Peptides. J Chem Theory Comput 2021; 17:599-604. [PMID: 33411518 DOI: 10.1021/acs.jctc.0c00948] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The challenges posed by intrinsically disordered proteins (IDPs) to atomistic and coarse-grained (CG) simulations are boosting efforts to develop and reparametrize current force fields. An assessment of the dynamical behavior of IDPs' and unstructured peptides with the CG SIRAH force field suggests that the current version achieves a fair description of IDPs' conformational flexibility. Moreover, we found a remarkable capability to capture the effect of point mutations in loosely structured peptides.
Collapse
Affiliation(s)
- Florencia Klein
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo, CP 11400, Uruguay.,Graduate Program in Chemistry, Facultad de Química, Universidad de la República, Montevideo 11800, Uruguay
| | - Exequiel E Barrera
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo, CP 11400, Uruguay.,Instituto de Histología y Embriología (IHEM) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), CC56, Universidad Nacional de Cuyo (UNCuyo), M5500 Mendoza, Argentina
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo, CP 11400, Uruguay.,Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China
| |
Collapse
|
56
|
Fatafta H, Samantray S, Sayyed-Ahmad A, Coskuner-Weber O, Strodel B. Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2021; 183:135-185. [PMID: 34656328 DOI: 10.1016/bs.pmbts.2021.06.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.
Collapse
Affiliation(s)
- Hebah Fatafta
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Suman Samantray
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; AICES Graduate School, RWTH Aachen University, Aachen, Germany
| | | | - Orkid Coskuner-Weber
- Molecular Biotechnology, Turkish-German University, Sahinkaya Caddesi, Istanbul, Turkey
| | - Birgit Strodel
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.
| |
Collapse
|
57
|
Abstract
Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.
Collapse
|
58
|
Regmi R, Srinivasan S, Latham AP, Kukshal V, Cui W, Zhang B, Bose R, Schlau-Cohen GS. Phosphorylation-Dependent Conformations of the Disordered Carboxyl-Terminus Domain in the Epidermal Growth Factor Receptor. J Phys Chem Lett 2020; 11:10037-10044. [PMID: 33179922 PMCID: PMC8063277 DOI: 10.1021/acs.jpclett.0c02327] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
The epidermal growth factor receptor (EGFR), a receptor tyrosine kinase, regulates basic cellular functions and is a major target for anticancer therapeutics. The carboxyl-terminus domain is a disordered region of EGFR that contains the tyrosine residues, which undergo autophosphorylation followed by docking of signaling proteins. Local phosphorylation-dependent secondary structure has been identified and is thought to be associated with the signaling cascade. Deciphering and distinguishing the overall conformations, however, have been challenging because of the disordered nature of the carboxyl-terminus domain and resultant lack of well-defined three-dimensional structure for most of the domain. We investigated the overall conformational states of the isolated EGFR carboxyl-terminus domain using single-molecule Förster resonance energy transfer and coarse-grained simulations. Our results suggest that electrostatic interactions between charged residues emerge within the disordered domain upon phosphorylation, producing a looplike conformation. This conformation may enable binding of downstream signaling proteins and potentially reflect a general mechanism in which electrostatics transiently generate functional architectures in disordered regions of a well-folded protein.
Collapse
Affiliation(s)
- Raju Regmi
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Shwetha Srinivasan
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Vandna Kukshal
- Department of Medicine and Department of Cell Biology and Physiology, Washington University School of Medicine, St. Louis, Missouri 63110, United States
| | - Weidong Cui
- Department of Chemistry, Washington University School of Medicine, St. Louis, Missouri 63110, United States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Ron Bose
- Department of Medicine and Department of Cell Biology and Physiology, Washington University School of Medicine, St. Louis, Missouri 63110, United States
| | - Gabriela S Schlau-Cohen
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
59
|
Ozgur B, Sayar M. Representation of the conformational ensemble of peptides in coarse grained simulations. J Chem Phys 2020; 153:054108. [DOI: 10.1063/5.0012391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Mehmet Sayar
- Chemical and Biological Engineering and Mechanical Engineering Departments, College of Engineering, Koç University, 34450 Istanbul, Turkey
| |
Collapse
|
60
|
Kilgore HR, Latham AP, Ressler VT, Zhang B, Raines RT. Structure and Dynamics of N-Glycosylated Human Ribonuclease 1. Biochemistry 2020; 59:3148-3156. [PMID: 32544330 DOI: 10.1021/acs.biochem.0c00191] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Glycosylation is a common modification that can endow proteins with altered physical and biological properties. Ribonuclease 1 (RNase 1), which is the human homologue of the archetypal enzyme RNase A, undergoes N-linked glycosylation at asparagine residues 34, 76, and 88. We have produced the three individual glycoforms that display the core heptasaccharide, Man5GlcNAc2, and analyzed the structure of each glycoform by using small-angle X-ray scattering along with molecular dynamics simulations. The glycan on Asn34 is relatively compact and rigid, donates hydrogen bonds that "cap" the carbonyl groups at the C-terminus of an α-helix, and enhances protein thermostability. In contrast, the glycan on Asn88 is flexible and can even enter the enzymic active site, hindering catalysis. The N-glycosylation of Asn76 has less pronounced consequences. These data highlight the diverse behaviors of Man5GlcNAc2 pendants and provide a structural underpinning to the functional consequences of protein glycosylation.
Collapse
Affiliation(s)
- Henry R Kilgore
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Valerie T Ressler
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Ronald T Raines
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
61
|
Mioduszewski Ł, Różycki B, Cieplak M. Pseudo-Improper-Dihedral Model for Intrinsically Disordered Proteins. J Chem Theory Comput 2020; 16:4726-4733. [PMID: 32436706 PMCID: PMC7588027 DOI: 10.1021/acs.jctc.0c00338] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
We present a new coarse-grained Cα-based protein model with a nonradial multibody pseudo-improper-dihedral potential that is transferable, time-independent, and suitable for molecular dynamics. It captures the nature of backbone and side-chain interactions between amino acid residues by adapting a simple improper dihedral term for a one-bead-per-residue model. It is parameterized for intrinsically disordered proteins and applicable to simulations of such proteins and their assemblies on millisecond time scales.
Collapse
Affiliation(s)
- Łukasz Mioduszewski
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Bartosz Różycki
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
62
|
Ferrie JJ, Petersson EJ. A Unified De Novo Approach for Predicting the Structures of Ordered and Disordered Proteins. J Phys Chem B 2020; 124:5538-5548. [PMID: 32525675 DOI: 10.1021/acs.jpcb.0c02924] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
As recognition of the abundance and relevance of intrinsically disordered proteins (IDPs) continues to grow, demand increases for methods that can rapidly predict the conformational ensembles populated by these proteins. To date, IDP simulations have largely been dominated by molecular dynamics (MD) simulations, which require significant compute times and/or complex hardware. Recent developments in MD have afforded methods capable of simulating both ordered and disordered proteins, yet to date, accurate fold prediction from a sequence has been dominated by Monte Carlo (MC)-based methods such as Rosetta. To overcome the limitations of current approaches in IDP simulation using Rosetta while maintaining its utility for modeling folded domains, we developed PyRosetta-based algorithms that allow for the accurate de novo prediction of proteins across all degrees of foldedness along with structural ensembles of disordered proteins. Our simulations have accuracy comparable to state-of-the-art MD with vastly reduced computational demands.
Collapse
Affiliation(s)
- John J Ferrie
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, United States
| | - E James Petersson
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, United States
| |
Collapse
|
63
|
Wäschenbach L, Gertzen CGW, Keitel V, Gohlke H. Dimerization energetics of the G-protein coupled bile acid receptor TGR5 from all-atom simulations. J Comput Chem 2019; 41:874-884. [PMID: 31880348 DOI: 10.1002/jcc.26135] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 12/07/2019] [Accepted: 12/09/2019] [Indexed: 12/20/2022]
Abstract
We describe the first extensive energetic evaluation of GPCR dimerization on the atomistic level by means of potential of mean force (PMF) computations and implicit solvent/implicit membrane end-point free energy calculations (MM-PBSA approach). Free energies of association computed from the PMFs show that the formation of both the 1/8 and 4/5 interface is energetically favorable for TGR5, the first GPCR known to be activated by hydrophobic bile acids and neurosteroids. Furthermore, formation of the 1/8 interface is favored over that of the 4/5 interface. Both results are in line with our previous FRET experiments in live cells. Differences in lipid-protein interactions are identified to contribute to the observed differences in free energies of association. A per-residue decomposition of the MM-PBSA effective binding energy reveals hot spot residues specific for both interfaces that form clusters. This knowledge may be used to guide the design of dimerization inhibitors or perform mutational studies to explore physiological consequences of distorted TGR5 association. Finally, we characterized the role of Y111, located in the conserved (D/E)RY motif, as a facilitator of TGR5 interactions. The types of computations performed here should be transferable to other transmembrane proteins that form dimers or higher oligomers as long as good structural models of the dimeric or oligomeric states are available. Such computations may help to overcome current restrictions due to an imperfect energetic representation of protein association at the coarse-grained level. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Lucas Wäschenbach
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Christoph G W Gertzen
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany.,Clinic for Gastroenterology, Hepatology and Infectious Diseases, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), and Institute for Complex Systems-Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| | - Verena Keitel
- Clinic for Gastroenterology, Hepatology and Infectious Diseases, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), and Institute for Complex Systems-Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| |
Collapse
|
64
|
Latham AP, Zhang B. Maximum Entropy Optimized Force Field for Intrinsically Disordered Proteins. J Chem Theory Comput 2019; 16:773-781. [PMID: 31756104 DOI: 10.1021/acs.jctc.9b00932] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Intrinsically disordered proteins (IDPs) constitute a significant fraction of eukaryotic proteomes. High-resolution characterization of IDP conformational ensembles can help elucidate their roles in a wide range of biological processes but remains challenging both experimentally and computationally. Here, we present a generic algorithm to improve the accuracy of coarse-grained IDP models using a diverse set of experimental measurements. It combines maximum entropy optimization and least-squares regression to systematically adjust model parameters and improve the agreement between simulation and experiment. We successfully applied the algorithm to derive a transferable force field, which we term the maximum entropy optimized force field (MOFF), for de novo prediction of IDP structures. Statistical analysis of force field parameters reveals features of amino acid interactions not captured by potentials designed to work well for folded proteins. We anticipate its combination of efficiency and accuracy will make MOFF useful for studying the phase separation of IDPs, which drives the formation of various biological compartments.
Collapse
Affiliation(s)
- Andrew P Latham
- Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Bin Zhang
- Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
65
|
Roche J, Potoyan DA. Disorder Mediated Oligomerization of DISC1 Proteins Revealed by Coarse-Grained Molecular Dynamics Simulations. J Phys Chem B 2019; 123:9567-9575. [PMID: 31614085 DOI: 10.1021/acs.jpcb.9b07467] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Disrupted-in-schizophrenia-1 (DISC1) is a scaffold protein of significant importance for neuro-development and a prominent candidate protein in the etiology of mental disorders. In this work, we investigate the role of conformational heterogeneity and local structural disorder in the oligomerization pathway of the full-length DISC1 and of two truncation variants. Through extensive coarse-grained molecular dynamics simulations with a predictive energy landscape-based model, we shed light on the interplay of local and global disorder which lead to different oligomerization pathways. We found that both global conformational heterogeneity and local structural disorder play an important role in shaping distinct oligomerization trends of DISC1. This study also sheds light on the differences in oligomerization pathways of the full-length protein compared to the truncated variants produced by a chromosomal translocation associated with schizophrenia. We report that oligomerization of full-length DISC1 sequence works in a nonadditive manner with respect to truncated fragments that do not mirror the conformational landscape or binding affinities of the full-length unit.
Collapse
Affiliation(s)
- Julien Roche
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology , Iowa State University , Ames , Iowa 50011 , United States
| | - Davit A Potoyan
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology , Iowa State University , Ames , Iowa 50011 , United States.,Department of Chemistry , Iowa State University , Ames , Iowa 50011 , United States.,Bioinformatics and Computational Biology Program , Iowa State University , Ames , Iowa 50011 , United States
| |
Collapse
|
66
|
Hermann MR, Hub JS. SAXS-Restrained Ensemble Simulations of Intrinsically Disordered Proteins with Commitment to the Principle of Maximum Entropy. J Chem Theory Comput 2019; 15:5103-5115. [DOI: 10.1021/acs.jctc.9b00338] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Markus R. Hermann
- Institute for Microbiology and Genetics, Georg-August-Universität Göttingen, 37077 Göttingen, Germany
| | - Jochen S. Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Campus E2 6, 66123 Saarbrücken, Germany
| |
Collapse
|
67
|
Baul U, Chakraborty D, Mugnai ML, Straub JE, Thirumalai D. Sequence Effects on Size, Shape, and Structural Heterogeneity in Intrinsically Disordered Proteins. J Phys Chem B 2019; 123:3462-3474. [PMID: 30913885 PMCID: PMC6920032 DOI: 10.1021/acs.jpcb.9b02575] [Citation(s) in RCA: 105] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Intrinsically disordered proteins (IDPs) lack well-defined three-dimensional structures, thus challenging the archetypal notion of structure-function relationships. Determining the ensemble of conformations that IDPs explore under physiological conditions is the first step toward understanding their diverse cellular functions. Here, we quantitatively characterize the structural features of IDPs as a function of sequence and length using coarse-grained simulations. For diverse IDP sequences, with the number of residues ( NT) ranging from 20 to 441, our simulations not only reproduce the radii of gyration ( Rg) obtained from experiments, but also predict the full scattering intensity profiles in excellent agreement with small-angle X-ray scattering experiments. The Rg values are well-described by the standard Flory scaling law, Rg = Rg0 NTν, with ν ≈ 0.588, making it tempting to assert that IDPs behave as polymers in a good solvent. However, clustering analysis reveals that the menagerie of structures explored by IDPs is diverse, with the extent of heterogeneity being highly sequence-dependent, even though ensemble-averaged properties, such as the dependence of Rg on chain length, may suggest synthetic polymer-like behavior in a good solvent. For example, we show that for the highly charged Prothymosin-α, a substantial fraction of conformations is highly compact. Even if the sequence compositions are similar, as is the case for α-Synuclein and a truncated construct from the Tau protein, there are substantial differences in the conformational heterogeneity. Taken together, these observations imply that metrics based on net charge or related quantities alone cannot be used to anticipate the phases of IDPs, either in isolation or in complex with partner IDPs or RNA. Our work sets the stage for probing the interactions of IDPs with each other, with folded protein domains, or with partner RNAs, which are critical for describing the structures of stress granules and biomolecular condensates with important cellular functions.
Collapse
Affiliation(s)
- Upayan Baul
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Mauro L. Mugnai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - D. Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
68
|
Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins. Biomolecules 2019; 9:biom9040146. [PMID: 30979035 PMCID: PMC6523529 DOI: 10.3390/biom9040146] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2019] [Revised: 04/09/2019] [Accepted: 04/10/2019] [Indexed: 01/09/2023] Open
Abstract
Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.
Collapse
|
69
|
Bignon E, Rizza S, Filomeni G, Papaleo E. Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase. Comput Struct Biotechnol J 2019; 17:415-429. [PMID: 30996821 PMCID: PMC6451115 DOI: 10.1016/j.csbj.2019.03.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 03/17/2019] [Accepted: 03/21/2019] [Indexed: 12/25/2022] Open
Abstract
Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.
Collapse
Affiliation(s)
- Emmanuelle Bignon
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Salvatore Rizza
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Giuseppe Filomeni
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Department of Biology, University of Rome Tor Vergata, Rome, Italy
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Translational Disease Systems Biology, Faculty of Health and Medical Sciences, Novo Nordisk Foundation Center for Protein Research University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
70
|
Structural and Dynamical Order of a Disordered Protein: Molecular Insights into Conformational Switching of PAGE4 at the Systems Level. Biomolecules 2019; 9:biom9020077. [PMID: 30813315 PMCID: PMC6406393 DOI: 10.3390/biom9020077] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Revised: 02/10/2019] [Accepted: 02/10/2019] [Indexed: 01/10/2023] Open
Abstract
Folded proteins show a high degree of structural order and undergo (fairly constrained) collective motions related to their functions. On the other hand, intrinsically disordered proteins (IDPs), while lacking a well-defined three-dimensional structure, do exhibit some structural and dynamical ordering, but are less constrained in their motions than folded proteins. The larger structural plasticity of IDPs emphasizes the importance of entropically driven motions. Many IDPs undergo function-related disorder-to-order transitions driven by their interaction with specific binding partners. As experimental techniques become more sensitive and become better integrated with computational simulations, we are beginning to see how the modest structural ordering and large amplitude collective motions of IDPs endow them with an ability to mediate multiple interactions with different partners in the cell. To illustrate these points, here, we use Prostate-associated gene 4 (PAGE4), an IDP implicated in prostate cancer (PCa) as an example. We first review our previous efforts using molecular dynamics simulations based on atomistic AWSEM to study the conformational dynamics of PAGE4 and how its motions change in its different physiologically relevant phosphorylated forms. Our simulations quantitatively reproduced experimental observations and revealed how structural and dynamical ordering are encoded in the sequence of PAGE4 and can be modulated by different extents of phosphorylation by the kinases HIPK1 and CLK2. This ordering is reflected in changing populations of certain secondary structural elements as well as in the regularity of its collective motions. These ordered features are directly correlated with the functional interactions of WT-PAGE4, HIPK1-PAGE4 and CLK2-PAGE4 with the AP-1 signaling axis. These interactions give rise to repeated transitions between (high HIPK1-PAGE4, low CLK2-PAGE4) and (low HIPK1-PAGE4, high CLK2-PAGE4) cell phenotypes, which possess differing sensitivities to the standard PCa therapies, such as androgen deprivation therapy (ADT). We argue that, although the structural plasticity of an IDP is important in promoting promiscuous interactions, the modulation of the structural ordering is important for sculpting its interactions so as to rewire with agility biomolecular interaction networks with significant functional consequences.
Collapse
|
71
|
Latham AP, Zhang B. Improving Coarse-Grained Protein Force Fields with Small-Angle X-ray Scattering Data. J Phys Chem B 2019; 123:1026-1034. [PMID: 30620594 DOI: 10.1021/acs.jpcb.8b10336] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Small-angle X-ray scattering (SAXS) experiments provide valuable structural data for biomolecules in solution. We develop a highly efficient maximum entropy approach to fit SAXS data by introducing minimal biases to a coarse-grained protein force field, the associative memory, water mediated, structure, and energy model (AWSEM). We demonstrate that the resulting force field, AWSEM-SAXS, succeeds in reproducing scattering profiles and models protein structures with shapes that are in much better agreement with experimental results. Quantitative metrics further reveal a modest, but consistent, improvement in the accuracy of modeled structures when SAXS data are incorporated into the force field. Additionally, when applied to a multiconformational protein, we find that AWSEM-SAXS is able to recover the population of different protein conformations from SAXS data alone. We, therefore, conclude that the maximum entropy approach is effective in fine-tuning the force field to better characterize both protein structure and conformational fluctuation.
Collapse
Affiliation(s)
- Andrew P Latham
- Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Bin Zhang
- Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
72
|
|