Waterdock 2.0: Water placement prediction for Holo-structures with a pymol plugin.

Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.GAR is supported by the Memorial Sloan Kettering Cancer Center, NIH grant P30 CA008748

Waterdock 2.0: Water placement prediction for Holo-structures with a pymol plugin.

Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.GAR is supported by the Memorial Sloan Kettering Cancer Center, NIH grant P30 CA008748

Dynamics based alignment of proteins: an alternative approach to quantify dynamic similarity

<p>Abstract</p> <p>Background</p> <p>The dynamic motions of many proteins are central to their function. It therefore follows that the dynamic requirements of a protein are evolutionary constrained. In order to assess and quantify this, one needs to compare the dynamic motions of different proteins. Comparing the dynamics of distinct proteins may also provide insight into how protein motions are modified by variations in sequence and, consequently, by structure. The optimal way of comparing complex molecular motions is, however, far from trivial. The majority of comparative molecular dynamics studies performed to date relied upon prior sequence or structural alignment to define which residues were equivalent in 3-dimensional space.</p> <p>Results</p> <p>Here we discuss an alternative methodology for comparative molecular dynamics that does not require any prior alignment information. We show it is possible to align proteins based solely on their dynamics and that we can use these dynamics-based alignments to quantify the dynamic similarity of proteins. Our method was tested on 10 representative members of the PDZ domain family.</p> <p>Conclusions</p> <p>As a result of creating pair-wise dynamics-based alignments of PDZ domains, we have found evolutionarily conserved patterns in their backbone dynamics. The dynamic similarity of PDZ domains is highly correlated with their structural similarity as calculated with Dali. However, significant differences in their dynamics can be detected indicating that sequence has a more refined role to play in protein dynamics than just dictating the overall fold. We suggest that the method should be generally applicable.</p

On the atomic structure of cocaine in solution

Cocaine is an amphiphilic drug which has the ability to cross the blood–brain barrier (BBB). Here, a combination of neutron diffraction and computation has been used to investigate the atomic scale structure of cocaine in aqueous solutions. Both the observed conformation and hydration of cocaine appear to contribute to its ability to cross hydrophobic layers afforded by the BBB, as the average conformation yields a structure which might allow cocaine to shield its hydrophilic regions from a lipophilic environment. Specifically, the carbonyl oxygens and amine group on cocaine, on average, form ~5 bonds with the water molecules in the surrounding solvent, and the top 30% of water molecules within 4 Å of cocaine are localized in the cavity formed by an internal hydrogen bond within the cocaine molecule. This water mediated internal hydrogen bonding suggests a mechanism of interaction between cocaine and the BBB that negates the need for deprotonation prior to interaction with the lipophilic portions of this barrier. This finding also has important implications for understanding how neurologically active molecules are able to interact with both the blood stream and BBB and emphasizes the use of structural measurements in solution in order to understand important biological function.Peer ReviewedPostprint (author's final draft

A7DB: a relational database for mutational, physiological and pharmacological data related to the α7 nicotinic acetylcholine receptor

BACKGROUND: Nicotinic acetylcholine receptors (nAChRs) are pentameric proteins that are important drug targets for a variety of diseases including Alzheimer's, schizophrenia and various forms of epilepsy. One of the most intensively studied nAChR subunits in recent years has been α7. This subunit can form functional homomeric pentamers (α7)(5), which can make interpretation of physiological and structural data much simpler. The growing amount of structural, pharmacological and physiological data for these receptors indicates the need for a dedicated and accurate database to provide a means to access this information in a coherent manner. DESCRIPTION: A7DB is a new relational database of manually curated experimental physiological data associated with the α7 nAChR. It aims to store as much of the pharmacology, physiology and structural data pertaining to the α7 nAChR. The data is accessed via web interface that allows a user to search the data in multiple ways: 1) a simple text query 2) an incremental query builder 3) an interactive query builder and 4) a file-based uploadable query. It currently holds more than 460 separately reported experiments on over 85 mutations. CONCLUSIONS: A7DB will be a useful tool to molecular biologists and bioinformaticians not only working on the α7 receptor family of proteins but also in the more general context of nicotinic receptor modelling. Furthermore it sets a precedent for expansion with the inclusion of all nicotinic receptor families and eventually all cys-loop receptor families

The Startle Disease Mutation E103K Impairs Activation of Human Homomeric α1 Glycine Receptors by Disrupting an Intersubunit Salt Bridge across the Agonist Binding Site

Glycine receptors (GlyR) belong to the pentameric ligand-gated ion channel (pLGIC) superfamily and mediate fast inhibitory transmission in the vertebrate CNS. Disruption of glycinergic transmission by inherited mutations produces startle disease in man. Many startle mutations are in GlyRs and provide useful clues to the function of the channel domains. E103K is one of few startle mutations found in the extracellular agonist binding site of the channel, in loop A of the principal side of the subunit interface. Homology modeling shows that the side chain of Glu-103 is close to that of Arg-131, in loop E of the complementary side of the binding site, and may form a salt bridge at the back of the binding site, constraining its size. We investigated this hypothesis in recombinant human α1 GlyR by site-directed mutagenesis and functional measurements of agonist efficacy and potency by whole cell patch clamp and single channel recording. Despite its position near the binding site, E103K causes hyperekplexia by impairing the efficacy of glycine, its ability to gate the channel once bound, which is very high in wild type GlyR. Mutating Glu-103 and Arg-131 caused various degrees of loss-of-function in the action of glycine, whereas mutations in Arg-131 enhanced the efficacy of the slightly bigger partial agonist sarcosine (N-methylglycine). The effects of the single charge-swapping mutations of these two residues were largely rescued in the double mutant, supporting the possibility that they interact via a salt bridge that normally constrains the efficacy of larger agonist molecules

Oxidation Resistance 1 Modulates Glycolytic Pathways in the Cerebellum via an Interaction with Glucose-6-Phosphate Isomerase

Glucose metabolism is essential for the brain: it not only provides the required energy for cellular function and communication but also participates in balancing the levels of oxidative stress in neurons. Defects in glucose metabolism have been described in neurodegenerative disease; however, it remains unclear how this fundamental process contributes to neuronal cell death in these disorders. Here, we investigated the molecular mechanisms driving the selective neurodegeneration in an ataxic mouse model lacking oxidation resistance 1 (Oxr1) and discovered an unexpected function for this protein as a regulator of the glycolytic enzyme, glucose-6-phosphate isomerase (GPI/Gpi1). Initially, we present a dysregulation of metabolites of glucose metabolism at the pre-symptomatic stage in the Oxr1 knockout cerebellum. We then demonstrate that Oxr1 and Gpi1 physically and functionally interact and that the level of Gpi1 oligomerisation is disrupted when Oxr1 is deleted in vivo. Furthermore, we show that Oxr1 modulates the additional and less well-understood roles of Gpi1 as a cytokine and neuroprotective factor. Overall, our data identify a new molecular function for Oxr1, establishing this protein as important player in neuronal survival, regulating both oxidative stress and glucose metabolism in the brain

JGromacs: A Java Package for Analyzing Protein Simulations

UNLABELLED: In this paper, we introduce JGromacs, a Java API (Application Programming Interface) that facilitates the development of cross-platform data analysis applications for Molecular Dynamics (MD) simulations. The API supports parsing and writing file formats applied by GROMACS (GROningen MAchine for Chemical Simulations), one of the most widely used MD simulation packages. JGromacs builds on the strengths of object-oriented programming in Java by providing a multilevel object-oriented representation of simulation data to integrate and interconvert sequence, structure, and dynamics information. The easy-to-learn, easy-to-use, and easy-to-extend framework is intended to simplify and accelerate the implementation and development of complex data analysis algorithms. Furthermore, a basic analysis toolkit is included in the package. The programmer is also provided with simple tools (e.g., XML-based configuration) to create applications with a user interface resembling the command-line interface of GROMACS applications. AVAILABILITY: JGromacs and detailed documentation is freely available from http://sbcb.bioch.ox.ac.uk/jgromacs under a GPLv3 license

Quantifying Water-Mediated Protein–Ligand Interactions in a Glutamate Receptor: A DFT Study

It is becoming increasingly clear that careful treatment of water molecules in ligand–protein interactions is required in many cases if the correct binding pose is to be identified in molecular docking. Water can form complex bridging networks and can play a critical role in dictating the binding mode of ligands. A particularly striking example of this can be found in the ionotropic glutamate receptors. Despite possessing similar chemical moieties, crystal structures of glutamate and α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) in complex with the ligand-binding core of the GluA2 ionotropic glutamate receptor revealed, contrary to all expectation, two distinct modes of binding. The difference appears to be related to the position of water molecules within the binding pocket. However, it is unclear exactly what governs the preference for water molecules to occupy a particular site in any one binding mode. In this work we use density functional theory (DFT) calculations to investigate the interaction energies and polarization effects of the various components of the binding pocket. Our results show (i) the energetics of a key water molecule are more favorable for the site found in the glutamate-bound mode compared to the alternative site observed in the AMPA-bound mode, (ii) polarization effects are important for glutamate but less so for AMPA, (iii) ligand–system interaction energies alone can predict the correct binding mode for glutamate, but for AMPA alternative modes of binding have similar interaction energies, and (iv) the internal energy is a significant factor for AMPA but not for glutamate. We discuss the results within the broader context of rational drug-design