Adaptive free energy sampling in multidimensional collective variable space using boxed molecular dynamics

The past decade has seen the development of a new class of rare event methods in which molecular configuration space is divided into a set of boundaries/interfaces, and then short trajectories are run between boundaries. For all these methods, an important concern is how to generate boundaries. In this paper, we outline an algorithm for adaptively generating boundaries along a free energy surface in multi-dimensional collective variable (CV) space, building on the boxed molecular dynamics (BXD) rare event algorithm. BXD is a simple technique for accelerating the simulation of rare events and free energy sampling which has proven useful for calculating kinetics and free energy profiles in reactive and non-reactive molecular dynamics (MD) simulations across a range of systems, in both NVT and NVE ensembles. Two key developments outlined in this paper make it possible to automate BXD, and to adaptively map free energy and kinetics in complex systems. First, we have generalized BXD to multidimensional CV space. Using strategies from rigid-body dynamics, we have derived a simple and general velocity-reflection procedure that conserves energy for arbitrary collective variable definitions in multiple dimensions, and show that it is straightforward to apply BXD to sampling in multidimensional CV space so long as the Cartesian gradients ∇CV are available. Second, we have modified BXD to undertake on-the-fly statistical analysis during a trajectory, harnessing the information content latent in the dynamics to automatically determine boundary locations. Such automation not only makes BXD considerably easier to use; it also guarantees optimal boundaries, speeding up convergence. We have tested the multidimensional adaptive BXD procedure by calculating the potential of mean force for a chemical reaction recently investigated using both experimental and computational approaches - i.e., F + CD3CN → DF + D2CN in both the gas phase and a strongly coupled explicit CD3CN solvent. The results obtained using multidimensional adaptive BXD agree well with previously published experimental and computational results, providing good evidence for its reliability

Early dynamic changes to monocytes following major surgery are associated with subsequent infections

Background: Post-operative infections are a common cause of morbidity following major surgery. Little is understood about how major surgery perturbs immune function leading to heightened risk of subsequent infection. Through analysis of paired blood samples obtained immediately before and 24 h following surgery, we evaluated changes in circulating immune cell phenotype and function across the first 24 h, to identify early immune changes associated with subsequent infection. / Methods: We conducted a prospective observational study of adult patients undergoing major elective gastrointestinal, gynecological, or maxillofacial surgery requiring planned admission to the post-anesthetic care unit. Patients were followed up to hospital discharge or death. Outcome data collected included mortality, length of stay, unplanned intensive care unit admission, and post-operative infections (using the standardized endpoints in perioperative medicine–core outcome measures for perioperative and anesthetic care criteria). Peripheral blood mononuclear cells were isolated prior to and 24 h following surgery from which cellular immune traits including activation and functional status were assessed by multi-parameter flow cytometry and serum immune analytes compared by enzyme-linked immunosorbent assay (ELISA). / Results: Forty-eight patients were recruited, 26 (54%) of whom developed a post-operative infection. We observed reduced baseline pre- and post-operative monocyte CXCR4 and CD80 expression (chemokine receptors and co-stimulation markers, respectively) in patients who subsequently developed an infection as well as a profound and selective post-operative increase in CD4+ lymphocyte IL-7 receptor expression in the infection group only. Higher post-operative monocyte count was significantly associated with the development of post-operative infection (false discovery rate < 1%; adjusted p-value = 0.001) with an area under the receiver operating characteristic curve of 0.84 (p < 0.0001). / Conclusion: Lower monocyte chemotaxis markers, higher post-operative circulating monocyte counts, and reduced co-stimulatory signals are associated with subsequent post-operative infections. Identifying the underlying mechanisms and therapeutics to reverse defects in immune cell function requires further exploration

Early dynamic changes to monocytes following major surgery are associated with subsequent infections

BackgroundPost-operative infections are a common cause of morbidity following major surgery. Little is understood about how major surgery perturbs immune function leading to heightened risk of subsequent infection. Through analysis of paired blood samples obtained immediately before and 24 h following surgery, we evaluated changes in circulating immune cell phenotype and function across the first 24 h, to identify early immune changes associated with subsequent infection.MethodsWe conducted a prospective observational study of adult patients undergoing major elective gastrointestinal, gynecological, or maxillofacial surgery requiring planned admission to the post-anesthetic care unit. Patients were followed up to hospital discharge or death. Outcome data collected included mortality, length of stay, unplanned intensive care unit admission, and post-operative infections (using the standardized endpoints in perioperative medicine–core outcome measures for perioperative and anesthetic care criteria). Peripheral blood mononuclear cells were isolated prior to and 24 h following surgery from which cellular immune traits including activation and functional status were assessed by multi-parameter flow cytometry and serum immune analytes compared by enzyme-linked immunosorbent assay (ELISA).ResultsForty-eight patients were recruited, 26 (54%) of whom developed a post-operative infection. We observed reduced baseline pre- and post-operative monocyte CXCR4 and CD80 expression (chemokine receptors and co-stimulation markers, respectively) in patients who subsequently developed an infection as well as a profound and selective post-operative increase in CD4+ lymphocyte IL-7 receptor expression in the infection group only. Higher post-operative monocyte count was significantly associated with the development of post-operative infection (false discovery rate &lt; 1%; adjusted p-value = 0.001) with an area under the receiver operating characteristic curve of 0.84 (p &lt; 0.0001).ConclusionLower monocyte chemotaxis markers, higher post-operative circulating monocyte counts, and reduced co-stimulatory signals are associated with subsequent post-operative infections. Identifying the underlying mechanisms and therapeutics to reverse defects in immune cell function requires further exploration

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA

Sitting at the edge: How biomolecules use hydrophobicity to tune their interactions and function

Water near hydrophobic surfaces is like that at a liquid-vapor interface, where fluctuations in water density are substantially enhanced compared to that in bulk water. Here we use molecular simulations with specialized sampling techniques to show that water density fluctuations are similarly enhanced, even near hydrophobic surfaces of complex biomolecules, situating them at the edge of a dewetting transition. Consequently, water near these surfaces is sensitive to subtle changes in surface conformation, topology, and chemistry, any of which can tip the balance towards or away from the wet state, and thus significantly alter biomolecular interactions and function. Our work also resolves the long-standing puzzle of why some biological surfaces dewet and other seemingly similar surfaces do not.Comment: 12 pages, 4 figure

Computational Analysis of Phosphopeptide Binding to the Polo-Box Domain of the Mitotic Kinase PLK1 Using Molecular Dynamics Simulation

The Polo-Like Kinase 1 (PLK1) acts as a central regulator of mitosis and is over-expressed in a wide range of human tumours where high levels of expression correlate with a poor prognosis. PLK1 comprises two structural elements, a kinase domain and a polo-box domain (PBD). The PBD binds phosphorylated substrates to control substrate phosphorylation by the kinase domain. Although the PBD preferentially binds to phosphopeptides, it has a relatively broad sequence specificity in comparison with other phosphopeptide binding domains. We analysed the molecular determinants of recognition by performing molecular dynamics simulations of the PBD with one of its natural substrates, CDC25c. Predicted binding free energies were calculated using a molecular mechanics, Poisson-Boltzmann surface area approach. We calculated the per-residue contributions to the binding free energy change, showing that the phosphothreonine residue and the mainchain account for the vast majority of the interaction energy. This explains the very broad sequence specificity with respect to other sidechain residues. Finally, we considered the key role of bridging water molecules at the binding interface. We employed inhomogeneous fluid solvation theory to consider the free energy of water molecules on the protein surface with respect to bulk water molecules. Such an analysis highlights binding hotspots created by elimination of water molecules from hydrophobic surfaces. It also predicts that a number of water molecules are stabilized by the presence of the charged phosphate group, and that this will have a significant effect on the binding affinity. Our findings suggest a molecular rationale for the promiscuous binding of the PBD and highlight a role for bridging water molecules at the interface. We expect that this method of analysis will be very useful for probing other protein surfaces to identify binding hotspots for natural binding partners and small molecule inhibitors

Instability of aquaglyceroporin (Aqp) 2 contributes to drug resistance in trypanosoma brucei

Defining mode of action is vital for both developing new drugs and predicting potential resistance mechanisms. Sensitivity of African trypanosomes to pentamidine and melarsoprol is predominantly mediated by aquaglyceroporin 2 (TbAQP2), a channel associated with water/glycerol transport. TbAQP2 is expressed at the flagellar pocket membrane and chimerisation with TbAQP3 renders parasites resistant to both drugs. Two models for how TbAQP2 mediates pentamidine sensitivity have emerged; that TbAQP2 mediates pentamidine translocation across the plasma membrane or via binding to TbAQP2, with subsequent endocytosis and presumably transport across the endosomal/lysosomal membrane, but as trafficking and regulation of TbAQPs is uncharacterised this remains unresolved. We demonstrate that TbAQP2 is organised as a high order complex, is ubiquitylated and is transported to the lysosome. Unexpectedly, mutation of potential ubiquitin conjugation sites, i.e. cytoplasmic-oriented lysine residues, reduced folding and tetramerization efficiency and triggered ER retention. Moreover, TbAQP2/TbAQP3 chimerisation, as observed in pentamidine-resistant parasites, also leads to impaired oligomerisation, mislocalisation and increased turnover. These data suggest that TbAQP2 stability is highly sensitive to mutation and that instability contributes towards the emergence of drug resistance

Accommodating a Non-Conservative Internal Mutation by WaterMediated Hydrogen-Bonding Between β-Sheet Strands: A Comparison of Human and Rat Type B (Mitochondrial) Cytochrome b5

Mammalian type B (mitochondrial) cytochromes b5 exhibit greater amino acid sequence diversity than their type A (microsomal) counterparts, as exemplified by the type B proteins from human (hCYB5B) and rat (rCYB5B). The comparison of X-ray crystal structures of hCYB5B and rCYB5B reported herein reveals a striking difference in packing involving the five-stranded β-sheet, attributable to fully buried residue 21 in strand β4. The greater bulk of Leu21 in hCYB5B in comparison to Thr21 in rCYB5B results in a substantial displacement of the first two residues in β5, and consequent loss of two of the three hydrogen bonds between β5 and β4. Hydrogen-bonding between the residues is instead mediated by two well-ordered, fully buried water molecules. In a 10 ns molecular dynamics simulation, one of the buried water molecules in the hCYB5B structure exchanged readily with solvent via intermediates having three water molecules sandwiched between β4 and β5. When the buried water molecules were removed prior to a second 10 ns simulation, β4 and β5 formed persistent hydrogen bonds identical to those in rCYB5B, but the Leu21 side chain was forced to adopt a rarely observed conformation. Despite the apparently greater ease of water access to the interior of hCYB5B than of rCYB5B suggested by these observations, the two proteins exhibit virtually identical stability, dynamic and redox properties. The results provide new insight into the factors stabilizing the cytochrome b5 fold