Equilibrium molecular thermodynamics from Kirkwood sampling.

We present two methods for barrierless equilibrium sampling of molecular systems based on the recently proposed Kirkwood method (J. Chem. Phys. 2009, 130, 134102). Kirkwood sampling employs low-order correlations among internal coordinates of a molecule for random (or non-Markovian) sampling of the high dimensional conformational space. This is a geometrical sampling method independent of the potential energy surface. The first method is a variant of biased Monte Carlo, where Kirkwood sampling is used for generating trial Monte Carlo moves. Using this method, equilibrium distributions corresponding to different temperatures and potential energy functions can be generated from a given set of low-order correlations. Since Kirkwood samples are generated independently, this method is ideally suited for massively parallel distributed computing. The second approach is a variant of reservoir replica exchange, where Kirkwood sampling is used to construct a reservoir of conformations, which exchanges conformations with the replicas performing equilibrium sampling corresponding to different thermodynamic states. Coupling with the Kirkwood reservoir enhances sampling by facilitating global jumps in the conformational space. The efficiency of both methods depends on the overlap of the Kirkwood distribution with the target equilibrium distribution. We present proof-of-concept results for a model nine-atom linear molecule and alanine dipeptide.This research was funded by the European Research Council and EPSRC grant EP/I001352/1. Y.O. was supported, in part, by the JSPS Grant-in-Aid for Scientific Research on Innovative Areas (“Dynamical Ordering and Integrated Functions”).This is the final published version. It first appeared at http://pubs.acs.org/doi/abs/10.1021/acs.jpcb.5b01800

Characterizing Loop Dynamics and Ligand Recognition in Human- and Avian-Type Influenza Neuraminidases via Generalized Born Molecular Dynamics and End-Point Free Energy Calculations

The comparative dynamics and inhibitor binding free energies of group-1 and group-2 pathogenic influenza A subtype neuraminidase (NA) enzymes are of fundamental biological interest and relevant to structure-based drug design studies for antiviral compounds. In this work, we present seven generalized Born molecular dynamics simulations of avian (N1)- and human (N9)-type NAs in order to probe the comparative flexibility of the two subtypes, both with and without the inhibitor oseltamivir bound. The enhanced sampling obtained through the implicit solvent treatment suggests several provocative insights into the dynamics of the two subtypes, including that the group-2 enzymes may exhibit similar motion in the 430-binding site regions but different 150-loop motion. End-point free energy calculations elucidate the contributions to inhibitor binding free energies and suggest that entropic considerations cannot be neglected when comparing across the subtypes. We anticipate the findings presented here will have broad implications for the development of novel antiviral compounds against both seasonal and pandemic influenza strains

c-Met activation leads to the establishment of a TGFβ-receptor regulatory network in bladder cancer progression

Treatment of muscle-invasive bladder cancer remains a major clinical challenge. Aberrant HGF/c-MET upregulation and activation is frequently observed in bladder cancer correlating with cancer progression and invasion. However, the mechanisms underlying HGF/c-MET-mediated invasion in bladder cancer remains unknown. As part of a negative feedback loop SMAD7 binds to SMURF2 targeting the TGFβ receptor for degradation. Under these conditions, SMAD7 acts as a SMURF2 agonist by disrupting the intramolecular interactions within SMURF2. We demonstrate that HGF stimulates TGFβ signalling through c-SRC-mediated phosphorylation of SMURF2 resulting in loss of SMAD7 binding and enhanced SMURF2 C2-HECT interaction, inhibiting SMURF2 and enhancing TGFβ receptor stabilisation. This upregulation of the TGFβ pathway by HGF leads to TGFβ-mediated EMT and invasion. In vivo we show that TGFβ receptor inhibition prevents bladder cancer invasion. Furthermore, we make a rationale for the use of combinatorial TGFβ and MEK inhibitors for treatment of high-grade non-muscle-invasive bladder cancers

Absolute Single-Molecule Entropies from Quasi-Harmonic Analysis of Microsecond Molecular Dynamics: Correction Terms and Convergence Properties

The convergence properties of the absolute single-molecule configurational entropy and the correction terms used to estimate it are investigated using microsecond molecular dynamics simulation of a peptide test system and an improved methodology. The results are compared with previous applications for systems of diverse chemical nature. It is shown that (i) the effect of anharmonicity is small, (ii) the effect of pairwise correlation is typically large, and (iii) the latter affects to a larger extent the entropy estimate of thermodynamic states characterized by a higher motional correlation. The causes of such deviations from a quasi-harmonic behavior are explained. This improved approach provides entropies also for molecular systems undergoing conformational transitions and characterized by highly frustrated energy surfaces, thus not limited to systems sampling a single quasi-harmonic basin. Overall, this study emphasizes the need for extensive phase-space sampling in order to obtain a reliable estimation of entropic contributions

A Comparison of Magnesium and Beryllium Acceptors in GaN Grown by rf-Plasma Assisted Molecular Beam Epitaxy

ABSTRACT Step-doped structures of both magnesium and beryllium were grown in GaN and analyzed using secondary ion mass spectrometry. Dopant incorporation was studied as a function of substrate temperature and dopant flux for Ga-polarity and N-polarity GaN. Incorporation is different for each polarity, with Mg incorporating by up to a factor of 20 times more (30 times more with atomic hydrogen) on the Ga-face, while Be incorporates more readily on the N-face. The effect of atomic hydrogen on the incorporation kinetics of both Mg and Be is also discussed. Mg and Be both undergo surface segregation during growth. Photoluminescence measurements suggest that Be is a p-type dopant with an optical activation energy of approximately 100 meV

The Distinct Conformational Dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species

Disordered Structural Ensembles of Vasopressin and Oxytocin and Their Mutants

Vasopressin and oxytocin are intrinsically disordered cyclic nonapeptides belonging to a family of neurohypophysial hormones. Although unique in their functions, these peptides differ only by two residues and both feature a tocin ring formed by the disulfide bridge between first and sixth cysteine residues. This sequence and structural similarity are experimentally linked to oxytocin agonism at vasopressin receptors and vasopressin antagonism at oxytocin receptors. Yet single- or double-residue mutations in both peptides have been shown to have drastic impacts on their activities at either receptor, and possibly the ability to bind to their neurophysin carrier protein. In this study we perform molecular dynamics simulations of the unbound native and mutant sequences of the oxytocin and vasopressin hormones to characterize their structural ensembles. We classify the subpopulations of these structural ensembles on the basis of the distributions of radius of gyration and secondary structure and hydrogen-bonding features of the canonical tocin ring and disordered tail region. We then relate the structural changes observed in the unbound form of the different hormone sequences to experimental information about peptide receptor binding, and more indirectly, carrier protein binding affinity, receptor activity, and protease degradation. This study supports the hypothesis that the structural characteristics of the unbound form of an IDP can be used to predict structural or functional preferences of its functional bound form

Subthalamic nucleus stimulation affects orbitofrontal cortex in facial emotion recognition: a pet study

Deep brain stimulation (DBS) of the bilateral subthalamic nucleus (STN) in Parkinson's disease is thought to produce adverse events such as emotional disorders, and in a recent study, we found fear recognition to be impaired as a result. These changes have been attributed to disturbance of the STN's limbic territory and would appear to confirm that the negative emotion recognition network passes through the STN. In addition, it is now widely acknowledged that damage to the orbitofrontal cortex (OFC), especially the right side, can result in impaired recognition of facial emotions (RFE). In this context, we hypothesized that this reduced recognition of fear is correlated with modifications in the cerebral glucose metabolism of the right OFC. The objective of the present study was first, to reinforce our previous results by demonstrating reduced fear recognition in our Parkinson's disease patient group following STN DBS and, second, to correlate these emotional performances with glucose metabolism using 18FDG-PET. The 18FDG-PET and RFE tasks were both performed by a cohort of 13 Parkinson's disease patients 3 months before and 3 months after surgery for STN DBS. As predicted, we observed a significant reduction in fear recognition following surgery and obtained a positive correlation between these neuropsychological results and changes in glucose metabolism, especially in the right OFC. These results confirm the role of the STN as a key basal ganglia structure in limbic circuits

VX Hydrolysis by Human Serum Paraoxonase 1: A Comparison of Experimental and Computational Results

Human Serum paraoxonase 1 (HuPON1) is an enzyme that has been shown to hydrolyze a variety of chemicals including the nerve agent VX. While wildtype HuPON1 does not exhibit sufficient activity against VX to be used as an in vivo countermeasure, it has been suggested that increasing HuPON1's organophosphorous hydrolase activity by one or two orders of magnitude would make the enzyme suitable for this purpose. The binding interaction between HuPON1 and VX has recently been modeled, but the mechanism for VX hydrolysis is still unknown. In this study, we created a transition state model for VX hydrolysis (VXts) in water using quantum mechanical/molecular mechanical simulations, and docked the transition state model to 22 experimentally characterized HuPON1 variants using AutoDock Vina. The HuPON1-VXts complexes were grouped by reaction mechanism using a novel clustering procedure. The average Vina interaction energies for different clusters were compared to the experimentally determined activities of HuPON1 variants to determine which computational procedures best predict how well HuPON1 variants will hydrolyze VX. The analysis showed that only conformations which have the attacking hydroxyl group of VXts coordinated by the sidechain oxygen of D269 have a significant correlation with experimental results. The results from this study can be used for further characterization of how HuPON1 hydrolyzes VX and design of HuPON1 variants with increased activity against VX.United States. Defense Threat Reduction Agenc

Molecular determinants of binding to the Plasmodium subtilisin-like protease 1.

PfSUB1, a subtilisin-like protease of the human malaria parasite Plasmodium falciparum, is known to play important roles during the life cycle of the parasite and has emerged as a promising antimalarial drug target. In order to provide a detailed understanding of the origin of binding determinants of PfSUB1 substrates, we performed molecular dynamics simulations in combination with MM-GBSA free energy calculations using a homology model of PfSUB1 in complex with different substrate peptides. Key interactions, as well as residues that potentially make a major contribution to the binding free energy, are identified at the prime and nonprime side of the scissile bond and comprise peptide residues P4 to P2'. This finding stresses the requirement for peptide substrates to interact with both prime and nonprime side residues of the PfSUB1 binding site. Analyzing the energetic contributions of individual amino acids within the peptide-PfSUB1 complexes indicated that van der Waals interactions and the nonpolar part of solvation energy dictate the binding strength of the peptides and that the most favorable interactions are formed by peptide residues P4 and P1. Hot spot residues identified in PfSUB1 are dispersed over the entire binding site, but clustered areas of hot spots also exist and suggest that either the S4-S2 or the S1-S2' binding site should be exploited in efforts to design small molecule inhibitors. The results are discussed with respect to which binding determinants are specific to PfSUB1 and, therefore, might allow binding selectivity to be obtained