skip to main content

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol—and consequently, its perceived accuracy—must be based primarily on the scientific question that the simulations are designed to address.

  more » « less
Award ID(s):
1817929
PAR ID:
10469786
Author(s) / Creator(s):
; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Membranes
Volume:
13
Issue:
7
ISSN:
2077-0375
Page Range / eLocation ID:
629
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Soft Matter)
    null (Ed.)
    The plasma membrane of eukaryotic cells is known to be compositionally asymmetric. Certain phospholipids, such as sphingomyelin and phosphatidylcholine species, are predominantly localized in the outer leaflet, while phosphatidylethanolamine and phosphatidylserine species primarily reside in the inner leaflet. While phospholipid asymmetry between the membrane leaflets is well established, there is no consensus about cholesterol distribution between the two leaflets. We have performed a systematic study, via molecular simulations, of how the spatial distribution of cholesterol molecules in different “asymmetric” lipid bilayers are affected by the lipids’ backbone, head-type, unsaturation, and chain-length by considering an asymmetric bilayer mimicking the plasma membrane lipids of red blood cells, as well as seventeen other asymmetric bilayers comprising of different lipid types. Our results reveal that the distribution of cholesterol in the leaflets is solely a function of the extent of ordering of the lipids within the leaflets. The ratio of the amount of cholesterol matches the ratio of lipid order in the two leaflets, thus providing a quantitative relationship between the two. These results are understood by the observation that asymmetric bilayers with equimolar amount of lipids in the two leaflets develop tensile and compressive stresses due to differences in the extent of lipid order. These stresses are alleviated by the transfer of cholesterol from the leaflet in compressive stress to the one in tensile stress. These findings are important in understanding the biology of the cell membrane, especially with regard to the composition of the membrane leaflets. 
    more » « less
  2. ; ; ( , Proteins: Structure, Function, and Bioinformatics)
    Abstract

    Protein structure refinement is the last step in protein structure prediction pipelines. Physics‐based refinement via molecular dynamics (MD) simulations has made significant progress during recent years. During CASP14, we tested a new refinement protocol based on an improved sampling strategy via MD simulations. MD simulations were carried out at an elevated temperature (360 K). An optimized use of biasing restraints and the use of multiple starting models led to enhanced sampling. The new protocol generally improved the model quality. In comparison with our previous protocols, the CASP14 protocol showed clear improvements. Our approach was successful with most initial models, many based on deep learning methods. However, we found that our approach was not able to refine machine‐learning models from the AlphaFold2 group, often decreasing already high initial qualities. To better understand the role of refinement given new types of models based on machine‐learning, a detailed analysis via MD simulations and Markov state modeling is presented here. We continue to find that MD‐based refinement has the potential to improve AI predictions. We also identified several practical issues that make it difficult to realize that potential. Increasingly important is the consideration of inter‐domain and oligomeric contacts in simulations; the presence of large kinetic barriers in refinement pathways also continues to present challenges. Finally, we provide a perspective on how physics‐based refinement could continue to play a role in the future for improving initial predictions based on machine learning‐based methods.

     
    more » « less
  3. ; ; ; ; ; ; ; ; ; ( , Nanoscale)
    Membrane undulations play a vital role in many biological processes, including the regulation of membrane protein activity. The asymmetric lipid composition of most biological membranes complicates theoretical description of these bending fluctuations, yet experimental data that would inform any such a theory is scarce. Here, we used neutron spin-echo (NSE) spectroscopy to measure the bending fluctuations of large unilamellar vesicles (LUV) having an asymmetric transbilayer distribution of high- and low-melting lipids. The asymmetric vesicles were prepared using cyclodextrin-mediated lipid exchange, and were composed of an outer leaflet enriched in egg sphingomyelin (ESM) and an inner leaflet enriched in 1-palmitoyl-2-oleoyl-phosphoethanolamine (POPE), which have main transition temperatures of 37 °C and 25 °C, respectively. The overall membrane bending rigidity was measured at three temperatures: 15 °C, where both lipids are in a gel state; 45 °C, where both lipids are in a fluid state; and 30 °C, where there is gel-fluid co-existence. Remarkably, the dynamics for the fluid asymmetric LUVs (aLUVs) at 30 °C and 45 °C do not follow trends predicted by their symmetric counterparts. At 30 °C, compositional asymmetry suppressed the bending fluctuations, with the asymmetric bilayer exhibiting a larger bending modulus than that of symmetric bilayers corresponding to either the outer or inner leaflet. We conclude that the compositional asymmetry and leaflet coupling influence the internal dissipation within the bilayer and result in membrane properties that cannot be directly predicted from corresponding symmetric bilayers. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Membranes)
    KCNE3 is a potassium channel accessory transmembrane protein that regulates the function of various voltage-gated potassium channels such as KCNQ1. KCNE3 plays an important role in the recycling of potassium ion by binding with KCNQ1. KCNE3 can be found in the small intestine, colon, and in the human heart. Despite its biological significance, there is little information on the structural dynamics of KCNE3 in native-like membrane environments. Molecular dynamics (MD) simulations are a widely used as a tool to study the conformational dynamics and interactions of proteins with lipid membranes. In this study, we have utilized all-atom molecular dynamics simulations to characterize the molecular motions and the interactions of KCNE3 in a bilayer composed of: a mixture of POPC and POPG lipids (3:1), POPC alone, and DMPC alone. Our MD simulation results suggested that the transmembrane domain (TMD) of KCNE3 is less flexible and more stable when compared to the N- and C-termini of KCNE3 in all three membrane environments. The conformational flexibility of N- and C-termini varies across these three lipid environments. The MD simulation results further suggested that the TMD of KCNE3 spans the membrane width, having residue A69 close to the center of the lipid bilayers and residues S57 and S82 close to the lipid bilayer membrane surfaces. These results are consistent with previous biophysical studies of KCNE3. The outcomes of these MD simulations will help design biophysical experiments and complement the experimental data obtained on KCNE3 to obtain a more detailed understanding of its structural dynamics in the native membrane environment. 
    more » « less
  5. ; ; ; ( , Advanced Biosystems)
    Abstract

    Cellular membranes play host to a wide variety of morphologically and chemically complex processes. Although model membranes, like liposomes, are already widely used to reconstitute and study these processes, better tools are needed for making model bilayers that faithfully mimic cellular membranes. Existing methods for fabricating cell‐sized (µm) or organelle‐sized (tens to hundreds of nanometers) lipid vesicles have distinctly different requirements. Of particular note for biology, it remains challenging for any technique to efficiently encapsulate fragile cargo molecules or to generate liposomes with stable, asymmetric lipid leaflets within the bilayer. Here a tunable microfluidic device and protocol for fabricating liposomes with desired diameters ranging from ≈10 µm to ≈100 nm are described. Lipid vesicle size is templated by the simple inclusion of a polycarbonate filter within the microfluidic system and tuned with flow rate. It is shown that the vesicles made with this device are stable, unilamellar, lipid asymmetric, and capable of supporting transmembrane protein assembly, peripheral membrane protein binding, as well as soluble cargo encapsulation (including designer nanocages for biotechnology applications). These fabricated vesicles provide a new platform for studying the biophysically rich processes found within lipid–lipid and lipid–protein systems typically associated with cellular membranes.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email