Computed protein-protein enthalpy signatures as a tool for identifying conformation sampling problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area

The KDEL trafficking receptor exploits pH to tune the strength of an unusual short hydrogen bond

The endoplasmic reticulum (ER) is the main site of protein synthesis in eukaryotic cells and requires a high concentration of luminal chaperones to function. During protein synthesis, ER luminal chaperones are swept along the secretory pathway and must be retrieved to maintain cell viability. ER protein retrieval is achieved by the KDEL receptor, which recognises a C-terminal Lys-Asp-Glu-Leu (KDEL) sequence. Recognition of ER proteins by the KDEL receptor is pH dependent, with binding occurring under acidic conditions in the Golgi and release under conditions of higher pH in the ER. Recent crystal structures of the KDEL receptor in the apo and peptide bound state suggested that peptide binding drives the formation of a short-hydrogen bond that locks the KDEL sequence in the receptor and activates the receptor for COPI binding in the cytoplasm. Using quantum mechanical calculations we demonstrate that the strength of this short hydrogen bond is reinforced following protonation of a nearby histidine, providing a conceptual link between receptor protonation and KDEL peptide binding. Protonation also controls the water networks adjacent to the peptide binding site, leading to a conformational change that ultimately allows the receptor-complex to be recognized by the COPI system

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

Limitations of non-polarizable force fields in describing anion binding poses in non-polar synthetic hosts

Transmembrane anion transport by synthetic ionophores has received increasing interest not only because of its relevance for understanding endogenous anion transport, but also because of potential implications for therapeutic routes in disease states where chloride transport is impaired. Computational studies can shed light on the binding recognition process and can deepen our mechanistic understanding of them. However, the ability of molecular mechanics methods to properly capture solvation and binding properties of anions is known to be challenging. Consequently, polarizable models have been suggested to improve the accuracy of such calculations. In this study, we calculate binding free energies for different anions to the synthetic ionophore, biotin[6]uril hexamethyl ester in acetonitrile and to biotin[6]uril hexaacid in water by employing non-polarizable and polarizable force fields. Anion binding shows strong solvent dependency consistent with experimental studies. In water, the binding strengths are iodide > bromide > chloride, and reversed in acetonitrile. These trends are well captured by both classes of force fields. However, the free energy profiles obtained from potential of mean force calculations and preferred binding positions of anions depend on the treatment of electrostatics. Results from simulations using the AMOEBA force-field, which recapitulate the observed binding positions, suggest strong effects from multipoles dominate with a smaller contribution from polarization. The oxidation status of the macrocycle was also found to influence anion recognition in water. Overall, these results have implications for the understanding of anion host interactions not just in synthetic ionophores, but also in narrow cavities of biological ion channels

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

Structural basis for antibiotic transport and inhibition in PepT2

The uptake and elimination of beta-lactam antibiotics in the human body are facilitated by the proton-coupled peptide transporters PepT1 (SLC15A1) and PepT2 (SLC15A2). The mechanism by which SLC15 family transporters recognize and discriminate between different drug classes and dietary peptides remains unclear, hampering efforts to improve antibiotic pharmacokinetics through targeted drug design and delivery. Here, we present cryo-EM structures of the proton-coupled peptide transporter, PepT2 from Rattus norvegicus, in complex with the widely used beta-lactam antibiotics cefadroxil, amoxicillin and cloxacillin. Our structures, combined with pharmacophore mapping, molecular dynamics simulations and biochemical assays, establish the mechanism of beta-lactam antibiotic recognition and the important role of protonation in drug binding and transport

Molecular basis for pH sensing in the KDEL trafficking receptor

Trafficking receptors control protein localization through the recognition of specific signal sequences that specify unique cellular locations. Differences in luminal pH are important for the vectorial trafficking of cargo receptors. The KDEL receptor is responsible for maintaining the integrity of the ER by retrieving luminally localized folding chaperones in a pH-dependent mechanism. Structural studies have revealed the end states of KDEL receptor activation and the mechanism of selective cargo binding. However, precisely how the KDEL receptor responds to changes in luminal pH remains unclear. To explain the mechanism of pH sensing, we combine analysis of X-ray crystal structures of the KDEL receptor at neutral and acidic pH with advanced computational methods and cell-based assays. We show a critical role for ordered water molecules that allows us to infer a direct connection between protonation in different cellular compartments and the consequent changes in the affinity of the receptor for cargo