CHARMM-GUI Membrane Builder, http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulation through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builderthat allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types such as phosphoinositides, cardiolipin, sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: cardiolipin lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From (1)H,(1)H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.

Gangliosides are a class of glycosphingolipids (GSLs) with amphiphilic character that are found at the outer leaflet of the cell membranes, where their ability to organize into special domains makes them vital cell membrane components. However, a molecular understanding of GSL-rich membranes in terms of their clustered organization, stability, and dynamics is still elusive. To gain molecular insight into the organization and dynamics of GSL-rich membranes, we performed all-atom molecular-dynamics simulations of bicomponent ganglioside GM1 in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) phospholipid bilayers with varying concentrations of GM1 (10%, 20%, and 30%). Overall, the simulations show very good agreement with available experimental data, including x-ray electron density profiles along the membrane normal, NMR carbohydrate proton-proton distances, and x-ray crystal structures. This validates the quality of our model systems for investigating GM1 clustering through an ordered-lipid-cluster analysis. The increase in GM1 concentration induces tighter lipid packing, driven mainly by inter-GM1 carbohydrate-carbohydrate interactions, leading to a greater preference for the positive curvature of GM1-containing membranes and larger cluster sizes of ordered-lipid clusters (with a composite of GM1 and POPC). These clusters tend to segregate and form a large percolated cluster at a 30% GM1 concentration at 293 K. At a higher temperature of 330 K, however, the segregation is not maintained.

The asymmetric outer membrane of Gram-negative bacteria is formed of the inner leaflet with phospholipids and the outer leaflet with lipopolysaccharides (LPS). Outer membrane protein F (OmpF) is a trimeric porin responsible for the passive transport of small molecules across the outer membrane of Escherichia coli. Here, we report the impact of different levels of heterogeneity in LPS environments on the structure and dynamics of OmpF using all-atom molecular dynamics simulations. The simulations provide insight into the flexibility and dynamics of LPS components that are highly dependent on local environments, with lipid A being the most rigid and O-antigen being the most flexible. Increased flexibility of O-antigen polysaccharides is observed in heterogeneous LPS systems, where the adjacent O-antigen repeating units are weakly interacting and thus more dynamic, compared to homogeneous LPS systems in which LPS interacts strongly with each other with limited overall flexibility due to dense packing. The model systems were validated by comparing molecular-level details of interactions between OmpF surface residues and LPS core sugars with experimental data, establishing the importance of LPS core oligosaccharides in shielding OmpF surface epitopes recognized by monoclonal antibodies. There are LPS environmental influences on the movement of bulk ions (K(+) and Cl(-)), but the ion selectivity of OmpF is mainly affected by bulk ion concentration.

Compstatin is a 13-residue disulfide-bridged peptide that inhibits a key step in the activation of the human complement system. Compstatin and its derivatives have shown great promise for the treatment of many clinical disorders associated with unbalanced complement activity. To obtain more potent compstatin analogues, we have now performed an N-methylation scan of the peptide backbone and amino acid substitutions at position 13. One analogue (Ac-I[CVW(Me)QDW-Sar-AHRC](NMe)I-NH2) displayed a 1,000-fold increase in both potency (IC50=62 nM) and binding affinity for C3b (KD=2.3 nM) over that of the original compstatin. Biophysical analysis using surface plasmon resonance and isothermal titration calorimetry suggests that the improved binding originates from more favorable free conformation and stronger hydrophobic interactions. This study provides a series of significantly improved drug leads for therapeutic applications in complement-related diseases, and offers new insights into the structure-activity relationships of compstatin analogues.

The outer membrane of Gram-negative bacteria is an asymmetric membrane with lipopolysaccharides on the external leaflet and phospholipids on the periplasmic leaflet. This outer membrane contains mainly β-barrel transmembrane proteins and lipidated periplasmic proteins (lipoproteins). The multisubunit protein β-barrel assembly machine (BAM) catalyzes the insertion and folding of the β-barrel proteins into this membrane. In Escherichia coli, the BAM complex consists of five subunits, a core transmembrane β-barrel with a long periplasmic domain (BamA) and four lipoproteins (BamB/C/D/E). The BamA periplasmic domain is composed of five globular subdomains in tandem called POTRA motifs that are key to BAM complex formation and interaction with the substrate β-barrel proteins. The BAM complex is believed to undergo conformational cycling while facilitating insertion of client proteins into the outer membrane. Reports describing variable conformations and dynamics of the periplasmic POTRA domain have been published. Therefore, elucidation of the conformational dynamics of the POTRA domain in full-length BamA is important to understand the function of this molecular complex. Using molecular dynamics simulations, we present evidence that the conformational flexibility of the POTRA domain is modulated by binding to the periplasmic surface of a native lipid membrane. Furthermore, membrane binding of the POTRA domain is compatible with both BamB and BamD binding, suggesting that conformational selection of different POTRA domain conformations may be involved in the mechanism of BAM-facilitated insertion of outer membrane β-barrel proteins.

Summary While membrane simulations are widely employed to study the structure and dynamics of various lipid bilayers and membrane proteins in the bilayers, simulations of lipopolysaccharides (LPS) in membrane environments have been limited due to its structural complexity, difficulties in building LPS-membrane systems, and lack of appropriate molecular force field. In this work, as a first step to extend CHARMM-GUI Membrane Builder to incorporate LPS molecules and to explore their structures and dynamics in membrane environments using molecular dynamics simulations, we describe step-by-step procedures to build LPS bilayer systems using CHARMM and the recently developed CHARMM carbohydrate and lipid force fields. Such procedures are illustrated by building various bilayers of Escherichia coli O6 LPS and their preliminary simulation results are given in terms of per-LPS area and density distributions of various components along the membrane normal.