CHARMM-GUI Supports Hydrogen Mass Repartitioning and Different Protonation States of Phosphates in Lipopolysaccharides

Abstract

Hydrogen mass repartitioning (HMR) that permits time steps of all-atom 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. Consistence 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 PO4 1or PO4 2group, but different phosphate protonation states make slight impacts on lipid packing and conformational properties of LPS acyl chains. Systems with PO4 1have slightly smaller area per lipid and thus slightly more ordered lipid A acyl chains compared to those with PO4 2-, due to more electrostatic repulsion between PO4 2even with neutralizing Ca 2+ 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.

Binding of Human ACE2 and RBD of Omicron Enhanced by Unique Interaction Patterns Among SARS-CoV-2 Variants of Concern

“…The 4 fs simulation time-step was used with the hydrogen mass repartitioning method. 25,26 Experimental Methods…”

CHARMM-GUI LBS Finder & Refiner for Ligand Binding Site Prediction and Refinement

“…56 In addition, we support hydrogen mass repartitioning (HMR) as the default option to accelerate MD simulations with 4 fs time steps. 69,70 ■ APPLICATIONS For our benchmark tests, we selected data sets from the directory of useful decoys-enhanced (DUD-E) and Astex. 71,72 For the computationally predicted structures benchmark test, we obtained protein sequences from DUD-E that have amino acid residues <300.…”

“…Then, users can generate simulation input files together with the LBS (Cα + SC-COM or Cα-only) refinement distance restraints for 5 different MD simulation programs (NAMD, AMBER, GROMACS, OpenMM, and GENESIS) . In addition, we support hydrogen mass repartitioning (HMR) as the default option to accelerate MD simulations with 4 fs time steps. , …”

Insights into the structural peculiarities of the N-terminal and receptor binding domains of the spike protein from the SARS-CoV-2 Omicron variant

“…Since cation-π interactions play an important role in biomolecular identification processes and stability of protein structure, WYF parameter was used to improve the MD simulation, which contains free energies, of cation-indole-choline and indole-trimethylated lysines interactions related to protein-protein interfaces [ 31 ]. Hydrogen mass repartitioning (HMR) technique permits all-atom MD simulation to employ time steps of up to 4 fs by devoting a heavy atom mass to hydrogen atom(s) linked to the heavy atom [ 32 ]. The linear-constraint-solving (LINCS) algorithms was used to limit bonds involving hydrogen atoms [ 33 , 34 ].…”