Hydrogen mass repartitioning (HMR) that permits time steps of allatom molecular dynamics simulation up to 4 fs by increasing the mass of hydrogen atoms has been used in protein and phospholipid bilayers simulations to improve conformational sampling. Molecular simulation input generation via CHARMM-GUI now supports HMR for diverse simulation programs. In addition, considering ambiguous pH at the bacterial outer membrane surface, different protonation states, either -2e or -1e, of phosphate groups in lipopolysaccharides (LPS) are also supported in CHARMM-GUI LPS Modeler. To examine the robustness of HMR and the influence of protonation states of phosphate groups on LPS bilayer properties, eight different LPS-type all-atom systems with two phosphate protonation states are modeled and simulated utilizing both OpenMM 2-fs (standard) and 4-fs (HMR) schemes. Consistency in the conformational space sampled by standard and HMR simulations shows the reliability of HMR even in LPS, one of the most complex biomolecules. For systems with different protonation states, similar conformations are sampled with a PO41- or PO(4)(-)(2)group, but different phosphate protonation states make slight impacts on lipid packing and conformational properties of LPS acyl chains. Systems with PO41- have a slightly smaller area per lipid and thus slightly more ordered lipid A acyl chains compared to those with PO42-, due to more electrostatic repulsion between PO42- even with neutralizing Ca2+ ions. HMR and different protonation states of phosphates of LPS available in CHARMM-GUI are expected to be useful for further investigations of biological systems of diverse origin.

Lipopolysaccharides (LPS) present in the outer leaflet of Gram-negative bacterial outer membranes protect the bacteria from external threats and influence antibiotic permeability as well as immune system recognition. The structure of lipid A, the anchor of an LPS molecule to the outer membrane, can make direct influences on membrane properties. Particularly, in Vibrio cholerae, a Gram-negative bacterium responsible for cholera, a severe diarrheal disease, modifications of lipid A structures grant antibiotic resistance and are a primary factor that led to the current cholera pandemic. However, the difference in structural properties incurred by such modifications has not been fully explored. In this work, five symmetric bilayer systems comprised of distinct lipid A structures of Vibrio cholerae LPS with O1 O-antigen were modeled and simulated to explore influences of different lipid A types on membrane properties. All-atom molecular dynamics simulations reveal that membrane properties such as hydrophobic thickness, acyl chain order parameter, and area per lipid are largely impacted by lipid A modifications due to differences in composition and acyl chain distortions. The modified lipid A is also less negatively charged, which possibly reveals a resistance mechanism to cationic antimicrobial peptide evasion. These findings present a possible explanation for Vibrio cholerae's immune system evasion properties and establish the differences between the lipid A types, which should be of use for any future study of the Gram-negative bacteria.