Glycolipids (such as glycoglycerolipids, glycosphingolipids, and glycosylphosphatidylinositol) and lipoglycans (such as lipopolysaccharides (LPS), lipooligosaccharides (LOS), mycobacterial lipoarabinomannan, and mycoplasma lipoglycans) are typically found on the surface of cell membranes and play crucial roles in various cellular functions. Characterizing their structure and dynamics at the molecular level is essential to understand their biological roles, but systematic generation of glycolipid and lipoglycan structures is challenging because of great variations in lipid structures and glycan sequences (i.e., carbohydrate types and their linkages). To facilitate the generation of all-atom glycolipid/LPS/LOS structures, we have developed Glycolipid Modeler and LPS Modeler in CHARMM-GUI (http://www.charmm-gui.org), a webbased interface that simplifies building of complex biological simulation systems. In addition, we have incorporated these modules into Membrane Builder, so that users can readily build a complex symmetric or asymmetric biological membrane system with various glycolipids and LPS/LOS. These tools are expected to be useful in innovative and novel glycolipid/LPS/LOS modeling and simulation research by easing tedious and intricate steps in modeling complex biological systems, and shall provide insight into structures, dynamics, and underlying mechanisms of complex glycolipid-/LPS-/LOS-containing biological membrane systems.

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the molecular details of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project.

A complex cell envelope, composed of a mixture of lipid types including complex lipopolysaccharides, protects bacteria from the external environment. Clearly, the proteins embedded within the various components of the cell envelope have an intricate relationship with their local environment. Therefore, to obtain meaningful results, molecular simulations need to mimic as far as possible this chemically heterogeneous system. However, setting up such systems for computational studies is far from trivial, and consequently the vast majority of simulations of outer membrane proteins still rely on oversimplified phospholipid membrane models. This work presents an update of CHARMM-GUI Martini Maker for coarse-grained modeling and simulation of complex bacterial membranes with lipopolysaccharides. The qualities of the outer membrane systems generated by Martini Maker are validated by simulating them in bilayer, vesicle, nanodisc, and micelle environments (with and without outer membrane proteins) using the Martini force field. We expect this new feature in Martini Maker to be a useful tool for modeling large, complicated bacterial outer membrane systems in a user-friendly manner.

Lipid A is the lipid anchor of a lipopolysaccharide in the outer leaflet of the outer membrane of Gram-negative bacteria. In general, lipid A consists of two phosphorylated N-acetyl glucosamine and several acyl chains that are directly linked to the two sugars. Depending on the bacterial species and environments, the acyl chain number and length vary, and lipid A can be chemically modified with phosphoethanolamine, aminoarabinose, or glycine residues, which are key to bacterial pathogenesis. In this work, homogeneous lipid bilayers of 21 distinct lipid A types from 12 bacterial species are modeled and simulated to investigate the differences and similarities of their membrane properties. In addition, different neutralizing ion types (Ca 2þ , K þ , and Na þ ) are considered to examine the ion's influence on the membrane properties. The trajectory analysis shows that (1) the area per lipid is mostly correlated to the acyl chain number, and the area per lipid increases as a function of the acyl chain number; (2) the hydrophobic thickness is mainly determined by the average acyl chain length with slight dependence on the acyl chain number, and the hydrophobic thickness generally increases with the average acyl chain length; (3) a good correlation is observed among the area per lipid, hydrophobic thickness, and acyl chain order; and (4) although the influence of neutralizing ion types on the area per lipid and hydrophobic thickness is minimal, Ca 2þ stays longer on the membrane surface than K þ or Na þ , consequently leading to lower lateral diffusion and a higher compressibility modulus, which agrees well with available experiments.

A polarizable empirical force field based on the classical Drude oscillator is presented for the hexopyranose form of selected monosaccharides. Parameter optimization targeted quantum mechanical (QM) dipole moments, solute - water interaction energies, vibrational frequencies, and conformational energies. Validation of the model was based on experimental data on crystals, densities of aqueous-sugar solutions, diffusion constants of glucose, and rotational preferences of the exocylic hydroxymethyl of D-glucose and D-galactose in aqueous solution as well as additional QM data. Notably, the final model involves a single electrostatic model for all sixteen diastereomers of the monosaccharides, indicating the transferability of the polarizable model. The presented parameters are anticipated to lay the foundation for a comprehensive polarizable force field for saccharides that will be compatible with the polarizable Drude parameters for lipids and proteins, allowing for simulations of glycolipids and glycoproteins.

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.

Conformational sampling for a set of ten α- or β-(1→6)-linked oligosaccharides has been studied using explicit solvent Hamiltonian Replica Exchange (HREX) simulations and NMR spectroscopy techniques. Validation of the force field and simulation methodology is done by comparing calculated transglycosidic J-coupling constants and proton-proton distances with the corresponding NMR data. Initial calculations showed poor agreement, for example, with > 3 Hz deviation of the calculated 3J(H5,H6R) values from the experimental data, prompting optimization of the ω torsion angle parameters associated with (1→6)-linkages. The resulting force field is in overall good agreement (i.e. within ~ 0.5 Hz deviation) from experimental 3J(H5,H6R) values, although some small limitations are evident. Detailed hydrogen bonding analysis indicates that most of the compounds lack direct intramolecular H-bonds between the two monosaccharides; however, minor sampling of the O6· · ·HO2′ hydrogen bond is present in three compounds. The results verify the role of the gauche effect between O5 and O6 atoms in gluco- and manno-configured pyranosides causing the ω torsion angle to sample an equilibrium between the gt and gg rotamers. Conversely, galacto-configured pyranosides sample a population distribution in equilibrium between gt and tg rotamers, while the gg rotamer populations are minor. Water radial distribution functions suggest decreased accessibility to the O6 atom in the (1→6)-linkage as compared to the O6' atom in the non-reducing sugar. The role of bridging water molecules between two sugar moieties on the distributions of ω torsion angles in oligosaccharides is also explored.

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.