1
|
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
|
2
|
Duran T, Naik S, Sharifi L, DiLuzio WR, Chanda A, Chaudhuri B. Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations. Int J Pharm 2024; 655:123985. [PMID: 38484860 DOI: 10.1016/j.ijpharm.2024.123985] [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] [Received: 11/11/2023] [Revised: 03/05/2024] [Accepted: 03/07/2024] [Indexed: 03/19/2024]
Abstract
The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.
Collapse
Affiliation(s)
- Tibo Duran
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Connecticut, Storrs, CT 06269, USA
| | - Shivangi Naik
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
| | - Leila Sharifi
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Connecticut, Storrs, CT 06269, USA
| | - Willow R DiLuzio
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
| | - Arani Chanda
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
| | - Bodhisattwa Chaudhuri
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Connecticut, Storrs, CT 06269, USA; Department of Chemical & Biomolecular Engineering, University of Connecticut, Storrs, CT 06269, USA; Institute of Material Sciences (IMS), University of Connecticut, Storrs, CT, USA.
| |
Collapse
|
3
|
Latham AP, Zhu L, Sharon DA, Ye S, Willard AP, Zhang X, Zhang B. Microphase Separation Produces Interfacial Environment within Diblock Biomolecular Condensates. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.03.30.534967. [PMID: 37034777 PMCID: PMC10081284 DOI: 10.1101/2023.03.30.534967] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvironments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in an interfacial environment, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.
Collapse
Affiliation(s)
- Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Longchen Zhu
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
- These authors contributed equally to this work
| | - Dina A Sharon
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
- These authors contributed equally to this work
| | - Songtao Ye
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
| | - Adam P Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Xin Zhang
- Department of Chemistry, School of Science and Research Center for Industries of the Future, Westlake University, 600 Dunyu Road, Hangzhou 310030, Zhejiang Province, China
- Westlake Laboratory of Life Sciences and Biomedicine, 18 Shilongshan Road, Hangzhou 310024, Zhejiang Province, China
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|
4
|
Mondal S, Cui Q. Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates. J Phys Chem B 2024; 128:2087-2099. [PMID: 38407041 DOI: 10.1021/acs.jpcb.3c08149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/27/2024]
Abstract
Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.
Collapse
Affiliation(s)
- Sayantan Mondal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
5
|
Zerze GH. Optimizing the Martini 3 Force Field Reveals the Effects of the Intricate Balance between Protein-Water Interaction Strength and Salt Concentration on Biomolecular Condensate Formation. J Chem Theory Comput 2024; 20:1646-1655. [PMID: 37043540 DOI: 10.1021/acs.jctc.2c01273] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2023]
Abstract
Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. The MARTINI force field offers a coarse-grained protein model with a resolution that preserves molecular details with an explicit (CG) solvent. Despite its relatively higher resolution, it can still achieve condensate formation in a reasonable computing time with explicit solvent and ionic species. Therefore, it is highly desirable to tune this force field to be able to reproduce the experimentally observed properties of the condensate formation. In this work, we studied the condensate formation of the low-sequence complexity domain of fused in sarcoma protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying the salt concentration. We found that the condensate formation is sensitive both to the protein-water interaction strength and the presence of salt. While the unmodified MARTINI force field yields a complete collapse of proteins into one dense phase (i.e., no dilute phase), we reported a range of modified protein-water interaction strength that is capable of capturing the experimentally found transfer free energy between dense and dilute phases. We also found that the condensates lose their spherical shape upon the addition of salt, especially when the protein-water interactions are weak. Interchain amino acid contact map analysis showed one explanation for this observation: the protein-protein contact fraction reduces as salt is added to systems (when the protein-water interactions are weak), consistent with electrostatic screening effects. This reduction might be responsible for the condensates becoming nonspherical upon the addition of salt by reducing the need for minimizing the interfacial area. However, as the protein-water interactions become stronger to the extent that makes the transfer free energy agree well with experimentally observed transfer free energy, we found an increase in the protein-protein contact fraction upon the addition of salt, consistent with the salting-out effects. Therefore, we concluded that there is an intricate balance between screening effects and salting-out effects upon the addition of salt and this balance is highly sensitive to the strength of protein-water interactions.
Collapse
Affiliation(s)
- Gül H Zerze
- William A. Brookshire Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| |
Collapse
|
6
|
Kharche S, Yadav M, Hande V, Prakash S, Sengupta D. Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model. J Chem Inf Model 2024; 64:837-850. [PMID: 38291973 DOI: 10.1021/acs.jcim.3c00802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.
Collapse
Affiliation(s)
- Shalmali Kharche
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Manjul Yadav
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Vrushali Hande
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Shikha Prakash
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
| | - Durba Sengupta
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
7
|
Hosseini AN, van der Spoel D. Martini on the Rocks: Can a Coarse-Grained Force Field Model Crystals? J Phys Chem Lett 2024; 15:1079-1088. [PMID: 38261634 PMCID: PMC10839907 DOI: 10.1021/acs.jpclett.4c00012] [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: 01/02/2024] [Revised: 01/18/2024] [Accepted: 01/19/2024] [Indexed: 01/25/2024]
Abstract
Computational chemistry is an important tool in numerous scientific disciplines, including drug discovery and structural biology. Coarse-grained models offer simple representations of molecular systems that enable simulations of large-scale systems. Because there has been an increase in the adoption of such models for simulations of biomolecular systems, critical evaluation is warranted. Here, the stability of the amyloid peptide and organic crystals is evaluated using the Martini 3 coarse-grained force field. The crystals change shape drastically during the simulations. Radial distribution functions show that the distance between backbone beads in β-sheets increases by ∼1 Å, breaking the crystals. The melting points of organic compounds are much too low in the Martini force field. This suggests that Martini 3 lacks the specific interactions needed to accurately simulate peptides or organic crystals without imposing artificial restraints. The problems may be exacerbated by the use of the 12-6 potential, suggesting that a softer potential could improve this model for crystal simulations.
Collapse
Affiliation(s)
- A. Najla Hosseini
- Department of Cell and Molecular
Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular
Biology, Uppsala University, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
8
|
Versini R, Sritharan S, Aykac Fas B, Tubiana T, Aimeur SZ, Henri J, Erard M, Nüsse O, Andreani J, Baaden M, Fuchs P, Galochkina T, Chatzigoulas A, Cournia Z, Santuz H, Sacquin-Mora S, Taly A. A Perspective on the Prospective Use of AI in Protein Structure Prediction. J Chem Inf Model 2024; 64:26-41. [PMID: 38124369 DOI: 10.1021/acs.jcim.3c01361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
AlphaFold2 (AF2) and RoseTTaFold (RF) have revolutionized structural biology, serving as highly reliable and effective methods for predicting protein structures. This article explores their impact and limitations, focusing on their integration into experimental pipelines and their application in diverse protein classes, including membrane proteins, intrinsically disordered proteins (IDPs), and oligomers. In experimental pipelines, AF2 models help X-ray crystallography in resolving the phase problem, while complementarity with mass spectrometry and NMR data enhances structure determination and protein flexibility prediction. Predicting the structure of membrane proteins remains challenging for both AF2 and RF due to difficulties in capturing conformational ensembles and interactions with the membrane. Improvements in incorporating membrane-specific features and predicting the structural effect of mutations are crucial. For intrinsically disordered proteins, AF2's confidence score (pLDDT) serves as a competitive disorder predictor, but integrative approaches including molecular dynamics (MD) simulations or hydrophobic cluster analyses are advocated for accurate dynamics representation. AF2 and RF show promising results for oligomeric models, outperforming traditional docking methods, with AlphaFold-Multimer showing improved performance. However, some caveats remain in particular for membrane proteins. Real-life examples demonstrate AF2's predictive capabilities in unknown protein structures, but models should be evaluated for their agreement with experimental data. Furthermore, AF2 models can be used complementarily with MD simulations. In this Perspective, we propose a "wish list" for improving deep-learning-based protein folding prediction models, including using experimental data as constraints and modifying models with binding partners or post-translational modifications. Additionally, a meta-tool for ranking and suggesting composite models is suggested, driving future advancements in this rapidly evolving field.
Collapse
Affiliation(s)
- Raphaelle Versini
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Sujith Sritharan
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Burcu Aykac Fas
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Thibault Tubiana
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
| | - Sana Zineb Aimeur
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Julien Henri
- Sorbonne Université, CNRS, Laboratoire de Biologie, Computationnelle et Quantitative UMR 7238, Institut de Biologie Paris-Seine, 4 Place Jussieu, F-75005 Paris, France
| | - Marie Erard
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Oliver Nüsse
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Jessica Andreani
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
| | - Marc Baaden
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Patrick Fuchs
- Sorbonne Université, École Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, 75005 Paris, France
- Université de Paris, UFR Sciences du Vivant, 75013 Paris, France
| | - Tatiana Galochkina
- Université Paris Cité and Université des Antilles and Université de la Réunion, INSERM, BIGR, F-75014 Paris, France
| | - Alexios Chatzigoulas
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
- Department of Informatics and Telecommunications, National and Kapodistrian University of Athens, 15784 Athens, Greece
| | - Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
- Department of Informatics and Telecommunications, National and Kapodistrian University of Athens, 15784 Athens, Greece
| | - Hubert Santuz
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Sophie Sacquin-Mora
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Antoine Taly
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| |
Collapse
|
9
|
Ghafouri H, Lazar T, Del Conte A, Tenorio Ku LG, Tompa P, Tosatto SCE, Monzon AM. PED in 2024: improving the community deposition of structural ensembles for intrinsically disordered proteins. Nucleic Acids Res 2024; 52:D536-D544. [PMID: 37904608 PMCID: PMC10767937 DOI: 10.1093/nar/gkad947] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 10/10/2023] [Accepted: 10/13/2023] [Indexed: 11/01/2023] Open
Abstract
The Protein Ensemble Database (PED) (URL: https://proteinensemble.org) is the primary resource for depositing structural ensembles of intrinsically disordered proteins. This updated version of PED reflects advancements in the field, denoting a continual expansion with a total of 461 entries and 538 ensembles, including those generated without explicit experimental data through novel machine learning (ML) techniques. With this significant increment in the number of ensembles, a few yet-unprecedented new entries entered the database, including those also determined or refined by electron paramagnetic resonance or circular dichroism data. In addition, PED was enriched with several new features, including a novel deposition service, improved user interface, new database cross-referencing options and integration with the 3D-Beacons network-all representing efforts to improve the FAIRness of the database. Foreseeably, PED will keep growing in size and expanding with new types of ensembles generated by accurate and fast ML-based generative models and coarse-grained simulations. Therefore, among future efforts, priority will be given to further develop the database to be compatible with ensembles modeled at a coarse-grained level.
Collapse
Affiliation(s)
| | - Tamas Lazar
- VIB-VUB Center for Structural Biology, Vlaams Instituut voor Biotechnologie (VIB), Brussels, Belgium
- Structural Biology Brussels, Department of Bioengineering, Vrije Universiteit Brussel (VUB), Brussels, Belgium
| | - Alessio Del Conte
- Department of Biomedical Sciences, University of Padova, Padova, Italy
| | | | - Peter Tompa
- VIB-VUB Center for Structural Biology, Vlaams Instituut voor Biotechnologie (VIB), Brussels, Belgium
- Structural Biology Brussels, Department of Bioengineering, Vrije Universiteit Brussel (VUB), Brussels, Belgium
- Institute of Enzymology, Research Centre for Natural Sciences (RCNS), Budapest, Hungary
| | | | | |
Collapse
|
10
|
Mehta N, Mondal S, Watson ET, Cui Q, Chapman ER. The juxtamembrane linker of synaptotagmin 1 regulates Ca 2+ binding via liquid-liquid phase separation. Nat Commun 2024; 15:262. [PMID: 38177243 PMCID: PMC10766989 DOI: 10.1038/s41467-023-44414-5] [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: 08/11/2023] [Accepted: 12/12/2023] [Indexed: 01/06/2024] Open
Abstract
Synaptotagmin (syt) 1, a Ca2+ sensor for synaptic vesicle exocytosis, functions in vivo as a multimer. Syt1 senses Ca2+ via tandem C2-domains that are connected to a single transmembrane domain via a juxtamembrane linker. Here, we show that this linker segment harbors a lysine-rich, intrinsically disordered region that is necessary and sufficient to mediate liquid-liquid phase separation (LLPS). Interestingly, condensate formation negatively regulates the Ca2+-sensitivity of syt1. Moreover, Ca2+ and anionic phospholipids facilitate the observed phase separation, and increases in [Ca2+]i promote the fusion of syt1 droplets in living cells. Together, these observations suggest a condensate-mediated feedback loop that serves to fine-tune the ability of syt1 to trigger release, via alterations in Ca2+ binding activity and potentially through the impact of LLPS on membrane curvature during fusion reactions. In summary, the juxtamembrane linker of syt1 emerges as a regulator of syt1 function by driving self-association via LLPS.
Collapse
Affiliation(s)
- Nikunj Mehta
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin-Madison, Madison, WI, 53705, USA
| | - Sayantan Mondal
- Department of Chemistry, Boston University, Boston, MA, 02215, USA
| | - Emma T Watson
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin-Madison, Madison, WI, 53705, USA
| | - Qiang Cui
- Department of Chemistry, Boston University, Boston, MA, 02215, USA
| | - Edwin R Chapman
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin-Madison, Madison, WI, 53705, USA.
| |
Collapse
|
11
|
Naglekar A, Chattopadhyay A, Sengupta D. Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment. J Phys Chem B 2023; 127:11000-11010. [PMID: 38111968 DOI: 10.1021/acs.jpcb.3c05930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.
Collapse
Affiliation(s)
- Amit Naglekar
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Amitabha Chattopadhyay
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
- CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad 500007, India
| | - Durba Sengupta
- CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
12
|
Mehta N, Mondal S, Watson ET, Cui Q, Chapman ER. The juxtamembrane linker of synaptotagmin 1 regulates Ca 2+ binding via liquid-liquid phase separation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.11.551903. [PMID: 37609296 PMCID: PMC10441399 DOI: 10.1101/2023.08.11.551903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2023]
Abstract
Synaptotagmin (syt) 1, a Ca2+ sensor for synaptic vesicle exocytosis, functions in vivo as a multimer. Syt1 senses Ca2+ via tandem C2-domains that are connected to a single transmembrane domain via a juxtamembrane linker. Here, we show that this linker segment harbors a lysine-rich, intrinsically disordered region that is necessary and sufficient to mediate liquid-liquid phase separation (LLPS). Interestingly, condensate formation negatively regulates the Ca2+-sensitivity of syt1. Moreover, Ca2+ and anionic phospholipids facilitate the observed phase separation, and increases in [Ca2+]i promote the fusion of syt1 droplets in living cells. Together, these observations suggest a condensate-mediated feedback loop that serves to fine-tune the ability of syt1 to trigger release, via alterations in Ca2+ binding activity and potentially through the impact of LLPS on membrane curvature during fusion reactions. In summary, the juxtamembrane linker of syt1 emerges as a regulator of syt1 function by driving self-association via LLPS.
Collapse
Affiliation(s)
- Nikunj Mehta
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin–Madison, Madison, WI 53705, United States
| | - Sayantan Mondal
- Department of Chemistry, Boston University, Boston, MA 02215, United States
| | - Emma T. Watson
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin–Madison, Madison, WI 53705, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, Boston, MA 02215, United States
| | - Edwin R. Chapman
- Howard Hughes Medical Institute, Department of Neuroscience, University of Wisconsin–Madison, Madison, WI 53705, United States
| |
Collapse
|
13
|
Yamada T, Miyazaki Y, Harada S, Kumar A, Vanni S, Shinoda W. Improved Protein Model in SPICA Force Field. J Chem Theory Comput 2023; 19:8967-8977. [PMID: 37989551 DOI: 10.1021/acs.jctc.3c01016] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.
Collapse
Affiliation(s)
- Teppei Yamada
- Graduate School of Natural Science and Technology, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Shogo Harada
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Ashutosh Kumar
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
14
|
Ingólfsson HI, Rizuan A, Liu X, Mohanty P, Souza PCT, Marrink SJ, Bowers MT, Mittal J, Berry J. Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation. Biophys J 2023; 122:4370-4381. [PMID: 37853696 PMCID: PMC10720261 DOI: 10.1016/j.bpj.2023.10.016] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 05/27/2023] [Accepted: 10/16/2023] [Indexed: 10/20/2023] Open
Abstract
The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent C⍺ model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.
Collapse
Affiliation(s)
- Helgi I Ingólfsson
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California.
| | - Azamat Rizuan
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas
| | - Xikun Liu
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California; Department of Chemistry & Biochemistry, University of California Santa Barbara, Santa Barbara, California
| | - Priyesh Mohanty
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS & University of Lyon, Lyon, France; Laboratory of Biology and Modeling of the Cell, École Normale Supérieure de Lyon, Université Claude Bernard Lyon 1, CNRS UMR 5239 and Inserm U1293, 46 Allée d'Italie, Lyon, France
| | - Siewert J Marrink
- Groningen Biomolecular Science and Biotechnology Institute, University of Groningen, Groningen, the Netherlands
| | - Michael T Bowers
- Department of Chemistry & Biochemistry, University of California Santa Barbara, Santa Barbara, California
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M College of Engineering, College Station, Texas; Department of Chemistry, Texas A&M University, College Station, Texas; Interdisciplinary Graduate Program in Genetics and Genomics, Texas A&M University, College Station, Texas
| | - Joel Berry
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California
| |
Collapse
|
15
|
Forsberg Z, Stepnov AA, Tesei G, Wang Y, Buchinger E, Kristiansen SK, Aachmann FL, Arleth L, Eijsink VGH, Lindorff-Larsen K, Courtade G. The effect of linker conformation on performance and stability of a two-domain lytic polysaccharide monooxygenase. J Biol Chem 2023; 299:105262. [PMID: 37734553 PMCID: PMC10598543 DOI: 10.1016/j.jbc.2023.105262] [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] [Received: 06/23/2023] [Revised: 09/05/2023] [Accepted: 09/14/2023] [Indexed: 09/23/2023] Open
Abstract
A considerable number of lytic polysaccharide monooxygenases (LPMOs) and other carbohydrate-active enzymes are modular, with catalytic domains being tethered to additional domains, such as carbohydrate-binding modules, by flexible linkers. While such linkers may affect the structure, function, and stability of the enzyme, their roles remain largely enigmatic, as do the reasons for natural variation in length and sequence. Here, we have explored linker functionality using the two-domain cellulose-active ScLPMO10C from Streptomyces coelicolor as a model system. In addition to investigating the WT enzyme, we engineered three linker variants to address the impact of both length and sequence and characterized these using small-angle X-ray scattering, NMR, molecular dynamics simulations, and functional assays. The resulting data revealed that, in the case of ScLPMO10C, linker length is the main determinant of linker conformation and enzyme performance. Both the WT and a serine-rich variant, which have the same linker length, demonstrated better performance compared with those with either a shorter linker or a longer linker. A highlight of our findings was the substantial thermostability observed in the serine-rich variant. Importantly, the linker affects thermal unfolding behavior and enzyme stability. In particular, unfolding studies show that the two domains unfold independently when mixed, whereas the full-length enzyme shows one cooperative unfolding transition, meaning that the impact of linkers in biomass-processing enzymes is more complex than mere structural tethering.
Collapse
Affiliation(s)
- Zarah Forsberg
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences (NMBU), Ås, Norway.
| | - Anton A Stepnov
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences (NMBU), Ås, Norway
| | - Giulio Tesei
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark
| | - Yong Wang
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark; College of Life Sciences, Zhejiang University, Hangzhou, China
| | - Edith Buchinger
- Vectron Biosolutions AS, Trondheim, Norway; Department of Biotechnology and Food Science, NTNU Norwegian University of Science and Technology, Trondheim, Norway
| | - Sandra K Kristiansen
- Department of Biotechnology and Food Science, NTNU Norwegian University of Science and Technology, Trondheim, Norway
| | - Finn L Aachmann
- Department of Biotechnology and Food Science, NTNU Norwegian University of Science and Technology, Trondheim, Norway
| | - Lise Arleth
- X-ray and Neutron Science, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
| | - Vincent G H Eijsink
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences (NMBU), Ås, Norway
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Department of Biology, Linderstrøm-Lang Centre for Protein Science, University of Copenhagen, Copenhagen, Denmark
| | - Gaston Courtade
- Department of Biotechnology and Food Science, NTNU Norwegian University of Science and Technology, Trondheim, Norway.
| |
Collapse
|
16
|
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
|
17
|
Álvarez D, Sapia J, Vanni S. Computational modeling of membrane trafficking processes: From large molecular assemblies to chemical specificity. Curr Opin Cell Biol 2023; 83:102205. [PMID: 37451175 DOI: 10.1016/j.ceb.2023.102205] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Revised: 06/16/2023] [Accepted: 06/16/2023] [Indexed: 07/18/2023]
Abstract
In the last decade, molecular dynamics (MD) simulations have become an essential tool to investigate the molecular properties of membrane trafficking processes, often in conjunction with experimental approaches. The combination of MD simulations with recent developments in structural biology, such as cryo-electron microscopy and artificial intelligence-based structure determination, opens new, exciting possibilities for future investigations. However, the full potential of MD simulations to provide a molecular view of the complex and dynamic processes involving membrane trafficking can only be realized if certain limitations are addressed, and especially those concerning the quality of coarse-grain models, which, despite recent successes in describing large-scale systems, still suffer from far-from-ideal chemical accuracy. In this review, we will highlight recent success stories of MD simulations in the investigation of membrane trafficking processes, their implications for future research, and the challenges that lie ahead in this specific research domain.
Collapse
Affiliation(s)
- Daniel Álvarez
- Department of Biology, University of Fribourg, Switzerland; Departamento de Química Física y Analítica, Universidad de Oviedo, Spain
| | - Jennifer Sapia
- Department of Biology, University of Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Switzerland; Université Côte d'Azur, Inserm, CNRS, Institut de Pharmacologie Moléculaire et Cellulaire, Valbonne, France.
| |
Collapse
|
18
|
Yu H, Qin L, Zhou J. Effect of Oil Polarity on the Protein Adsorption at Oil-Water Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:10701-10710. [PMID: 37470337 DOI: 10.1021/acs.langmuir.3c01541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/21/2023]
Abstract
Protein adsorption at oil-water interfaces has received much attention in applications of food emulsion and biocatalysis. The protein activity is influenced by the protein orientation and conformation. The oil polarity is expected to influence the orientation and conformation of adsorbed proteins by modulating intermolecular interactions. Hence, it is possible to tune the protein emulsion stability and activity by varying the oil polarity. Martini v3.0-based coarse-grained molecular dynamics (CGMD) simulations were employed to investigate the effect of oil polarity on the orientation and conformation of hydrophobin (HFBI) and Candida antarctica lipase B (CALB) adsorbed at triolein-water, hexadecane-water, and octanol-water interfaces for the first time. The protein adsorption orientation was predicted through the hydrophobic dipole, indicating that protein adsorption exists in preferred orientations at hydrophobic oil interfaces. The conformation of the adsorbed HFBI is well conserved, whereas relatively larger conformational changes occur during the CALB adsorption as the oil hydrophobicity increases. Comparisons on the adsorption interaction energy of proteins with oils confirm the relationship between the oil polarity and the interaction strength of proteins with oils. In addition, CGMD simulations allow longer time scale simulations of the behaviors of protein adsorption at oil-water interfaces.
Collapse
Affiliation(s)
- Hai Yu
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Lanlan Qin
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Jian Zhou
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| |
Collapse
|
19
|
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
|
20
|
Fagerberg E, Skepö M. Comparative Performance of Computer Simulation Models of Intrinsically Disordered Proteins at Different Levels of Coarse-Graining. J Chem Inf Model 2023. [PMID: 37339604 DOI: 10.1021/acs.jcim.3c00113] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2023]
Abstract
Coarse-graining is commonly used to decrease the computational cost of simulations. However, coarse-grained models are also considered to have lower transferability, with lower accuracy for systems outside the original scope of parametrization. Here, we benchmark a bead-necklace model and a modified Martini 2 model, both coarse-grained models, for a set of intrinsically disordered proteins, with the different models having different degrees of coarse-graining. The SOP-IDP model has earlier been used for this set of proteins; thus, those results are included in this study to compare how models with different levels of coarse-graining compare. The sometimes naive expectation of the least coarse-grained model performing best does not hold true for the experimental pool of proteins used here. Instead, it showed the least good agreement, indicating that one should not necessarily trust the otherwise intuitive notion of a more advanced model inherently being better in model choice.
Collapse
Affiliation(s)
- Eric Fagerberg
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
| | - Marie Skepö
- Theoretical Chemistry, Lund University, POB 124, SE-221 00 Lund, Sweden
- LINXS - Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
21
|
Paladino A, Vitagliano L, Graziano G. The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins. BIOLOGY 2023; 12:biology12050754. [PMID: 37237566 DOI: 10.3390/biology12050754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 05/15/2023] [Accepted: 05/18/2023] [Indexed: 05/28/2023]
Abstract
Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.
Collapse
Affiliation(s)
- Antonella Paladino
- Institute of Biostructures and Bioimaging, CNR, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Luigi Vitagliano
- Institute of Biostructures and Bioimaging, CNR, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Giuseppe Graziano
- Department of Science and Technology, University of Sannio, via Francesco de Sanctis snc, 82100 Benevento, Italy
| |
Collapse
|
22
|
Di Bartolo AL, Caparotta M, Masone D. Intrinsic Disorder in α-Synuclein Regulates the Exocytotic Fusion Pore Transition. ACS Chem Neurosci 2023. [PMID: 37192400 DOI: 10.1021/acschemneuro.3c00040] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023] Open
Abstract
Today, it is widely accepted that intrinsic disorder is strongly related to the cell cycle, during mitosis, differentiation, and apoptosis. Of particular interest are hybrid proteins possessing both structured and unstructured domains that are critical in human health and disease, such as α-synuclein. In this work, we describe how α-synuclein interacts with the nascent fusion pore as it evolves toward expansion. We unveil the key role played by its intrinsically disordered region as a thermodynamic regulator of the nucleation-expansion energy barrier. By analyzing a truncated variant of α-synuclein that lacks the disordered region, we find that the landscape of protein interactions with PIP2 and POPS lipids is highly altered, ultimately increasing the energy cost for the fusion pore to transit from nucleation to expansion. We conclude that the intrinsically disordered region in full-length α-synuclein recognizes and allocates pivotal protein:lipid interactions during membrane remodeling in the first stages of the fusion pore.
Collapse
Affiliation(s)
- Ary Lautaro Di Bartolo
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| | - Marcelo Caparotta
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Diego Masone
- Instituto de Histología y Embriología de Mendoza (IHEM) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
- Facultad de Ingeniería, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| |
Collapse
|
23
|
Jussupow A, Kaila VRI. Effective Molecular Dynamics from Neural Network-Based Structure Prediction Models. J Chem Theory Comput 2023; 19:1965-1975. [PMID: 36961997 PMCID: PMC11181330 DOI: 10.1021/acs.jctc.2c01027] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Indexed: 03/26/2023]
Abstract
Recent breakthroughs in neural network-based structure prediction methods, such as AlphaFold2 and RoseTTAFold, have dramatically improved the quality of computational protein structure prediction. These models also provide statistical confidence scores that can estimate uncertainties in the predicted structures, but it remains unclear to what extent these scores are related to the intrinsic conformational dynamics of proteins. Here, we compare AlphaFold2 prediction scores with explicit large-scale molecular dynamics simulations of 28 one- and two-domain proteins with varying degrees of flexibility. We demonstrate a strong correlation between the statistical prediction scores and the explicit motion derived from extensive atomistic molecular dynamics simulations and further derive an elastic network model based on the statistical scores of AlphFold2 (AF-ENM), which we benchmark in combination with coarse-grained molecular dynamics simulations. We show that our AF-ENM method reproduces the global protein dynamics with improved accuracy, providing a powerful way to derive effective molecular dynamics using neural network-based structure prediction models.
Collapse
Affiliation(s)
- Alexander Jussupow
- Department of Biochemistry
and Biophysics, Stockholm University, 10691 Stockholm, Sweden
| | - Ville R. I. Kaila
- Department of Biochemistry
and Biophysics, Stockholm University, 10691 Stockholm, Sweden
| |
Collapse
|
24
|
Saurabh S, Nadendla K, Purohit SS, Sivakumar PM, Cetinel S. Fuzzy Drug Targets: Disordered Proteins in the Drug-Discovery Realm. ACS OMEGA 2023; 8:9729-9747. [PMID: 36969402 PMCID: PMC10034788 DOI: 10.1021/acsomega.2c07708] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 02/17/2023] [Indexed: 06/18/2023]
Abstract
Intrinsically disordered proteins (IDPs) and regions (IDRs) form a large part of the eukaryotic proteome. Contrary to the structure-function paradigm, the disordered proteins perform a myriad of functions in vivo. Consequently, they are involved in various disease pathways and are plausible drug targets. Unlike folded proteins, that have a defined structure and well carved out drug-binding pockets that can guide lead molecule selection, the disordered proteins require alternative drug-development methodologies that are based on an acceptable picture of their conformational ensemble. In this review, we discuss various experimental and computational techniques that contribute toward understanding IDP "structure" and describe representative pursuances toward IDP-targeting drug development. We also discuss ideas on developing rational drug design protocols targeting IDPs.
Collapse
Affiliation(s)
- Suman Saurabh
- Molecular
Sciences Research Hub, Department of Chemistry, Imperial College London, London W12 0BZ, U.K.
| | - Karthik Nadendla
- Center
for Misfolding Diseases, Yusuf Hamied Department of Chemistry, Lensfield
Road, University of Cambridge, Cambridge CB2 1EW, U.K.
| | - Shubh Sanket Purohit
- Department
of Clinical Haematology, Sahyadri Superspeciality
Hospital, Pune, Maharashtra 411038, India
| | - Ponnurengam Malliappan Sivakumar
- Institute
of Research and Development, Duy Tan University, Da Nang 550000, Vietnam
- School
of Medicine and Pharmacy, Duy Tan University, Da Nang 550000, Vietnam
- Nanotechnology
Research and Application Center (SUNUM), Sabanci University, Istanbul 34956, Turkey
| | - Sibel Cetinel
- Nanotechnology
Research and Application Center (SUNUM), Sabanci University, Istanbul 34956, Turkey
- Faculty of
Engineering and Natural Sciences, Molecular Biology, Genetics and
Bioengineering Program, Sabanci University, Istanbul 34956, Turkey
| |
Collapse
|
25
|
Thomasen FE, Cuneo MJ, Mittag T, Lindorff-Larsen K. Conformational and oligomeric states of SPOP from small-angle X-ray scattering and molecular dynamics simulations. eLife 2023; 12:e84147. [PMID: 36856266 PMCID: PMC9998093 DOI: 10.7554/elife.84147] [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: 10/12/2022] [Accepted: 02/20/2023] [Indexed: 03/02/2023] Open
Abstract
Speckle-type POZ protein (SPOP) is a substrate adaptor in the ubiquitin proteasome system, and plays important roles in cell-cycle control, development, and cancer pathogenesis. SPOP forms linear higher-order oligomers following an isodesmic self-association model. Oligomerization is essential for SPOP's multivalent interactions with substrates, which facilitate phase separation and localization to biomolecular condensates. Structural characterization of SPOP in its oligomeric state and in solution is, however, challenging due to the inherent conformational and compositional heterogeneity of the oligomeric species. Here, we develop an approach to simultaneously and self-consistently characterize the conformational ensemble and the distribution of oligomeric states of SPOP by combining small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. We build initial conformational ensembles of SPOP oligomers using coarse-grained molecular dynamics simulations, and use a Bayesian/maximum entropy approach to refine the ensembles, along with the distribution of oligomeric states, against a concentration series of SAXS experiments. Our results suggest that SPOP oligomers behave as rigid, helical structures in solution, and that a flexible linker region allows SPOP's substrate-binding domains to extend away from the core of the oligomers. Additionally, our results are in good agreement with previous characterization of the isodesmic self-association of SPOP. In the future, the approach presented here can be extended to other systems to simultaneously characterize structural heterogeneity and self-assembly.
Collapse
Affiliation(s)
- F Emil Thomasen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of CopenhagenCopenhagenDenmark
| | - Matthew J Cuneo
- Department of Structural Biology, St. Jude Children’s Research HospitalMemphisUnited States
| | - Tanja Mittag
- Department of Structural Biology, St. Jude Children’s Research HospitalMemphisUnited States
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of CopenhagenCopenhagenDenmark
| |
Collapse
|
26
|
Claveras Cabezudo A, Athanasiou C, Tsengenes A, Wade RC. Scaling Protein-Water Interactions in the Martini 3 Coarse-Grained Force Field to Simulate Transmembrane Helix Dimers in Different Lipid Environments. J Chem Theory Comput 2023; 19:2109-2119. [PMID: 36821400 DOI: 10.1021/acs.jctc.2c00950] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Abstract
Martini 3, the latest version of the widely used Martini force field for coarse-grained molecular dynamics simulations, is a promising tool to investigate proteins in phospholipid bilayers. However, simulating other lipid environments, such as detergent micelles, presents challenges due to the absence of validated parameters for their constituent molecules. Here, we propose parameters for the micelle-forming surfactant, dodecylphosphocholine (DPC). These result in micelle assembly with aggregation numbers in agreement with the experimental values. However, we identified a lack of hydrophobic interactions between transmembrane helix protein dimers and the tails of DPC molecules, preventing insertion and stabilization of the protein in the micelles. This problem was also observed for protein insertion by self-assembling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or dipalmitoylphosphatidylcholine (DPPC) bilayers. We propose the reduction of the nonbonded interactions between protein and water beads by 10% as a simple and effective solution to this problem that enables protein encapsulation in phospholipid micelles and bilayers without altering protein dimerization or the bilayer structure.
Collapse
Affiliation(s)
- Ainara Claveras Cabezudo
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany
| | - Christina Athanasiou
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Heidelberg Biosciences International Graduate School, Heidelberg University, Im Neuenheimer Feld 501, 69120 Heidelberg, Germany
| | - Alexandros Tsengenes
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Heidelberg Biosciences International Graduate School, Heidelberg University, Im Neuenheimer Feld 501, 69120 Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
27
|
Zhang M, Xue B, Li Q, Shi R, Cao Y, Wang W, Li J. Sequence Tendency for the Interaction between Low-Complexity Intrinsically Disordered Proteins. JACS AU 2023; 3:93-104. [PMID: 36711093 PMCID: PMC9875249 DOI: 10.1021/jacsau.2c00414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 12/16/2022] [Accepted: 12/16/2022] [Indexed: 06/18/2023]
Abstract
Reversible interaction between intrinsically disordered proteins (IDPs) is considered as the driving force for liquid-liquid phase separation (LLPS), while the detailed description of such a transient interaction process still remains a challenge. And the mechanisms underlying the behavior of IDP interaction, for example, the possible relationship with its inherent conformational fluctuations and sequence features, remain elusive. Here, we use atomistic molecular dynamics (MD) simulation to investigate the reversible association of the LAF-1 RGG domain, the IDP with ultra-low LLPS concentration (0.06 mM). We find that the duration of the association between two RGG domains is highly heterogeneous, and the sustained associations essentially dominate the IDP interaction. More interestingly, such sustained associations are mediated by a finite region, that is, the C-terminal region 138-168 (denoted as a contact-prone region). We noticed that such sequence tendency is attributed to the extended conformation of the RGG domain during its inherent conformational fluctuations. Hence, our results suggest that there is a certain region in this low-complexity IDP which can essentially dominate their interaction and should be also important to the LLPS. And the inherent conformational fluctuations are actually essential for the emergence of such a hot region of IDP interaction. The importance of this hot region to LLPS is verified by experiment.
Collapse
Affiliation(s)
- Moxin Zhang
- Zhejiang
Province Key Laboratory of Quantum Technology and Device, School of
Physics, Zhejiang University, Hangzhou310058, China
| | - Bin Xue
- Collaborative
Innovation Center of Advanced Microstructures, National Laboratory
of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing210093, China
| | - Qingtai Li
- Collaborative
Innovation Center of Advanced Microstructures, National Laboratory
of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing210093, China
| | - Rui Shi
- Zhejiang
Province Key Laboratory of Quantum Technology and Device, School of
Physics, Zhejiang University, Hangzhou310058, China
| | - Yi Cao
- Collaborative
Innovation Center of Advanced Microstructures, National Laboratory
of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing210093, China
| | - Wei Wang
- Collaborative
Innovation Center of Advanced Microstructures, National Laboratory
of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing210093, China
| | - Jingyuan Li
- Zhejiang
Province Key Laboratory of Quantum Technology and Device, School of
Physics, Zhejiang University, Hangzhou310058, China
| |
Collapse
|
28
|
He X, Man VH, Gao J, Wang J. Investigation of the Structure of Full-Length Tau Proteins with Coarse-Grained and All-Atom Molecular Dynamics Simulations. ACS Chem Neurosci 2023; 14:209-217. [PMID: 36563129 DOI: 10.1021/acschemneuro.2c00381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Tau proteins not only have many important biological functions but also are associated with several neurodegenerative diseases, such as Parkinson's disease and Alzheimer's disease (AD). However, it is still a challenge to identify the atomic structure of full-length tau proteins due to their lengthy and disordered characteristics and the factor that there are no crystal structures of full-length tau proteins available. We performed multi- and large-scale molecular dynamics simulations of the full-length tau monomer (the 2N4R isoform and 441 residues) in aqueous solution under biological conditions with coarse-grained and all-atom force fields. The obtained atomic structures produced radii of gyration and chemical shifts that are in excellent agreement with those of experiment. The generated monomer structure ensemble would be very useful for further studying the oligomerization mechanism and discovering tau oligomerization inhibitors, which are important events in AD drug development.
Collapse
Affiliation(s)
- Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Jie Gao
- Department of Neuroscience, The Ohio State University Wexner Medical Center, Columbus, Ohio 43210, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
29
|
Pesce F, Newcombe EA, Seiffert P, Tranchant EE, Olsen JG, Grace CR, Kragelund BB, Lindorff-Larsen K. Assessment of models for calculating the hydrodynamic radius of intrinsically disordered proteins. Biophys J 2023; 122:310-321. [PMID: 36518077 PMCID: PMC9892621 DOI: 10.1016/j.bpj.2022.12.013] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 11/18/2022] [Accepted: 12/09/2022] [Indexed: 12/15/2022] Open
Abstract
Diffusion measurements by pulsed-field gradient NMR and fluorescence correlation spectroscopy can be used to probe the hydrodynamic radius of proteins, which contains information about the overall dimension of a protein in solution. The comparison of this value with structural models of intrinsically disordered proteins is nonetheless impaired by the uncertainty of the accuracy of the methods for computing the hydrodynamic radius from atomic coordinates. To tackle this issue, we here build conformational ensembles of 11 intrinsically disordered proteins that we ensure are in agreement with measurements of compaction by small-angle x-ray scattering. We then use these ensembles to identify the forward model that more closely fits the radii derived from pulsed-field gradient NMR diffusion experiments. Of the models we examined, we find that the Kirkwood-Riseman equation provides the best description of the hydrodynamic radius probed by pulsed-field gradient NMR experiments. While some minor discrepancies remain, our results enable better use of measurements of the hydrodynamic radius in integrative modeling and for force field benchmarking and parameterization.
Collapse
Affiliation(s)
- Francesco Pesce
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Estella A Newcombe
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Pernille Seiffert
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Emil E Tranchant
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Johan G Olsen
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Christy R Grace
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee
| | - Birthe B Kragelund
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
30
|
Sun B, Kekenes-Huskey PM. Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling. Q Rev Biophys 2023; 56:e2. [PMID: 36628457 PMCID: PMC11070111 DOI: 10.1017/s003358352300001x] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.
Collapse
Affiliation(s)
- Bin Sun
- Research Center for Pharmacoinformatics (The State-Province Key Laboratories of Biomedicine-Pharmaceutics of China), Department of Medicinal Chemistry and Natural Medicine Chemistry, College of Pharmacy, Harbin Medical University, Harbin 150081, China
| | | |
Collapse
|
31
|
Mohanty P, Kapoor U, Sundaravadivelu Devarajan D, Phan TM, Rizuan A, Mittal J. Principles Governing the Phase Separation of Multidomain Proteins. Biochemistry 2022; 61:2443-2455. [PMID: 35802394 PMCID: PMC9669140 DOI: 10.1021/acs.biochem.2c00210] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
A variety of membraneless organelles, often termed "biological condensates", play an important role in the regulation of cellular processes such as gene transcription, translation, and protein quality control. On the basis of experimental and theoretical investigations, liquid-liquid phase separation (LLPS) has been proposed as a possible mechanism for the origin of biological condensates. LLPS requires multivalent macromolecules that template the formation of long-range, intermolecular interaction networks and results in the formation of condensates with defined composition and material properties. Multivalent interactions driving LLPS exhibit a wide range of modes from highly stereospecific to nonspecific and involve both folded and disordered regions. Multidomain proteins serve as suitable macromolecules for promoting phase separation and achieving disparate functions due to their potential for multivalent interactions and regulation. Here, we aim to highlight the influence of the domain architecture and interdomain interactions on the phase separation of multidomain protein condensates. First, the general principles underlying these interactions are illustrated on the basis of examples of multidomain proteins that are predominantly associated with nucleic acid binding and protein quality control and contain both folded and disordered regions. Next, the examples showcase how LLPS properties of folded and disordered regions can be leveraged to engineer multidomain constructs that form condensates with the desired assembly and functional properties. Finally, we highlight the need for improvements in coarse-grained computational models that can provide molecular-level insights into multidomain protein condensates in conjunction with experimental efforts.
Collapse
Affiliation(s)
- Priyesh Mohanty
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | - Utkarsh Kapoor
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | | | - Tien Minh Phan
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | - Azamat Rizuan
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, TX 77843
| |
Collapse
|
32
|
Salvi N, Zapletal V, Jaseňáková Z, Zachrdla M, Padrta P, Narasimhan S, Marquardsen T, Tyburn JM, Žídek L, Blackledge M, Ferrage F, Kadeřávek P. Convergent views on disordered protein dynamics from NMR and computational approaches. Biophys J 2022; 121:3785-3794. [PMID: 36131545 PMCID: PMC9674986 DOI: 10.1016/j.bpj.2022.09.016] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Revised: 09/07/2022] [Accepted: 09/15/2022] [Indexed: 11/02/2022] Open
Abstract
Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.
Collapse
Affiliation(s)
- Nicola Salvi
- Institut de Biologie Structurale (IBS), CEA, CNRS, University Grenoble Alpes, Grenoble, France
| | - Vojtěch Zapletal
- National Centre for Biomolecular Research, Faculty of Science and Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Zuzana Jaseňáková
- National Centre for Biomolecular Research, Faculty of Science and Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Milan Zachrdla
- Laboratoire des Biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, Paris, France
| | - Petr Padrta
- Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Subhash Narasimhan
- National Centre for Biomolecular Research, Faculty of Science and Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | | | | | - Lukáš Žídek
- National Centre for Biomolecular Research, Faculty of Science and Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Martin Blackledge
- Institut de Biologie Structurale (IBS), CEA, CNRS, University Grenoble Alpes, Grenoble, France.
| | - Fabien Ferrage
- Laboratoire des Biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, Paris, France.
| | - Pavel Kadeřávek
- Central European Institute of Technology, Masaryk University, Brno, Czech Republic.
| |
Collapse
|
33
|
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
|
34
|
Marrink SJ, Monticelli L, Melo MN, Alessandri R, Tieleman DP, Souza PCT. Two decades of Martini: Better beads, broader scope. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1620] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials University of Groningen Groningen The Netherlands
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| | - Manuel N. Melo
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa Oeiras Portugal
| | - Riccardo Alessandri
- Pritzker School of Molecular Engineering University of Chicago Chicago Illinois USA
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences University of Calgary Alberta Canada
| | - Paulo C. T. Souza
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| |
Collapse
|
35
|
Jussupow A, Lopez A, Baumgart M, Mader SL, Sattler M, Kaila VRI. Extended conformational states dominate the Hsp90 chaperone dynamics. J Biol Chem 2022; 298:102101. [PMID: 35667441 PMCID: PMC9251789 DOI: 10.1016/j.jbc.2022.102101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Revised: 05/29/2022] [Accepted: 05/30/2022] [Indexed: 11/16/2022] Open
Abstract
The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.
Collapse
Affiliation(s)
- Alexander Jussupow
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | - Abraham Lopez
- Center of Integrated Protein Science, Department Chemie, Technische Universität München Lichtenbergstr. 4, 85747 Garching (Germany); Institute of Structural Biology, Helmholtz Zentrum Mu¨nchen, Ingolstädter Landstrasse 1, Neuherberg 85764, Germany
| | - Mona Baumgart
- Center of Integrated Protein Science, Department Chemie, Technische Universität München Lichtenbergstr. 4, 85747 Garching (Germany)
| | - Sophie L Mader
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | - Michael Sattler
- Center of Integrated Protein Science, Department Chemie, Technische Universität München Lichtenbergstr. 4, 85747 Garching (Germany); Institute of Structural Biology, Helmholtz Zentrum Mu¨nchen, Ingolstädter Landstrasse 1, Neuherberg 85764, Germany
| | - Ville R I Kaila
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden.
| |
Collapse
|