1
|
Thomasen FE, Skaalum T, Kumar A, Srinivasan S, Vanni S, Lindorff-Larsen K. Rescaling protein-protein interactions improves Martini 3 for flexible proteins in solution. Nat Commun 2024; 15:6645. [PMID: 39103332 DOI: 10.1038/s41467-024-50647-9] [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: 05/29/2023] [Accepted: 07/15/2024] [Indexed: 08/07/2024] Open
Abstract
Multidomain proteins with flexible linkers and disordered regions play important roles in many cellular processes, but characterizing their conformational ensembles is difficult. We have previously shown that the coarse-grained model, Martini 3, produces too compact ensembles in solution, that may in part be remedied by strengthening protein-water interactions. Here, we show that decreasing the strength of protein-protein interactions leads to improved agreement with experimental data on a wide set of systems. We show that the 'symmetry' between rescaling protein-water and protein-protein interactions breaks down when studying interactions with or within membranes; rescaling protein-protein interactions better preserves the binding specificity of proteins with lipid membranes, whereas rescaling protein-water interactions preserves oligomerization of transmembrane helices. We conclude that decreasing the strength of protein-protein interactions improves the accuracy of Martini 3 for IDPs and multidomain proteins, both in solution and in the presence of a lipid membrane.
Collapse
Affiliation(s)
- F Emil Thomasen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200, Copenhagen N, Denmark.
| | - Tórur Skaalum
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200, Copenhagen N, Denmark
| | - Ashutosh Kumar
- Department of Biology, University of Fribourg, Fribourg, Switzerland
- Swiss National Center for Competence in Research (NCCR) Bio-inspired Materials, University of Fribourg, Chemin des Verdiers 4, CH-1700, Fribourg, Switzerland
| | | | - Stefano Vanni
- Department of Biology, University of Fribourg, Fribourg, Switzerland.
- Swiss National Center for Competence in Research (NCCR) Bio-inspired Materials, University of Fribourg, Chemin des Verdiers 4, CH-1700, Fribourg, Switzerland.
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200, Copenhagen N, Denmark.
| |
Collapse
|
2
|
Nemchinova M, Schuurman-Wolters GK, Whittaker JJ, Arkhipova V, Marrink SJ, Poolman B, Guskov A. Exploring the Ligand Binding and Conformational Dynamics of the Substrate-Binding Domain 1 of the ABC Transporter GlnPQ. J Phys Chem B 2024. [PMID: 39090964 DOI: 10.1021/acs.jpcb.4c02662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
The adenosine triphosphate (ATP)-binding cassette (ABC) importer GlnPQ from Lactococcus lactis has two sequential covalently linked substrate-binding domains (SBDs), which capture the substrates and deliver them to the translocon. The two SBDs differ in their ligand specificities, binding affinities and the distance to the transmembrane domain; interestingly, both SBDs can bind their ligands simultaneously without affecting each other. In this work, we studied the binding of ligands to both SBDs using X-ray crystallography and molecular dynamics simulations. We report three high-resolution structures of SBD1, namely, the wild-type SBD1 with bound asparagine or arginine, and E184D SBD1 with glutamine bound. Molecular dynamics (MD) simulations provide a detailed insight into the dynamics associated with open-closed transitions of the SBDs.
Collapse
Affiliation(s)
- Mariia Nemchinova
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Gea K Schuurman-Wolters
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Jacob J Whittaker
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Valentina Arkhipova
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Siewert J Marrink
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Bert Poolman
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Albert Guskov
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| |
Collapse
|
3
|
Chiariello MG, Zarmiento-Garcia R, Marrink SJ. Martini 3 Coarse-Grained Model for the Cofactors Involved in Photosynthesis. Int J Mol Sci 2024; 25:7947. [PMID: 39063190 PMCID: PMC11277265 DOI: 10.3390/ijms25147947] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2024] [Revised: 07/15/2024] [Accepted: 07/18/2024] [Indexed: 07/28/2024] Open
Abstract
As a critical step in advancing the simulation of photosynthetic complexes, we present the Martini 3 coarse-grained (CG) models of key cofactors associated with light harvesting (LHCII) proteins and the photosystem II (PSII) core complex. Our work focuses on the parametrization of beta-carotene, plastoquinone/quinol, violaxanthin, lutein, neoxanthin, chlorophyll A, chlorophyll B, and heme. We derived the CG parameters to match the all-atom reference simulations, while structural and thermodynamic properties of the cofactors were compared to experimental values when available. To further assess the reliability of the parameterization, we tested the behavior of these cofactors within their physiological environments, specifically in a lipid bilayer and bound to photosynthetic complexes. The results demonstrate that our CG models maintain the essential features required for realistic simulations. This work lays the groundwork for detailed simulations of the PSII-LHCII super-complex, providing a robust parameter set for future studies.
Collapse
Affiliation(s)
| | | | - Siewert-Jan Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands; (M.G.C.); (R.Z.-G.)
| |
Collapse
|
4
|
Saurabh S, Lei L, Li Z, Seddon JM, Lu JR, Kalonia C, Bresme F. Adsorption of monoclonal antibody fragments at the water-oil interface: A coarse-grained molecular dynamics study. APL Bioeng 2024; 8:026128. [PMID: 38948350 PMCID: PMC11211994 DOI: 10.1063/5.0207959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 06/06/2024] [Indexed: 07/02/2024] Open
Abstract
Monoclonal antibodies (mAbs) can undergo structural changes due to interaction with oil-water interfaces during storage. Such changes can lead to aggregation, resulting in a loss of therapeutic efficacy. Therefore, understanding the microscopic mechanism controlling mAb adsorption is crucial to developing strategies that can minimize the impact of interfaces on the therapeutic properties of mAbs. In this study, we used MARTINI coarse-grained molecular dynamics simulations to investigate the adsorption of the Fab and Fc domains of the monoclonal antibody COE3 at the oil-water interface. Our aim was to determine the regions on the protein surface that drive mAb adsorption. We also investigate the role of protein concentration on protein orientation and protrusion to the oil phase. While our structural analyses compare favorably with recent neutron reflectivity measurements, we observe some differences. Unlike the monolayer at the interface predicted by neutron reflectivity experiments, our simulations indicate the presence of a secondary diffused layer near the interface. We also find that under certain conditions, protein-oil interaction can lead to a considerable distortion in the protein structure, resulting in enhanced adsorption behavior.
Collapse
Affiliation(s)
- Suman Saurabh
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College, W12 0BZ London, United Kingdom
| | - Li Lei
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College, W12 0BZ London, United Kingdom
| | - Zongyi Li
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, Oxford Road, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - John M. Seddon
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College, W12 0BZ London, United Kingdom
| | - Jian R. Lu
- Biological Physics Group, School of Physics and Astronomy, Faculty of Science and Engineering, Oxford Road, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Cavan Kalonia
- Dosage Form Design and Development, BioPharmaceutical Development, BioPharmaceuticals R&D, AstraZeneca, Gaithersburg, Maryland 20878, USA
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College, W12 0BZ London, United Kingdom
| |
Collapse
|
5
|
Hu M, Dong X, Shi Q, Sun Y. Identification of a broad-spectrum high-affinity peptide ligand for the purification of spike proteins. J Chromatogr A 2024; 1723:464912. [PMID: 38643740 DOI: 10.1016/j.chroma.2024.464912] [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: 03/06/2024] [Revised: 04/07/2024] [Accepted: 04/11/2024] [Indexed: 04/23/2024]
Abstract
Since the outbreak of coronavirus disease 2019, the global demand for vaccines has increased rapidly to prevent infection and protect high-risk populations. However, identifying viral mutations poses an additional challenge for chromatographic purification of vaccines and subunit vaccines. In this study, a new affinity peptide model, X1VX2GLNX3WX4RYSK, was established, and a library of 612 peptides was generated for ligand screening. Based on a multistep strategy of ligand screening, 18 candidate peptides were obtained. The top ranking peptide, LP14 (YVYGLNIWLRYSK), and two other representative peptides, LP02 and LP06, with lower rankings were compared via molecular dynamics simulation. The results revealed that peptide binding to the receptor binding domain (RBD) was driven by hydrophobic interactions and the key residues involved in the binding were identified. Surface plasmon resonance analysis further confirmed that LP14 had the highest affinity for the wild RBD (Kd=0.520 μmol/L), and viral mutation had little influence on the affinity of LP14, demonstrating its great potential as a broad-spectrum ligand for RBD purification. Finally, chromatographic performance of LP14-coupled gel-packed column verified that both wild and omicron RBDs could be purified and were eluted by 0.1 mol/L Gly-HCl buffer (pH 3.0). This research identified a broad-spectrum peptide for RBD purification based on rational design and demonstrated its potential application in the purification of RBDs from complex feedstock.
Collapse
Affiliation(s)
- Mengke Hu
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China
| | - Xiaoyan Dong
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Key Laboratory of Systems Bioengineering and Frontiers Science Center for Synthetic Biology (Ministry of Education), Tianjin University, Tianjin 300350, China
| | - Qinghong Shi
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Key Laboratory of Systems Bioengineering and Frontiers Science Center for Synthetic Biology (Ministry of Education), Tianjin University, Tianjin 300350, China.
| | - Yan Sun
- Department of Biochemical Engineering, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China; Key Laboratory of Systems Bioengineering and Frontiers Science Center for Synthetic Biology (Ministry of Education), Tianjin University, Tianjin 300350, China.
| |
Collapse
|
6
|
Yang HW, Ju SP, Hsieh YT, Yang YC. Design single-stranded DNA aptamer of cluster of differentiation 47 protein by stochastic tunnelling-basin hopping-discrete molecular dynamics method. J Biomol Struct Dyn 2024; 42:3969-3982. [PMID: 37261868 DOI: 10.1080/07391102.2023.2217511] [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: 01/27/2023] [Accepted: 05/18/2023] [Indexed: 06/02/2023]
Abstract
The formation of the Cluster of Differentiation 47 (CD47, PDB code: 2JJT)/signal regulatory protein α (SIRPα) complex is very important as it protects healthy cells from immune clearance while promoting macrophage phagocytosis for tumour elimination. Although several antibodies have been developed for cancer therapy, new function-blocking aptamers are still under development. This study aims to design the aptamer AptCD47, which can block the formation of the CD47/SIRPα complex. This study employs the MARTINI coarse-grained (CG) force field and the stochastic tunnelling-basin hopping-discrete molecular dynamics (STUN-BH-DMD) method to identify the most stable AptCD47/CD47 complexes. Coarse-grained molecular dynamics (CGMD) simulations were used to obtain root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analyses. The results demonstrate that the formation of AptCD47/CD47 complexes renders the CD47 structure more stable than the single CD47 molecule in a water environment. The minimum energy pathway (MEP) obtained by the nudged elastic band (NEB) method indicates that the binding processes of 5'-ATTCAATTCC-3' and 5'-AGTGCAATCT-3' to CD47 are barrierless, which is much lower than the binding barrier of SIRPα to CD47 of about 14.23 kcal/mol. Therefore, these two AptCD47/CD47 complexes can create a high spatial binding barrier for SIRPα, preventing the formation of a stable CD47/SIRPα complex. The proposed numerical process with the MARTINI CG force field can be used to design CD47 aptamers that efficiently block SIRPα from binding to CD47.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hung-Wei Yang
- Department of Biomedical Engineering, National Cheng Kung University, Tainan City, Taiwan
- Medical Device Innovation Center, National Cheng Kung University, Tainan City, Taiwan
| | - Shin-Pon Ju
- Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan
- Department of Medicinal and Applied Chemistry, Kaohsiung Medical University, Kaohsiung, Taiwan
| | - Yun-Te Hsieh
- Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan
| | - Yung-Cheng Yang
- Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan
| |
Collapse
|
7
|
Schuhmann F, Ramsay JL, Kattnig DR, Solov’yov IA. Structural Rearrangements of Pigeon Cryptochrome 4 Undergoing a Complete Redox Cycle. J Phys Chem B 2024; 128:3844-3855. [PMID: 38568745 PMCID: PMC11056986 DOI: 10.1021/acs.jpcb.4c00424] [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/20/2024] [Revised: 03/04/2024] [Accepted: 03/11/2024] [Indexed: 04/26/2024]
Abstract
Cryptochrome is currently the major contender of a protein to underpin magnetoreception, the ability to sense the Earth's magnetic field. Among various types of cryptochromes, cryptochrome 4 has been identified as the likely magnetoreceptor in migratory birds. All-atom molecular dynamics (MD) studies have offered first insights into the structural dynamics of cryptochrome but are limited to a short time scale due to large computational demands. Here, we employ coarse-grained MD simulations to investigate the emergence of long-lived states and conformational changes in pigeon cryptochrome 4. Our coarse-grained simulations complete the picture by permitting observation on a significantly longer time scale. We observe conformational transitions in the phosphate-binding loop of pigeon cryptochrome 4 upon activation and identify prominent motions in residues 440-460, suggesting a possible role as a signaling state of the protein or as a gated interaction site for forming protein complexes that might facilitate downstream processes. The findings highlight the importance of considering longer time scales in studying cryptochrome dynamics and magnetoreception. Coarse-grained MD simulations offer a valuable tool to unravel the complex behavior of cryptochrome proteins and shed new light on the mechanisms underlying their role in magnetoreception. Further exploration of these conformational changes and their functional implications may contribute to a deeper understanding of the molecular mechanisms of magnetoreception in birds.
Collapse
Affiliation(s)
- Fabian Schuhmann
- Institute
of Physics, Carl von Ossietzky Universität
Oldenburg, Carl-von-Ossietzky Str. 9-11, Oldenburg 26129, Germany
- Niels
Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, Copenhagen 2100, Denmark
| | - Jessica L. Ramsay
- Living
Systems Institute and Department of Physics, University of Exeter, Stocker Rd., Exeter EX4
4QD, U.K.
| | - Daniel R. Kattnig
- Living
Systems Institute and Department of Physics, University of Exeter, Stocker Rd., Exeter EX4
4QD, U.K.
| | - Ilia A. Solov’yov
- Institute
of Physics, Carl von Ossietzky Universität
Oldenburg, Carl-von-Ossietzky Str. 9-11, Oldenburg 26129, Germany
- Research
Centre for Neurosensory Science, Carl von
Ossietzky Universität Oldenburg, Carl-von-Ossietzky-Str. 9-11, Oldenburg 26129, Germany
- Center
for Nanoscale Dynamics (CENAD), Carl von
Ossietzky Universität Oldenburg, Ammerländer Heerstr. 114-118, Oldenburg 26129, Germany
| |
Collapse
|
8
|
Rogoża NH, Krupa MA, Krupa P, Sieradzan AK. Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies. Molecules 2024; 29:1919. [PMID: 38731411 PMCID: PMC11085604 DOI: 10.3390/molecules29091919] [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: 03/26/2024] [Revised: 04/16/2024] [Accepted: 04/19/2024] [Indexed: 05/13/2024] Open
Abstract
Fullerenes, particularly C60, exhibit unique properties that make them promising candidates for various applications, including drug delivery and nanomedicine. However, their interactions with biomolecules, especially proteins, remain not fully understood. This study implements both explicit and implicit C60 models into the UNRES coarse-grained force field, enabling the investigation of fullerene-protein interactions without the need for restraints to stabilize protein structures. The UNRES force field offers computational efficiency, allowing for longer timescale simulations while maintaining accuracy. Five model proteins were studied: FK506 binding protein, HIV-1 protease, intestinal fatty acid binding protein, PCB-binding protein, and hen egg-white lysozyme. Molecular dynamics simulations were performed with and without C60 to assess protein stability and investigate the impact of fullerene interactions. Analysis of contact probabilities reveals distinct interaction patterns for each protein. FK506 binding protein (1FKF) shows specific binding sites, while intestinal fatty acid binding protein (1ICN) and uteroglobin (1UTR) exhibit more generalized interactions. The explicit C60 model shows good agreement with all-atom simulations in predicting protein flexibility, the position of C60 in the binding pocket, and the estimation of effective binding energies. The integration of explicit and implicit C60 models into the UNRES force field, coupled with recent advances in coarse-grained modeling and multiscale approaches, provides a powerful framework for investigating protein-nanoparticle interactions at biologically relevant scales without the need to use restraints stabilizing the protein, thus allowing for large conformational changes to occur. These computational tools, in synergy with experimental techniques, can aid in understanding the mechanisms and consequences of nanoparticle-biomolecule interactions, guiding the design of nanomaterials for biomedical applications.
Collapse
Affiliation(s)
- Natalia H. Rogoża
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Bażyńskiego 8, 80-309 Gdańsk, Poland; (N.H.R.); (M.A.K.); (A.K.S.)
| | - Magdalena A. Krupa
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Bażyńskiego 8, 80-309 Gdańsk, Poland; (N.H.R.); (M.A.K.); (A.K.S.)
| | - Pawel Krupa
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Bażyńskiego 8, 80-309 Gdańsk, Poland; (N.H.R.); (M.A.K.); (A.K.S.)
| |
Collapse
|
9
|
Dash R, Jabbari E. A Structure Independent Molecular Fragment Interfuse Model for Mesoscale Dissipative Particle Dynamics Simulation of Peptides. ACS OMEGA 2024; 9:18001-18022. [PMID: 38680324 PMCID: PMC11044228 DOI: 10.1021/acsomega.3c09534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 03/07/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
There is a need to develop robust computational models for mesoscale simulation of the structure of peptides over large length scales toward the discovery of novel peptides for medical applications to address the issues of peptide aggregation, enzymatic degradation, and short half-life. The primary objective was to predict the structure and conformation of peptides whose native structures are not known. This work presents a new model for computation of interaction parameters between the beads in coarse-grained dissipative particle dynamics (DPD) simulation that is properly calibrated for amino acids, supports compressibility requirement of water molecules, and accounts for subtle differences in the structure of amino acids and the charge in the side chain of charged amino acids. This new model is referred to as Structure Independent Molecular Fragment Interfuse Model, abbreviated as SIMFIM, because it accounts for specific interactions between different beads, which represent molecular fragments of the amino acids, in calculating nonbonded interaction parameters in the absence of knowing the actual peptide structure. The electrostatic interactions are incorporated in this model by using a normal distribution of charges around the center of the beads to prevent the collapse of oppositely charged soft beads. The uniquely parameterized DPD force field in the SIMFIM model is optimized for a given peptide with respect to the degree of coarse-grained graining for simulating the peptide over long times and length scales. The SIMFIM model was tested in this work using four peptides, namely, TrpZip2, Rubrivinodin, Lihuanodin, and IC3-CB1/Gai peptides, whose structures were sourced from the Protein Data Bank. The SIMFIM model predicted radius of gyration (Rg) values for the peptides closer to the actual structures as compared to the conventional model, and there was less deviation between the predicted and actual structures of the peptides.
Collapse
Affiliation(s)
- Ricky
Anshuman Dash
- Biomimetic Materials and
Tissue Engineering Laboratory, Chemical Engineering Department, University of South Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Esmaiel Jabbari
- Biomimetic Materials and
Tissue Engineering Laboratory, Chemical Engineering Department, University of South Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| |
Collapse
|
10
|
Grassmann G, Miotto M, Desantis F, Di Rienzo L, Tartaglia GG, Pastore A, Ruocco G, Monti M, Milanetti E. Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments. Chem Rev 2024; 124:3932-3977. [PMID: 38535831 PMCID: PMC11009965 DOI: 10.1021/acs.chemrev.3c00550] [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: 07/31/2023] [Revised: 02/20/2024] [Accepted: 02/21/2024] [Indexed: 04/11/2024]
Abstract
Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.
Collapse
Affiliation(s)
- Greta Grassmann
- Department
of Biochemical Sciences “Alessandro Rossi Fanelli”, Sapienza University of Rome, Rome 00185, Italy
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Mattia Miotto
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Fausta Desantis
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- The
Open University Affiliated Research Centre at Istituto Italiano di
Tecnologia, Genoa 16163, Italy
| | - Lorenzo Di Rienzo
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
| | - Gian Gaetano Tartaglia
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
- Center
for Human Technologies, Genoa 16152, Italy
| | - Annalisa Pastore
- Experiment
Division, European Synchrotron Radiation
Facility, Grenoble 38043, France
| | - Giancarlo Ruocco
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| | - Michele Monti
- RNA
System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
| | - Edoardo Milanetti
- Center
for Life Nano & Neuro Science, Istituto
Italiano di Tecnologia, Rome 00161, Italy
- Department
of Physics, Sapienza University, Rome 00185, Italy
| |
Collapse
|
11
|
Andreev G, Kovalenko M, Bozdaganyan ME, Orekhov PS. Colabind: A Cloud-Based Approach for Prediction of Binding Sites Using Coarse-Grained Simulations with Molecular Probes. J Phys Chem B 2024; 128:3211-3219. [PMID: 38514440 DOI: 10.1021/acs.jpcb.3c07853] [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: 03/23/2024]
Abstract
Binding site prediction is a crucial step in understanding protein-ligand and protein-protein interactions (PPIs) with broad implications in drug discovery and bioinformatics. This study introduces Colabind, a robust, versatile, and user-friendly cloud-based approach that employs coarse-grained molecular dynamics simulations in the presence of molecular probes, mimicking fragments of drug-like compounds. Our method has demonstrated high effectiveness when validated across a diverse range of biological targets spanning various protein classes, successfully identifying orthosteric binding sites, as well as known druggable allosteric or PPI sites, in both experimentally determined and AI-predicted protein structures, consistently placing them among the top-ranked sites. Furthermore, we suggest that careful inspection of the identified regions with a high affinity for specific probes can provide valuable insights for the development of pharmacophore hypotheses. The approach is available at https://github.com/porekhov/CG_probeMD.
Collapse
Affiliation(s)
- Georgy Andreev
- Insilico Medicine AI Ltd., Masdar City 145748, United Arab Emirates
| | - Max Kovalenko
- Division of Scientific Computing, Department of Information Technology, Uppsala University, Uppsala 752 37, Sweden
| | | | - Philipp S Orekhov
- Faculty of Biology, Shenzhen MSU-BIT University, Shenzhen 518172, China
| |
Collapse
|
12
|
Yang S, Song C. Switching Go̅ -Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions. J Chem Theory Comput 2024; 20:2618-2629. [PMID: 38447049 DOI: 10.1021/acs.jctc.3c01222] [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: 03/08/2024]
Abstract
Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and β2-adrenergic receptor (β2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.
Collapse
Affiliation(s)
- Song Yang
- Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Chen Song
- Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| |
Collapse
|
13
|
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
|
14
|
Plazinski W, Lutsyk V, Plazinska A. Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field. J Chem Theory Comput 2024; 20:2273-2283. [PMID: 38427574 PMCID: PMC10938637 DOI: 10.1021/acs.jctc.3c01155] [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: 10/19/2023] [Revised: 02/16/2024] [Accepted: 02/16/2024] [Indexed: 03/03/2024]
Abstract
Coarse-grained (CG) level molecular dynamics simulations are routinely used to study various biomolecular processes. The Martini force field is currently the most widely adopted parameter set for such simulations. The functional form of this and several other CG force fields enforces secondary protein structure support by employing a variety of harmonic potentials or restraints that favor the protein's native conformation. We propose a straightforward method to calculate the energetic consequences of transitions between predefined conformational states in systems in which multiple factors can affect protein conformational equilibria. This method is designed for use within the Martini force field and involves imposing conformational transitions by linking a Martini-inherent elastic network to the coupling parameter λ. We demonstrate the applicability of our method using the example of five biomolecular systems that undergo experimentally characterized conformational transitions between well-defined structures (Staphylococcal nuclease, C-terminal segment of surfactant protein B, LAH4 peptide, and β2-adrenergic receptor) as well as between folded and unfolded states (GCN4 leucine zipper protein). The results show that the relative free energy changes associated with protein conformational transitions, which are affected by various factors, such as pH, mutations, solvent, and lipid membrane composition, are correctly reproduced. The proposed method may be a valuable tool for understanding how different conditions and modifications affect conformational equilibria in proteins.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Jerzy
Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
- Department
of Biopharmacy, Medical University of Lublin, Chodzki 4a, Lublin 20-093, Poland
| | - Valery Lutsyk
- Jerzy
Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
| | - Anita Plazinska
- Department
of Biopharmacy, Medical University of Lublin, Chodzki 4a, Lublin 20-093, Poland
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
Smith ER, Theodorakis PE. Multiscale simulation of fluids: coupling molecular and continuum. Phys Chem Chem Phys 2024; 26:724-744. [PMID: 38113114 DOI: 10.1039/d3cp03579d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.
Collapse
Affiliation(s)
- Edward R Smith
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge, Middlesex UB8 3PH, UK.
| | | |
Collapse
|
18
|
Patel KN, Chavda D, Manna M. Molecular Docking of Intrinsically Disordered Proteins: Challenges and Strategies. Methods Mol Biol 2024; 2780:165-201. [PMID: 38987470 DOI: 10.1007/978-1-0716-3985-6_11] [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: 07/12/2024]
Abstract
Intrinsically disordered proteins (IDPs) are a novel class of proteins that have established a significant importance and attention within a very short period of time. These proteins are essentially characterized by their inherent structural disorder, encoded mainly by their amino acid sequences. The profound abundance of IDPs and intrinsically disordered regions (IDRs) in the biological world delineates their deep-rooted functionality. IDPs and IDRs convey such extensive functionality through their unique dynamic nature, which enables them to carry out huge number of multifaceted biomolecular interactions and make them "interaction hub" of the cellular systems. Additionally, with such widespread functions, their misfunctioning is also intimately associated with multiple diseases. Thus, understanding the dynamic heterogeneity of various IDPs along with their interactions with respective binding partners is an important field with immense potentials in biomolecular research. In this context, molecular docking-based computational approaches have proven to be remarkable in case of ordered proteins. Molecular docking methods essentially model the biomolecular interactions in both structural and energetic terms and use this information to characterize the putative interactions between the two participant molecules. However, direct applications of the conventional docking methods to study IDPs are largely limited by their structural heterogeneity and demands for unique IDP-centric strategies. Thus, in this chapter, we have presented an overview of current methodologies for successful docking operations involving IDPs and IDRs. These specialized methods majorly include the ensemble-based and fragment-based approaches with their own benefits and limitations. More recently, artificial intelligence and machine learning-assisted approaches are also used to significantly reduce the complexity and computational burden associated with various docking applications. Thus, this chapter aims to provide a comprehensive summary of major challenges and recent advancements of molecular docking approaches in the IDP field for their better utilization and greater applicability.Asp (D).
Collapse
Affiliation(s)
- Keyur N Patel
- Applied Phycology and Biotechnology Division, CSIR Central Salt and Marine Chemicals Research Institute, Bhavnagar, Gujarat, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, Uttar Pradesh, India
| | - Dhruvil Chavda
- Applied Phycology and Biotechnology Division, CSIR Central Salt and Marine Chemicals Research Institute, Bhavnagar, Gujarat, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, Uttar Pradesh, India
| | - Moutusi Manna
- Applied Phycology and Biotechnology Division, CSIR Central Salt and Marine Chemicals Research Institute, Bhavnagar, Gujarat, India.
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, Uttar Pradesh, India.
| |
Collapse
|
19
|
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
|
20
|
Santo KP, Neimark AV. Adsorption of pulmonary and exogeneous surfactants on SARS-CoV-2 spike protein. J Colloid Interface Sci 2023; 650:28-39. [PMID: 37392497 PMCID: PMC10279468 DOI: 10.1016/j.jcis.2023.06.121] [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: 03/23/2023] [Revised: 06/06/2023] [Accepted: 06/17/2023] [Indexed: 07/03/2023]
Abstract
COVID-19 is transmitted by airborne particles containing virions of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Coronavirus virions represent nanoparticles enveloped by a lipid bilayer decorated by a "crown" of Spike protein protrusions. Virus transmission into the cells is induced by binding of Spike proteins with ACE2 receptors of alveolar epithelial cells. Active clinical search is ongoing for exogenous surfactants and biologically active chemicals capable of hindering virion-receptor binding. Here, we explore by using coarse-grained molecular dynamics simulations the physico-chemical mechanisms of adsorption of selected pulmonary surfactants, zwitterionic dipalmitoyl phosphatidyl choline and cholesterol, and exogeneous anionic surfactant, sodium dodecyl sulfate, on the S1-domain of the Spike protein. We show that surfactants form micellar aggregates that selectively adhere to the specific regions of the S1-domain that are responsible for binding with ACE2 receptors. We find distinctly higher cholesterol adsorption and stronger cholesterol-S1 interactions in comparison with other surfactants, that is consistent with the experimental observations of the effects of cholesterol on COVID-19 infection. Distribution of adsorbed surfactant along the protein residue chain is highly specific and inhomogeneous with preferential adsorption around specific amino acid sequences. We observe preferential adsorption of surfactants on cationic arginine and lysine residues in the receptor-binding domain (RBD) that play an important role in ACE2 binding and are present in higher amounts in Delta and Omicron variants, which may lead to blocking direct Spike-ACE2 interactions. Our findings of strong selective adhesion of surfactant aggregates to Spike proteins have important implications for informing clinical search for therapeutic surfactants for curing and preventing COVID-19 caused by SARS-CoV-2 and its variants.
Collapse
Affiliation(s)
- Kolattukudy P Santo
- Department of Chemical and Biochemical Engineering, Rutgers, the State University of New Jersey, Piscataway, NJ 08854, USA
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers, the State University of New Jersey, Piscataway, NJ 08854, USA.
| |
Collapse
|
21
|
Laurent H, Hughes MDG, Walko M, Brockwell DJ, Mahmoudi N, Youngs TGA, Headen TF, Dougan L. Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering. Biomacromolecules 2023; 24:4869-4879. [PMID: 37874935 PMCID: PMC10646990 DOI: 10.1021/acs.biomac.3c00583] [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: 06/14/2023] [Revised: 10/03/2023] [Indexed: 10/26/2023]
Abstract
Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.
Collapse
Affiliation(s)
- Harrison Laurent
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
| | - Matt D. G. Hughes
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| | - Martin Walko
- School
of Chemistry, University of Leeds, Leeds, United
Kingdom, LS2 9JT
| | - David J. Brockwell
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| | - Najet Mahmoudi
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Tristan G. A. Youngs
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Thomas F. Headen
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Lorna Dougan
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| |
Collapse
|
22
|
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
|
23
|
Mishra S, van Aalst EJ, Wylie BJ, Brady LJ. Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes. Front Mol Biosci 2023; 10:1264454. [PMID: 37867558 PMCID: PMC10588454 DOI: 10.3389/fmolb.2023.1264454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 09/22/2023] [Indexed: 10/24/2023] Open
Abstract
YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.
Collapse
Affiliation(s)
- Surabhi Mishra
- Department of Oral Biology, University of Florida, Gainesville, FL, United States
| | - Evan J. van Aalst
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, United States
| | - Benjamin J. Wylie
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, United States
| | - L. Jeannine Brady
- Department of Oral Biology, University of Florida, Gainesville, FL, United States
| |
Collapse
|
24
|
Morstein J, Shrestha R, Van QN, López CA, Arora N, Tonelli M, Liang H, Chen D, Zhou Y, Hancock JF, Stephen AG, Turbyville TJ, Shokat KM. Direct Modulators of K-Ras-Membrane Interactions. ACS Chem Biol 2023; 18:2082-2093. [PMID: 37579045 PMCID: PMC10510109 DOI: 10.1021/acschembio.3c00413] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Accepted: 07/31/2023] [Indexed: 08/16/2023]
Abstract
Protein-membrane interactions (PMIs) are ubiquitous in cellular signaling. Initial steps of signal transduction cascades often rely on transient and dynamic interactions with the inner plasma membrane leaflet to populate and regulate signaling hotspots. Methods to target and modulate these interactions could yield attractive tool compounds and drug candidates. Here, we demonstrate that the conjugation of a medium-chain lipid tail to the covalent K-Ras(G12C) binder MRTX849 at a solvent-exposed site enables such direct modulation of PMIs. The conjugated lipid tail interacts with the tethered membrane and changes the relative membrane orientation and conformation of K-Ras(G12C), as shown by molecular dynamics (MD) simulation-supported NMR studies. In cells, this PMI modulation restricts the lateral mobility of K-Ras(G12C) and disrupts nanoclusters. The described strategy could be broadly applicable to selectively modulate transient PMIs.
Collapse
Affiliation(s)
- Johannes Morstein
- Department
of Cellular and Molecular Pharmacology and Howard Hughes Medical Institute, University of California, San Francisco, California 94158, United States
| | - Rebika Shrestha
- NCI
RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Que N. Van
- NCI
RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - César A. López
- Theoretical
Biology and Biophysics Group, Los Alamos
National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Neha Arora
- Department
of Integrative Biology and Pharmacology, McGovern Medical School, University of Texas Health Science Center, Houston, Texas 77030, United States
| | - Marco Tonelli
- National
Magnetic Resonance Facility at Madison, Biochemistry Department, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Hong Liang
- Department
of Integrative Biology and Pharmacology, McGovern Medical School, University of Texas Health Science Center, Houston, Texas 77030, United States
| | - De Chen
- NCI
RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Yong Zhou
- Department
of Integrative Biology and Pharmacology, McGovern Medical School, University of Texas Health Science Center, Houston, Texas 77030, United States
| | - John F. Hancock
- Department
of Integrative Biology and Pharmacology, McGovern Medical School, University of Texas Health Science Center, Houston, Texas 77030, United States
| | - Andrew G. Stephen
- NCI
RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Thomas J. Turbyville
- NCI
RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Kevan M. Shokat
- Department
of Cellular and Molecular Pharmacology and Howard Hughes Medical Institute, University of California, San Francisco, California 94158, United States
| |
Collapse
|
25
|
Wellawatte GP, Hocky GM, White AD. Neural potentials of proteins extrapolate beyond training data. J Chem Phys 2023; 159:085103. [PMID: 37642255 PMCID: PMC10474891 DOI: 10.1063/5.0147240] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Accepted: 07/31/2023] [Indexed: 08/31/2023] Open
Abstract
We evaluate neural network (NN) coarse-grained (CG) force fields compared to traditional CG molecular mechanics force fields. We conclude that NN force fields are able to extrapolate and sample from unseen regions of the free energy surface when trained with limited data. Our results come from 88 NN force fields trained on different combinations of clustered free energy surfaces from four protein mapped trajectories. We used a statistical measure named total variation similarity to assess the agreement between reference free energy surfaces from mapped atomistic simulations and CG simulations from trained NN force fields. Our conclusions support the hypothesis that NN CG force fields trained with samples from one region of the proteins' free energy surface can, indeed, extrapolate to unseen regions. Additionally, the force matching error was found to only be weakly correlated with a force field's ability to reconstruct the correct free energy surface.
Collapse
Affiliation(s)
- Geemi P. Wellawatte
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Glen M. Hocky
- Department of Chemistry, Simons Center for Computational Physical Chemistry, New York University, New York, New York 10003, USA
| | - Andrew D. White
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
26
|
Vishnoi S, Bhattacharya S, Walsh EM, Okoh GI, Thompson D. Computational Peptide Design Cotargeting Glucagon and Glucagon-like Peptide-1 Receptors. J Chem Inf Model 2023; 63:4934-4947. [PMID: 37523325 PMCID: PMC10428222 DOI: 10.1021/acs.jcim.3c00752] [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: 05/17/2023] [Indexed: 08/02/2023]
Abstract
Peptides are sustainable alternatives to conventional therapeutics for G protein-coupled receptor (GPCR) linked disorders, promising biocompatible and tailorable next-generation therapeutics for metabolic disorders including type-2 diabetes, as agonists of the glucagon receptor (GCGR) and the glucagon-like peptide-1 receptor (GLP-1R). However, single agonist peptides activating GLP-1R to stimulate insulin secretion also suppress obesity-linked glucagon release. Hence, bioactive peptides cotargeting GCGR and GLP-1R may remediate the blood glucose and fatty acid metabolism imbalance, tackling both diabetes and obesity to supersede current monoagonist therapy. Here, we design and model optimized peptide sequences starting from peptide sequences derived from earlier phage-displayed library screening, identifying those with predicted molecular binding profiles for dual agonism of GCGR and GLP-1R. We derive design rules from extensive molecular dynamics simulations based on peptide-receptor binding. Our newly designed coagonist peptide exhibits improved predicted coupled binding affinity for GCGR and GLP-1R relative to endogenous ligands and could in the future be tested experimentally, which may provide superior glycemic and weight loss control.
Collapse
Affiliation(s)
- Shubham Vishnoi
- Department
of Physics, Bernal Institute, University
of Limerick, Limerick V94T9PX, Ireland
| | - Shayon Bhattacharya
- Department
of Physics, Bernal Institute, University
of Limerick, Limerick V94T9PX, Ireland
| | | | | | - Damien Thompson
- Department
of Physics, Bernal Institute, University
of Limerick, Limerick V94T9PX, Ireland
| |
Collapse
|
27
|
Thallmair V, Schultz L, Evers S, Jolie T, Goecke C, Leitner MG, Thallmair S, Oliver D. Localization of the tubby domain, a PI(4,5)P2 biosensor, to E-Syt3-rich endoplasmic reticulum-plasma membrane junctions. J Cell Sci 2023; 136:jcs260848. [PMID: 37401342 PMCID: PMC10445746 DOI: 10.1242/jcs.260848] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 06/26/2023] [Indexed: 07/05/2023] Open
Abstract
The phospholipid phosphatidylinositol (4,5)-bisphosphate [PI(4,5)P2] acts as a signaling lipid at the plasma membrane (PM) with pleiotropic regulatory actions on multiple cellular processes. Signaling specificity might result from spatiotemporal compartmentalization of the lipid and from combinatorial binding of PI(4,5)P2 effector proteins to additional membrane components. Here, we analyzed the spatial distribution of tubbyCT, a paradigmatic PI(4,5)P2-binding domain, in live mammalian cells by total internal reflection fluorescence (TIRF) microscopy and molecular dynamics simulations. We found that unlike other well-characterized PI(4,5)P2 recognition domains, tubbyCT segregates into distinct domains within the PM. TubbyCT enrichment occurred at contact sites between PM and endoplasmic reticulum (ER) (i.e. at ER-PM junctions) as shown by colocalization with ER-PM markers. Localization to these sites was mediated in a combinatorial manner by binding to PI(4,5)P2 and by interaction with a cytosolic domain of extended synaptotagmin 3 (E-Syt3), but not other E-Syt isoforms. Selective localization to these structures suggests that tubbyCT is a novel selective reporter for a ER-PM junctional pool of PI(4,5)P2. Finally, we found that association with ER-PM junctions is a conserved feature of tubby-like proteins (TULPs), suggesting an as-yet-unknown function of TULPs.
Collapse
Affiliation(s)
- Veronika Thallmair
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps UniversityMarburg, 35037 Marburg, Germany
| | - Lea Schultz
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
| | - Saskia Evers
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
| | - Theresa Jolie
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
| | - Christian Goecke
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
| | - Michael G. Leitner
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
- Institute of Physiology, Department of Physiology and Medical Physics, Medical University of Innsbruck, 6020 Innsbruck, Austria
- Drug Discovery Sciences, Boehringer Ingelheim Pharma GmbH&Co.KG, Birkendorfer Str. 65, 88400 Biberach an der Riß, Germany
| | - Sebastian Thallmair
- Frankfurt Institute for Advanced Studies, 60438 Frankfurt am Main, Germany
- Groningen Biomolecular Sciences and Biotechnology Institute and The Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Dominik Oliver
- Department of Neurophysiology, Institute of Physiology and Pathophysiology, Philipps University Marburg, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps UniversityMarburg, 35037 Marburg, Germany
- Center for Mind, Brain and Behavior (CMBB), Universities of Marburg and Giessen, 35032 Marburg, Germany
| |
Collapse
|
28
|
MacCallum JL, Hu S, Lenz S, Souza PCT, Corradi V, Tieleman DP. An implementation of the Martini coarse-grained force field in OpenMM. Biophys J 2023; 122:2864-2870. [PMID: 37050876 PMCID: PMC10398343 DOI: 10.1016/j.bpj.2023.04.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 03/05/2023] [Accepted: 04/06/2023] [Indexed: 04/14/2023] Open
Abstract
We describe a complete implementation of Martini 2 and Martini 3 in the OpenMM molecular dynamics software package. Martini is a widely used coarse-grained force field with applications in biomolecular simulation, materials, and broader areas of chemistry. It is implemented as a force field but makes extensive use of facilities unique to the GROMACS software, including virtual sites and bonded terms that are not commonly used in standard atomistic force fields. OpenMM is a flexible molecular dynamics package widely used for methods development and is competitive in speed on GPUs with other commonly used packages. OpenMM has facilities to easily implement new force field terms, external forces and fields, and other nonstandard features, which we use to implement all force field terms used in Martini 2 and Martini 3. This allows Martini simulations, starting with GROMACS topology files that are processed by custom scripts, with all the added flexibility of OpenMM. We provide a GitHub repository with test cases, compare accuracy and performance between GROMACS and OpenMM, and discuss the limitations of our implementation in terms of direct comparison with GROMACS. We describe a use case that implements the Modeling Employing Limited Data method to apply experimental constraints in a Martini simulation to efficiently determine the structure of a protein complex. We also discuss issues and a potential solution with the Martini 2 topology for cholesterol.
Collapse
Affiliation(s)
- Justin L MacCallum
- Department of Chemistry, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada.
| | - Shangnong Hu
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - Stefan Lenz
- Department of Chemistry, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB - UMR 5086), CNRS & University of Lyon, Lyon, France
| | - Valentina Corradi
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - D Peter Tieleman
- Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada; Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada.
| |
Collapse
|
29
|
Su R, Zeng J, Marcink TC, Porotto M, Moscona A, O’Shaughnessy B. Host Cell Membrane Capture by the SARS-CoV-2 Spike Protein Fusion Intermediate. ACS CENTRAL SCIENCE 2023; 9:1213-1228. [PMID: 37396856 PMCID: PMC10255576 DOI: 10.1021/acscentsci.3c00158] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Indexed: 07/04/2023]
Abstract
Cell entry by SARS-CoV-2 is accomplished by the S2 subunit of the spike S protein on the virion surface by capture of the host cell membrane and fusion with the viral envelope. Capture and fusion require the prefusion S2 to transit to its potent fusogenic form, the fusion intermediate (FI). However, the FI structure is unknown, detailed computational models of the FI are unavailable, and the mechanisms and timing of membrane capture and fusion are not established. Here, we constructed a full-length model of the SARS-CoV-2 FI by extrapolating from known SARS-CoV-2 pre- and postfusion structures. In atomistic and coarse-grained molecular dynamics simulations the FI was remarkably flexible and executed giant bending and extensional fluctuations due to three hinges in the C-terminal base. The simulated configurations and their giant fluctuations are quantitatively consistent with SARS-CoV-2 FI configurations measured recently using cryo-electron tomography. Simulations suggested a host cell membrane capture time of ∼2 ms. Isolated fusion peptide simulations identified an N-terminal helix that directed and maintained binding to the membrane but grossly underestimated the binding time, showing that the fusion peptide environment is radically altered when attached to its host fusion protein. The large configurational fluctuations of the FI generated a substantial exploration volume that aided capture of the target membrane, and may set the waiting time for fluctuation-triggered refolding of the FI that draws the viral envelope and host cell membrane together for fusion. These results describe the FI as machinery that uses massive configurational fluctuations for efficient membrane capture and suggest novel potential drug targets.
Collapse
Affiliation(s)
- Rui Su
- Department
of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Jin Zeng
- Department
of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Tara C. Marcink
- Department
of Pediatrics, Columbia University Vagelos
College of Physicians & Surgeons, New York, New York 10032, United States
- Center
for Host−Pathogen Interaction, Columbia
University Vagelos College of Physicians & Surgeons, New York, New York 10032, United States
| | - Matteo Porotto
- Department
of Pediatrics, Columbia University Vagelos
College of Physicians & Surgeons, New York, New York 10032, United States
- Center
for Host−Pathogen Interaction, Columbia
University Vagelos College of Physicians & Surgeons, New York, New York 10032, United States
- Department
of Experimental Medicine, University of
Campania “Luigi Vanvitelli”, 81100 Caserta, Italy
| | - Anne Moscona
- Department
of Pediatrics, Columbia University Vagelos
College of Physicians & Surgeons, New York, New York 10032, United States
- Center
for Host−Pathogen Interaction, Columbia
University Vagelos College of Physicians & Surgeons, New York, New York 10032, United States
- Department
of Microbiology & Immunology, Columbia
University Vagelos College of Physicians & Surgeons, New York, New York 10032, United States
- Department
of Physiology, Columbia University Vagelos
College of Physicians & Surgeons, New York, New York 10032, United States
| | - Ben O’Shaughnessy
- Department
of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|
30
|
Sami S, Marrink SJ. Reactive Martini: Chemical Reactions in Coarse-Grained Molecular Dynamics Simulations. J Chem Theory Comput 2023. [PMID: 37327401 DOI: 10.1021/acs.jctc.2c01186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Chemical reactions are ubiquitous in both materials and the biophysical sciences. While coarse-grained (CG) molecular dynamics simulations are often needed to study the spatiotemporal scales present in these fields, chemical reactivity has not been explored thoroughly in CG models. In this work, a new approach to model chemical reactivity is presented for the widely used Martini CG Martini model. Employing tabulated potentials with a single extra particle for the angle dependence, the model provides a generic framework for capturing bonded topology changes using nonbonded interactions. As a first example application, the reactive model is used to study the macrocycle formation of benzene-1,3-dithiol molecules through the formation of disulfide bonds. We show that starting from monomers, macrocycles with sizes in agreement with experimental results are obtained using reactive Martini. Overall, our reactive Martini framework is general and can be easily extended to other systems. All of the required scripts and tutorials to explain its use are provided online.
Collapse
Affiliation(s)
- Selim Sami
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| |
Collapse
|
31
|
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: 10] [Impact Index Per Article: 10.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
|
32
|
Frangos ZJ, Wilson KA, Aitken HM, Cantwell Chater R, Vandenberg RJ, O'Mara ML. Membrane cholesterol regulates inhibition and substrate transport by the glycine transporter, GlyT2. Life Sci Alliance 2023; 6:e202201708. [PMID: 36690444 PMCID: PMC9873984 DOI: 10.26508/lsa.202201708] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2022] [Revised: 01/10/2023] [Accepted: 01/11/2023] [Indexed: 01/25/2023] Open
Abstract
Membrane cholesterol binds to and modulates the function of various SLC6 neurotransmitter transporters, including stabilizing the outward-facing conformation of the dopamine and serotonin transporters. Here, we investigate how cholesterol binds to GlyT2 (SLC6A5), modulates glycine transport rate, and influences bioactive lipid inhibition of GlyT2. Bioactive lipid inhibitors are analgesics that bind to an allosteric site accessible from the extracellular solution when GlyT2 adopts an outward-facing conformation. Using molecular dynamics simulations, mutagenesis, and cholesterol depletion experiments, we show that bioactive lipid inhibition of glycine transport is modulated by the recruitment of membrane cholesterol to a binding site formed by transmembrane helices 1, 5, and 7. Recruitment involves cholesterol flipping from its membrane orientation, and insertion of the 3' hydroxyl group into the cholesterol binding cavity, close to the allosteric site. The synergy between cholesterol and allosteric inhibitors provides a novel mechanism of inhibition and a potential avenue for the development of potent GlyT2 inhibitors as alternative therapeutics for the treatment of neuropathic pain and therapeutics that target other SLC6 transporters.
Collapse
Affiliation(s)
- Zachary J Frangos
- Molecular Biomedicine Theme, School of Medical Sciences, University of Sydney, Sydney, Australia
| | - Katie A Wilson
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australia
- Department of Biochemistry, Memorial University of Newfoundland, St. John's, Canada
| | - Heather M Aitken
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australia
- Australian Institute of Bioengineering and Nanotechnology, The University of Queensland, St Lucia, Australia
| | - Ryan Cantwell Chater
- Molecular Biomedicine Theme, School of Medical Sciences, University of Sydney, Sydney, Australia
| | - Robert J Vandenberg
- Molecular Biomedicine Theme, School of Medical Sciences, University of Sydney, Sydney, Australia
| | - Megan L O'Mara
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australia
- Australian Institute of Bioengineering and Nanotechnology, The University of Queensland, St Lucia, Australia
| |
Collapse
|
33
|
van Teijlingen A, Smith MC, Tuttle T. Short Peptide Self-Assembly in the Martini Coarse-Grain Force Field Family. Acc Chem Res 2023; 56:644-654. [PMID: 36866851 PMCID: PMC10035038 DOI: 10.1021/acs.accounts.2c00810] [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/04/2023]
Abstract
ConspectusPivotal to the success of any computational experiment is the ability to make reliable predictions about the system under study and the time required to yield these results. Biomolecular interactions is one area of research that sits in every camp of resolution vs the time required, from the quantum mechanical level to in vivo studies. At an approximate midpoint, there is coarse-grained molecular dynamics, for which the Martini force fields have become the most widely used, fast enough to simulate the entire membrane of a mitochondrion though lacking atom-specific precision. While many force fields have been parametrized to account for a specific system under study, the Martini force field has aimed at casting a wider net with more generalized bead types that have demonstrated suitability for broad use and reuse in applications from protein-graphene oxide coassembly to polysaccharides interactions.In this Account, the progressive (Martini versions 1 through 3) and peripheral (Sour Martini, constant pH, Martini Straight, Dry Martini, etc.) developmental trajectory of the Martini force field will be analyzed in terms of self-assembling systems with a focus on short (two to three amino acids) peptide self-assembly in aqueous environments. In particular, this will focus on the effects of the Martini solvent model and compare how changes in bead definitions and mapping have effects on different systems. Considerable effort in the development of Martini has been expended to reduce the "stickiness" of amino acids to better simulate proteins in bilayers. We have included in this Account a short study of dipeptide self-assembly in water, using all mainstream Martini force fields, to examine their ability to reproduce this behavior. The three most recently released versions of Martini and variations in their solvents are used to simulate in triplicate all 400 dipeptides of the 20 gene-encoded amino acids. The ability of the force fields to model the self-assembly of the dipeptides in aqueoues environments is determined by the measurement of the aggregation propensity, and additional descriptors are used to gain further insight into the dipeptide aggregates.
Collapse
Affiliation(s)
| | - Melissa C Smith
- Pure & Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, U.K
| | - Tell Tuttle
- Pure & Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, U.K
| |
Collapse
|
34
|
Duncan AL, Pezeshkian W. Mesoscale simulations: An indispensable approach to understand biomembranes. Biophys J 2023:S0006-3495(23)00123-6. [PMID: 36809878 DOI: 10.1016/j.bpj.2023.02.017] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 12/10/2022] [Accepted: 02/13/2023] [Indexed: 02/23/2023] Open
Abstract
Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.
Collapse
Affiliation(s)
- Anna L Duncan
- Department of Chemistry, Aarhus University, Aarhus C, Denmark.
| | - Weria Pezeshkian
- Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
35
|
Hilpert C, Beranger L, Souza PCT, Vainikka PA, Nieto V, Marrink SJ, Monticelli L, Launay G. Facilitating CG Simulations with MAD: The MArtini Database Server. J Chem Inf Model 2023; 63:702-710. [PMID: 36656159 DOI: 10.1021/acs.jcim.2c01375] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The MArtini Database (MAD - https://mad.ibcp.fr) is a web server designed for the sharing of structures and topologies of molecules parametrized with the Martini coarse-grained (CG) force field. MAD can also convert atomistic structures into CG structures and prepare complex systems (including proteins, lipids, etc.) for molecular dynamics (MD) simulations at the CG level. It is dedicated to the generation of input files for Martini 3, the most recent version of this popular CG force field. Specifically, the MAD server currently includes tools to submit or retrieve CG models of a wide range of molecules (lipids, carbohydrates, nanoparticles, etc.), transform atomistic protein structures into CG structures and topologies, with fine control on the process and assemble biomolecules into large systems, and deliver all files necessary to start simulations in the GROMACS MD engine.
Collapse
Affiliation(s)
- Cécile Hilpert
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Louis Beranger
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Paulo C T Souza
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Petteri A Vainikka
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Vincent Nieto
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Luca Monticelli
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| | - Guillaume Launay
- Microbiologie Moléculaire et Biochimie Structurale (MMSB), UMR 5086 CNRS & University of Lyon. 7 passage du Vercors, 69367 Lyon, France
| |
Collapse
|
36
|
Amyot R, Kodera N, Flechsig H. BioAFMviewer software for simulation atomic force microscopy of molecular structures and conformational dynamics. J Struct Biol X 2023; 7:100086. [PMID: 36865763 PMCID: PMC9972558 DOI: 10.1016/j.yjsbx.2023.100086] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Revised: 02/09/2023] [Accepted: 02/11/2023] [Indexed: 02/15/2023] Open
Abstract
Atomic force microscopy (AFM) and high-speed scanning have significantly advanced real time observation of biomolecular dynamics, with applications ranging from single molecules to the cellular level. To facilitate the interpretation of resolution-limited imaging, post-experimental computational analysis plays an increasingly important role to understand AFM measurements. Data-driven simulation of AFM, computationally emulating experimental scanning, and automatized fitting has recently elevated the understanding of measured AFM topographies by inferring the underlying full 3D atomistic structures. Providing an interactive user-friendly interface for simulation AFM, the BioAFMviewer software has become an established tool within the Bio-AFM community, with a plethora of applications demonstrating how the obtained full atomistic information advances molecular understanding beyond topographic imaging. This graphical review illustrates the BioAFMviewer capacities and further emphasizes the importance of simulation AFM to complement experimental observations.
Collapse
|
37
|
Gomes PSFC, Forrester M, Pace M, Gomes DEB, Bernardi RC. May the force be with you: The role of hyper-mechanostability of the bone sialoprotein binding protein during early stages of Staphylococci infections. Front Chem 2023; 11:1107427. [PMID: 36846849 PMCID: PMC9944720 DOI: 10.3389/fchem.2023.1107427] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Accepted: 01/25/2023] [Indexed: 02/10/2023] Open
Abstract
The bone sialoprotein-binding protein (Bbp) is a mechanoactive MSCRAMM protein expressed on the surface of Staphylococcus aureus that mediates adherence of the bacterium to fibrinogen-α (Fgα), a component of the bone and dentine extracellular matrix of the host cell. Mechanoactive proteins like Bbp have key roles in several physiological and pathological processes. Particularly, the Bbp: Fgα interaction is important in the formation of biofilms, an important virulence factor of pathogenic bacteria. Here, we investigated the mechanostability of the Bbp: Fgα complex using in silico single-molecule force spectroscopy (SMFS), in an approach that combines results from all-atom and coarse-grained steered molecular dynamics (SMD) simulations. Our results show that Bbp is the most mechanostable MSCRAMM investigated thus far, reaching rupture forces beyond the 2 nN range in typical experimental SMFS pulling rates. Our results show that high force-loads, which are common during initial stages of bacterial infection, stabilize the interconnection between the protein's amino acids, making the protein more "rigid". Our data offer new insights that are crucial on the development of novel anti-adhesion strategies.
Collapse
Affiliation(s)
- Priscila S. F. C. Gomes
- Department of Physics, College of Sciences and Mathematics, Auburn University, Auburn, AL, United States
| | - Meredith Forrester
- Department of Physics, College of Sciences and Mathematics, Auburn University, Auburn, AL, United States
| | - Margaret Pace
- Department of Physics, College of Sciences and Mathematics, Auburn University, Auburn, AL, United States
| | - Diego E. B. Gomes
- Department of Physics, College of Sciences and Mathematics, Auburn University, Auburn, AL, United States
| | | |
Collapse
|
38
|
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
|
39
|
Prakash S, Krishna A, Sengupta D. Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding. Chemphyschem 2023; 24:e202200509. [PMID: 36200760 DOI: 10.1002/cphc.202200509] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 10/06/2022] [Indexed: 02/04/2023]
Abstract
The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP2 and PIP3 lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.
Collapse
Affiliation(s)
- Shikha Prakash
- CSIR - National Chemical Laboratory, Dr. Homi Bhabha Road, Pune, 411008, India
| | - Anjali Krishna
- CSIR - National Chemical Laboratory, Dr. Homi Bhabha Road, Pune, 411008, India.,Current Address: School of Biological Sciences, University of Auckland, Auckland, 1010, New Zealand
| | - Durba Sengupta
- CSIR - National Chemical Laboratory, Dr. Homi Bhabha Road, Pune, 411008, India
| |
Collapse
|
40
|
Stevens JA, Grünewald F, van Tilburg PAM, König M, Gilbert BR, Brier TA, Thornburg ZR, Luthey-Schulten Z, Marrink SJ. Molecular dynamics simulation of an entire cell. Front Chem 2023; 11:1106495. [PMID: 36742032 PMCID: PMC9889929 DOI: 10.3389/fchem.2023.1106495] [Citation(s) in RCA: 26] [Impact Index Per Article: 26.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 01/09/2023] [Indexed: 01/19/2023] Open
Abstract
The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell's components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell's spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future.
Collapse
Affiliation(s)
- Jan A. Stevens
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Fabian Grünewald
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - P. A. Marco van Tilburg
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Melanie König
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Benjamin R. Gilbert
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Troy A. Brier
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Zane R. Thornburg
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Zaida Luthey-Schulten
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Champaign, IL, United States
| | - Siewert J. Marrink
- Molecular Dynamics Group, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands,*Correspondence: Siewert J. Marrink,
| |
Collapse
|
41
|
Vuillemot R, Mirzaei A, Harastani M, Hamitouche I, Fréchin L, Klaholz BP, Miyashita O, Tama F, Rouiller I, Jonic S. MDSPACE: Extracting Continuous Conformational Landscapes from Cryo-EM Single Particle Datasets Using 3D-to-2D Flexible Fitting based on Molecular Dynamics Simulation. J Mol Biol 2023; 435:167951. [PMID: 36638910 DOI: 10.1016/j.jmb.2023.167951] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2022] [Revised: 12/08/2022] [Accepted: 01/03/2023] [Indexed: 01/12/2023]
Abstract
This article presents an original approach for extracting atomic-resolution landscapes of continuous conformational variability of biomolecular complexes from cryo electron microscopy (cryo-EM) single particle images. This approach is based on a new 3D-to-2D flexible fitting method, which uses molecular dynamics (MD) simulation and is embedded in an iterative conformational-landscape refinement scheme. This new approach is referred to as MDSPACE, which stands for Molecular Dynamics simulation for Single Particle Analysis of Continuous Conformational hEterogeneity. The article describes the MDSPACE approach and shows its performance using synthetic and experimental datasets.
Collapse
Affiliation(s)
- Rémi Vuillemot
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France; Department of Biochemistry & Pharmacology and Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Victoria, Australia
| | - Alex Mirzaei
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Mohamad Harastani
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Ilyes Hamitouche
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Léo Fréchin
- Centre for Integrative Biology, Department of Integrated Structural Biology, IGBMC-UMR 7104 CNRS, U964 Inserm, Université de Strasbourg, Strasbourg, France
| | - Bruno P Klaholz
- Centre for Integrative Biology, Department of Integrated Structural Biology, IGBMC-UMR 7104 CNRS, U964 Inserm, Université de Strasbourg, Strasbourg, France
| | | | - Florence Tama
- RIKEN Center for Computational Science, Kobe, Japan; Institute of Transformative Biomolecules, Graduate School of Science, Nagoya University, Nagoya, Japan; Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
| | - Isabelle Rouiller
- Department of Biochemistry & Pharmacology and Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Victoria, Australia
| | - Slavica Jonic
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France.
| |
Collapse
|
42
|
Waclawiková B, Cesar Telles de Souza P, Schwalbe M, Neochoritis CG, Hoornenborg W, Nelemans SA, Marrink SJ, El Aidy S. Potential binding modes of the gut bacterial metabolite, 5-hydroxyindole, to the intestinal L-type calcium channels and its impact on the microbiota in rats. Gut Microbes 2023; 15:2154544. [PMID: 36511640 PMCID: PMC9754111 DOI: 10.1080/19490976.2022.2154544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Intestinal microbiota and microbiota-derived metabolites play a key role in regulating the host physiology. Recently, we have identified a gut-bacterial metabolite, namely 5-hydroxyindole, as a potent stimulant of intestinal motility via its modulation of L-type voltage-gated calcium channels located on the intestinal smooth muscle cells. Dysregulation of L-type voltage-gated calcium channels is associated with various gastrointestinal motility disorders, including constipation, making L-type voltage-gated calcium channels an important target for drug development. Nonetheless, the majority of currently available drugs are associated with alteration of the gut microbiota. Using 16S rRNA sequencing this study shows that, when administered orally, 5-hydroxyindole has only marginal effects on the rat cecal microbiota. Molecular dynamics simulations propose potential-binding pockets of 5-hydroxyindole in the α1 subunit of the L-type voltage-gated calcium channels and when its stimulatory effect on the rat colonic contractility was compared to 16 different analogues, ex-vivo, 5-hydroxyindole stood as the most potent enhancer of the intestinal contractility. Overall, the present findings imply a potential role of microbiota-derived metabolites as candidate therapeutics for targeted treatment of slow intestinal motility-related disorders including constipation.
Collapse
Affiliation(s)
- Barbora Waclawiková
- Host-Microbe Metabolic Interactions, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Groningen, The Netherlands
| | - Paulo Cesar Telles de Souza
- Molecular Microbiology and Structural Biochemistry (MMSB - UMR 5086), CNRS & University of Lyon, Lyon, France
| | - Markus Schwalbe
- Host-Microbe Metabolic Interactions, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Groningen, The Netherlands
| | | | - Warner Hoornenborg
- Department of Behavioral Neurosciences, Cluster Neurobiology, Groningen Institute of for Evolutionary Life Sciences (GELIFES), University of Groningen, Groningen, The Netherlands
| | - Sieger A. Nelemans
- Department of Molecular Neurobiology, Groningen Institute for Evolutionary Life Sciences (GELIFES), University of Groningen, Groningen, The Netherlands
| | - Siewert J. Marrink
- Molecular Dynamics, Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - Sahar El Aidy
- Host-Microbe Metabolic Interactions, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Groningen, The Netherlands,CONTACT Sahar El Aidy Host-Microbe Metabolic Interactions, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, Groningen, The Netherlands
| |
Collapse
|
43
|
Jensen LE, Rao S, Schuschnig M, Cada AK, Martens S, Hummer G, Hurley JH. Membrane curvature sensing and stabilization by the autophagic LC3 lipidation machinery. SCIENCE ADVANCES 2022; 8:eadd1436. [PMID: 36516251 PMCID: PMC9750143 DOI: 10.1126/sciadv.add1436] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Accepted: 11/10/2022] [Indexed: 05/28/2023]
Abstract
How the highly curved phagophore membrane is stabilized during autophagy initiation is a major open question in autophagosome biogenesis. Here, we use in vitro reconstitution on membrane nanotubes and molecular dynamics simulations to investigate how core autophagy proteins in the LC3 (Microtubule-associated proteins 1A/1B light chain 3) lipidation cascade interact with curved membranes, providing insight into their possible roles in regulating membrane shape during autophagosome biogenesis. ATG12(Autophagy-related 12)-ATG5-ATG16L1 was up to 100-fold enriched on highly curved nanotubes relative to flat membranes. At high surface density, ATG12-ATG5-ATG16L1 binding increased the curvature of the nanotubes. While WIPI2 (WD repeat domain phosphoinositide-interacting protein 2) binding directs membrane recruitment, the amphipathic helix α2 of ATG16L1 is responsible for curvature sensitivity. Molecular dynamics simulations revealed that helix α2 of ATG16L1 inserts shallowly into the membrane, explaining its curvature-sensitive binding to the membrane. These observations show how the binding of the ATG12-ATG5-ATG16L1 complex to the early phagophore rim could stabilize membrane curvature and facilitate autophagosome growth.
Collapse
Affiliation(s)
- Liv E. Jensen
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA 94720, USA
- California Institute for Quantitative Biosciences, University of California, Berkeley, CA 94720, USA
- Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD 20815, USA
| | - Shanlin Rao
- Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD 20815, USA
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main, Germany
| | - Martina Schuschnig
- Department of Biochemistry and Cell Biology, Max Perutz Labs, University of Vienna, Vienna BioCenter, Vienna, Austria
| | - A. King Cada
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA 94720, USA
- California Institute for Quantitative Biosciences, University of California, Berkeley, CA 94720, USA
| | - Sascha Martens
- Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD 20815, USA
- Department of Biochemistry and Cell Biology, Max Perutz Labs, University of Vienna, Vienna BioCenter, Vienna, Austria
| | - Gerhard Hummer
- Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD 20815, USA
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main, Germany
- Institute of Biophysics, Goethe University Frankfurt, Frankfurt am Main 60438, Germany
| | - James H. Hurley
- Department of Molecular and Cell Biology, University of California Berkeley, Berkeley, CA 94720, USA
- California Institute for Quantitative Biosciences, University of California, Berkeley, CA 94720, USA
- Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD 20815, USA
- Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, CA 94720, USA
| |
Collapse
|
44
|
Ahmadpour-Samani P, Zahedi P. An investigation on nematic-isotropic phase transition, viscosity and diffusion coefficient of liquid crystalline elastomers at different temperatures using molecular dynamics simulation. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
45
|
Sang P, Chen YQ, Liu MT, Wang YT, Yue T, Li Y, Yin YR, Yang LQ. Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants. Int J Mol Sci 2022; 23:ijms232314796. [PMID: 36499120 PMCID: PMC9740405 DOI: 10.3390/ijms232314796] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 11/18/2022] [Accepted: 11/23/2022] [Indexed: 12/02/2022] Open
Abstract
To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBDWT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBDOMIs) feature increased global structural fluctuations when compared with RBDWT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBDOMIs than RBDWT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBDBA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBDBA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.
Collapse
Affiliation(s)
- Peng Sang
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
- Key Laboratory of Bioinformatics and Computational Biology, Department of Education of Yunnan Province, Dali University, Dali 671000, China
- Yunnan Key Laboratory of Screening and Research on Anti-Pathogenic Plant Resources from West Yunnan, Dali University, Dali 671000, China
| | - Yong-Qin Chen
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
| | - Meng-Ting Liu
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
| | - Yu-Ting Wang
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
| | - Ting Yue
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
| | - Yi Li
- College of Mathematics and Computer Science, Dali University, Dali 671000, China
| | - Yi-Rui Yin
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
| | - Li-Quan Yang
- College of Agriculture and Biological Science, Dali University, Dali 671000, China
- Key Laboratory of Bioinformatics and Computational Biology, Department of Education of Yunnan Province, Dali University, Dali 671000, China
- Yunnan Key Laboratory of Screening and Research on Anti-Pathogenic Plant Resources from West Yunnan, Dali University, Dali 671000, China
- Correspondence:
| |
Collapse
|
46
|
Thallmair V, Schultz L, Zhao W, Marrink SJ, Oliver D, Thallmair S. Two cooperative binding sites sensitize PI(4,5)P 2 recognition by the tubby domain. SCIENCE ADVANCES 2022; 8:eabp9471. [PMID: 36070381 PMCID: PMC9451155 DOI: 10.1126/sciadv.abp9471] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Accepted: 07/20/2022] [Indexed: 05/22/2023]
Abstract
Phosphoinositides (PIs) are lipid signaling molecules that operate by recruiting proteins to cellular membranes via PI recognition domains. The dominant PI of the plasma membrane is phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2]. One of only two PI(4,5)P2 recognition domains characterized in detail is the tubby domain. It is essential for targeting proteins into cilia involving reversible membrane association. However, the PI(4,5)P2 binding properties of tubby domains have remained enigmatic. Here, we used coarse-grained molecular dynamics simulations to explore PI(4,5)P2 binding by the prototypic tubby domain. The comparatively low PI(4,5)P2 affinity of the previously described canonical binding site is underpinned in a cooperative manner by a previously unknown, adjacent second binding site. Mutations in the previously unknown site impaired PI(4,5)P2-dependent plasma membrane localization in living cells and PI(4,5)P2 interaction in silico, emphasizing its importance for PI(4,5)P2 affinity. The two-ligand binding mode may serve to sharpen the membrane association-dissociation cycle of tubby-like proteins that underlies delivery of ciliary cargo.
Collapse
Affiliation(s)
- Veronika Thallmair
- Institute for Physiology and Pathophysiology, Philipps University Marburg, Deutschhausstr. 1-2, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps University Marburg, Marburg, Germany
| | - Lea Schultz
- Institute for Physiology and Pathophysiology, Philipps University Marburg, Deutschhausstr. 1-2, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps University Marburg, Marburg, Germany
| | - Wencai Zhao
- Institute for Physiology and Pathophysiology, Philipps University Marburg, Deutschhausstr. 1-2, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps University Marburg, Marburg, Germany
| | - Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Dominik Oliver
- Institute for Physiology and Pathophysiology, Philipps University Marburg, Deutschhausstr. 1-2, 35037 Marburg, Germany
- DFG Research Training Group, Membrane Plasticity in Tissue Development and Remodeling, GRK 2213, Philipps University Marburg, Marburg, Germany
- Corresponding author. (S.T.); (D.O.)
| | - Sebastian Thallmair
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
- Frankfurt Institute for Advanced Studies, Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany
- Corresponding author. (S.T.); (D.O.)
| |
Collapse
|
47
|
Towards Molecular Understanding of the Functional Role of UbiJ-UbiK2 Complex in Ubiquinone Biosynthesis by Multiscale Molecular Modelling Studies. Int J Mol Sci 2022; 23:ijms231810323. [PMID: 36142227 PMCID: PMC9499169 DOI: 10.3390/ijms231810323] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/30/2022] [Accepted: 09/01/2022] [Indexed: 11/17/2022] Open
Abstract
Ubiquinone (UQ) is a polyisoprenoid lipid found in the membranes of bacteria and eukaryotes. UQ has important roles, notably in respiratory metabolisms which sustain cellular bioenergetics. Most steps of UQ biosynthesis take place in the cytosol of E. coli within a multiprotein complex called the Ubi metabolon, that contains five enzymes and two accessory proteins, UbiJ and UbiK. The SCP2 domain of UbiJ was proposed to bind the hydrophobic polyisoprenoid tail of UQ biosynthetic intermediates in the Ubi metabolon. How the newly synthesised UQ might be released in the membrane is currently unknown. In this paper, we focused on better understanding the role of the UbiJ-UbiK2 heterotrimer forming part of the metabolon. Given the difficulties to gain functional insights using biophysical techniques, we applied a multiscale molecular modelling approach to study the UbiJ-UbiK2 heterotrimer. Our data show that UbiJ-UbiK2 interacts closely with the membrane and suggests possible pathways to enable the release of UQ into the membrane. This study highlights the UbiJ-UbiK2 complex as the likely interface between the membrane and the enzymes of the Ubi metabolon and supports that the heterotrimer is key to the biosynthesis of UQ8 and its release into the membrane of E. coli.
Collapse
|
48
|
López CA, Zhang X, Aydin F, Shrestha R, Van QN, Stanley CB, Carpenter TS, Nguyen K, Patel LA, Chen D, Burns V, Hengartner NW, Reddy TJE, Bhatia H, Di Natale F, Tran TH, Chan AH, Simanshu DK, Nissley DV, Streitz FH, Stephen AG, Turbyville TJ, Lightstone FC, Gnanakaran S, Ingólfsson HI, Neale C. Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework. J Chem Theory Comput 2022; 18:5025-5045. [PMID: 35866871 DOI: 10.1021/acs.jctc.2c00168] [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/28/2022]
Abstract
The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.
Collapse
Affiliation(s)
- Cesar A López
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Xiaohua Zhang
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Fikret Aydin
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Rebika Shrestha
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Que N Van
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Christopher B Stanley
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Timothy S Carpenter
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Kien Nguyen
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Lara A Patel
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - De Chen
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Violetta Burns
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nicolas W Hengartner
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tyler J E Reddy
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Harsh Bhatia
- Computing Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Francesco Di Natale
- Computing Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Timothy H Tran
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Albert H Chan
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Dhirendra K Simanshu
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Dwight V Nissley
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Frederick H Streitz
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Andrew G Stephen
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Thomas J Turbyville
- NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States
| | - Felice C Lightstone
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Sandrasegaram Gnanakaran
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Helgi I Ingólfsson
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Chris Neale
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
49
|
Magi Meconi G, Sasselli IR, Bianco V, Onuchic JN, Coluzza I. Key aspects of the past 30 years of protein design. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 85:086601. [PMID: 35704983 DOI: 10.1088/1361-6633/ac78ef] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 06/15/2022] [Indexed: 06/15/2023]
Abstract
Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.
Collapse
Affiliation(s)
- Giulia Magi Meconi
- Computational Biophysics Lab, Center for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramon 182, 20014, Donostia-San Sebastián, Spain
| | - Ivan R Sasselli
- Computational Biophysics Lab, Center for Cooperative Research in Biomaterials (CIC biomaGUNE), Basque Research and Technology Alliance (BRTA), Paseo de Miramon 182, 20014, Donostia-San Sebastián, Spain
| | | | - Jose N Onuchic
- Center for Theoretical Biological Physics, Department of Physics & Astronomy, Department of Chemistry, Department of Biosciences, Rice University, Houston, TX 77251, United States of America
| | - Ivan Coluzza
- BCMaterials, Basque Center for Materials, Applications and Nanostructures, Bld. Martina Casiano, UPV/EHU Science Park, Barrio Sarriena s/n, 48940 Leioa, Spain
- Basque Foundation for Science, Ikerbasque, 48009, Bilbao, Spain
| |
Collapse
|
50
|
Sahoo A, Lee PY, Matysiak S. Transferable and Polarizable Coarse Grained Model for Proteins─ProMPT. J Chem Theory Comput 2022; 18:5046-5055. [PMID: 35793442 DOI: 10.1021/acs.jctc.2c00269] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The application of classical molecular dynamics (MD) simulations at atomic resolution (fine-grained level, FG), to most biomolecular processes, remains limited because of the associated computational complexity of representing all the atoms. This problem is magnified in the presence of protein-based biomolecular systems that have a very large conformational space, and MD simulations with fine-grained resolution have slow dynamics to explore this space. Current transferable coarse grained (CG) force fields in literature are either limited to only peptides with the environment encoded in an implicit form or cannot capture transitions into secondary/tertiary peptide structures from a primary sequence of amino acids. In this work, we present a transferable CG force field with an explicit representation of the environment for accurate simulations with proteins. The force field consists of a set of pseudoatoms representing different chemical groups that can be joined/associated together to create different biomolecular systems. This preserves the transferability of the force field to multiple environments and simulation conditions. We have added electronic polarization that can respond to environmental heterogeneity/fluctuations and couple it to protein's structural transitions. The nonbonded interactions are parametrized with physics-based features such as solvation and partitioning free energies determined by thermodynamic calculations and matched with experiments and/or atomistic simulations. The bonded potentials are inferred from corresponding distributions in nonredundant protein structure databases. We present validations of the CG model with simulations of well-studied aqueous protein systems with specific protein fold types─Trp-cage, Trpzip4, villin, WW-domain, and β-α-β. We also explore the applications of the force field to study aqueous aggregation of Aβ 16-22 peptides.
Collapse
Affiliation(s)
- Abhilash Sahoo
- Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Pei-Yin Lee
- Chemical Physics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Silvina Matysiak
- Fischell Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|