scholarly journals Implicit micelle model for membrane proteins using super-ellipsoid approximation

2019 ◽  
Author(s):  
Takaharu Mori ◽  
Yuji Sugita
Keyword(s):  
Membrane Proteins ◽  
Computational Cost ◽  
Surface Effects ◽  
Md Simulations ◽  
Curved Surface ◽  
Hydrophobic Core ◽  
Surfactant Micelles ◽  
Membrane Mimetic ◽  
Explicit Solvent ◽  
Protein Gp41

AbstractSurfactant micelles are often utilized as membrane mimetics for structure determination and functional analysis of membrane proteins. Although curved-surface effects of the micelle can perturb their structure, it is difficult to assess such effects and membrane mimetic artifacts by experimental and theoretical methods. Here, we propose an implicit micelle model (IMIC) to be used in molecular dynamics (MD) simulations of membrane proteins. IMIC is an extension of the IMM1 implicit membrane model by introducing a super-ellipsoid approximation to represent the curved-surface effects. Most of the parameters for IMIC are obtained from all-atom explicit solvent MD simulations of twelve membrane proteins in various micelles. In simulations of the HIV envelop protein gp41, M13 major coat protein gp8, and amyloid precursor protein (APP) dimer, curved-surface and compact hydrophobic-core effects are exhibited. The MD simulations with IMIC provide accurate structure predictions of membrane proteins in various micelle environments quickly with smaller computational cost than that necessary for explicit solvent/micelle model.

scholarly journals Conformational Ensembles by NMR and MD Simulations in Model Heptapeptides with Select Tri-Peptide Motifs

2021 ◽  
Vol 22 (3) ◽  
pp. 1364
Author(s):  
V. V. Krishnan ◽  
Timothy Bentley ◽  
Alina Xiong ◽  
Kalyani Maitra
Keyword(s):  
Principal Component ◽  
Md Simulations ◽  
Conformational Space ◽  
Random Coil ◽  
Radius Of Gyration ◽  
Small Peptides ◽  
Conformational Ensembles ◽  
Solution State ◽  
Peptide Motifs ◽  
Explicit Solvent

Both nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations are routinely used in understanding the conformational space sampled by peptides in the solution state. To investigate the role of single-residue change in the ensemble of conformations sampled by a set of heptapeptides, AEVXEVG with X = L, F, A, or G, comprehensive NMR, and MD simulations were performed. The rationale for selecting the particular model peptides is based on the high variability in the occurrence of tri-peptide E*L between the transmembrane β-barrel (TMB) than in globular proteins. The ensemble of conformations sampled by E*L was compared between the three sets of ensembles derived from NMR spectroscopy, MD simulations with explicit solvent, and the random coil conformations. In addition to the estimation of global determinants such as the radius of gyration of a large sample of structures, the ensembles were analyzed using principal component analysis (PCA). In general, the results suggest that the -EVL- peptide indeed adopts a conformational preference that is distinctly different not only from a random distribution but also from other peptides studied here. The relatively straightforward approach presented herein could help understand the conformational preferences of small peptides in the solution state.

scholarly journals The role of BamA in the biogenesis of beta-barrel membrane proteins

2014 ◽  
Vol 70 (a1) ◽  
pp. C578-C578
Author(s):  
Nicholas Noinaj ◽  
Adam Kuszak ◽  
Curtis Balusek ◽  
JC Gumbart ◽  
Petra Lukacik ◽  
...  
Keyword(s):  
Membrane Proteins ◽  
Outer Membrane ◽  
Crystal Structures ◽  
Outer Membrane Proteins ◽  
Hydrophobic Surface ◽  
Md Simulations ◽  
Structural Features ◽  
Bam Complex ◽  
Beta Barrel

Beta-barrel membrane proteins are essential for nutrient import, signaling, motility, and survival. In Gram-negative bacteria, the beta-barrel assembly machinery (BAM) complex is responsible for the biogenesis of beta-barrel outer membrane proteins (OMPs), with homologous complexes found in mitochondria and chloroplasts. Despite their essential roles, exactly how these OMPs are formed remains unknown. The BAM complex consists of a central and essential component called BamA (an OMP itself) and four lipoproteins called BamB-E. While the structure of the lipoproteins have been reported, the structure of full length BamA has been elusive. Recently though, we described the structure of BamA from two species of bacteria: Neisseria gonorrhoeae and Haemophilus ducreyi. BamA consists of a large periplasmic domain attached to a 16-strand transmembrane beta-barrel domain. Together, our crystal structures and molecule dynamics (MD) simulations revealed several structural features which gave clues to the mechanism by which BamA catalyzes beta-barrel assembly. The first is that the interior cavity is accessible in one BamA structure and conformationally closed in the other. Second, an exterior rim of the beta-barrel has a distinctly narrowed hydrophobic surface, locally destabilizing the outer membrane. Third, the beta-barrel can undergo lateral opening, suggesting a route from the interior cavity in BamA into the outer membrane. And fourth, a surface exposed exit pore positioned above the lateral opening site which may play a role in the biogenesis of extracellular loops. In this presentation, the crystal structures and MD simulations of BamA will be presented along with our work looking at the role of these four structural features in the role of BamA within the BAM complex.

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 A hierarchical Bayesian framework for force field selection in molecular dynamics simulations

2016 ◽  
Vol 374 (2060) ◽  
pp. 20150032 ◽  
Author(s):  
S. Wu ◽  
P. Angelikopoulos ◽  
C. Papadimitriou ◽  
R. Moser ◽  
P. Koumoutsakos
Keyword(s):  
Experimental Data ◽  
Molecular Dynamics ◽  
Force Field ◽  
Computational Cost ◽  
Md Simulations ◽  
Bayesian Framework ◽  
Underlying Structure ◽  
Hierarchical Bayesian ◽  
Wide Range ◽  
Selection Of

We present a hierarchical Bayesian framework for the selection of force fields in molecular dynamics (MD) simulations. The framework associates the variability of the optimal parameters of the MD potentials under different environmental conditions with the corresponding variability in experimental data. The high computational cost associated with the hierarchical Bayesian framework is reduced by orders of magnitude through a parallelized Transitional Markov Chain Monte Carlo method combined with the Laplace Asymptotic Approximation. The suitability of the hierarchical approach is demonstrated by performing MD simulations with prescribed parameters to obtain data for transport coefficients under different conditions, which are then used to infer and evaluate the parameters of the MD model. We demonstrate the selection of MD models based on experimental data and verify that the hierarchical model can accurately quantify the uncertainty across experiments; improve the posterior probability density function estimation of the parameters, thus, improve predictions on future experiments; identify the most plausible force field to describe the underlying structure of a given dataset. The framework and associated software are applicable to a wide range of nanoscale simulations associated with experimental data with a hierarchical structure.

Mechanical Properties of α-Helices Estimated Using Molecular Dynamics and Finite Element Simulations

2008 ◽  
Author(s):  
Peyman Honarmandi ◽  
Philip Bransford ◽  
Roger D. Kamm
Keyword(s):  
Mechanical Properties ◽  
Molecular Dynamics ◽  
Finite Element ◽  
Solid Mechanics ◽  
Normal Mode Analysis ◽  
Computational Cost ◽  
Md Simulations ◽  
Elastic Stiffness ◽  
Mode Analysis ◽  
Molecular Volume

Mechanical properties of biomolecules and their response to mechanical forces may be studied using Molecular Dynamics (MD) simulations. However, high computational cost is a primary drawback of MD simulations. This paper presents a computational framework based on the integration of the Finite Element Method (FEM) with MD simulations to calculate the mechanical properties of polyalanine α-helix proteins. In this method, proteins are treated as continuum elastic solids with molecular volume defined exclusively by their atomic surface. Therefore, all solid mechanics theories would be applicable for the presumed elastic media. All-atom normal mode analysis is used to calculate protein’s elastic stiffness as input to the FEM. In addition, constant force molecular dynamics (CFMD) simulations can be used to predict other effective mechanical properties, such as the Poisson’s Ratio. Force versus strain data help elucidate the mechanical behavior of α-helices upon application of constant load. The proposed method may be useful in identifying the mechanical properties of any protein or protein assembly with known atomic structure.

scholarly journals Molecular dynamics simulations suggest stabilizing mutations in a de novo designed α/β protein

2019 ◽  
Vol 32 (7) ◽  
pp. 317-329
Author(s):  
Matthew Gill ◽  
Michelle E McCully
Keyword(s):  
Molecular Dynamics ◽  
Secondary Structure ◽  
Surface Area ◽  
De Novo ◽  
Md Simulations ◽  
Solvent Accessible Surface Area ◽  
Hydrophobic Core ◽  
Design Strategies ◽  
Major Interest ◽  
Core Packing

Abstract Designing functional proteins that can withstand extreme heat is beneficial for industrial and protein therapeutic applications. Thus, elucidating the atomic-level determinants of thermostability is a major interest for rational protein design. To that end, we compared the structure and dynamics of a set of previously designed, thermostable proteins based on the activation domain of human procarboxypeptidase A2 (AYEwt). The mutations in these designed proteins were intended to increase hydrophobic core packing and inter-secondary-structure interactions. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations of AYEwt and three designed variants at both 25 and 100°C. Our MD simulations agreed with the relative experimental stabilities of the designs based on their secondary structure content, Cα root-mean-square deviation/fluctuation, and buried-residue solvent accessible surface area. Using a contact analysis, we found that the designs stabilize inter-secondary structure interactions and buried hydrophobic surface area, as intended. Based on our analysis, we designed three additional variants to test the role of helix stabilization, core packing, and a Phe → Met mutation on thermostability. We performed the additional MD simulations and analysis on these variants, and these data supported our predictions.

scholarly journals An Adaptable Phospholipid Membrane Mimetic System for Solution NMR Studies of Membrane Proteins

2017 ◽  
Vol 139 (42) ◽  
pp. 14829-14832 ◽  
Author(s):  
Chih-Ta Henry Chien ◽  
Lukas R. Helfinger ◽  
Mark J. Bostock ◽  
Andras Solt ◽  
Yi Lei Tan ◽  
...  
Keyword(s):  
Membrane Proteins ◽  
Solution Nmr ◽  
Phospholipid Membrane ◽  
Nmr Studies ◽  
Membrane Mimetic
Sign in / Sign up

Export Citation Format

Share Document