Learning Force Field Parameters from Differentiable Particle-Field Molecular Dynamics.

We develop ∂-HylleraasMD (∂-HyMD), a fully end-to-end differentiable molecular dynamics software based on the Hamiltonian hybrid particle-field formalism, and use it to establish a protocol for automated optimization of force field parameters. ∂-HyMD is templated on the recently released HylleraaasMD software, while using the JAX autodiff framework as the main engine for the differentiable dynamics. ∂-HyMD exploits an embarrassingly parallel optimization algorithm by spawning independent simulations, whose trajectories are simultaneously processed by reverse mode automatic differentiation to calculate the gradient of the loss function, which is in turn used for iterative optimization of the force-field parameters. We show that parallel organization facilitates the convergence of the minimization procedure, avoiding the known memory and numerical stability issues of differentiable molecular dynamics approaches. We showcase the effectiveness of our implementation by producing a library of force field parameters for standard phospholipids, with either zwitterionic or anionic heads and with saturated or unsaturated tails. Compared to the all-atom reference, the force field obtained by ∂-HyMD yields better density profiles than the parameters derived from previously utilized gradient-free optimization procedures. Moreover, ∂-HyMD models can predict with good accuracy properties not included in the learning objective, such as lateral pressure profiles, and are transferable to other systems, including triglycerides.

Accelerating Lipid Flip-Flop at Low Concentrations: A General Mechanism for Membrane Binding Peptides.

We report a physicochemical investigation of the lipid transport properties of model lipid membranes in the presence of the antimicrobial peptide indolicidin through comparisons of experimental SANS/SAXS scattering techniques to fully atomistic molecular dynamics simulations. In agreement with the experiment, we show that upon peripheral binding of the peptides, even at low concentrations, lipid flip-flop dynamics is greatly accelerated. Computer modeling elucidates the interplay between structural changes and lipid dynamics induced by peptides and proposes a mechanism for the mode of action of antimicrobial peptides, assessing the major role of entropy for the catalysis of the flipping events. The mechanism introduced here is universal for all peptides with preferential peripheral binding to the membrane as it does not depend on the specific amino acid sequence.

The structural basis of the multi-step allosteric activation of Aurora B kinase.

Aurora B, together with IN-box, the C-terminal part of INCENP, forms an enzymatic complex that ensures faithful cell division. The [Aurora B/IN-box] complex is activated by autophosphorylation in the Aurora B activation loop and in IN-box, but it is not clear how these phosphorylations activate the enzyme. We used a combination of experimental and computational studies to investigate the effects of phosphorylation on the molecular dynamics and structure of [Aurora B/IN-box]. In addition, we generated partially phosphorylated intermediates to analyze the contribution of each phosphorylation independently. We found that the dynamics of Aurora and IN-box are interconnected, and IN-box plays both positive and negative regulatory roles depending on the phosphorylation status of the enzyme complex. Phosphorylation in the activation loop of Aurora B occurs intramolecularly and prepares the enzyme complex for activation, but two phosphorylated sites are synergistically responsible for full enzyme activity.

HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter.

We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.

What Drives Chorismate Mutase to Top Performance? Insights from a Combined In Silico and In Vitro Study.

Unlike typical chorismate mutases, the enzyme from Mycobacterium tuberculosis (MtCM) has only low activity on its own. Remarkably, its catalytic efficiency kcat/Km can be boosted more than 100-fold by complex formation with a partner enzyme. Recently, an autonomously fully active MtCM variant was generated using directed evolution, and its structure was solved by X-ray crystallography. However, key residues were involved in crystal contacts, challenging the functional interpretation of the structural changes. Here, we address these challenges by microsecond molecular dynamics simulations, followed up by additional kinetic and structural analyses of selected sets of specifically engineered enzyme variants. A comparison of wild-type MtCM with naturally and artificially activated MtCMs revealed the overall dynamic profiles of these enzymes as well as key interactions between the C-terminus and the active site loop. In the artificially evolved variant of this model enzyme, this loop is preorganized and stabilized by Pro52 and Asp55, two highly conserved residues in typical, highly active chorismate mutases. Asp55 stretches across the active site and helps to appropriately position active site residues Arg18 and Arg46 for catalysis. The role of Asp55 can be taken over by another acidic residue, if introduced at position 88 close to the C-terminus of MtCM, as suggested by molecular dynamics simulations and confirmed by kinetic investigations of engineered variants.

Can the roles of polar and non-polar moieties be reversed in non-polar solvents?

Using thermodynamic integration, we study the solvation free energy of 18 amino acid side chain equivalents in solvents with different polarities, ranging from the most polar water to the most non-polar cyclohexane. The amino acid side chain equivalents are obtained from the 20 natural amino acids by replacing the backbone part with a hydrogen atom, and discarding proline and glycine that have special properties. A detailed analysis of the relative solvation free energies suggests how it is possible to achieve a robust and unambiguous hydrophobic scale for the amino acids. By discriminating the relative contributions of the entropic and enthalpic terms, we find strong negative correlations in water and ethanol, associated with the well-known entropy-enthalpy compensation, and a much reduced correlation in cyclohexane. This shows that in general the role of the polar and non-polar moieties cannot be reversed in a non-polar solvent. Our findings are compared with past experimental as well as numerical results, and may shed additional light on the unique role of water as a biological solvent.

Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents?

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse-grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states coexist in solution. Our simulations also provide indications as to how to escape this region to produce thermodynamically stable micellar aggregates.