1
|
Li S, Wu B, Luo YL, Han W. Simulations of Functional Motions of Super Large Biomolecules with a Mixed-Resolution Model. J Chem Theory Comput 2024; 20:2228-2245. [PMID: 38374639 PMCID: PMC10938502 DOI: 10.1021/acs.jctc.3c01046] [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: 09/22/2023] [Revised: 01/18/2024] [Accepted: 01/29/2024] [Indexed: 02/21/2024]
Abstract
Many large protein machines function through an interplay between large-scale movements and intricate conformational changes. Understanding functional motions of these proteins through simulations becomes challenging for both all-atom and coarse-grained (CG) modeling techniques because neither approach alone can readily capture the full details of these motions. In this study, we develop a multiscale model by employing the popular MARTINI CG model to represent a heterogeneous environment and structurally stable proteins and using the united-atom (UA) model PACE to describe proteins undergoing subtle conformational changes. PACE was previously developed to be compatible with the MARTINI solvent and membrane. Here, we couple the protein descriptions of the two models by directly mixing UA and CG interaction parameters to greatly simplify parameter determination. Through extensive validations with diverse protein systems in solution or membrane, we demonstrate that only additional parameter rescaling is needed to enable the resulting model to recover the stability of native structures of proteins under mixed representation. Moreover, we identify the optimal scaling factors that can be applied to various protein systems, rendering the model potentially transferable. To further demonstrate its applicability for realistic systems, we apply the model to a mechanosensitive ion channel Piezo1 that has peripheral arms for sensing membrane tension and a central pore for ion conductance. The model can reproduce the coupling between Piezo1's large-scale arm movement and subtle pore opening in response to membrane stress while consuming much less computational costs than all-atom models. Therefore, our model shows promise for studying functional motions of large protein machines.
Collapse
Affiliation(s)
- Shu Li
- Centre
for Artificial Intelligence Driven Drug Discovery, Faculty of Applied
Sciences, Macao Polytechnic University, Macao 999078, China
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Bohua Wu
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun Lyna Luo
- Department
of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Wei Han
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Department
of Chemistry, Faculty of Science, Hong Kong
Baptist University, Hong Kong SAR 999077, China
- Shenzhen
Bay Laboratory, Institute of Chemical Biology, Shenzhen 518132, China
| |
Collapse
|
2
|
Sarter T, Friess W. Molecular Dynamics Study of Protein Aggregation at Moving Interfaces. Mol Pharm 2024; 21:1214-1221. [PMID: 38321750 DOI: 10.1021/acs.molpharmaceut.3c00865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
Repeated compression and dilation of a protein film adsorbed to an interface lead to aggregation and entry of film fragments into the bulk. This is a major mechanism for protein aggregate formation in drug products upon mechanical stress, such as shaking or pumping. To gain a better understanding of these events, we developed a molecular dynamics (MD) setup, which would, in a later stage, allow for in silico formulation optimization. In contrast to previous approaches, the molecules of our model protein human growth hormone displayed realistic shapes, surfaces, and interactions with each other and the interface. This enabled quantitative assessment of protein cluster formation. Simulation outcomes aligned with experimental data on subvisible particles and turbidity, thereby validating the model. Computational and experimental results indicated that compression speed does not affect the aggregation behavior of preformed protein films but rather their regeneration. Protein clusters that formed during compression disassembled upon relaxation, suggesting that the particles originate from a partly compressed state. Desorption studies via steered MD revealed that proteins from compressed systems are more likely to detach as clusters, implying that compression effects at the interface translate into aggregates present in the bulk solution. With the possibility of studying the impact of different variables upon compression and dilation at the interface on a molecular level, our model contributes to the understanding of the mechanisms of protein aggregation at moving interfaces. It also enables further studies to change formulation parameters, interfaces, or proteins.
Collapse
Affiliation(s)
- Tim Sarter
- Department of Pharmacy, Pharmaceutical Technology and Biopharmaceutics, Ludwig-Maximilians-Universität München, 81377 Munich, Germany
| | - Wolfgang Friess
- Department of Pharmacy, Pharmaceutical Technology and Biopharmaceutics, Ludwig-Maximilians-Universität München, 81377 Munich, Germany
| |
Collapse
|
3
|
Bergh C, Rovšnik U, Howard R, Lindahl E. Discovery of lipid binding sites in a ligand-gated ion channel by integrating simulations and cryo-EM. eLife 2024; 12:RP86016. [PMID: 38289224 PMCID: PMC10945520 DOI: 10.7554/elife.86016] [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] [Indexed: 02/01/2024] Open
Abstract
Ligand-gated ion channels transduce electrochemical signals in neurons and other excitable cells. Aside from canonical ligands, phospholipids are thought to bind specifically to the transmembrane domain of several ion channels. However, structural details of such lipid contacts remain elusive, partly due to limited resolution of these regions in experimental structures. Here, we discovered multiple lipid interactions in the channel GLIC by integrating cryo-electron microscopy and large-scale molecular simulations. We identified 25 bound lipids in the GLIC closed state, a conformation where none, to our knowledge, were previously known. Three lipids were associated with each subunit in the inner leaflet, including a buried interaction disrupted in mutant simulations. In the outer leaflet, two intrasubunit sites were evident in both closed and open states, while a putative intersubunit site was preferred in open-state simulations. This work offers molecular details of GLIC-lipid contacts particularly in the ill-characterized closed state, testable hypotheses for state-dependent binding, and a multidisciplinary strategy for modeling protein-lipid interactions.
Collapse
Affiliation(s)
- Cathrine Bergh
- Science for Life Laboratory & Swedish e-Science Research Center, Department of Applied Physics, KTH Royal Institute of TechnologyStockholmSweden
| | - Urška Rovšnik
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm UniversityStockholmSweden
| | - Rebecca Howard
- Science for Life Laboratory & Swedish e-Science Research Center, Department of Applied Physics, KTH Royal Institute of TechnologyStockholmSweden
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm UniversityStockholmSweden
| | - Erik Lindahl
- Science for Life Laboratory & Swedish e-Science Research Center, Department of Applied Physics, KTH Royal Institute of TechnologyStockholmSweden
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm UniversityStockholmSweden
| |
Collapse
|
4
|
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
|
5
|
Gasparello J, Verona M, Chilin A, Gambari R, Marzaro G. Assessing the interaction between hemoglobin and the receptor binding domain of SARS-CoV-2 spike protein through MARTINI coarse-grained molecular dynamics. Int J Biol Macromol 2023; 253:127088. [PMID: 37774812 DOI: 10.1016/j.ijbiomac.2023.127088] [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: 07/10/2023] [Revised: 09/22/2023] [Accepted: 09/24/2023] [Indexed: 10/01/2023]
Abstract
The emergence of different coronavirus-related diseases in the 2000's (SARS, MERS, and Covid-19) warrants the need of a complete understanding of the pathological, biological, and biochemical behavior of this class of pathogens. Great attention has been paid to the SARS-CoV-2 Spike protein, and its interaction with the human ACE2 has been thoroughly investigated. Recent findings suggested that the SARS-CoV-2 components may interact with different human proteins, and hemoglobin has very recently been demonstrated as a potential target for the Spike protein. Here we have investigated the interaction between either adult or fetal hemoglobin and the receptor binding domain of the Spike protein at molecular level through advanced molecular dynamics techniques and proposed rational binding modes and energy estimations. Our results agree with biochemical data previously reported in literature. We also demonstrated that co-incubation of pulmonary epithelial cells with hemoglobin strongly reduces the pro-inflammatory effects exerted by the concomitant administration of Spike protein.
Collapse
Affiliation(s)
- Jessica Gasparello
- Department of Life Sciences and Biotechnology, University of Ferrara, via Fossato di Mortara 74, 44121 Ferrara, Italy
| | - Marco Verona
- Department of Pharmaceutical and Pharmacological Sciences, University of Padova, via Marzolo 5, 35313 Padova, Italy
| | - Adriana Chilin
- Department of Pharmaceutical and Pharmacological Sciences, University of Padova, via Marzolo 5, 35313 Padova, Italy
| | - Roberto Gambari
- Department of Life Sciences and Biotechnology, University of Ferrara, via Fossato di Mortara 74, 44121 Ferrara, Italy
| | - Giovanni Marzaro
- Department of Pharmaceutical and Pharmacological Sciences, University of Padova, via Marzolo 5, 35313 Padova, Italy.
| |
Collapse
|
6
|
Koukos PI, Dehghani-Ghahnaviyeh S, Velez-Vega C, Manchester J, Tieleman DP, Duca JS, Souza PCT, Cournia Z. Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications. J Chem Theory Comput 2023; 19:8901-8918. [PMID: 38019969 DOI: 10.1021/acs.jctc.3c00604] [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/01/2023]
Abstract
Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.
Collapse
Affiliation(s)
- Panagiotis I Koukos
- Biomedical Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Sepehr Dehghani-Ghahnaviyeh
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Camilo Velez-Vega
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - John Manchester
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - D Peter Tieleman
- Department of Biological Sciences, University of Calgary, Calgary T2N 1N4 Alberta, Canada
- Centre for Molecular Simulation, University of Calgary, Calgary T2N 1N4 Alberta, Canada
| | - José S Duca
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry, (MMSB, UMR 5086), CNRS & University of Lyon, 69367 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, 69364 Lyon, France
| | - Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Di Marino D, Conflitti P, Motta S, Limongelli V. Structural basis of dimerization of chemokine receptors CCR5 and CXCR4. Nat Commun 2023; 14:6439. [PMID: 37833254 PMCID: PMC10575954 DOI: 10.1038/s41467-023-42082-z] [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] [Received: 10/10/2021] [Accepted: 09/28/2023] [Indexed: 10/15/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.
Collapse
Affiliation(s)
- Daniele Di Marino
- Department of Life and Environmental Sciences - New York-Marche Structural Biology Centre (NY-MaSBiC), Polytechnic University of Marche, Via Brecce Bianche, 60131, Ancona, Italy
- Neuronal Death and Neuroprotection Unit, Department of Neuroscience, Mario Negri Institute for Pharmacological Research-IRCCS, Via Mario Negri 2, 20156, Milan, Italy
- National Biodiversity Future Center (NBFC), Palermo, Italy
| | - Paolo Conflitti
- Università della Svizzera italiana (USI), Faculty of Biomedical Sciences, Euler Institute, Via G. Buffi 13, CH-6900, Lugano, Switzerland
| | - Stefano Motta
- Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 1, 20126, Milan, Italy
| | - Vittorio Limongelli
- Università della Svizzera italiana (USI), Faculty of Biomedical Sciences, Euler Institute, Via G. Buffi 13, CH-6900, Lugano, Switzerland.
| |
Collapse
|
10
|
Gasparello J, Marzaro G, Papi C, Gentili V, Rizzo R, Zurlo M, Scapoli C, Finotti A, Gambari R. Effects of Sulforaphane on SARS‑CoV‑2 infection and NF‑κB dependent expression of genes involved in the COVID‑19 'cytokine storm'. Int J Mol Med 2023; 52:76. [PMID: 37477130 PMCID: PMC10555481 DOI: 10.3892/ijmm.2023.5279] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Accepted: 06/09/2023] [Indexed: 07/22/2023] Open
Abstract
Since its spread at the beginning of 2020, the coronavirus disease 2019 (COVID‑19) pandemic represents one of the major health problems. Despite the approval, testing, and worldwide distribution of anti‑severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) vaccines, the development of specific antiviral agents targeting the SARS‑CoV‑2 life cycle with high efficiency, and/or interfering with the associated 'cytokine storm', is highly required. A recent study, conducted by the authors' group indicated that sulforaphane (SFN) inhibits the expression of IL‑6 and IL‑8 genes induced by the treatment of IB3‑1 bronchial cells with a recombinant spike protein of SARS‑CoV‑2. In the present study, the ability of SFN to inhibit SARS‑CoV‑2 replication and the expression of pro‑inflammatory genes encoding proteins of the COVID‑19 'cytokine storm' was evaluated. SARS‑CoV‑2 replication was assessed in bronchial epithelial Calu‑3 cells. Moreover, SARS‑CoV‑2 replication and expression of pro‑inflammatory genes was evaluated by reverse transcription quantitative droplet digital PCR. The effects on the expression levels of NF‑κB were assessed by western blotting. Molecular dynamics simulations of NF‑kB/SFN interactions were conducted with Gromacs 2021.1 software under the Martini 2 CG force field. Computational studies indicated that i) SFN was stably bound with the NF‑κB monomer; ii) a ternary NF‑kB/SFN/DNA complex was formed; iii) the SFN interacted with both the protein and the nucleic acid molecules modifying the binding mode of the latter, and impairing the full interaction between the NF‑κB protein and the DNA molecule. This finally stabilized the inactive complex. Molecular studies demonstrated that SFN i) inhibits the SARS‑CoV‑2 replication in infected Calu‑3 cells, decreasing the production of the N‑protein coding RNA sequences, ii) decreased NF‑κB content in SARS‑CoV‑2 infected cells and inhibited the expression of NF‑kB‑dependent IL‑1β and IL‑8 gene expression. The data obtained in the present study demonstrated inhibitory effects of SFN on the SARS‑CoV‑2 life cycle and on the expression levels of the pro‑inflammatory genes, sustaining the possible use of SFN in the management of patients with COVID‑19.
Collapse
Affiliation(s)
- Jessica Gasparello
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| | - Giovanni Marzaro
- Department of Pharmaceutical and Pharmacological Sciences, University of Padova, I-35131 Padova
| | - Chiara Papi
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| | - Valentina Gentili
- Department of Chemical and Pharmaceutical Sciences, University of Ferrara, I-44121 Ferrara, Italy
| | - Roberta Rizzo
- Department of Chemical and Pharmaceutical Sciences, University of Ferrara, I-44121 Ferrara, Italy
| | - Matteo Zurlo
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| | - Chiara Scapoli
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| | - Alessia Finotti
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| | - Roberto Gambari
- Department of Life Sciences and Biotechnology, University of Ferrara, I-44121 Ferrara
| |
Collapse
|
11
|
Spinti JK, Neiva Nunes F, Melo MN. Room for improvement in the initial martini 3 parameterization of peptide interactions. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/01/2023]
|
12
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
13
|
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
|
14
|
Pezeshkian W, Grünewald F, Narykov O, Lu S, Arkhipova V, Solodovnikov A, Wassenaar TA, Marrink SJ, Korkin D. Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling. Structure 2023; 31:492-503.e7. [PMID: 36870335 DOI: 10.1016/j.str.2023.02.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 11/15/2022] [Accepted: 02/07/2023] [Indexed: 03/06/2023]
Abstract
Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.
Collapse
Affiliation(s)
- Weria Pezeshkian
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
| | - Fabian Grünewald
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands
| | - Oleksandr Narykov
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | - Senbao Lu
- Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | | | | | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Institute for Life Science and Technology, Hanze University of Applied Sciences, 9747AS Groningen, the Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands.
| | - Dmitry Korkin
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA; Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA.
| |
Collapse
|
15
|
Chiariello M, Grünewald F, Zarmiento-Garcia R, Marrink SJ. pH-Dependent Conformational Switch Impacts Stability of the PsbS Dimer. J Phys Chem Lett 2023; 14:905-911. [PMID: 36662680 PMCID: PMC9900633 DOI: 10.1021/acs.jpclett.2c03760] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2022] [Accepted: 01/13/2023] [Indexed: 06/17/2023]
Abstract
The photosystem II PsbS protein triggers the photoprotective mechanism of plants by sensing the acidification of the thylakoid lumen. Despite the mechanism of action of PsbS would require a pH-dependent monomerization of the dimeric form, a clear connection between the pH-induced structural changes and the dimer stability is missing. Here, by applying constant pH coarse-grained and all-atom molecular dynamics simulations, we investigate the pH-dependent structural response of the PsbS dimer. We find that the pH variation leads to structural changes in the lumen-exposed helices, located at the dimeric interface, providing an effective switch between PsbS inactive and active form. Moreover, the monomerization free energies reveal that in the neutral pH conformation, where the network of H-bond interactions at the dimeric interface is destroyed, the protein-protein interaction is weaker. Our results show how the pH-dependent conformations of PsbS affect their dimerization propensity, which is at the basis of the photoprotective mechanism.
Collapse
Affiliation(s)
| | | | - Rubi Zarmiento-Garcia
- Groningen
Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Nijenborgh 4, 9747
AG Groningen, The Netherlands
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Nijenborgh 4, 9747
AG Groningen, The Netherlands
| |
Collapse
|
16
|
Keller F, Alavizargar A, Wedlich-Söldner R, Heuer A. The impact of bilayer composition on the dimerization properties of the Slg1 stress sensor TMD from a multiscale analysis. Phys Chem Chem Phys 2023; 25:1299-1309. [PMID: 36533706 DOI: 10.1039/d2cp03497b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The impact of mutual interactions between the transmembrane domains of membrane proteins and lipids on bilayer properties has gained major attraction. Most simulation studies of membranes rely on the Martini force field, which has proven extremely helpful in providing molecular insights into realistic systems. Accordingly, an evaluation of the accuracy of the Martini force field is crucial to be able to correctly interpret the reported data. In this study, we combine atomistic and coarse-grained Martini simulations to investigate the properties of transmembrane domains (TMDs) in a model yeast membrane. The results show that the TMD binding state (monomeric and dimeric with positive or negative crossing angle) and the membrane composition significantly influence the properties around the TMDs and change TMD-TMD and TMD-lipid affinities. Furthermore, ergosterol (ERG) exhibits a strong affinity to TMD dimers. Importantly, the right-handed TMD dimer configuration is stabilized via TMD-TMD contacts by the addition of asymmetric anionic phosphatidylserine (PS). The coarse-grained simulations corroborate many of these findings, with two notable exceptions: a systematic overestimation of TMD-ERG interaction and lack of stabilization of the right-handed TMD dimers with the addition of PS.
Collapse
Affiliation(s)
- Fabian Keller
- Institut für Physikalische Chemie, Corrensstraße 28, Münster, Germany.
| | | | | | - Andreas Heuer
- Institut für Physikalische Chemie, Corrensstraße 28, Münster, Germany.
| |
Collapse
|
17
|
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
|
18
|
Darmawan KK, Karagiannis TC, Hughes JG, Small DM, Hung A. Molecular modeling of lactoferrin for food and nutraceutical applications: insights from in silico techniques. Crit Rev Food Sci Nutr 2022; 63:9074-9097. [PMID: 35503258 DOI: 10.1080/10408398.2022.2067824] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Lactoferrin is a protein, primarily found in milk that has attracted the interest of the food industries due to its health properties. Nevertheless, the instability of lactoferrin has limited its commercial application. Recent studies have focused on encapsulation to enhance the stability of lactoferrin. However, the molecular insights underlying the changes of structural properties of lactoferrin and the interaction with protectants remain poorly understood. Computational approaches have proven useful in understanding the structural properties of molecules and the key binding with other constituents. In this review, comprehensive information on the structure and function of lactoferrin and the binding with various molecules for food purposes are reviewed, with a special emphasis on the use of molecular dynamics simulations. The results demonstrate the application of modeling and simulations to determine key residues of lactoferrin responsible for its stability and interactions with other biomolecular components under various conditions, which are also associated with its functional benefits. These have also been extended into the potential creation of enhanced lactoferrin for commercial purposes. This review provides valuable strategies in designing novel nutraceuticals for food science practitioners and those who have interests in acquiring familiarity with the application of computational modeling for food and health purposes.
Collapse
Affiliation(s)
- Kevion K Darmawan
- School of Science, STEM College, RMIT University, Melbourne, Australia
| | - Tom C Karagiannis
- Epigenomic Medicine, Department of Diabetes, Central Clinical School, Monash University, Melbourne, Australia
- Department of Clinical Pathology, The University of Melbourne, Melbourne, Australia
| | - Jeff G Hughes
- School of Science, STEM College, RMIT University, Melbourne, Australia
| | - Darryl M Small
- School of Science, STEM College, RMIT University, Melbourne, Australia
| | - Andrew Hung
- School of Science, STEM College, RMIT University, Melbourne, Australia
| |
Collapse
|
19
|
Alavizargar A, Elting A, Wedlich-Söldner R, Heuer A. Lipid-Mediated Association of the Slg1 Transmembrane Domains in Yeast Plasma Membranes. J Phys Chem B 2022; 126:3240-3256. [PMID: 35446028 DOI: 10.1021/acs.jpcb.2c00192] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Clustering of transmembrane proteins underlies a multitude of fundamental biological processes at the plasma membrane (PM) such as receptor activation, lateral domain formation, and mechanotransduction. The self-association of the respective transmembrane domains (TMDs) has also been suggested to be responsible for the micron-scaled patterns seen for integral membrane proteins in the budding yeast PM. However, the underlying interplay between the local lipid composition and the TMD identity is still not mechanistically understood. In this work, we combined coarse-grained molecular dynamics simulations of simplified bilayer systems with high-resolution live-cell microscopy to analyze the distribution of a representative helical yeast TMD from the PM sensor Slg1 within different lipid environments. In our simulations, we specifically evaluated the effects of acyl chain saturation and anionic lipid head groups on the association of two TMDs. We found that weak lipid-protein interactions significantly affect the configuration of TMD dimers and the free energy of association. Increased amounts of unsaturated phospholipids (PLs) strongly reduced the helix-helix interaction, while the presence of anionic phosphatidylserine (PS) hardly affected the dimer formation. We could experimentally confirm this surprising lack of effect of PS using the network factor, a mesoscopic measure of PM pattern formation in yeast cells. Simulations also showed that the formation of TMD dimers in turn increased the order parameter of the surrounding lipids and induced long-range perturbations in lipid organization. In summary, our results shed new light on the mechanisms of lipid-mediated dimerization of TMDs in complex lipid mixtures.
Collapse
Affiliation(s)
- Azadeh Alavizargar
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| | - Annegret Elting
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Roland Wedlich-Söldner
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Andreas Heuer
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| |
Collapse
|
20
|
Salahub DR. Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond. Phys Chem Chem Phys 2022; 24:9051-9081. [PMID: 35389399 DOI: 10.1039/d1cp05928a] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.
Collapse
Affiliation(s)
- Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS-Centre for Molecular Simulation, IQST-Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
21
|
Thomasen FE, Pesce F, Roesgaard MA, Tesei G, Lindorff-Larsen K. Improving Martini 3 for Disordered and Multidomain Proteins. J Chem Theory Comput 2022; 18:2033-2041. [PMID: 35377637 DOI: 10.1021/acs.jctc.1c01042] [Citation(s) in RCA: 35] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.
Collapse
Affiliation(s)
- F Emil Thomasen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Francesco Pesce
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Mette Ahrensback Roesgaard
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Giulio Tesei
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
22
|
Wang J, Miao Y. Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics. J Chem Theory Comput 2022; 18:1275-1285. [PMID: 35099970 DOI: 10.1021/acs.jctc.1c00974] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
23
|
Challenges and frontiers of computational modelling of biomolecular recognition. QRB DISCOVERY 2022. [DOI: 10.1017/qrd.2022.11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Abstract
Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modelling. Here, we review the challenges and computational approaches developed to characterise biomolecular binding, including molecular docking, molecular dynamics simulations (especially enhanced sampling) and machine learning. Further improvements are still needed in order to accurately and efficiently characterise binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.
Collapse
|
24
|
Sawade K, Peter C. Multiscale simulations of protein and membrane systems. Curr Opin Struct Biol 2021; 72:203-208. [PMID: 34953308 DOI: 10.1016/j.sbi.2021.11.010] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 11/01/2021] [Accepted: 11/10/2021] [Indexed: 02/07/2023]
Abstract
Classical multiscale simulations are perfectly suited to investigate biological soft matter systems. Owing to the bridging between microscopically realistic and lower-resolution models or the integration of a hierarchy of subsystems, one gets access to biologically relevant system sizes and timescales. In recent years, increasingly complex systems and processes have come into focus such as multidomain proteins, phase separation processes in biopolymer solutions, multicomponent biomembranes, or multiprotein complexes up to entire viruses. The review shows factors that have contributed to this progress - from improved models to machine-learning-based analysis and scale-bridging methods.
Collapse
Affiliation(s)
- Kevin Sawade
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78 457, Konstanz, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78 457, Konstanz, Germany.
| |
Collapse
|
25
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|