scholarly journals Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

2021 ◽  
Author(s):  
Andreas Haahr Larsen ◽  
Lilya Tata ◽  
Laura John ◽  
Mark S.P. Sansom
Keyword(s):  
Molecular Dynamics ◽  
Membrane Binding ◽  
Coiled Coil ◽  
Binding Mode ◽  
Coarse Grained ◽  
Total Binding Energy ◽  
Endosomal Trafficking ◽  
Fyve Domain ◽  
C Terminus ◽  
Endosomal Membranes

AbstractEarly Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The total binding energy is 70 kJ/mol, of which 50-60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem, allowing binding via a mechanism resembling that of a suction cup connected to a locally rigid stem via a flexible hinge.

scholarly journals Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

2021 ◽  
Vol 17 (9) ◽  
pp. e1008807
Author(s):  
Andreas Haahr Larsen ◽  
Lilya Tata ◽  
Laura H. John ◽  
Mark S. P. Sansom
Keyword(s):  
Molecular Dynamics ◽  
Membrane Binding ◽  
Coiled Coil ◽  
Binding Mode ◽  
The Other ◽  
Coarse Grained ◽  
Total Binding Energy ◽  
Fyve Domain ◽  
C Terminus ◽  
Endosomal Membranes

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50–60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 μm2 in the search for a second endosome with which to interact.

scholarly journals Chiral twisting in cytoskeletal polymers regulates filament size and orientation

2018 ◽  
Author(s):  
Handuo Shi ◽  
David Quint ◽  
Ajay Gopinathan ◽  
Kerwyn Casey Huang
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Caulobacter Crescentus ◽  
Membrane Binding ◽  
Membrane Curvature ◽  
Cytoskeletal Proteins ◽  
Coarse Grained ◽  
Biophysical Model ◽  
Dynamics Simulations

AbstractWhile cytoskeletal proteins in the actin family are structurally similar, as filaments they act as critical components of diverse cellular processes across all kingdoms of life. In many rod-shaped bacteria, the actin homolog MreB directs cell-wall insertion and maintains cell shape, but it remains unclear how structural changes to MreB affect its physiological function. To bridge this gap, we performed molecular dynamics simulations forCaulobacter crescentusMreB and then utilized a coarse-grained biophysical model to successfully predict MreB filament propertiesin vivo.We discovered that MreB double protofilaments exhibit left-handed twisting that is dependent on the bound nucleotide and membrane binding; the degree of twisting determines the limit length and orientation of MreB filamentsin vivo.Membrane binding of MreB also induces a stable membrane curvature that is physiologically relevant. Together, our data empower the prediction of cytoskeletal filament size from molecular dynamics simulations, providing a paradigm for connecting protein filament structure and mechanics to cellular functions.

scholarly journals Binding Mode of SARS-CoV2 Fusion Peptide to Human Cellular Membrane

2020 ◽  
Author(s):  
D. Gorgun ◽  
M. Lihan ◽  
K. Kapoor ◽  
E. Tajkhorshid
Keyword(s):  
Molecular Dynamics ◽  
Host Cell ◽  
Membrane Binding ◽  
Binding Mode ◽  
Fusion Peptide ◽  
Atomic Level ◽  
Human Cells ◽  
Spike Protein ◽  
Cellular Membranes ◽  
Membrane Bound

ABSTRACTInfection of human cells by the SARS-CoV2 relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics (MD) simulations taking advantage of the Highly Mobile Membrane Mimetic (HMMM) model, to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level, and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas. In 60% of the simulations, the FP reaches a stable, membrane-bound configuration where the peptide deeply penetrated into the membrane. Clustering of the results reveals two major membrane binding modes, the helix-binding mode and the loop-binding mode. Taken into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, we propose that the helix-binding mode represents more closely the biologically relevant form. In the helix-binding mode, the helix is stabilized in an oblique angle with respect to the membrane with its N-terminus tilted towards the membrane core. Analysis of the FP-lipid interactions shows the involvement of specific residues of the helix in membrane binding previously described as the fusion active core residues. Taken together, the results shed light on a key step involved in SARS-CoV2 infection with potential implications in designing novel inhibitors.SIGNIFICANCEA key step in cellular infection by the SARS-CoV2 virus is its attachment to and penetration into the plasma membrane of human cells. These processes hinge upon the membrane interaction of the viral fusion peptide, a segment exposed by the spike protein upon its conformational changes after encountering the host cell. In this study, using molecular dynamics simulations, we describe how the fusion peptide from the SARS-CoV2 virus binds human cellular membranes and characterize, at an atomic level, lipid-protein interactions important for the stability of the bound state.

Investigation of the Microscopic Viscoelastic Property for Cross-linked Polymer Network by Molecular Dynamics Simulation

2011 ◽  
Vol 39 (1) ◽  
pp. 44-58 ◽  
Author(s):  
Y. Masumoto ◽  
Y. Iida
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulation ◽  
Analytical Method ◽  
Viscoelastic Properties ◽  
Polymer Network ◽  
Dynamics Simulation ◽  
Coarse Grained ◽  
The Cross ◽  
Microscopic Dynamic ◽  
Coarse Grained Molecular Dynamics

Abstract The purpose of this work is to develop a new analytical method for simulating the microscopic mechanical property of the cross-linked polymer system using the coarse-grained molecular dynamics simulation. This new analytical method will be utilized for the molecular designing of the tire rubber compound to improve the tire performances such as rolling resistance and wet traction. First, we evaluate the microscopic dynamic viscoelastic properties of the cross-linked polymer using coarse-grained molecular dynamics simulation. This simulation has been conducted by the coarse-grained molecular dynamics program in the OCTA) (http://octa.jp/). To simplify the problem, we employ the bead-spring model, in which a sequence of beads connected by springs denotes a polymer chain. The linear polymer chains that are cross-linked by the cross-linking agents express the three-dimensional cross-linked polymer network. In order to obtain the microscopic dynamic viscoelastic properties, oscillatory deformation is applied to the simulation cell. By applying the time-temperature reduction law to this simulation result, we can evaluate the dynamic viscoelastic properties in the wide deformational frequency range including the rubbery state. Then, the stress is separated into the nonbonding stress and the bonding stress. We confirm that the contribution of the nonbonding stress is larger at lower temperatures. On the other hand, the contribution of the bonding stress is larger at higher temperatures. Finally, analyzing a change of microscopic structure in dynamic oscillatory deformation, we determine that the temperature/frequency dependence of bond stress response to a dynamic oscillatory deformation depends on the temperature dependence of the average bond length in the equilibrium structure and the temperature/frequency dependence of bond orientation. We show that our simulation is a useful tool for studying the microscopic properties of a cross-linked polymer.

Parameterization of Divalent Cations for Coarse-Grained Simulations

2020 ◽  
Author(s):  
Florencia Klein ◽  
Daniela Cáceres-Rojas ◽  
Monica Carrasco ◽  
Juan Carlos Tapia ◽  
Julio Caballero ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Metal Ions ◽  
Molecular Dynamics Simulations ◽  
Divalent Cations ◽  
Computational Cost ◽  
Data Bank ◽  
Coarse Grained ◽  
Interaction Parameters ◽  
Dynamics Simulations ◽  
Dynamical Description

<p>Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges to achieve an accurate structural and dynamical description of many biological assemblies. This is particularly the case for coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or simply not available for the vast majority of CG- force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations. Using this simple approach, we provide a set of interaction parameters for Calcium, Magnesium, and Zinc ions, which cover more than 80% of the metal-bound structures reported on the PDB. Simulations performed using the SIRAH force field on several proteins and DNA systems show that using the present approach it is possible to obtain non-bonded interaction parameters that obviate the use of topological constraints. </p>

scholarly journals Computational Investigations on the Binding Mode of Ligands for the Cannabinoid-Activated G Protein-Coupled Receptor GPR18

Biomolecules ◽  
2020 ◽  
Vol 10 (5) ◽  
pp. 686 ◽  
Author(s):  
Alexander Neumann ◽  
Viktor Engel ◽  
Andhika B. Mahardhika ◽  
Clara T. Schoeder ◽  
Vigneshwaran Namasivayam ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulation ◽  
G Protein ◽  
Phenyl Ring ◽  
Binding Mode ◽  
Dynamics Simulation ◽  
Simulation Studies ◽  
Binding Modes ◽  
G Protein Coupled Receptor ◽  
G Protein Coupled

GPR18 is an orphan G protein-coupled receptor (GPCR) expressed in cells of the immune system. It is activated by the cannabinoid receptor (CB) agonist ∆9-tetrahydrocannabinol (THC). Several further lipids have been proposed to act as GPR18 agonists, but these results still require unambiguous confirmation. In the present study, we constructed a homology model of the human GPR18 based on an ensemble of three GPCR crystal structures to investigate the binding modes of the agonist THC and the recently reported antagonists which feature an imidazothiazinone core to which a (substituted) phenyl ring is connected via a lipophilic linker. Docking and molecular dynamics simulation studies were performed. As a result, a hydrophobic binding pocket is predicted to accommodate the imidazothiazinone core, while the terminal phenyl ring projects towards an aromatic pocket. Hydrophobic interaction of Cys251 with substituents on the phenyl ring could explain the high potency of the most potent derivatives. Molecular dynamics simulation studies suggest that the binding of imidazothiazinone antagonists stabilizes transmembrane regions TM1, TM6 and TM7 of the receptor through a salt bridge between Asp118 and Lys133. The agonist THC is presumed to bind differently to GPR18 than to the distantly related CB receptors. This study provides insights into the binding mode of GPR18 agonists and antagonists which will facilitate future drug design for this promising potential drug target.

scholarly journals Study of Endogen Substrates, Drug Substrates and Inhibitors Binding Conformations on MRP4 and Its Variants by Molecular Docking and Molecular Dynamics

Molecules ◽  
2021 ◽  
Vol 26 (4) ◽  
pp. 1051
Author(s):  
Edgardo Becerra ◽  
Giovanny Aguilera-Durán ◽  
Laura Berumen ◽  
Antonio Romo-Mancillas ◽  
Guadalupe García-Alcocer
Keyword(s):  
Molecular Dynamics ◽  
Molecular Docking ◽  
Binding Energy ◽  
Binding Site ◽  
Binding Sites ◽  
Adenosine Monophosphate ◽  
Competitive Inhibitor ◽  
Umbrella Sampling ◽  
Coarse Grained ◽  
Nucleotide Polymorphisms

Multidrug resistance protein-4 (MRP4) belongs to the ABC transporter superfamily and promotes the transport of xenobiotics including drugs. A non-synonymous single nucleotide polymorphisms (nsSNPs) in the ABCC4 gene can promote changes in the structure and function of MRP4. In this work, the interaction of certain endogen substrates, drug substrates, and inhibitors with wild type-MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using protein threading modeling, molecular docking, all-atom, coarse grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively). The results showed that the three MRP4 structures had significantly different conformations at given sites, leading to differences in the docking scores (DS) and binding sites of three different groups of molecules. Folic acid (FA) had the highest variation in DS on G187W concerning WT-MRP4. WT-MRP4, G187W, Y556C, and FA had different conformations through 25 ns AA-MD. Umbrella sampling simulations indicated that the Y556C-FA complex was the most stable one with or without ATP. In Y556C, the cyclic adenosine monophosphate (cAMP) and ceefourin-1 binding sites are located out of the entrance of the inner cavity, which suggests that both cAMP and ceefourin-1 may not be transported. The binding site for cAMP and ceefourin-1 is quite similar and the affinity (binding energy) of ceefourin-1 to WT-MRP4, G187W, and Y556C is greater than the affinity of cAMP, which may suggest that ceefourin-1 works as a competitive inhibitor. In conclusion, the nsSNPs G187W and Y556C lead to changes in protein conformation, which modifies the ligand binding site, DS, and binding energy.

Sign in / Sign up

Export Citation Format

Share Document