Open AccessLipopolysaccharide Simulations Are Sensitive to Phosphate Charge and Ion Parameterization
Author(s) -
Amy Rice,
Mary Rooney,
Alexander I. Greenwood,
Myriam Cotten,
Jeff Wereszczynski
Publication year - 2020
Publication title -
journal of chemical theory and computation
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.001
H-Index - 185
eISSN - 1549-9626
pISSN - 1549-9618
DOI - 10.1021/acs.jctc.9b00868
Subject(s) - protonation , lipid bilayer , chemistry , deprotonation , lipid a , chemical physics , lipopolysaccharide , molecular dynamics , bacterial outer membrane , bilayer , membrane , molecule , hydrogen bond , ion , static electricity , computational chemistry , potential of mean force , biophysics , organic chemistry , biochemistry , physics , biology , escherichia coli , gene , endocrinology , quantum mechanics
The high proportion of lipopolysaccharide (LPS) molecules in the outer membrane of Gram-negative bacteria makes it a highly effective barrier to small molecules, antibiotic drugs, and other antimicrobial agents. Given this vital role in protecting bacteria from potentially hostile environments, simulations of LPS bilayers and outer membrane systems represent a critical tool for understanding the mechanisms of bacterial resistance and the development of new antibiotic compounds that circumvent these defenses. The basis of these simulations is parameterizations of LPS, which have been developed for all major molecular dynamics force fields. However, these parameterizations differ in both the protonation state of LPS and how LPS membranes behave in the presence of various ion species. To address these discrepancies and understand the effects of phosphate charge on bilayer properties, simulations were performed for multiple distinct LPS chemotypes with different ion parameterizations in both protonated or deprotonated lipid A states. These simulations show that bilayer properties, such as the area per lipid and inter-lipid hydrogen bonding, are highly influenced by the choice of phosphate group charges, cation type, and ion parameterization, with protonated LPS and monovalent cations with modified nonbonded parameters providing the best match to the experiments. Additionally, alchemical free energy simulations were performed to determine theoretical p K a values for LPS and subsequently validated by 31 P solid-state nuclear magnetic resonance experiments. Results from these complementary computational and experimental studies demonstrate that the protonated state dominates at physiological pH, contrary to the deprotonated form modeled by many LPS force fields. Overall, these results highlight the sensitivity of LPS simulations to phosphate charge and ion parameters while offering recommendations for how existing models should be updated for consistency between force fields as well as to best match experiments.