1
|
Rizuan A, Shenoy J, Mohanty P, dos Passos PMS, Mercado Ortiz JF, Bai L, Viswanathan R, Wang SH, Johnson V, Mamede LD, Ayala YM, Ghirlando R, Mittal J, Fawzi NL. Structural details of helix-mediated TDP-43 C-terminal domain multimerization. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.05.602258. [PMID: 39005345 PMCID: PMC11245101 DOI: 10.1101/2024.07.05.602258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
The primarily disordered C-terminal domain (CTD) of TAR DNA binding protein-43 (TDP-43), a key nuclear protein in RNA metabolism, forms neuronal inclusions in several neurodegenerative diseases. A conserved region (CR, spanning residues 319-341) in CTD forms transient helix-helix contacts important for its higher-order oligomerization and function that are disrupted by ALS-associated mutations. However, the structural details of CR assembly and the explanation for several ALS-associated variants' impact on phase separation and function remain unclear due to challenges in analyzing the dynamic association of TDP-43 CTD using traditional structural biology approaches. By employing an integrative approach, combining biophysical experiments, biochemical assays, AlphaFold2-Multimer (AF2-Multimer), and atomistic simulations, we generated structural models of helical oligomerization of TDP-43 CR. Using NMR, we first established that the native state of TDP-43 CR under physiological conditions is α-helical. Next, alanine scanning mutagenesis revealed that while hydrophobic residues in the CR are important for CR assembly, phase separation and TDP-43 nuclear retention function, polar residues down regulate these processes. Finally, pairing AF2-Multimer modeling with AAMD simulations indicated that dynamic, oligomeric assemblies of TDP-43 that are stabilized by a methionine-rich core with specific contributions from a tryptophan/leucine pair. In conclusion, our results advance the structural understanding of the mechanisms driving TDP-43 function and provide a window into the initial stages of its conversion into pathogenic aggregates.
Collapse
Affiliation(s)
- Azamat Rizuan
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | - Jayakrishna Shenoy
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Priyesh Mohanty
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | - Patricia M. S. dos Passos
- Edward Doisy Department of Biochemistry and Molecular Biology, Saint Louis University School of Medicine, St. Louis, Missouri
| | - José F. Mercado Ortiz
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Leanna Bai
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Renjith Viswanathan
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Szu-Huan Wang
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Victoria Johnson
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| | - Lohany D. Mamede
- Edward Doisy Department of Biochemistry and Molecular Biology, Saint Louis University School of Medicine, St. Louis, Missouri
| | - Yuna M. Ayala
- Edward Doisy Department of Biochemistry and Molecular Biology, Saint Louis University School of Medicine, St. Louis, Missouri
| | - Rodolfo Ghirlando
- Laboratory of Molecular Biology, National Institute of Diabetes, Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
- Department of Chemistry, Texas A&M University, College Station, TX 77843
- Interdisciplinary Graduate Program in Genetics and Genomics, Texas A&M University, College Station, TX 77843
| | - Nicolas L. Fawzi
- Department of Molecular Biology, Cell Biology & Biochemistry, Brown University, Providence, RI 02912
| |
Collapse
|
2
|
Borges-Araújo L, Pereira GP, Valério M, Souza PCT. Assessing the Martini 3 protein model: A review of its path and potential. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:141014. [PMID: 38670324 DOI: 10.1016/j.bbapap.2024.141014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/13/2024] [Accepted: 04/09/2024] [Indexed: 04/28/2024]
Abstract
Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Gilberto P Pereira
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Mariana Valério
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Paulo C T Souza
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France.
| |
Collapse
|
3
|
Mohanty P, Phan TM, Mittal J. Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.03.592468. [PMID: 38766024 PMCID: PMC11100600 DOI: 10.1101/2024.05.03.592468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.
Collapse
|
4
|
Holehouse AS, Kragelund BB. The molecular basis for cellular function of intrinsically disordered protein regions. Nat Rev Mol Cell Biol 2024; 25:187-211. [PMID: 37957331 DOI: 10.1038/s41580-023-00673-0] [Citation(s) in RCA: 42] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/26/2023] [Indexed: 11/15/2023]
Abstract
Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.
Collapse
Affiliation(s)
- Alex S Holehouse
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St Louis, MO, USA.
- Center for Biomolecular Condensates, Washington University in St Louis, St Louis, MO, USA.
| | - Birthe B Kragelund
- REPIN, Structural Biology and NMR Laboratory, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
5
|
Galano‐Frutos JJ, Sancho J. Energy, water, and protein folding: A molecular dynamics-based quantitative inventory of molecular interactions and forces that make proteins stable. Protein Sci 2024; 33:e4905. [PMID: 38284492 PMCID: PMC10804899 DOI: 10.1002/pro.4905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 12/12/2023] [Accepted: 01/05/2024] [Indexed: 01/30/2024]
Abstract
Protein folding energetics can be determined experimentally on a case-by-case basis but it is not understood in sufficient detail to provide deep control in protein design. The fundamentals of protein stability have been outlined by calorimetry, protein engineering, and biophysical modeling, but these approaches still face great difficulty in elucidating the specific contributions of the intervening molecules and physical interactions. Recently, we have shown that the enthalpy and heat capacity changes associated to the protein folding reaction can be calculated within experimental error using molecular dynamics simulations of native protein structures and their corresponding unfolded ensembles. Analyzing in depth molecular dynamics simulations of four model proteins (CI2, barnase, SNase, and apoflavodoxin), we dissect here the energy contributions to ΔH (a key component of protein stability) made by the molecular players (polypeptide and solvent molecules) and physical interactions (electrostatic, van der Waals, and bonded) involved. Although the proteins analyzed differ in length, isoelectric point and fold class, their folding energetics is governed by the same quantitative pattern. Relative to the unfolded ensemble, the native conformations are enthalpically stabilized by comparable contributions from protein-protein and solvent-solvent interactions, and almost equally destabilized by interactions between protein and solvent molecules. The native protein surface seems to interact better with water than the unfolded one, but this is outweighed by the unfolded surface being larger. From the perspective of physical interactions, the native conformations are stabilized by van de Waals and Coulomb interactions and destabilized by conformational strain arising from bonded interactions. Also common to the four proteins, the sign of the heat capacity change is set by interactions between protein and solvent molecules or, from the alternative perspective, by Coulomb interactions.
Collapse
Affiliation(s)
- Juan José Galano‐Frutos
- Biocomputation and Complex Systems Physics Institute (BIFI)‐Joint Unit GBsC‐CSICUniversity of ZaragozaZaragozaSpain
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de CienciasUniversity of ZaragozaZaragozaSpain
| | - Javier Sancho
- Biocomputation and Complex Systems Physics Institute (BIFI)‐Joint Unit GBsC‐CSICUniversity of ZaragozaZaragozaSpain
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de CienciasUniversity of ZaragozaZaragozaSpain
- Aragon Health Research Institute (IIS Aragón)ZaragozaSpain
| |
Collapse
|
6
|
Mohanty P, Rizuan A, Kim YC, Fawzi NL, Mittal J. A complex network of interdomain interactions underlies the conformational ensemble of monomeric TDP-43 and modulates its phase behavior. Protein Sci 2024; 33:e4891. [PMID: 38160320 PMCID: PMC10804676 DOI: 10.1002/pro.4891] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 12/07/2023] [Accepted: 12/28/2023] [Indexed: 01/03/2024]
Abstract
TAR DNA-binding protein 43 (TDP-43) is a multidomain protein involved in the regulation of RNA metabolism, and its aggregates have been observed in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Numerous studies indicate TDP-43 can undergo liquid-liquid phase separation (LLPS) in vitro and is a component of biological condensates. Homo-oligomerization via the folded N-terminal domain (aa:1-77) and the conserved helical region (aa:319-341) of the disordered, C-terminal domain is found to be an important driver of TDP-43 phase separation. However, a comprehensive molecular view of TDP-43 phase separation, particularly regarding the nature of heterodomain interactions, is lacking due to the challenges associated with its stability and purification. Here, we utilize all-atom and coarse-grained (CG) molecular dynamics (MD) simulations to uncover the network of interdomain interactions implicated in TDP-43 phase separation. All-atom simulations uncovered the presence of transient, interdomain interactions involving flexible linkers, RNA-recognition motif (RRM) domains and a charged segment of disordered C-terminal domain (CTD). CG simulations indicate these inter-domain interactions which affect the conformational landscape of TDP-43 in the dilute phase are also prevalent in the condensed phase. Finally, sequence and surface charge distribution analysis coupled with all-atom simulations (at high salt) confirmed that the transient interdomain contacts are predominantly electrostatic in nature. Overall, our findings from multiscale simulations lead to a greater appreciation of the complex interaction network underlying the structural landscape and phase separation of TDP-43.
Collapse
Affiliation(s)
- Priyesh Mohanty
- Artie McFerrin Department of Chemical EngineeringTexas A&M UniversityCollege StationTexasUSA
| | - Azamat Rizuan
- Artie McFerrin Department of Chemical EngineeringTexas A&M UniversityCollege StationTexasUSA
| | - Young C. Kim
- Naval Research LaboratoryCenter for Materials Physics and TechnologyWashingtonDistrict of ColumbiaUSA
| | - Nicolas L. Fawzi
- Department of Molecular Biology, Cell Biology and BiochemistryProvidenceRhode IslandUSA
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical EngineeringTexas A&M UniversityCollege StationTexasUSA
- Department of ChemistryTexas A&M UniversityCollege StationTexasUSA
- Interdisciplinary Graduate Program in Genetics and GenomicsTexas A&M UniversityCollege StationTexasUSA
| |
Collapse
|
7
|
AlRawashdeh S, Barakat KH. Applications of Molecular Dynamics Simulations in Drug Discovery. Methods Mol Biol 2024; 2714:127-141. [PMID: 37676596 DOI: 10.1007/978-1-0716-3441-7_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.
Collapse
Affiliation(s)
- Sara AlRawashdeh
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB, Canada
| | - Khaled H Barakat
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB, Canada.
| |
Collapse
|
8
|
Aduriz-Arrizabalaga J, Lopez X, De Sancho D. Atomistic molecular simulations of Aβ-Zn conformational ensembles. Proteins 2024; 92:134-144. [PMID: 37746887 DOI: 10.1002/prot.26590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 09/01/2023] [Accepted: 09/07/2023] [Indexed: 09/26/2023]
Abstract
The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.
Collapse
Affiliation(s)
- Julen Aduriz-Arrizabalaga
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, UPV/EHU & Donostia International Physics Center (DIPC), Donostia-San Sebastian, Euskadi, Spain
| | - Xabier Lopez
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, UPV/EHU & Donostia International Physics Center (DIPC), Donostia-San Sebastian, Euskadi, Spain
| | - David De Sancho
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, UPV/EHU & Donostia International Physics Center (DIPC), Donostia-San Sebastian, Euskadi, Spain
| |
Collapse
|
9
|
Balasubramanian S, Maharana S, Srivastava A. "Boundary residues" between the folded RNA recognition motif and disordered RGG domains are critical for FUS-RNA binding. J Biol Chem 2023; 299:105392. [PMID: 37890778 PMCID: PMC10687056 DOI: 10.1016/j.jbc.2023.105392] [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: 02/21/2023] [Revised: 09/19/2023] [Accepted: 10/19/2023] [Indexed: 10/29/2023] Open
Abstract
Fused in sarcoma (FUS) is an abundant RNA-binding protein, which drives phase separation of cellular condensates and plays multiple roles in RNA regulation. The RNA-binding ability of FUS protein is crucial to its cellular function. Here, our molecular simulation study on the FUS-RNA complex provides atomic resolution insights into the observations from biochemical studies and also illuminates our understanding of molecular driving forces that mediate the structure, stability, and interaction of the RNA recognition motif (RRM) and RGG domains of FUS with a stem-loop junction RNA. We observe clear cooperativity and division of labor among the ordered (RRM) and disordered domains (RGG1 and RGG2) of FUS that leads to an organized and tighter RNA binding. Irrespective of the length of RGG2, the RGG2-RNA interaction is confined to the stem-loop junction and the proximal stem regions. On the other hand, the RGG1 interactions are primarily with the longer RNA stem. We find that the C terminus of RRM, which make up the "boundary residues" that connect the folded RRM with the long disordered RGG2 stretch of the protein, plays a critical role in FUS-RNA binding. Our study provides high-resolution molecular insights into the FUS-RNA interactions and forms the basis for understanding the molecular origins of full-length FUS interaction with RNA.
Collapse
Affiliation(s)
| | - Shovamayee Maharana
- Department of Molecular and Cell Biology, Indian Institute of Science Bangalore, Bangalore, Karnataka, India
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, Bangalore, Karnataka, India.
| |
Collapse
|
10
|
Vural D, Shrestha UR, Petridis L, Smith JC. Water molecule ordering on the surface of an intrinsically disordered protein. Biophys J 2023; 122:4326-4335. [PMID: 37838830 PMCID: PMC10722392 DOI: 10.1016/j.bpj.2023.10.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 09/26/2023] [Accepted: 10/10/2023] [Indexed: 10/16/2023] Open
Abstract
The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.
Collapse
Affiliation(s)
- Derya Vural
- Department of Physics, Marmara University, Istanbul, Türkiye; Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, Tennessee; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee.
| | - Utsab R Shrestha
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, Tennessee; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee
| | - Loukas Petridis
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, Tennessee; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee
| | - Jeremy C Smith
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, Tennessee; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee
| |
Collapse
|
11
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
12
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python Package for the Analysis of Simulations of Intrinsically Disordered Proteins. J Chem Theory Comput 2023; 19:5609-5620. [PMID: 37463458 PMCID: PMC11188088 DOI: 10.1021/acs.jctc.3c00190] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite simulation analysis of unfolded regions of proteins (SOURSOP), an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
Affiliation(s)
- Jared M. Lalmansingh
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex T. Keeley
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana-Champaign, IL 61801, USA
| | - Kiersten M. Ruff
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Rohit V. Pappu
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex S. Holehouse
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, 63110, USA
| |
Collapse
|
13
|
Saar KL, Qian D, Good LL, Morgunov AS, Collepardo-Guevara R, Best RB, Knowles TPJ. Theoretical and Data-Driven Approaches for Biomolecular Condensates. Chem Rev 2023; 123:8988-9009. [PMID: 37171907 PMCID: PMC10375482 DOI: 10.1021/acs.chemrev.2c00586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Indexed: 05/14/2023]
Abstract
Biomolecular condensation processes are increasingly recognized as a fundamental mechanism that living cells use to organize biomolecules in time and space. These processes can lead to the formation of membraneless organelles that enable cells to perform distinct biochemical processes in controlled local environments, thereby supplying them with an additional degree of spatial control relative to that achieved by membrane-bound organelles. This fundamental importance of biomolecular condensation has motivated a quest to discover and understand the molecular mechanisms and determinants that drive and control this process. Within this molecular viewpoint, computational methods can provide a unique angle to studying biomolecular condensation processes by contributing the resolution and scale that are challenging to reach with experimental techniques alone. In this Review, we focus on three types of dry-lab approaches: theoretical methods, physics-driven simulations and data-driven machine learning methods. We review recent progress in using these tools for probing biomolecular condensation across all three fields and outline the key advantages and limitations of each of the approaches. We further discuss some of the key outstanding challenges that we foresee the community addressing next in order to develop a more complete picture of the molecular driving forces behind biomolecular condensation processes and their biological roles in health and disease.
Collapse
Affiliation(s)
- Kadi L. Saar
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
- Transition
Bio Ltd., Cambridge, United Kingdom
| | - Daoyuan Qian
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Lydia L. Good
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
- Laboratory
of Chemical Physics, National Institute of Diabetes and Digestive
and Kidney Diseases, National Institutes
of Health, Bethesda, Maryland 20892, United States
| | - Alexey S. Morgunov
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Rosana Collepardo-Guevara
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
- Department
of Genetics, University of Cambridge, Cambridge CB2 3EH, United Kingdom
| | - Robert B. Best
- Laboratory
of Chemical Physics, National Institute of Diabetes and Digestive
and Kidney Diseases, National Institutes
of Health, Bethesda, Maryland 20892, United States
| | - Tuomas P. J. Knowles
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Cambridge CB2 1EW, United Kingdom
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
14
|
Appadurai R, Koneru JK, Bonomi M, Robustelli P, Srivastava A. Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding. J Chem Theory Comput 2023; 19:4711-4727. [PMID: 37338049 PMCID: PMC11108026 DOI: 10.1021/acs.jctc.3c00224] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.
Collapse
Affiliation(s)
- Rajeswari Appadurai
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | | | - Massimiliano Bonomi
- Structural Bioinformatics Unit, Department of Structural Biology and Chemistry. CNRS UMR 3528, C3BI, CNRS USR 3756, Institut Pasteur, Paris, France
| | - Paul Robustelli
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755, USA
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| |
Collapse
|
15
|
Koren G, Meir S, Holschuh L, Mertens HDT, Ehm T, Yahalom N, Golombek A, Schwartz T, Svergun DI, Saleh OA, Dzubiella J, Beck R. Intramolecular structural heterogeneity altered by long-range contacts in an intrinsically disordered protein. Proc Natl Acad Sci U S A 2023; 120:e2220180120. [PMID: 37459524 PMCID: PMC10372579 DOI: 10.1073/pnas.2220180120] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2022] [Accepted: 06/02/2023] [Indexed: 07/20/2023] Open
Abstract
Short-range interactions and long-range contacts drive the 3D folding of structured proteins. The proteins' structure has a direct impact on their biological function. However, nearly 40% of the eukaryotes proteome is composed of intrinsically disordered proteins (IDPs) and protein regions that fluctuate between ensembles of numerous conformations. Therefore, to understand their biological function, it is critical to depict how the structural ensemble statistics correlate to the IDPs' amino acid sequence. Here, using small-angle X-ray scattering and time-resolved Förster resonance energy transfer (trFRET), we study the intramolecular structural heterogeneity of the neurofilament low intrinsically disordered tail domain (NFLt). Using theoretical results of polymer physics, we find that the Flory scaling exponent of NFLt subsegments correlates linearly with their net charge, ranging from statistics of ideal to self-avoiding chains. Surprisingly, measuring the same segments in the context of the whole NFLt protein, we find that regardless of the peptide sequence, the segments' structural statistics are more expanded than when measured independently. Our findings show that while polymer physics can, to some level, relate the IDP's sequence to its ensemble conformations, long-range contacts between distant amino acids play a crucial role in determining intramolecular structures. This emphasizes the necessity of advanced polymer theories to fully describe IDPs ensembles with the hope that it will allow us to model their biological function.
Collapse
Affiliation(s)
- Gil Koren
- The School of Physics and Astronomy, Department of Condensed Matter, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
| | - Sagi Meir
- The School of Physics and Astronomy, Department of Condensed Matter, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
| | - Lennard Holschuh
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universit Freiburg, FreiburgD-79104, Germany
| | | | - Tamara Ehm
- The School of Physics and Astronomy, Department of Condensed Matter, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
- Faculty of Physics and Center for NanoScience, Ludwig-Maximilians-Universität, MünchenD-80539, Germany
| | - Nadav Yahalom
- The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
- School of Chemistry, Raymond and Beverly Sackler Faculty of Exact Sciences and Tel Aviv University Center for Light–Matter Interaction, Tel Aviv University, Tel Aviv6997801, Israel
| | - Adina Golombek
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
- School of Chemistry, Raymond and Beverly Sackler Faculty of Exact Sciences and Tel Aviv University Center for Light–Matter Interaction, Tel Aviv University, Tel Aviv6997801, Israel
| | - Tal Schwartz
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
- School of Chemistry, Raymond and Beverly Sackler Faculty of Exact Sciences and Tel Aviv University Center for Light–Matter Interaction, Tel Aviv University, Tel Aviv6997801, Israel
| | - Dmitri I. Svergun
- European Molecular Biology Laboratory, Hamburg Unit, Hamburg22607, Germany
| | - Omar A. Saleh
- BMSE Program, University of California, Santa Barbara, CA93110
- Materials Department, University of California, Santa Barbara, CA93110
| | - Joachim Dzubiella
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universit Freiburg, FreiburgD-79104, Germany
- Cluster of Excellence livMatS @ FIT–Freiburg Center for Interactive Materials and Bioinspired Technologies, Albert-Ludwigs-Universit Freiburg, FreiburgD-79104, Germany
| | - Roy Beck
- The School of Physics and Astronomy, Department of Condensed Matter, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv69978, Israel
- The Center for Nanoscience and Nanotechnology, Tel Aviv University, Tel Aviv69978, Israel
| |
Collapse
|
16
|
Chowdhury UD, Paul A, Bhargava BL. Interaction of the tau fibrils with the neuronal membrane. Biophys Chem 2023; 298:107024. [PMID: 37104971 DOI: 10.1016/j.bpc.2023.107024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Revised: 04/17/2023] [Accepted: 04/17/2023] [Indexed: 04/29/2023]
Abstract
Tau proteins are recently gaining a lot of interest due to their active role in causing a range of tauopathies. Molecular mechanisms underlying the tau interaction with the neuronal membrane are hitherto unknown and difficult to characterize using experimental methods. Using the cryo-EM structure of the tau-fibrils we have used atomistic molecular dynamics simulation to model the tau fibril and neuronal membrane interaction using explicit solvation. The dynamics and structural characteristics of the tau fibril with the neuronal membrane are compared to the tau fibril in the aqueous phase to corroborate the effect of the neuronal membrane in the tau structure. Tau fibrils have been modelled using CHARMM-36m force field and the six component neuronal membrane composition is taken from the earlier simulation results. The timescale conceivable in our molecular dynamics simulations is of the order of microseconds which captures the onset of the interaction of the tau fibrils with the neuronal membrane. This interaction is found to impact the tau pathogenesis that finally causes neuronal toxicity. Our study initiates the understanding of tau conformational ensemble in the presence of neuronal membrane and sheds the light on the significant tau-membrane interactions.
Collapse
Affiliation(s)
- Unmesh D Chowdhury
- School of Chemical Sciences, National Institute of Science Education & Research-Bhubaneswar, An OCC of Homi Bhabha National Institute, P.O. Jatni, Khurda, Odisha 752050, India
| | - Arnav Paul
- School of Chemical Sciences, National Institute of Science Education & Research-Bhubaneswar, An OCC of Homi Bhabha National Institute, P.O. Jatni, Khurda, Odisha 752050, India
| | - B L Bhargava
- School of Chemical Sciences, National Institute of Science Education & Research-Bhubaneswar, An OCC of Homi Bhabha National Institute, P.O. Jatni, Khurda, Odisha 752050, India.
| |
Collapse
|
17
|
Tan C, Niitsu A, Sugita Y. Highly Charged Proteins and Their Repulsive Interactions Antagonize Biomolecular Condensation. JACS AU 2023; 3:834-848. [PMID: 37006777 PMCID: PMC10052238 DOI: 10.1021/jacsau.2c00646] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 02/09/2023] [Accepted: 02/09/2023] [Indexed: 06/19/2023]
Abstract
Biomolecular condensation is involved in various cellular processes; therefore, regulation of condensation is crucial to prevent deleterious protein aggregation and maintain a stable cellular environment. Recently, a class of highly charged proteins, known as heat-resistant obscure (Hero) proteins, was shown to protect other client proteins from pathological aggregation. However, the molecular mechanisms by which Hero proteins protect other proteins from aggregation remain unknown. In this study, we performed multiscale molecular dynamics (MD) simulations of Hero11, a Hero protein, and the C-terminal low-complexity domain (LCD) of the transactive response DNA-binding protein 43 (TDP-43), a client protein of Hero11, under various conditions to examine their interactions with each other. We found that Hero11 permeates into the condensate formed by the LCD of TDP-43 (TDP-43-LCD) and induces changes in conformation, intermolecular interactions, and dynamics of TDP-43-LCD. We also examined possible Hero11 structures in atomistic and coarse-grained MD simulations and found that Hero11 with a higher fraction of disordered region tends to assemble on the surface of the condensates. Based on the simulation results, we have proposed three possible mechanisms for Hero11's regulatory function: (i) In the dense phase, TDP-43-LCD reduces contact with each other and shows faster diffusion and decondensation due to the repulsive Hero11-Hero11 interactions. (ii) In the dilute phase, the saturation concentration of TDP-43-LCD is increased, and its conformation is relatively more extended and variant, induced by the attractive Hero11-TDP-43-LCD interactions. (iii) Hero11 on the surface of small TDP-43-LCD condensates can contribute to avoiding their fusion due to repulsive interactions. The proposed mechanisms provide new insights into the regulation of biomolecular condensation in cells under various conditions.
Collapse
Affiliation(s)
- Cheng Tan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Ai Niitsu
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
18
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python package for the analysis of simulations of intrinsically disordered proteins. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.16.528879. [PMID: 36824878 PMCID: PMC9949127 DOI: 10.1101/2023.02.16.528879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite SOURSOP, an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
|
19
|
Khatua P, Mondal S, Gupta M, Bandyopadhyay S. In Silico Studies to Predict the Role of Solvent in Guiding the Conformations of Intrinsically Disordered Peptides and Their Aggregated Protofilaments. ACS OMEGA 2022; 7:43337-43345. [PMID: 36506131 PMCID: PMC9730305 DOI: 10.1021/acsomega.2c06235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Accepted: 11/08/2022] [Indexed: 06/17/2023]
Abstract
The formation of amyloids due to the self-assembly of intrinsically disordered proteins or peptides is a hallmark for different neurodegenerative diseases. For example, amyloids formed by the amyloid beta (Aβ) peptides are responsible for the most devastating neuropathological disease, namely, Alzheimer's disease, while aggregation of α-synuclein peptides causes the etiology of another neuropathological disease, Parkinson's disease. Characterization of the intermediates and the final amyloid formed during the aggregation process is, therefore, crucial for microscopic understanding of the origin behind such diseases, as well as for the development of proper therapeutics to combat those. However, most of the research activities reported in this area have been directed toward examining the early stages of the aggregation process, including probing the conformational characteristics of the responsible protein/peptide in the monomeric state or in small oligomeric forms. This is because the small soluble oligomers have been found to be more deleterious than the final insoluble amyloids. This review discusses some of the recent findings obtained from our simulation studies on Aβ and α-synuclein monomers and small preformed Aβ aggregates. A molecular-level insight of the aggregation process with a special emphasis on the role of water in inducing the aggregation process has been provided.
Collapse
Affiliation(s)
- Prabir Khatua
- Molecular
Modeling Laboratory, Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Souvik Mondal
- Molecular
Modeling Laboratory, Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Madhulika Gupta
- Department
of Chemistry and Chemical Biology, Indian
Institute of Technology (Indian School of Mines) Dhanbad, Jharkhand 826004, India
| | - Sanjoy Bandyopadhyay
- Molecular
Modeling Laboratory, Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
20
|
McIvor JAP, Larsen DS, Mercadante D. Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM. J Phys Chem B 2022; 126:7833-7846. [PMID: 36125334 DOI: 10.1021/acs.jpcb.2c03837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the αR region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.
Collapse
Affiliation(s)
- Jordan A P McIvor
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Danaé S Larsen
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| |
Collapse
|
21
|
Matsubara D, Kasahara K, Dokainish HM, Oshima H, Sugita Y. Modified Protein-Water Interactions in CHARMM36m for Thermodynamics and Kinetics of Proteins in Dilute and Crowded Solutions. Molecules 2022; 27:molecules27175726. [PMID: 36080494 PMCID: PMC9457699 DOI: 10.3390/molecules27175726] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 08/30/2022] [Accepted: 08/30/2022] [Indexed: 11/16/2022] Open
Abstract
Proper balance between protein-protein and protein-water interactions is vital for atomistic molecular dynamics (MD) simulations of globular proteins as well as intrinsically disordered proteins (IDPs). The overestimation of protein-protein interactions tends to make IDPs more compact than those in experiments. Likewise, multiple proteins in crowded solutions are aggregated with each other too strongly. To optimize the balance, Lennard-Jones (LJ) interactions between protein and water are often increased about 10% (with a scaling parameter, λ = 1.1) from the existing force fields. Here, we explore the optimal scaling parameter of protein-water LJ interactions for CHARMM36m in conjunction with the modified TIP3P water model, by performing enhanced sampling MD simulations of several peptides in dilute solutions and conventional MD simulations of globular proteins in dilute and crowded solutions. In our simulations, 10% increase of protein-water LJ interaction for the CHARMM36m cannot maintain stability of a small helical peptide, (AAQAA)3 in a dilute solution and only a small modification of protein-water LJ interaction up to the 3% increase (λ = 1.03) is allowed. The modified protein-water interactions are applicable to other peptides and globular proteins in dilute solutions without changing thermodynamic properties from the original CHARMM36m. However, it has a great impact on the diffusive properties of proteins in crowded solutions, avoiding the formation of too sticky protein-protein interactions.
Collapse
Affiliation(s)
- Daiki Matsubara
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
| | - Kento Kasahara
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Osaka, Japan
| | - Hisham M. Dokainish
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako 351-0198, Saitama, Japan
| | - Hiraku Oshima
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe 650-0047, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako 351-0198, Saitama, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe 650-0047, Hyogo, Japan
- Correspondence: ; Tel.: +81-48-462-1407
| |
Collapse
|
22
|
Li L, Casalini T, Arosio P, Salvalaglio M. Modeling the Structure and Interactions of Intrinsically Disordered Peptides with Multiple Replica, Metadynamics-Based Sampling Methods and Force-Field Combinations. J Chem Theory Comput 2022; 18:1915-1928. [PMID: 35174713 PMCID: PMC9097291 DOI: 10.1021/acs.jctc.1c00889] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Indexed: 01/08/2023]
Abstract
Intrinsically disordered proteins play a key role in many biological processes, including the formation of biomolecular condensates within cells. A detailed characterization of their configurational ensemble and structure-function paradigm is crucial for understanding their biological activity and for exploiting them as building blocks in material sciences. In this work, we incorporate bias-exchange metadynamics and parallel-tempering well-tempered metadynamics with CHARMM36m and CHARMM22* to explore the structural and thermodynamic characteristics of a short archetypal disordered sequence derived from a DEAD-box protein. The conformational landscapes emerging from our simulations are largely congruent across methods and force fields. Nevertheless, differences in fine details emerge from varying combinations of force-fields and sampling methods. For this protein, our analysis identifies features that help to explain the low propensity of this sequence to undergo self-association in vitro, which are common to all force-field/sampling method combinations. Overall, our work demonstrates the importance of using multiple force-field and sampling method combinations for accurate structural and thermodynamic information in the study of disordered proteins.
Collapse
Affiliation(s)
- Lunna Li
- Thomas
Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, U.K.
| | - Tommaso Casalini
- Department
of Chemistry and Applied Biosciences, ETH
Zurich, Zurich 8093, Switzerland
| | - Paolo Arosio
- Department
of Chemistry and Applied Biosciences, ETH
Zurich, Zurich 8093, Switzerland
| | - Matteo Salvalaglio
- Thomas
Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, U.K.
| |
Collapse
|
23
|
Conformational ensembles of intrinsically disordered proteins and flexible multidomain proteins. Biochem Soc Trans 2022; 50:541-554. [PMID: 35129612 DOI: 10.1042/bst20210499] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 01/13/2022] [Accepted: 01/17/2022] [Indexed: 12/29/2022]
Abstract
Intrinsically disordered proteins (IDPs) and multidomain proteins with flexible linkers show a high level of structural heterogeneity and are best described by ensembles consisting of multiple conformations with associated thermodynamic weights. Determining conformational ensembles usually involves the integration of biophysical experiments and computational models. In this review, we discuss current approaches to determine conformational ensembles of IDPs and multidomain proteins, including the choice of biophysical experiments, computational models used to sample protein conformations, models to calculate experimental observables from protein structure, and methods to refine ensembles against experimental data. We also provide examples of recent applications of integrative conformational ensemble determination to study IDPs and multidomain proteins and suggest future directions for research in the field.
Collapse
|
24
|
Murthy AC, Tang WS, Jovic N, Janke AM, Seo DH, Perdikari TM, Mittal J, Fawzi NL. Molecular interactions contributing to FUS SYGQ LC-RGG phase separation and co-partitioning with RNA polymerase II heptads. Nat Struct Mol Biol 2021; 28:923-935. [PMID: 34759379 PMCID: PMC8654040 DOI: 10.1038/s41594-021-00677-4] [Citation(s) in RCA: 57] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Accepted: 09/28/2021] [Indexed: 01/23/2023]
Abstract
The RNA-binding protein FUS (Fused in Sarcoma) mediates phase separation in biomolecular condensates and functions in transcription by clustering with RNA polymerase II. Specific contact residues and interaction modes formed by FUS and the C-terminal heptad repeats of RNA polymerase II (CTD) have been suggested but not probed directly. Here we show how RGG domains contribute to phase separation with the FUS N-terminal low-complexity domain (SYGQ LC) and RNA polymerase II CTD. Using NMR spectroscopy and molecular simulations, we demonstrate that many residue types, not solely arginine-tyrosine pairs, form condensed-phase contacts via several interaction modes including, but not only sp2-π and cation-π interactions. In phases also containing RNA polymerase II CTD, many residue types form contacts, including both cation-π and hydrogen-bonding interactions formed by the conserved human CTD lysines. Hence, our data suggest a surprisingly broad array of residue types and modes explain co-phase separation of FUS and RNA polymerase II.
Collapse
Affiliation(s)
- Anastasia C Murthy
- Graduate Program in Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI, USA
| | - Wai Shing Tang
- Graduate Program in Physics, Brown University, Providence, RI, USA
| | - Nina Jovic
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX, USA
| | - Abigail M Janke
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI, USA
| | - Da Hee Seo
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI, USA
| | | | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX, USA.
| | - Nicolas L Fawzi
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI, USA.
| |
Collapse
|
25
|
Jephthah S, Pesce F, Lindorff-Larsen K, Skepö M. Force Field Effects in Simulations of Flexible Peptides with Varying Polyproline II Propensity. J Chem Theory Comput 2021; 17:6634-6646. [PMID: 34524800 PMCID: PMC8515809 DOI: 10.1021/acs.jctc.1c00408] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Five peptides previously suggested to possess polyproline II (PPII) structure have here been investigated by using atomistic molecular dynamics simulations to compare how well four different force fields known for simulating intrinsically disordered proteins relatively well (Amber ff99SB-disp, Amber ff99SB-ILDN, CHARM36IDPSFF, and CHARMM36m) can capture this secondary structure element. The results revealed that all force fields sample PPII structures but to different extents and with different propensities toward other secondary structure elements, in particular, the β-sheet and "random coils". A cluster analysis of the simulations of histatin 5 also revealed that the conformational ensembles of the force fields are quite different. We compared the simulations to circular dichroism and nuclear magnetic resonance spectroscopy experiments and conclude that further experiments and methods for interpreting them are needed to assess the accuracy of force fields in determining PPII structure.
Collapse
Affiliation(s)
- Stéphanie Jephthah
- Division of Theoretical Chemistry, Lund University, SE-221 00 Lund, Sweden
| | - Francesco Pesce
- Structural Biology and NMR Laboratory & the Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory & the Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Marie Skepö
- Division of Theoretical Chemistry, Lund University, SE-221 00 Lund, Sweden
| |
Collapse
|
26
|
Alston JJ, Soranno A, Holehouse AS. Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins. Methods 2021; 193:116-135. [PMID: 33831596 PMCID: PMC8713295 DOI: 10.1016/j.ymeth.2021.03.018] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Revised: 03/25/2021] [Accepted: 03/31/2021] [Indexed: 12/21/2022] Open
Abstract
Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.
Collapse
Affiliation(s)
- Jhullian J Alston
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA
| | - Andrea Soranno
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA.
| | - Alex S Holehouse
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis 63110, MO, USA; Center for Science and Engineering of Living Systems (CSELS), Washington University in St. Louis, St. Louis 63130, MO, USA.
| |
Collapse
|
27
|
Lindorff-Larsen K, Kragelund BB. On the potential of machine learning to examine the relationship between sequence, structure, dynamics and function of intrinsically disordered proteins. J Mol Biol 2021; 433:167196. [PMID: 34390736 DOI: 10.1016/j.jmb.2021.167196] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Revised: 08/03/2021] [Accepted: 08/04/2021] [Indexed: 11/29/2022]
Abstract
Intrinsically disordered proteins (IDPs) constitute a broad set of proteins with few uniting and many diverging properties. IDPs-and intrinsically disordered regions (IDRs) interspersed between folded domains-are generally characterized as having no persistent tertiary structure; instead they interconvert between a large number of different and often expanded structures. IDPs and IDRs are involved in an enormously wide range of biological functions and reveal novel mechanisms of interactions, and while they defy the common structure-function paradigm of folded proteins, their structural preferences and dynamics are important for their function. We here discuss open questions in the field of IDPs and IDRs, focusing on areas where machine learning and other computational methods play a role. We discuss computational methods aimed to predict transiently formed local and long-range structure, including methods for integrative structural biology. We discuss the many different ways in which IDPs and IDRs can bind to other molecules, both via short linear motifs, as well as in the formation of larger dynamic complexes such as biomolecular condensates. We discuss how experiments are providing insight into such complexes and may enable more accurate predictions. Finally, we discuss the role of IDPs in disease and how new methods are needed to interpret the mechanistic effects of genomic variants in IDPs.
Collapse
Affiliation(s)
- Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen. Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark.
| | - Birthe B Kragelund
- Structural Biology and NMR Laboratory & Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen. Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark.
| |
Collapse
|
28
|
Kauffmann C, Ceccolini I, Kontaxis G, Konrat R. Detecting anisotropic segmental dynamics in disordered proteins by cross-correlated spin relaxation. MAGNETIC RESONANCE (GOTTINGEN, GERMANY) 2021; 2:557-569. [PMID: 37905226 PMCID: PMC10539831 DOI: 10.5194/mr-2-557-2021] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/23/2021] [Accepted: 06/02/2021] [Indexed: 11/01/2023]
Abstract
Among the numerous contributions of Geoffrey Bodenhausen to NMR spectroscopy, his developments in the field of spin-relaxation methodology and theory will definitely have a long lasting impact. Starting with his seminal contributions to the excitation of multiple-quantum coherences, he and his group thoroughly investigated the intricate relaxation properties of these "forbidden fruits" and developed experimental techniques to reveal the relevance of previously largely ignored cross-correlated relaxation (CCR) effects, as "the essential is invisible to the eyes". Here we consider CCR within the challenging context of intrinsically disordered proteins (IDPs) and emphasize its potential and relevance for the studies of structural dynamics of IDPs in the future years to come. Conventionally, dynamics of globularly folded proteins are modeled and understood as deviations from otherwise rigid structures tumbling in solution. However, with increasing protein flexibility, as observed for IDPs, this apparent dichotomy between structure and dynamics becomes blurred. Although complex dynamics and ensemble averaging might impair the extraction of mechanistic details even further, spin relaxation uniquely encodes a protein's structural memory. Due to significant methodological developments, such as high-dimensional non-uniform sampling techniques, spin relaxation in IDPs can now be monitored in unprecedented resolution. Not embedded within a rigid globular fold, conventional 15 N spin probes might not suffice to capture the inherently local nature of IDP dynamics. To better describe and understand possible segmental motions of IDPs, we propose an experimental approach to detect the signature of anisotropic segmental dynamics by quantifying cross-correlated spin relaxation of individual 15 N 1 H N and 13 C ' 13 C α spin pairs. By adapting Geoffrey Bodenhausen's symmetrical reconversion principle to obtain zero frequency spectral density values, we can define and demonstrate more sensitive means to characterize anisotropic dynamics in IDPs.
Collapse
Affiliation(s)
- Clemens Kauffmann
- Department of Structural and Computational Biology, Max Perutz Laboratories, University of Vienna, Campus-Vienna-Biocenter 5, 1030 Vienna, Austria
| | - Irene Ceccolini
- Department of Structural and Computational Biology, Max Perutz Laboratories, University of Vienna, Campus-Vienna-Biocenter 5, 1030 Vienna, Austria
| | - Georg Kontaxis
- Department of Structural and Computational Biology, Max Perutz Laboratories, University of Vienna, Campus-Vienna-Biocenter 5, 1030 Vienna, Austria
| | - Robert Konrat
- Department of Structural and Computational Biology, Max Perutz Laboratories, University of Vienna, Campus-Vienna-Biocenter 5, 1030 Vienna, Austria
| |
Collapse
|
29
|
Feng JJ, Chen JN, Kang W, Wu YD. Accurate Structure Prediction for Protein Loops Based on Molecular Dynamics Simulations with RSFF2C. J Chem Theory Comput 2021; 17:4614-4628. [PMID: 34170125 DOI: 10.1021/acs.jctc.1c00341] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein loops, connecting the α-helices and β-strands, are involved in many important biological processes. However, due to their conformational flexibility, it is still challenging to accurately determine three-dimensional (3D) structures of long loops experimentally and computationally. Herein, we present a systematic study of the protein loop structure prediction via a total of ∼850 μs molecular dynamics (MD) simulations. For a set of 15 long (10-16 residues) and solvent-exposed loops, we first evaluated the performance of four state-of-the-art loop modeling algorithms, DaReUS-Loop, Sphinx, Rosetta-NGK, and MODELLER, on each loop, and none of them could accurately predict the structures for most loops. Then, temperature replica exchange molecular dynamics (REMD) simulations were conducted with three recent force fields, RSFF2C with TIP3P water model, CHARMM36m with CHARMM-modified TIP3P, and AMBER ff19SB with OPC. We found that our recently developed residue-specific force field RSFF2C performed the best and successfully predicted 12 out of 15 loops with a root-mean-square deviation (RMSD) < 1.5 Å. As an alternative with lower computational cost, normal MD simulations at high temperatures (380, 500, and 620 K) were investigated. Temperature-dependent performance was observed for each force field, and, for RSFF2C+TIP3P, we found that three independent 100-ns MD simulations at 500 K gave comparable results with REMD simulations. These results suggest that MD simulations, especially with enhanced sampling techniques such as replica exchange, with the RSFF2C force field could be useful for accurate loop structure prediction.
Collapse
Affiliation(s)
- Jia-Jie Feng
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Jia-Nan Chen
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Wei Kang
- Pingshan Translational Medicine Center, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Yun-Dong Wu
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China.,College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Shenzhen Bay Laboratory, Shenzhen 518132, China
| |
Collapse
|
30
|
Ahmed MC, Skaanning LK, Jussupow A, Newcombe EA, Kragelund BB, Camilloni C, Langkilde AE, Lindorff-Larsen K. Refinement of α-Synuclein Ensembles Against SAXS Data: Comparison of Force Fields and Methods. Front Mol Biosci 2021; 8:654333. [PMID: 33968988 PMCID: PMC8100456 DOI: 10.3389/fmolb.2021.654333] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 03/12/2021] [Indexed: 12/22/2022] Open
Abstract
The inherent flexibility of intrinsically disordered proteins (IDPs) makes it difficult to interpret experimental data using structural models. On the other hand, molecular dynamics simulations of IDPs often suffer from force-field inaccuracies, and long simulation times or enhanced sampling methods are needed to obtain converged ensembles. Here, we apply metainference and Bayesian/Maximum Entropy reweighting approaches to integrate prior knowledge of the system with experimental data, while also dealing with various sources of errors and the inherent conformational heterogeneity of IDPs. We have measured new SAXS data on the protein α-synuclein, and integrate this with simulations performed using different force fields. We find that if the force field gives rise to ensembles that are much more compact than what is implied by the SAXS data it is difficult to recover a reasonable ensemble. On the other hand, we show that when the simulated ensemble is reasonable, we can obtain an ensemble that is consistent with the SAXS data, but also with NMR diffusion and paramagnetic relaxation enhancement data.
Collapse
Affiliation(s)
- Mustapha Carab Ahmed
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark
| | - Line K Skaanning
- Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Alexander Jussupow
- Department of Chemistry, Institute for Advanced Study, Technical University of Munich, Munich, Germany
| | - Estella A Newcombe
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark.,Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Birthe B Kragelund
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark
| | - Carlo Camilloni
- Department of Chemistry, Institute for Advanced Study, Technical University of Munich, Munich, Germany.,Dipartimento di Bioscienze, Università degli Studi di Milano, Milan, Italy
| | - Annette E Langkilde
- Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
31
|
Ploetz EA, Karunaweera S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Applications of KBFF20. J Chem Theory Comput 2021; 17:2991-3009. [PMID: 33878264 DOI: 10.1021/acs.jctc.1c00076] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native β-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
32
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
33
|
Shea JE, Best RB, Mittal J. Physics-based computational and theoretical approaches to intrinsically disordered proteins. Curr Opin Struct Biol 2021; 67:219-225. [PMID: 33545530 PMCID: PMC8150118 DOI: 10.1016/j.sbi.2020.12.012] [Citation(s) in RCA: 73] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 12/24/2020] [Accepted: 12/28/2020] [Indexed: 02/06/2023]
Abstract
Intrinsically disordered proteins (IDPs) are an important class of proteins that do not fold to a well-defined three-dimensional shape but rather adopt an ensemble of inter-converting conformations. This feature makes their experimental characterization challenging and invites a theoretical and computational approach to complement experimental studies. In this review, we highlight the recent progress in developing new computational and theoretical approaches to study the structure and dynamics of monomeric and order higher assemblies of IDPs, with a particular emphasis on their phase separation into protein-rich condensates.
Collapse
Affiliation(s)
- Joan-Emma Shea
- Department of Chemistry & Biochemistry, University of California, Santa Barbara, CA 93106, United States; Department of Physics, University of California, Santa Barbara, CA 93106, United States.
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, 111 Research Drive, Bethlehem, PA 18015, United States.
| |
Collapse
|
34
|
Dong X, Bera S, Qiao Q, Tang Y, Lao Z, Luo Y, Gazit E, Wei G. Liquid-Liquid Phase Separation of Tau Protein Is Encoded at the Monomeric Level. J Phys Chem Lett 2021; 12:2576-2586. [PMID: 33686854 DOI: 10.1021/acs.jpclett.1c00208] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Liquid-liquid phase separation (LLPS) is involved in both physiological and pathological processes. The intrinsically disordered protein Tau and its K18 construct can undergo LLPS in a distinct temperature-dependent manner, and the LLPS of Tau protein can initiate Tau aggregation. However, the underlying mechanism driving Tau LLPS remains largely elusive. To understand the temperature-dependent LLPS behavior of Tau at the monomeric level, we explored the conformational ensemble of Tau at different temperatures by performing all-atom replica-exchange molecular dynamic simulation on K18 monomer with an accumulated simulation time of 26.4 μs. Our simulation demonstrates that the compactness, β-structure propensity, and intramolecular interaction of K18 monomer exhibit nonlinear temperature-dependent behavior. 295DNIKHV300/326GNIHHK331/337VEVKSE342 make significant contributions to the temperature dependence of the β propensity of K18 monomer, while the two fibril-nucleating cores display relatively high β propensity at all temperatures. At a specific temperature, K18 monomer adopts the most collapsed state with exposed sites for both persistent and transient interactions. Given that more collapsed polypeptide chains were reported to be more prone to phase separate, our results suggest that K18 monomer inherently possesses conformational characteristics favoring LLPS. Our simulation predicts the importance of 295DNIKHV300/326GNIHHK331/337VEVKSE342 to the temperature-dependent conformational properties of K18, which is corroborated by CD spectra, turbidity assays, and DIC microscopy. Taken together, we offer a computational and experimental approach to comprehend the structural basis for LLPS by amyloidal building blocks.
Collapse
Affiliation(s)
- Xuewei Dong
- Department of Physics, State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200433, People's Republic of China
| | - Santu Bera
- Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Qin Qiao
- Digital Medical Research Center, School of Basic Medical Sciences, Shanghai Key Laboratory of Medical Imaging Computing and Computer Assisted Intervention, Fudan University, Shanghai 200032, People's Republic of China
| | - Yiming Tang
- Department of Physics, State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200433, People's Republic of China
| | - Zenghui Lao
- Department of Physics, State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200433, People's Republic of China
| | - Yin Luo
- Department of Physics, State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200433, People's Republic of China
| | - Ehud Gazit
- Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Guanghong Wei
- Department of Physics, State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200433, People's Republic of China
| |
Collapse
|
35
|
Shrestha UR, Smith JC, Petridis L. Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations. Commun Biol 2021; 4:243. [PMID: 33623120 PMCID: PMC7902620 DOI: 10.1038/s42003-021-01759-1] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2020] [Accepted: 01/07/2021] [Indexed: 12/13/2022] Open
Abstract
Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.
Collapse
Affiliation(s)
- Utsab R Shrestha
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, TN, USA
| | - Jeremy C Smith
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, TN, USA
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN, USA
| | - Loukas Petridis
- Oak Ridge National Laboratory, Biosciences Division, UT/ORNL Center for Molecular Biophysics, Oak Ridge, TN, USA.
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN, USA.
| |
Collapse
|
36
|
Perdikari TM, Jovic N, Dignon GL, Kim YC, Fawzi NL, Mittal J. A predictive coarse-grained model for position-specific effects of post-translational modifications. Biophys J 2021; 120:1187-1197. [PMID: 33582133 DOI: 10.1016/j.bpj.2021.01.034] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Revised: 12/29/2020] [Accepted: 01/19/2021] [Indexed: 02/06/2023] Open
Abstract
Biomolecules undergo liquid-liquid phase separation (LLPS), resulting in the formation of multicomponent protein-RNA membraneless organelles in cells. However, the physiological and pathological role of post-translational modifications (PTMs) on the biophysics of phase behavior is only beginning to be probed. To study the effect of PTMs on LLPS in silico, we extend our transferable coarse-grained model of intrinsically disordered proteins to include phosphorylated and acetylated amino acids. Using the parameters for modified amino acids available for fixed-charge atomistic force fields, we parameterize the size and atomistic hydropathy of the coarse-grained-modified amino acid beads and, hence, the interactions between the modified and natural amino acids. We then elucidate how the number and position of phosphorylated and acetylated residues alter the protein's single-chain compactness and its propensity to phase separate. We show that both the number and the position of phosphorylated threonines/serines or acetylated lysines can serve as a molecular on/off switch for phase separation in the well-studied disordered regions of Fused in Sarcoma (FUS) and DDX3X, respectively. We also compare modified residues to their commonly used PTM mimics for their impact on chain properties. Importantly, we show that the model can predict and capture experimentally measured differences in the phase behavior for position-specific modifications, showing that the position of modifications can dictate phase separation. In sum, this model will be useful for studying LLPS of post-translationally modified intrinsically disordered proteins and predicting how modifications control phase behavior with position-specific resolution.
Collapse
Affiliation(s)
| | - Nina Jovic
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania
| | - Gregory L Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania; Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York
| | - Young C Kim
- Center for Materials Physics and Technology, U.S. Naval Research Laboratory, Washington D.C
| | - Nicolas L Fawzi
- Department of Molecular Pharmacology, Physiology, and Biotechnology; Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island.
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania.
| |
Collapse
|
37
|
Wang W. Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins. Phys Chem Chem Phys 2021; 23:777-784. [PMID: 33355572 DOI: 10.1039/d0cp05818a] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.
Collapse
Affiliation(s)
- Wenning Wang
- Department of Chemistry, Multiscale Research Institute of Complex Systems and Institute of Biomedical Sciences, Fudan University, Shanghai 200438, China.
| |
Collapse
|
38
|
Ramanathan A, Ma H, Parvatikar A, Chennubhotla SC. Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins. Curr Opin Struct Biol 2021; 66:216-224. [PMID: 33421906 DOI: 10.1016/j.sbi.2020.12.001] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 12/01/2020] [Accepted: 12/03/2020] [Indexed: 12/16/2022]
Abstract
We outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships-however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.
Collapse
Affiliation(s)
- Arvind Ramanathan
- Data Science & Learning Division, Argonne National Laboratory, Lemont, IL 60439, United States; Consortium for Advanced Science and Engineering (CASE), University of Chicago, Hyde Park, IL, United States.
| | - Heng Ma
- Data Science & Learning Division, Argonne National Laboratory, Lemont, IL 60439, United States
| | - Akash Parvatikar
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15260, United States
| | - S Chakra Chennubhotla
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15260, United States
| |
Collapse
|
39
|
Zheng W, Dignon GL, Jovic N, Xu X, Regy RM, Fawzi NL, Kim YC, Best RB, Mittal J. Molecular Details of Protein Condensates Probed by Microsecond Long Atomistic Simulations. J Phys Chem B 2020; 124:11671-11679. [PMID: 33302617 PMCID: PMC7879053 DOI: 10.1021/acs.jpcb.0c10489] [Citation(s) in RCA: 102] [Impact Index Per Article: 25.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The formation of membraneless organelles in cells commonly occurs via liquid-liquid phase separation (LLPS) and is in many cases driven by multivalent interactions between intrinsically disordered proteins (IDPs). Investigating the nature of these interactions, and their effect on dynamics within the condensed phase, is therefore of critical importance but very challenging for either simulation or experiment. Here, we study these interactions and their dynamics by pairing a novel multiscale simulation strategy with microsecond all-atom MD simulations of a condensed, IDP-rich phase. We simulate two IDPs this way, the low complexity domain of FUS and the N-terminal disordered domain of LAF-1, and find good agreement with experimental information about average density, water content, and residue-residue contacts. We go significantly beyond what is known from experiments by showing that ion partitioning within the condensed phase is largely driven by the charge distribution of the proteins and-in the cases considered-shows little evidence of preferential interactions of the ions with the proteins. Furthermore, we can probe the microscopic diffusive dynamics within the condensed phase, showing that water and ions are in dynamic equilibrium between dense and dilute phases, and their diffusion is reduced in the dense phase. Despite their high concentration in the condensate, the protein molecules also remain mobile, explaining the observed liquid-like properties of this phase. We finally show that IDP self-association is driven by a combination of nonspecific hydrophobic interactions as well as hydrogen bonds, salt bridges, and π-π and cation-π interactions. The simulation approach presented here allows the structural and dynamical properties of biomolecular condensates to be studied in microscopic detail and is generally applicable to single- and multicomponent systems of proteins and nucleic acids involved in LLPS.
Collapse
Affiliation(s)
- Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, Arizona 85281, United States
| | - Gregory L Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Nina Jovic
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Xichen Xu
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Roshan M Regy
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Nicolas L Fawzi
- Department of Molecular Pharmacology, Physiology and Biotechnology, Brown University, Providence, Rhode Island 02912, United States
| | - Young C Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, D.C. 20375, United States
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
40
|
Best RB. A Tale of Two Tyrosines. Biophys J 2020; 119:1927-1928. [PMID: 33120016 DOI: 10.1016/j.bpj.2020.09.036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Accepted: 09/29/2020] [Indexed: 10/23/2022] Open
Affiliation(s)
- Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland.
| |
Collapse
|
41
|
Development of Charge-Augmented Three-Point Water Model (CAIPi3P) for Accurate Simulations of Intrinsically Disordered Proteins. Int J Mol Sci 2020; 21:ijms21176166. [PMID: 32859072 PMCID: PMC7504337 DOI: 10.3390/ijms21176166] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 08/17/2020] [Accepted: 08/21/2020] [Indexed: 01/15/2023] Open
Abstract
Intrinsically disordered proteins (IDPs) are molecules without a fixed tertiary structure, exerting crucial roles in cellular signalling, growth and molecular recognition events. Due to their high plasticity, IDPs are very challenging in experimental and computational structural studies. To provide detailed atomic insight in IDPs' dynamics governing their functional mechanisms, all-atom molecular dynamics (MD) simulations are widely employed. However, the current generalist force fields and solvent models are unable to generate satisfactory ensembles for IDPs when compared to existing experimental data. In this work, we present a new solvation model, denoted as the Charge-Augmented Three-Point Water Model for Intrinsically Disordered Proteins (CAIPi3P). CAIPi3P has been generated by performing a systematic scan of atomic partial charges assigned to the widely popular molecular scaffold of the three-point TIP3P water model. We found that explicit solvent MD simulations employing CAIPi3P solvation considerably improved the small-angle X-ray scattering (SAXS) scattering profiles for three different IDPs. Not surprisingly, this improvement was further enhanced by using CAIPi3P water in combination with the protein force field parametrized for IDPs. We also demonstrated the applicability of CAIPi3P to molecular systems containing structured as well as intrinsically disordered regions/domains. Our results highlight the crucial importance of solvent effects for generating molecular ensembles of IDPs which reproduce the experimental data available. Hence, we conclude that our newly developed CAIPi3P solvation model is a valuable tool for molecular simulations of intrinsically disordered proteins and assessing their molecular dynamics.
Collapse
|
42
|
Jafarinia H, van der Giessen E, Onck PR. Phase Separation of Toxic Dipeptide Repeat Proteins Related to C9orf72 ALS/FTD. Biophys J 2020; 119:843-851. [PMID: 32730793 DOI: 10.1016/j.bpj.2020.07.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Revised: 06/28/2020] [Accepted: 07/06/2020] [Indexed: 10/23/2022] Open
Abstract
The expansion mutation in the C9orf72 gene is the most common known genetic cause for amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). This mutation can produce five dipeptide repeat proteins (DPRs), of which three are known to be toxic: poly-PR, poly-GR, and poly-GA. The toxicity of poly-GA is attributed to its aggregation in the cytoplasm, whereas for poly-PR and poly-GR, several toxicity pathways have been proposed. The toxicity of the DPRs has been shown to depend on their length, but the underlying molecular mechanism of this length dependence is not well understood. To address the possible role of phase separation in DPR toxicity, a one-bead-per-amino-acid (1BPA) coarse-grained molecular dynamics model is used to study the single-molecule and phase-separation properties of the DPRs. We find a strong dependence of the phase-separation behavior on both DPR length and concentration, with longer DPRs having a higher propensity to phase separate and form condensed phases with higher concentrations. The critical lengths required for phase separation (25 for poly-PR and 50 for poly-GA) are comparable to the toxicity threshold limit of 30 repeats found for the expansion mutation in patient cells, suggesting that phase separation could play an important role in DPR toxicity.
Collapse
Affiliation(s)
- Hamidreza Jafarinia
- Zernike Institute for Advanced Materials, University of Groningen, Groningen, the Netherlands
| | - Erik van der Giessen
- Zernike Institute for Advanced Materials, University of Groningen, Groningen, the Netherlands
| | - Patrick R Onck
- Zernike Institute for Advanced Materials, University of Groningen, Groningen, the Netherlands.
| |
Collapse
|
43
|
Cohan MC, Ruff KM, Pappu RV. Information theoretic measures for quantifying sequence-ensemble relationships of intrinsically disordered proteins. Protein Eng Des Sel 2020; 32:191-202. [PMID: 31375817 PMCID: PMC7462041 DOI: 10.1093/protein/gzz014] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 06/19/2019] [Indexed: 01/26/2023] Open
Abstract
Intrinsically disordered proteins (IDPs) contribute to a multitude of functions. De novo design of IDPs should open the door to modulating functions and phenotypes controlled by these systems. Recent design efforts have focused on compositional biases and specific sequence patterns as the design features. Analysis of the impact of these designs on sequence-function relationships indicates that individual sequence/compositional parameters are insufficient for describing sequence-function relationships in IDPs. To remedy this problem, we have developed information theoretic measures for sequence–ensemble relationships (SERs) of IDPs. These measures rely on prior availability of statistically robust conformational ensembles derived from all atom simulations. We show that the measures we have developed are useful for comparing sequence-ensemble relationships even when sequence is poorly conserved. Based on our results, we propose that de novo designs of IDPs, guided by knowledge of their SERs, should provide improved insights into their sequence–ensemble–function relationships.
Collapse
Affiliation(s)
- Megan C Cohan
- Department of Biomedical Engineering and Center for Science & Engineering of Living Systems (CSELS) Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis MO, USA
| | - Kiersten M Ruff
- Department of Biomedical Engineering and Center for Science & Engineering of Living Systems (CSELS) Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis MO, USA
| | - Rohit V Pappu
- Department of Biomedical Engineering and Center for Science & Engineering of Living Systems (CSELS) Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis MO, USA
| |
Collapse
|
44
|
Zheng W, Dignon G, Brown M, Kim YC, Mittal J. Hydropathy Patterning Complements Charge Patterning to Describe Conformational Preferences of Disordered Proteins. J Phys Chem Lett 2020; 11:3408-3415. [PMID: 32227994 PMCID: PMC7450210 DOI: 10.1021/acs.jpclett.0c00288] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Understanding the conformational ensemble of an intrinsically disordered protein (IDP) is of great interest due to its relevance to critical intracellular functions and diseases. It is now well established that the polymer scaling behavior can provide a great deal of information about the conformational properties as well as liquid-liquid phase separation of an IDP. It is, therefore, extremely desirable to be able to predict an IDP's scaling behavior from the protein sequence itself. The work in this direction so far has focused on highly charged proteins and how charge patterning can perturb their structural properties. As naturally occurring IDPs are composed of a significant fraction of uncharged amino acids, the rules based on charge content and patterning are only partially helpful in solving the problem. Here, we propose a new order parameter, sequence hydropathy decoration, which can provide a near-quantitative understanding of scaling and structural properties of IDPs devoid of charged residues. We combine this with a charge patterning parameter, sequence charge decoration, to obtain a general equation, parametrized from extensive coarse-grained simulation data, for predicting protein dimensions from the sequence. We finally test this equation against available experimental data and find a semiquantitative match in predicting the scaling behavior. We also provide guidance on how to extend this approach to experimental data, which should be feasible in the near future.
Collapse
Affiliation(s)
- Wenwei Zheng
- College of Integrative Sciences and Arts, Arizona State University, Mesa, Arizona 85212, United States
| | - Gregory Dignon
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Matthew Brown
- College of Integrative Sciences and Arts, Arizona State University, Mesa, Arizona 85212, United States
| | - Young C Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, D.C. 20375, United States
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
45
|
Horx P, Geyer A. Defining the mobility range of a hinge-type connection using molecular dynamics and metadynamics. PLoS One 2020; 15:e0230962. [PMID: 32282813 PMCID: PMC7153902 DOI: 10.1371/journal.pone.0230962] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Accepted: 03/12/2020] [Indexed: 01/29/2023] Open
Abstract
A designed disulfide-rich β-hairpin peptide that dimerizes spontaneously served as a hinge-type connection between proteins. Here, we analyze the range of dynamics of this hinge dimer with the aim of proposing new applications for the DNA-encodable peptide and establishing guidelines for the computational analysis of other disulfide hinges. A recent structural analysis based on nuclear magnetic resonance spectroscopy and ion mobility spectrometry revealed an averaged conformation in the hinge region which motivated us to investigate the dynamic behavior using a combination of molecular dynamics simulation, metadynamics and free energy surface analysis to characterize the conformational space available to the hinge. Principal component analysis uncovered two slow modes of the peptide, namely, the opening and closing motion and twisting of the two β-hairpins assembling the hinge. Applying a collective variable (CV) that mimics the first dominating mode, led to a major expansion of the conformational space. The description of the dynamics could be achieved by analysis of the opening angle and the twisting of the β-hairpins and, thus, offers a methodology that can also be transferred to other derivatives. It has been demonstrated that the hinge peptide’s lowest energy conformation consists of a large opening angle and strong twist but is separated by small energy barriers and can, thus, adopt a closed and untwisted structure. With the aim of proposing further applications for the hinge peptide, we simulated its behavior in the sterically congested environment of a four-helix bundle. Preliminary investigations show that one helix is pushed out and a three-helix bundle forms. The insights gained into the dynamics of the tetra-disulfide peptide and analytical guidelines developed in this study may contribute to the understanding of the structure and function of more complex hinge-type proteins, such as the IgG antibody family.
Collapse
Affiliation(s)
- Philip Horx
- Department of Chemistry, Philipps-Universität Marburg, Marburg, Germany
| | - Armin Geyer
- Department of Chemistry, Philipps-Universität Marburg, Marburg, Germany
- * E-mail:
| |
Collapse
|
46
|
Batys P, Morga M, Bonarek P, Sammalkorpi M. pH-Induced Changes in Polypeptide Conformation: Force-Field Comparison with Experimental Validation. J Phys Chem B 2020; 124:2961-2972. [PMID: 32182068 PMCID: PMC7590956 DOI: 10.1021/acs.jpcb.0c01475] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Indexed: 12/17/2022]
Abstract
Microsecond-long all-atom molecular dynamics (MD) simulations, circular dichroism, laser Doppler velocimetry, and dynamic light-scattering techniques have been used to investigate pH-induced changes in the secondary structure, charge, and conformation of poly l-lysine (PLL) and poly l-glutamic acid (PGA). The employed combination of the experimental methods reveals for both PLL and PGA a narrow pH range at which they are charged enough to form stable colloidal suspensions, maintaining their α-helix content above 60%; an elevated charge state of the peptides required for colloidal stability promotes the peptide solvation as a random coil. To obtain a more microscopic view on the conformations and to verify the modeling performance, peptide secondary structure and conformations rising in MD simulations are also examined using three different force fields, i.e., OPLS-AA, CHARMM27, and AMBER99SB*-ILDNP. Ramachandran plots reveal that in the examined setup the α-helix content is systematically overestimated in CHARMM27, while OPLS-AA overestimates the β-sheet fraction at lower ionization degrees. At high ionization degrees, the OPLS-AA force-field-predicted secondary structure fractions match the experimentally measured distribution most closely. However, the pH-induced changes in PLL and PGA secondary structure are reasonably captured only by the AMBER99SB*-ILDNP force field, with the exception of the fully charged PGA in which the α-helix content is overestimated. The comparison to simulations results shows that the examined force fields involve significant deviations in their predictions for charged homopolypeptides. The detailed mapping of secondary structure dependency on pH for the polypeptides, especially finding the stable colloidal α-helical regime for both examined peptides, has significant potential for practical applications of the charged homopolypeptides. The findings raise attention especially to the pH fine tuning as an underappreciated control factor in surface modification and self-assembly.
Collapse
Affiliation(s)
- Piotr Batys
- Jerzy
Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, PL-30239 Krakow, Poland
| | - Maria Morga
- Jerzy
Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, PL-30239 Krakow, Poland
| | - Piotr Bonarek
- Department
of Physical Biochemistry, Faculty of Biochemistry, Biophysics and
Biotechnology, Jagiellonian University, Krakow, Poland
| | - Maria Sammalkorpi
- Department of Chemistry and Materials Science and Department of Bioproducts and Biosystems, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| |
Collapse
|
47
|
Rieloff E, Skepö M. Phosphorylation of a Disordered Peptide-Structural Effects and Force Field Inconsistencies. J Chem Theory Comput 2020; 16:1924-1935. [PMID: 32050065 DOI: 10.1021/acs.jctc.9b01190] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Phosphorylation is one of the most abundant types of post-translational modifications of intrinsically disordered proteins (IDPs). This study examines the conformational changes in the 15-residue-long N-terminal fragment of the IDP statherin upon phosphorylation, using computer simulations with two different force fields: AMBER ff99SB-ILDN and CHARMM36m. The results from the simulations are compared with experimental small-angle X-ray scattering (SAXS) and circular dichroism data. In the unphosphorylated state, the two force fields are in excellent agreement regarding global structural properties such as size and shape. However, they exhibit some differences in the extent and type of the secondary structure. In the phosphorylated state, neither of the force fields performs well compared to the experimental data. Both force fields show a compaction of the peptide upon phosphorylation, greater than what is seen in SAXS experiments, although they differ in the local structure. While the CHARMM force field increases the fraction of bends in the peptide as a response to strong interactions between the phosphorylated residues and arginines, the AMBER force field shows an increase of the helical content in the N-terminal part of the peptide, where the phosphorylated residues reside, in better agreement with circular dichroism results.
Collapse
Affiliation(s)
- Ellen Rieloff
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Marie Skepö
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden.,LINXS-Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
48
|
Lemkul JA. Pairwise-additive and polarizable atomistic force fields for molecular dynamics simulations of proteins. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:1-71. [PMID: 32145943 DOI: 10.1016/bs.pmbts.2019.12.009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Protein force fields have been undergoing continual development since the first complete parameter sets were introduced nearly four decades ago. The functional forms that underlie these models have many common elements for the treatment of bonded and nonbonded forces, which are reviewed here. The most widely used force fields to date use a fixed-charge convention in which electronic polarization effects are treated via a mean-field approximation during partial charge assignment. Despite success in modeling folded proteins over many years, the fixed-charge assumption has limitations that cannot necessarily be overcome within their potential energy equations. To overcome these limitations, several force fields have recently been derived that explicitly treat electronic polarization effects with straightforward extensions of the potential energy functions used by nonpolarizable force fields. Here, we review the history of the most popular nonpolarizable force fields (AMBER, CHARMM, OPLS, and GROMOS) as well as studies that have validated them and applied them to studies of protein folding and misfolding. Building upon these force fields are more recent polarizable interaction potentials, including fluctuating charge models, POSSIM, AMOEBA, and the classical Drude oscillator. These force fields differ in their implementations but all attempt to model electronic polarization in a computationally tractable manner. Despite their recent emergence in the field of protein folding, several studies have already applied these polarizable models to challenging problems in this domain, including the role of polarization in folding free energies and sequence-specific effects on the stability of α-helical structures.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA, United States.
| |
Collapse
|
49
|
Abstract
Intrinsically disordered proteins (IDPs) can adopt a range of conformations from globules to swollen coils. This large range of conformational preferences for different IDPs raises the question of how conformational preferences are encoded by sequence. Global compositional features of a sequence such as the fraction of charged residues and the net charge per residue engender certain conformational biases. However, more specific sequence features such as the patterning of oppositely charged residues, expansion driving residues, or residues that can undergo posttranslational modifications can also influence the conformational ensembles of an IDP. Here, we outline how to calculate important global compositional features and patterning metrics that can be used to classify IDPs into different conformational classes and predict relative changes in conformation for sequences with the same amino acid composition. Although increased effort has been devoted to determining conformational properties of IDPs in recent years, quantitative predictions of conformation directly from sequence remain difficult and often inaccurate. Thus, if quantitative predictions of conformational properties are desired, then sequence-specific simulations must be performed.
Collapse
Affiliation(s)
- Kiersten M Ruff
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA.
| |
Collapse
|
50
|
Fagerberg E, Lenton S, Skepö M. Evaluating Models of Varying Complexity of Crowded Intrinsically Disordered Protein Solutions Against SAXS. J Chem Theory Comput 2019; 15:6968-6983. [DOI: 10.1021/acs.jctc.9b00723] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Eric Fagerberg
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
| | - Samuel Lenton
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
- LINXS - Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, 223 70 Lund, Sweden
| | - Marie Skepö
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
- LINXS - Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, 223 70 Lund, Sweden
| |
Collapse
|