1
|
Majumder A, Straub JE. Machine Learning Derived Collective Variables for the Study of Protein Homodimerization in Membrane. J Chem Theory Comput 2024. [PMID: 38918177 DOI: 10.1021/acs.jctc.4c00454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
The accurate calculation of equilibrium constants for protein-protein association is of fundamental importance to quantitative biology and remains an outstanding challenge for computational biophysics. Traditionally, equilibrium constants have been computed from one-dimensional free energy surfaces derived from sampling along a single collective variable. Importantly, recent advances in enhanced sampling methodology have facilitated the characterization of multidimensional free energy landscapes, often exposing multiple thermodynamically important minima missed by more restrictive sampling methods. A key to the effectiveness of this multidimensional sampling approach is the identification of collective variables that effectively define the configurational space of dissociated and associated states. Here we present the application of two machine learning methods for the unbiased determination of collective variables for enhanced sampling and analysis of protein-protein association. Our results both validate prior work, based on intuition derived collective variables, and demonstrate the effectiveness of the machine learning methods for the identification of collective variables for association reactions in complex biomolecular systems.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
2
|
Haque Pial T, Li Y, Olvera de la Cruz M. Microscopically segregated ligand distribution in co-assembled peptide-amphiphile nanofibers. SOFT MATTER 2024; 20:4640-4647. [PMID: 38819791 DOI: 10.1039/d4sm00315b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Peptide amphiphiles (PAs) self-assemble into cylindrical nanofibers with applications in protein purification, tissue engineering, and regenerative medicine. For these applications, functionalized PAs are often co-assembled with oppositely charged filler PAs. Finding the conditions at which these fibers are homogeneously mixed or segregated is crucial for the required application. We co-assemble negative C12VVEE fillers and positive C12VVKK-OEG4-Z33 ligands, which are important for antibody purifications. Our results show that the ligands tend to cluster and locally segregate in the fiber surfaces. The Z33s are overall neutral and form large aggregates in bulk solution due to short range attractions. However, full segregation of the C12VVKK-OEG4-Z33 is not observed in the cylindrical surface due to the electrostatic penalty of forming large domains of similarly charged molecules. This is commensurate with previous theoretical predictions, showing that the competition between short-range attractive interactions and long-range electrostatic repulsions leads to pattern formation in cylindrical surfaces. This work offers valuable insight into the design of functionalized nanofibers for various biomedical and chemical applications.
Collapse
Affiliation(s)
- Turash Haque Pial
- Department of Materials Science and Engineering, Northwestern University, 2220 Campus Drive, Evanston, IL 60208, USA.
- Center of Computation and Theory of Soft Materials, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
| | - Yang Li
- Center of Computation and Theory of Soft Materials, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
| | - Monica Olvera de la Cruz
- Department of Materials Science and Engineering, Northwestern University, 2220 Campus Drive, Evanston, IL 60208, USA.
- Center of Computation and Theory of Soft Materials, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
- Department of Physics and Astronomy, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA
| |
Collapse
|
3
|
Ito S, Sugita Y. Free-energy landscapes of transmembrane homodimers by bias-exchange adaptively biased molecular dynamics. Biophys Chem 2024; 307:107190. [PMID: 38290241 DOI: 10.1016/j.bpc.2024.107190] [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: 12/03/2023] [Revised: 01/21/2024] [Accepted: 01/21/2024] [Indexed: 02/01/2024]
Abstract
Membrane proteins play essential roles in various biological functions within the cell. One of the most common functional regulations involves the dimerization of two single-pass transmembrane (TM) helices. Glycophorin A (GpA) and amyloid precursor protein (APP) form TM homodimers in the membrane, which have been investigated both experimentally and computationally. The homodimer structures are well characterized using only four collective variables (CVs) when each TM helix is stable. The CVs are the interhelical distance, the crossing angle, and the Crick angles for two TM helices. However, conformational sampling with multi-dimensional replica-exchange umbrella sampling (REUS) requires too many replicas to sample all the CVs for exploring the conformational landscapes. Here, we show that the bias-exchange adaptively biased molecular dynamics (BE-ABMD) with the four CVs effectively explores the free-energy landscapes of the TM helix dimers of GpA, wild-type APP and its mutants in the IMM1 implicit membrane. Compared to the original ABMD, the bias-exchange algorithm in BE-ABMD can provide a more rapidly converged conformational landscape. The BE-ABMD simulations could also reveal TM packing interfaces of the membrane proteins and the dependence of the free-energy landscapes on the membrane thickness. This approach is valuable for numerous other applications, including those involving explicit solvent and a lipid bilayer in all-atom force fields or Martini coarse-grained models, and enhances our understanding of protein-protein interactions in biological membranes.
Collapse
Affiliation(s)
- Shingo Ito
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| |
Collapse
|
4
|
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
|
5
|
Majumder A, Straub JE. The role of structural heterogeneity in the homodimerization of transmembrane proteins. J Chem Phys 2023; 159:134101. [PMID: 37782254 PMCID: PMC10547497 DOI: 10.1063/5.0159801] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Accepted: 07/27/2023] [Indexed: 10/03/2023] Open
Abstract
The equilibrium association of transmembrane proteins plays a fundamental role in membrane protein function and cellular signaling. While the study of the equilibrium binding of single pass transmembrane proteins has received significant attention in experiment and simulation, the accurate assessment of equilibrium association constants remains a challenge to experiment and simulation. In experiment, there remain wide variations in association constants derived from experimental studies of the most widely studied transmembrane proteins. In simulation, state-of-the art methods have failed to adequately sample the thermodynamically relevant structures of the dimer state ensembles using coarse-grained models. In addition, all-atom force fields often fail to accurately assess the relative free energies of the dimer and monomer states. Given the importance of this fundamental biophysical process, it is essential to address these shortcomings. In this work, we establish an effective computational protocol for the calculation of equilibrium association constants for transmembrane homodimer formation. A set of transmembrane protein homodimers, used in the parameterization of the MARTINI v3 force field, are simulated using metadynamics, based on three collective variables. The method is found to be accurate and computationally efficient, providing a standard to be used in the future simulation studies using coarse-grained or all-atom models.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - John E. Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|
6
|
张 民, 龚 铭, 王 进, 陈 振, 周 良. [Research progress of coarse-grained molecular dynamics in drug carrier materials]. SHENG WU YI XUE GONG CHENG XUE ZA ZHI = JOURNAL OF BIOMEDICAL ENGINEERING = SHENGWU YIXUE GONGCHENGXUE ZAZHI 2023; 40:799-804. [PMID: 37666772 PMCID: PMC10477382 DOI: 10.7507/1001-5515.202303008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Received: 03/04/2023] [Revised: 07/13/2023] [Indexed: 09/06/2023]
Abstract
As one of the traditional computer simulation techniques, molecular simulation can intuitively display and quantify molecular structure and explain experimental phenomena from the microscopic molecular level. When the simulation system increases, the amount of calculation will also increase, which will cause a great burden on the simulation system. Coarse-grained molecular dynamics is a method of mesoscopic molecular simulation, which can simplify the molecular structure and improve computational efficiency, as a result, coarse-grained molecular dynamics is often used when simulating macromolecular systems such as drug carrier materials. In this article, we reviewed the recent research results of using coarse-grained molecular dynamics to simulate drug carriers, in order to provide a reference for future pharmaceutical preparation research and accelerate the entry of drug research into the era of precision drug design.
Collapse
Affiliation(s)
- 民权 张
- 江西科技师范大学 药学院 江西省药物分子设计与评价重点实验室(南昌 330013)Jiangxi Provincial Key Laboratory of Drug Design and Evaluation, School of Pharmacy, Jiangxi Science and Technology Normal University, Nanchang 330013, P. R. China
| | - 铭城 龚
- 江西科技师范大学 药学院 江西省药物分子设计与评价重点实验室(南昌 330013)Jiangxi Provincial Key Laboratory of Drug Design and Evaluation, School of Pharmacy, Jiangxi Science and Technology Normal University, Nanchang 330013, P. R. China
| | - 进 王
- 江西科技师范大学 药学院 江西省药物分子设计与评价重点实验室(南昌 330013)Jiangxi Provincial Key Laboratory of Drug Design and Evaluation, School of Pharmacy, Jiangxi Science and Technology Normal University, Nanchang 330013, P. R. China
| | - 振华 陈
- 江西科技师范大学 药学院 江西省药物分子设计与评价重点实验室(南昌 330013)Jiangxi Provincial Key Laboratory of Drug Design and Evaluation, School of Pharmacy, Jiangxi Science and Technology Normal University, Nanchang 330013, P. R. China
| | - 良良 周
- 江西科技师范大学 药学院 江西省药物分子设计与评价重点实验室(南昌 330013)Jiangxi Provincial Key Laboratory of Drug Design and Evaluation, School of Pharmacy, Jiangxi Science and Technology Normal University, Nanchang 330013, P. R. China
| |
Collapse
|
7
|
Sahoo AR, Souza PCT, Meng Z, Buck M. Transmembrane dimers of type 1 receptors sample alternate configurations: MD simulations using coarse grain Martini 3 versus AlphaFold2 Multimer. Structure 2023; 31:735-745.e2. [PMID: 37075749 PMCID: PMC10833135 DOI: 10.1016/j.str.2023.03.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 02/07/2023] [Accepted: 03/23/2023] [Indexed: 04/21/2023]
Abstract
Structures and dynamics of transmembrane (TM) receptor regions are key to understanding their signaling mechanism across membranes. Here we examine configurations of TM region dimers, assembled using the recent Martini 3 force field for coarse-grain (CG) molecular dynamics simulations. At first glance, our results show only a reasonable agreement with ab initio predictions using PREDDIMER and AlphaFold2 Multimer and with nuclear magnetic resonance (NMR)-derived structures. 5 of 11 CG TM structures are similar to the NMR structures (within <3.5 Å root-mean-square deviation [RMSD]) compared with 10 and 9 using PREDDIMER and AlphaFold2, respectively (with 8 structures of the later within 1.5 Å). Surprisingly, AlphaFold2 predictions are closer to NMR structures when the 2001 instead of 2020 database is used for training. The CG simulations reveal that alternative configurations of TM dimers readily interconvert with a predominant population. The implications for transmembrane signaling are discussed, including for the development of peptide-based pharmaceuticals.
Collapse
Affiliation(s)
- Amita R Sahoo
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS & University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Zhiyuan Meng
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA.
| |
Collapse
|
8
|
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
|
9
|
Claveras Cabezudo A, Athanasiou C, Tsengenes A, Wade RC. Scaling Protein-Water Interactions in the Martini 3 Coarse-Grained Force Field to Simulate Transmembrane Helix Dimers in Different Lipid Environments. J Chem Theory Comput 2023; 19:2109-2119. [PMID: 36821400 DOI: 10.1021/acs.jctc.2c00950] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Abstract
Martini 3, the latest version of the widely used Martini force field for coarse-grained molecular dynamics simulations, is a promising tool to investigate proteins in phospholipid bilayers. However, simulating other lipid environments, such as detergent micelles, presents challenges due to the absence of validated parameters for their constituent molecules. Here, we propose parameters for the micelle-forming surfactant, dodecylphosphocholine (DPC). These result in micelle assembly with aggregation numbers in agreement with the experimental values. However, we identified a lack of hydrophobic interactions between transmembrane helix protein dimers and the tails of DPC molecules, preventing insertion and stabilization of the protein in the micelles. This problem was also observed for protein insertion by self-assembling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or dipalmitoylphosphatidylcholine (DPPC) bilayers. We propose the reduction of the nonbonded interactions between protein and water beads by 10% as a simple and effective solution to this problem that enables protein encapsulation in phospholipid micelles and bilayers without altering protein dimerization or the bilayer structure.
Collapse
Affiliation(s)
- Ainara Claveras Cabezudo
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany
| | - Christina Athanasiou
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Heidelberg Biosciences International Graduate School, Heidelberg University, Im Neuenheimer Feld 501, 69120 Heidelberg, Germany
| | - Alexandros Tsengenes
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Heidelberg Biosciences International Graduate School, Heidelberg University, Im Neuenheimer Feld 501, 69120 Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
10
|
Jin J, Pak AJ, Durumeric AEP, Loose TD, Voth GA. Bottom-up Coarse-Graining: Principles and Perspectives. J Chem Theory Comput 2022; 18:5759-5791. [PMID: 36070494 PMCID: PMC9558379 DOI: 10.1021/acs.jctc.2c00643] [Citation(s) in RCA: 72] [Impact Index Per Article: 36.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Indexed: 01/14/2023]
Abstract
Large-scale computational molecular models provide scientists a means to investigate the effect of microscopic details on emergent mesoscopic behavior. Elucidating the relationship between variations on the molecular scale and macroscopic observable properties facilitates an understanding of the molecular interactions driving the properties of real world materials and complex systems (e.g., those found in biology, chemistry, and materials science). As a result, discovering an explicit, systematic connection between microscopic nature and emergent mesoscopic behavior is a fundamental goal for this type of investigation. The molecular forces critical to driving the behavior of complex heterogeneous systems are often unclear. More problematically, simulations of representative model systems are often prohibitively expensive from both spatial and temporal perspectives, impeding straightforward investigations over possible hypotheses characterizing molecular behavior. While the reduction in resolution of a study, such as moving from an atomistic simulation to that of the resolution of large coarse-grained (CG) groups of atoms, can partially ameliorate the cost of individual simulations, the relationship between the proposed microscopic details and this intermediate resolution is nontrivial and presents new obstacles to study. Small portions of these complex systems can be realistically simulated. Alone, these smaller simulations likely do not provide insight into collectively emergent behavior. However, by proposing that the driving forces in both smaller and larger systems (containing many related copies of the smaller system) have an explicit connection, systematic bottom-up CG techniques can be used to transfer CG hypotheses discovered using a smaller scale system to a larger system of primary interest. The proposed connection between different CG systems is prescribed by (i) the CG representation (mapping) and (ii) the functional form and parameters used to represent the CG energetics, which approximate potentials of mean force (PMFs). As a result, the design of CG methods that facilitate a variety of physically relevant representations, approximations, and force fields is critical to moving the frontier of systematic CG forward. Crucially, the proposed connection between the system used for parametrization and the system of interest is orthogonal to the optimization used to approximate the potential of mean force present in all systematic CG methods. The empirical efficacy of machine learning techniques on a variety of tasks provides strong motivation to consider these approaches for approximating the PMF and analyzing these approximations.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Alexander J. Pak
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Aleksander E. P. Durumeric
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy D. Loose
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
11
|
Miszta P, Pasznik P, Niewieczerzał S, Młynarczyk K, Filipek S. COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments. J Chem Theory Comput 2022; 18:5145-5156. [PMID: 35998323 PMCID: PMC9476660 DOI: 10.1021/acs.jctc.2c00140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The presented methodology is based on coarse-grained
representation
of biomolecules in implicit environments and is designed for the molecular
dynamics simulations of membrane proteins and their complexes. The
membrane proteins are not only found in the cell membrane but also
in all membranous compartments of the cell: Golgi apparatus, mitochondria,
endosomes and lysosomes, and they usually form large complexes. To
investigate such systems the methodology is proposed based on two
independent approaches combining the coarse-grained MARTINI model
for proteins and the effective energy function to mimic the water/membrane
environments. The latter is based on the implicit environment developed
for all-atom simulations in the IMM1 method. The force field solvation
parameters for COGRIMEN were initially calculated from IMM1 all-atom
parameters and then optimized using Genetic Algorithms. The new methodology
was tested on membrane proteins, their complexes and oligomers. COGRIMEN
method is implemented as a patch for NAMD program and can be useful
for fast and brief studies of large membrane protein complexes.
Collapse
Affiliation(s)
- Przemysław Miszta
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw 02-093, Poland
| | - Paweł Pasznik
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw 02-093, Poland
| | - Szymon Niewieczerzał
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw 02-093, Poland
| | - Krzysztof Młynarczyk
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw 02-093, Poland
| | - Sławomir Filipek
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw 02-093, Poland
| |
Collapse
|
12
|
Marrink SJ, Monticelli L, Melo MN, Alessandri R, Tieleman DP, Souza PCT. Two decades of Martini: Better beads, broader scope. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1620] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials University of Groningen Groningen The Netherlands
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| | - Manuel N. Melo
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa Oeiras Portugal
| | - Riccardo Alessandri
- Pritzker School of Molecular Engineering University of Chicago Chicago Illinois USA
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences University of Calgary Alberta Canada
| | - Paulo C. T. Souza
- Molecular Microbiology and Structural Biochemistry (MMSB ‐ UMR 5086) CNRS & University of Lyon Lyon France
| |
Collapse
|
13
|
Majumder A, Kwon S, Straub JE. On Computing Equilibrium Binding Constants for Protein-Protein Association in Membranes. J Chem Theory Comput 2022; 18:3961-3971. [PMID: 35580264 DOI: 10.1021/acs.jctc.2c00106] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein association in lipid membranes is fundamental to membrane protein function and of great biomedical relevance. All-atom and coarse-grained models have been extensively used to understand the protein-protein interactions in the membrane and to compute equilibrium association constants. However, slow translational and rotational diffusion of protein in membrane presents challenges to the effective sampling of conformations defining the ensembles of free and bound states contributing to the association equilibrium and the free energy of dimerization. We revisit the homodimerization equilibrium of the TM region of glycophorin A. Conformational sampling is performed using umbrella sampling along previously proposed one-dimensional collective variables and compared with sampling over a two-dimensional collective variable space using the MARTINI v2.2 force field. We demonstrate that the one-dimensional collective variables suffer from restricted sampling of the native homodimer conformations leading to a biased free energy landscape. Conversely, simulations along the two-dimensional collective variable effectively characterize the thermodynamically relevant native and non-native interactions contributing to the association equilibrium. These results demonstrate the challenges associated with accurately characterizing binding equilibria when multiple poses contribute to the bound state ensemble.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Seulki Kwon
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
14
|
Alavizargar A, Elting A, Wedlich-Söldner R, Heuer A. Lipid-Mediated Association of the Slg1 Transmembrane Domains in Yeast Plasma Membranes. J Phys Chem B 2022; 126:3240-3256. [PMID: 35446028 DOI: 10.1021/acs.jpcb.2c00192] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Clustering of transmembrane proteins underlies a multitude of fundamental biological processes at the plasma membrane (PM) such as receptor activation, lateral domain formation, and mechanotransduction. The self-association of the respective transmembrane domains (TMDs) has also been suggested to be responsible for the micron-scaled patterns seen for integral membrane proteins in the budding yeast PM. However, the underlying interplay between the local lipid composition and the TMD identity is still not mechanistically understood. In this work, we combined coarse-grained molecular dynamics simulations of simplified bilayer systems with high-resolution live-cell microscopy to analyze the distribution of a representative helical yeast TMD from the PM sensor Slg1 within different lipid environments. In our simulations, we specifically evaluated the effects of acyl chain saturation and anionic lipid head groups on the association of two TMDs. We found that weak lipid-protein interactions significantly affect the configuration of TMD dimers and the free energy of association. Increased amounts of unsaturated phospholipids (PLs) strongly reduced the helix-helix interaction, while the presence of anionic phosphatidylserine (PS) hardly affected the dimer formation. We could experimentally confirm this surprising lack of effect of PS using the network factor, a mesoscopic measure of PM pattern formation in yeast cells. Simulations also showed that the formation of TMD dimers in turn increased the order parameter of the surrounding lipids and induced long-range perturbations in lipid organization. In summary, our results shed new light on the mechanisms of lipid-mediated dimerization of TMDs in complex lipid mixtures.
Collapse
Affiliation(s)
- Azadeh Alavizargar
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| | - Annegret Elting
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Roland Wedlich-Söldner
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Andreas Heuer
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| |
Collapse
|
15
|
Thomasen FE, Pesce F, Roesgaard MA, Tesei G, Lindorff-Larsen K. Improving Martini 3 for Disordered and Multidomain Proteins. J Chem Theory Comput 2022; 18:2033-2041. [PMID: 35377637 DOI: 10.1021/acs.jctc.1c01042] [Citation(s) in RCA: 35] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.
Collapse
Affiliation(s)
- F Emil Thomasen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Francesco Pesce
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Mette Ahrensback Roesgaard
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Giulio Tesei
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
16
|
Borges-Araújo L, Souza PCT, Fernandes F, Melo MN. Improved Parameterization of Phosphatidylinositide Lipid Headgroups for the Martini 3 Coarse-Grain Force Field. J Chem Theory Comput 2021; 18:357-373. [PMID: 34962393 DOI: 10.1021/acs.jctc.1c00615] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Phosphoinositides are a family of membrane phospholipids that play crucial roles in membrane regulatory events. As such, these lipids are often a key part of molecular dynamics simulation studies of biological membranes, in particular of those employing coarse-grain models because of the potential long times and sizes of the involved membrane processes. Version 3 of the widely used Martini coarse-grain force field has been recently published, greatly refining many aspects of biomolecular interactions. In order to properly use it for lipid membrane simulations with phosphoinositides, we put forth the Martini 3-specific parameterization of inositol, phosphatidylinositol, and seven physiologically relevant phosphorylated derivatives of phosphatidylinositol. Compared to parameterizations for earlier Martini versions, focus was put on a more accurate reproduction of the behavior seen in both atomistic simulations and experimental studies, including the signaling-relevant phosphoinositide interaction with divalent cations. The models that we develop improve upon the conformational dynamics of phosphoinositides in the Martini force field and provide stable topologies at typical Martini time steps. They are able to reproduce experimentally known protein-binding poses as well as phosphoinositide aggregation tendencies. The latter was tested both in the presence and absence of calcium and included correct behavior of PI(4,5)P2 calcium-induced clusters, which can be of relevance for regulation.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- iBB-Institute for Bioengineering and Biosciences, Instituto Superior Técnico, Universidade de Lisboa, Lisbon 1049-001, Portugal.,Associate Laboratory i4HB─Institute for Health and Bioeconomy, at Instituto Superior Técnico, Universidade de Lisboa, Lisbon 1049-001, Portugal.,Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, Oeiras 2780-157, Portugal
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry, UMR 5086 CNRS & University of Lyon, 7 Passage du Vercors, Lyon F-69367, France
| | - Fábio Fernandes
- iBB-Institute for Bioengineering and Biosciences, Instituto Superior Técnico, Universidade de Lisboa, Lisbon 1049-001, Portugal.,Associate Laboratory i4HB─Institute for Health and Bioeconomy, at Instituto Superior Técnico, Universidade de Lisboa, Lisbon 1049-001, Portugal.,Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon 1049-001, Portugal
| | - Manuel N Melo
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, Oeiras 2780-157, Portugal
| |
Collapse
|
17
|
Sawade K, Peter C. Multiscale simulations of protein and membrane systems. Curr Opin Struct Biol 2021; 72:203-208. [PMID: 34953308 DOI: 10.1016/j.sbi.2021.11.010] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 11/01/2021] [Accepted: 11/10/2021] [Indexed: 02/07/2023]
Abstract
Classical multiscale simulations are perfectly suited to investigate biological soft matter systems. Owing to the bridging between microscopically realistic and lower-resolution models or the integration of a hierarchy of subsystems, one gets access to biologically relevant system sizes and timescales. In recent years, increasingly complex systems and processes have come into focus such as multidomain proteins, phase separation processes in biopolymer solutions, multicomponent biomembranes, or multiprotein complexes up to entire viruses. The review shows factors that have contributed to this progress - from improved models to machine-learning-based analysis and scale-bridging methods.
Collapse
Affiliation(s)
- Kevin Sawade
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78 457, Konstanz, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78 457, Konstanz, Germany.
| |
Collapse
|
18
|
Ugarte La Torre D, Takada S. Modeling lipid-protein interactions for coarse-grained lipid and Cα protein models. J Chem Phys 2021; 155:155101. [PMID: 34686048 DOI: 10.1063/5.0057278] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Biological membranes that play major roles in diverse functions are composed of numerous lipids and proteins, making them an important target for coarse-grained (CG) molecular dynamics (MD) simulations. Recently, we have developed the CG implicit solvent lipid force field (iSoLF) that has a resolution compatible with the widely used Cα protein representation [D. Ugarte La Torre and S. Takada, J. Chem. Phys. 153, 205101 (2020)]. In this study, we extended it and developed a lipid-protein interaction model that allows the combination of the iSoLF and the Cα protein force field, AICG2+. The hydrophobic-hydrophilic interaction is modeled as a modified Lennard-Jones potential in which parameters were tuned partly to reproduce the experimental transfer free energy and partly based on the free energy profile normal to the membrane surface from previous all-atom MD simulations. Then, the obtained lipid-protein interaction is tested for the configuration and placement of transmembrane proteins, water-soluble proteins, and peripheral proteins, showing good agreement with prior knowledge. The interaction is generally applicable and is implemented in the publicly available software, CafeMol.
Collapse
Affiliation(s)
- Diego Ugarte La Torre
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| |
Collapse
|
19
|
Sahoo AR, Buck M. Structural and Functional Insights into the Transmembrane Domain Association of Eph Receptors. Int J Mol Sci 2021; 22:ijms22168593. [PMID: 34445298 PMCID: PMC8395321 DOI: 10.3390/ijms22168593] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2021] [Revised: 08/02/2021] [Accepted: 08/06/2021] [Indexed: 12/04/2022] Open
Abstract
Eph receptors are the largest family of receptor tyrosine kinases and by interactions with ephrin ligands mediate a myriad of processes from embryonic development to adult tissue homeostasis. The interaction of Eph receptors, especially at their transmembrane (TM) domains is key to understanding their mechanism of signal transduction across cellular membranes. We review the structural and functional aspects of EphA1/A2 association and the techniques used to investigate their TM domains: NMR, molecular modelling/dynamics simulations and fluorescence. We also introduce transmembrane peptides, which can be used to alter Eph receptor signaling and we provide a perspective for future studies.
Collapse
Affiliation(s)
- Amita R. Sahoo
- Department of Physiology and Biophysics, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA;
| | - Matthias Buck
- Department of Physiology and Biophysics, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA;
- Department of Neurosciences, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA
- Department of Pharmacology, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA
- Case Comprehensive Cancer Center, School of Medicine, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA
- Correspondence:
| |
Collapse
|