Comparing an All-Atom and a Coarse-Grained Description of Lipid Bilayers in Terms of Enthalpies and Entropies: From MD Simulations to 2D Lattice Models

2019 ◽  
Vol 15 (11) ◽  
pp. 6393-6402
Author(s):  
Davit Hakobyan ◽  
Andreas Heuer
Keyword(s):  
Lipid Bilayers ◽  
Lattice Models ◽  
Md Simulations ◽  
Coarse Grained

scholarly journals In Silico Study of the Mechanism of Binding of the N-Terminal Region of α Synuclein to Synaptic-Like Membranes

Life ◽  
2020 ◽  
Vol 10 (6) ◽  
pp. 98 ◽  
Author(s):  
Carlos Navarro-Paya ◽  
Maximo Sanz-Hernandez ◽  
Alfonso De Simone
Keyword(s):  
Lipid Bilayers ◽  
Conformational Transition ◽  
Md Simulations ◽  
Bound State ◽  
Terminal Region ◽  
Coarse Grained ◽  
Biological Behavior ◽  
In Silico Study ◽  
Membrane Bound ◽  
Presynaptic Protein

The membrane binding by α-synuclein (αS), a presynaptic protein whose aggregation is strongly linked with Parkinson’s disease, influences its biological behavior under functional and pathological conditions. This interaction requires a conformational transition from a disordered-unbound to a partially helical membrane-bound state of the protein. In the present study, we used enhanced coarse-grained MD simulations to characterize the sequence and conformational determinants of the binding to synaptic-like vesicles by the N-terminal region of αS. This region is the membrane anchor and is of crucial importance for the properties of the physiological monomeric state of αS as well as for its aberrant aggregates. These results identify the key factors that play a role in the binding of αS with synaptic lipid bilayers in both membrane-tethered and membrane-locked conformational states.

Interactions of peripheral proteins with model membranes as viewed by molecular dynamics simulations

2014 ◽  
Vol 42 (5) ◽  
pp. 1418-1424 ◽  
Author(s):  
Antreas C. Kalli ◽  
Mark S. P. Sansom
Keyword(s):  
Lipid Bilayers ◽  
Molecular Mechanisms ◽  
Phospholipid Composition ◽  
Md Simulations ◽  
Multiscale Simulation ◽  
Lipid Binding ◽  
Coarse Grained ◽  
Peripheral Proteins ◽  
Dynamics Simulations ◽  
Phosphatidylinositol Phosphates

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein–lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

scholarly journals Binding of Ca2+-Independent C2 Domains to Lipid Membranes: a Multi-Scale Molecular Dynamics Study

2020 ◽  
Author(s):  
Andreas Haahr Larsen ◽  
Mark S.P. Sansom
Keyword(s):  
Molecular Dynamics ◽  
Lipid Bilayers ◽  
Md Simulations ◽  
Multiscale Simulation ◽  
Coarse Grained ◽  
Type I ◽  
Type Ii ◽  
Binding Modes ◽  
C2 Domains ◽  
Anionic Lipids

AbstractC2 domains facilitate protein-lipid interaction in cellular recognition and signalling processes. They possess a β-sandwich structure, with either type I or type II topology. C2 domains can interact with anionic lipid bilayers in either a Ca2+-dependent or a Ca2+-independent manner. The mechanism of recognition of anionic lipids by Ca2+-independent C2 domains is incompletely understood. We have used molecular dynamics (MD) simulations to explore the membrane interactions of six Ca2+– independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse grained MD simulations these C2 domains bound to lipid bilayers, forming transient interactions with zwitterionic (phosphatidylcholine, PC) bilayers compared to long lived interactions with anionic bilayers also containing either phosphatidylserine (PS) or PS and phosphatidylinositol bisphosphate (PIP2). Type I C2 domains bound non-canonically via the front, back or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops (as is typically the case for Ca2+-dependent C2 domains). C2 domains interacted strongly (up to 120 kJ/mol) with membranes containing PIP2 causing the bound anionic lipids to clustered around the protein. The C2 domains bound less strongly to anionic membranes without PIP2 (<50 kJ/mol), and most weakly to neutral membranes (<33 kJ/mol). Productive binding modes were identified and further analysed in atomistic simulations. For PTEN and SHIP2, CG simulations were also performed of the intact enzymes (i.e. phosphatase domain plus C2 domain) with PIP2-contating bilayers and the roles of the two domains in membrane localization were compared. From a methodological perspective, these studies establish a multiscale simulation protocol for studying membrane binding/recognition proteins, capable of revealing binding modes alongside details of lipid binding affinity and specificity.

scholarly journals Formation of aggregates, icosahedral structures and percolation clusters of fullerenes in lipids bilayers: The key role of lipid saturation

2020 ◽  
Author(s):  
Nililla Nisoh ◽  
Viwan Jarerattanachat ◽  
Mikko Karttunen ◽  
Jirasak Wong-ekkabut
Keyword(s):  
Lipid Bilayers ◽  
Carbon Nanoparticles ◽  
Md Simulations ◽  
Coarse Grained ◽  
Percolation Clusters ◽  
High Fullerene ◽  
Key Factor ◽  
Chain Lengths ◽  
Lateral Area

AbstractCarbon nanoparticles (CNPs) are attractive materials for a great number of applications but there are serious concerns regarding their influence on health and environment. Here, our focus is on the behavior of fullerenes in lipid bilayers with varying lipid saturations, chain lengths and fullerene concentrations using coarse-grained molecular dynamics (CG-MD) simulations. Our findings show that the lipid saturation level is a key factor in determining how fullerenes behave and where the fullerenes are located inside a lipid bilayer. In saturated and monounsaturated bilayers fullerenes aggregated and formed clusters with some of them showing icosahedral structures. In polyunsaturated lipid bilayers, no such structures were observed: In polyunsaturated lipid bilayers at high fullerene concentrations, connected percolation-like networks of fullerenes spanning the whole lateral area emerged at the bilayer center. In other systems only separate isolated aggregates were observed. The effects of fullerenes on lipid bilayers depend strongly on fullerene aggregation. When fullerenes aggregate, their interactions with the lipid tails change.

scholarly journals Evaluating coarse-grained MARTINI force-fields for capturing the ripple phase of lipid membranes

2020 ◽  
Author(s):  
Pradyumn Sharma ◽  
Rajat Desikan ◽  
K. Ganapathy Ayappa
Keyword(s):  
Lipid Bilayers ◽  
Lipid Membranes ◽  
Liquid Crystalline ◽  
Intermediate Phase ◽  
Finite Size ◽  
Fluid Phase ◽  
Md Simulations ◽  
Coarse Grained ◽  
Lamellar Phases ◽  
Ripple Phase

AbstractPhospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (Pβ′) membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained resolution. Here, we systematically assess the ability of five coarse-grained MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid force-field (FF) variants, parametrized to reproduce the DPPC gel and fluid phases, for their ability to capture the Pβ′ phase. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (> 2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated Pβ′ and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the Pβ′ phase, similar to previous studies, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically-trapped structures caused by finite size effects rather than being representative of true Pβ′ phases. We showed that even a MARTINI FF variant that could capture the tilted Lβ′ gel phase, a prerequisite for stabilizing the Pβ′ phase, could not capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs may require specific re-parametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.

scholarly journals Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

Molecules ◽  
2021 ◽  
Vol 26 (11) ◽  
pp. 3293
Author(s):  
Mateusz Zalewski ◽  
Sebastian Kmiecik ◽  
Michał Koliński
Keyword(s):  
Molecular Dynamics ◽  
Geometry Optimization ◽  
Computational Prediction ◽  
Md Simulations ◽  
Coarse Grained ◽  
Ligand Interaction ◽  
Vast Number ◽  
Docking Simulations ◽  
Peptide Models ◽  
Peptide Complexes

One of the major challenges in the computational prediction of protein–peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein–peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor–ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein–peptide complexes show that the proposed MD-based scoring approach can be used to identify protein–peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

scholarly journals Dynamical Behavior and Conformational Selection Mechanism of the Intrinsically Disordered Sic1 Kinase-Inhibitor Domain

Life ◽  
2020 ◽  
Vol 10 (7) ◽  
pp. 110 ◽  
Author(s):  
Davide Sala ◽  
Ugo Cosentino ◽  
Anna Ranaudo ◽  
Claudio Greco ◽  
Giorgio Moro
Keyword(s):  
Kinase Inhibitor ◽  
Dynamical Behavior ◽  
Md Simulations ◽  
Molecular Shape ◽  
Conformational Space ◽  
Coarse Grained ◽  
Conformational Ensemble ◽  
Selection Mechanism ◽  
Conformational Selection ◽  
Intrinsically Disordered

Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners.

scholarly journals Biophysical characterization of the SARS-CoV-2 E protein

2021 ◽  
Author(s):  
Aujan Mehregan ◽  
Sergio Perez-Conesa ◽  
Yuxuan Zhuang ◽  
Ahmad Elbahnsi ◽  
Diletta Pasini ◽  
...  
Keyword(s):  
Recombinant Expression ◽  
Structural Data ◽  
Md Simulations ◽  
Full Length ◽  
Coarse Grained ◽  
Biophysical Characterization ◽  
Viral Budding ◽  
Nmr Structures ◽  
E Protein

SARS-CoV-2 is the virus responsible for the COVID-19 pandemic which continues to wreak havoc across the world, over a year and a half after its effects were first reported in the general media. Current fundamental research efforts largely focus on the SARS-CoV-2 Spike protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length E protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening.

scholarly journals Molecular simulations of lipid membrane partitioning and translocation by bacterial quorum sensing modulators

PLoS ONE ◽  
2021 ◽  
Vol 16 (2) ◽  
pp. e0246187
Author(s):  
Tianyi Jin ◽  
Samarthaben J. Patel ◽  
Reid C. Van Lehn
Keyword(s):  
Quorum Sensing ◽  
Lipid Bilayers ◽  
Lipid Membrane ◽  
Computational Cost ◽  
Md Simulations ◽  
Umbrella Sampling ◽  
Communication Process ◽  
Free Energies ◽  
Membrane Partitioning ◽  
Solvent Model

Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.

scholarly journals Residue-Level Contact Reveals Modular Domain Interactions of PICK1 Are Driven by Both Electrostatic and Hydrophobic Forces

2021 ◽  
Vol 7 ◽  
Author(s):  
Amy O. Stevens ◽  
Yi He
Keyword(s):  
Hydrophobic Interactions ◽  
Pdz Domain ◽  
Md Simulations ◽  
Intramolecular Interactions ◽  
Coarse Grained ◽  
Scaffolding Protein ◽  
Concave Surface ◽  
Stable Complex ◽  
Bar Domain ◽  
C Terminus

PICK1 is a multi-domain scaffolding protein that is uniquely comprised of both a PDZ domain and a BAR domain. While previous experiments have shown that the PDZ domain and the linker positively regulate the BAR domain and the C-terminus negatively regulates the BAR domain, the details of internal regulation mechanisms are unknown. Molecular dynamics (MD) simulations have been proven to be a useful tool in revealing the intramolecular interactions at atomic-level resolution. PICK1 performs its biological functions in a dimeric form which is extremely computationally demanding to simulate with an all-atom force field. Here, we use coarse-grained MD simulations to expose the key residues and driving forces in the internal regulations of PICK1. While the PDZ and BAR domains do not form a stable complex, our simulations show the PDZ domain preferentially interacting with the concave surface of the BAR domain over other BAR domain regions. Furthermore, our simulations show that the short helix in the linker region can form interactions with the PDZ domain. Our results reveal that the surface of the βB-βC loop, βC strand, and αA-βD loop of the PDZ domain can form a group of hydrophobic interactions surrounding the linker helix. These interactions are driven by hydrophobic forces. In contrast, our simulations reveal a very dynamic C-terminus that most often resides on the convex surface of the BAR domain rather than the previously suspected concave surface. These interactions are driven by a combination of electrostatic and hydrophobic interactions.

Sign in / Sign up

Export Citation Format

Share Document