The quality of biomolecular dynamics simulations relies critically on the force field that is used to describe the interactions between particles in the system. Force fields, which are generally parameterized using experimental data on small molecules, can only prove themselves in realistic simulations of relevant biomolecular systems. In this work, we begin the validation of the new 53A6 GROMOS parameter set by examining three test cases. Simulations of the well-studied 129 residue protein hen egg-white lysozyme, of the DNA dodecamer d(CGCGAATTCGCG)(2), and a proteinogenic beta(3)-dodecapeptide were performed and analysed. It was found that the new parameter set performs as well as the previous parameter sets in terms of protein (45A3) and DNA (45A4) stability and that it is better at describing the folding-unfolding balance of the peptide. The latter is a property that is directly associated with the free enthalpy of hydration, to which the 53A6 parameter set was parameterized.

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Over the past decades, the GROMOS force field for biomolecular simulation has primarily been developed for performing molecular dynamics (MD) simulations of polypeptides and, to a lesser extent, sugars. When applied to DNA, the 43A1 and 45A3 parameter sets of the years 1996 and 2001 produced rather flexible double-helical structures, in which the Watson-Crick hydrogen-bonding content was more limited than expected. To improve on the currently available parameter sets, the nucleotide backbone torsional-angle parameters and the charge distribution of the nucleotide bases are reconsidered based on quantum-chemical data. The new 45A4 parameter set resulting from this refinement appears to perform well in terms of reproducing solution NMR data and canonical hydrogen bonding. The deviation between simulated and experimental observables is now of the same order of magnitude as the uncertainty in the experimental values themselves.

No abstract

Lipopolysaccharides (LPS) comprise the outermost layer of the Gram-negative bacteria cell envelope. Packed onto a lipid layer, the outer membrane displays remarkable physical-chemical differences compared to cell membranes. The carbohydrate-rich region confers a membrane asymmetry that underlies many biological processes such as endotoxicity, antibiotic resistance, and cell adhesion. Furthermore, unlike membrane proteins from other sources, integral outer-membrane proteins do not consist of transmembrane α helices; instead they consist of antiparallel β-barrels, which highlights the importance of the LPS membrane as a medium. In this work, we present an extension of the GLYCAM06 force field that has been specifically developed for LPS membranes using our Wolf2Pack program. This new set of parameters for lipopolysaccharide molecules expands the GLYCAM06 repertoire of monosaccharides to include phosphorylated N- and O-acetylglucosamine, 3-deoxy-d-manno-oct-2-ulosonic acid, l-glycero-D-manno-heptose and its O-carbamoylated variant, and N-alanine-d-galactosamine. A total of 1 μs of molecular dynamics simulations of the rough LPS membrane of Pseudomonas aeruginosa PA01 is used to showcase the added parameter set. The equilibration of the LPS membrane is shown to be significantly slower compared to phospholipid membranes, on the order of 500 ns. It is further shown that water molecules penetrate the hydrocarbon region up to the terminal methyl groups, much deeper than commonly observed for phospholipid bilayers, and in agreement with neutron diffraction measurements. A comparison of simulated structural, dynamical, and electrostatic properties against corresponding experimentally available data shows that the present parameter set reproduces well the overall structure and the permeability of LPS membranes in the liquid-crystalline phase.