skip to main content


Title: DRoP: Automated detection of conserved solvent‐binding sites on proteins
Abstract

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X‐ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein‐protein or protein‐ligand interactions, as well as conserved water‐binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.

 
more » « less
NSF-PAR ID:
10459766
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Proteins: Structure, Function, and Bioinformatics
Volume:
88
Issue:
1
ISSN:
0887-3585
Page Range / eLocation ID:
p. 152-165
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Protein‐protein interactions are either through direct contacts between two binding partners or mediated by structural waters. Both direct contacts and water‐mediated interactions are crucial to the formation of a protein‐protein complex. During the recent CAPRI rounds, a novel parallel searching strategy for predicting water‐mediated interactions is introduced into our protein‐protein docking method, MDockPP. Briefly, a FFT‐based docking algorithm is employed in generating putative binding modes, and an iteratively derived statistical potential‐based scoring function, ITScorePP, in conjunction with biological information is used to assess and rank the binding modes. Up to 10 binding modes are selected as the initial protein‐protein complex structures for MD simulations in explicit solvent. Water molecules near the interface are clustered based on the snapshots extracted from independent equilibrated trajectories. Then, protein‐ligand docking is employed for a parallel search for water molecules near the protein‐protein interface. The water molecules generated by ligand docking and the clustered water molecules generated by MD simulations are merged, referred to as the predicted structural water molecules. Here, we report the performance of this protocol for CAPRI rounds 28–29 and 31–35 containing 20 valid docking targets and 11 scoring targets. In the docking experiments, we predicted correct binding modes for nine targets, including one high‐accuracy, two medium‐accuracy, and six acceptable predictions. Regarding the two targets for the prediction of water‐mediated interactions, we achieved models ranked as “excellent” in accordance with the CAPRI evaluation criteria; one of these two targets is considered as a difficult target for structural water prediction. Proteins 2017; 85:424–434. © 2016 Wiley Periodicals, Inc.

     
    more » « less
  2. We propose a free energy calculation method for receptor–ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site. This modification leads to more accurate free energy calculation without performing separate free energy simulations for each binding pose. We applied our modification to simple model host–guest systems as a test, which have only two binding poses, by using a single decoupling method (SDM) in implicit solvent. The results showed that the binding free energies obtained from our method without knowing the two binding poses were in good agreement with the benchmark results obtained by explicit enumeration of the binding poses. Our method is applicable to other alchemical binding free energy calculation methods such as the double decoupling method (DDM) in explicit solvent. We performed a calculation for a protein–ligand system with explicit solvent using our modified thermodynamic path. The results of the free energy simulation along our modified path were in good agreement with the results of conventional DDM, which requires a separate binding free energy calculation for each of the binding poses of the example of phenol binding to T4 lysozyme in explicit solvent. © 2019 Wiley Periodicals, Inc.

     
    more » « less
  3. null (Ed.)
    Some methane-oxidizing bacteria use the ribosomally synthesized, posttranslationally modified natural product methanobactin (Mbn) to acquire copper for their primary metabolic enzyme, particulate methane monooxygenase. The operons encoding the machinery to biosynthesize and transport Mbns typically include genes for two proteins, MbnH and MbnP, which are also found as a pair in other genomic contexts related to copper homeostasis. While the MbnH protein, a member of the bacterial diheme cytochrome c peroxidase (bC c P)/MauG superfamily, has been characterized, the structure and function of MbnP, the relationship between the two proteins, and their role in copper homeostasis remain unclear. Biochemical characterization of MbnP from the methanotroph Methylosinus trichosporium OB3b now reveals that MbnP binds a single copper ion, present in the +1 oxidation state, with high affinity. Copper binding to MbnP in vivo is dependent on oxidation of the first tryptophan in a conserved WxW motif to a kynurenine, a transformation that occurs through an interaction of MbnH with MbnP. The 2.04-Å-resolution crystal structure of MbnP reveals a unique fold and an unusual copper-binding site involving a histidine, a methionine, a solvent ligand, and the kynurenine. Although the kynurenine residue may not serve as a Cu I primary-sphere ligand, being positioned ∼2.9 Å away from the Cu I ion, its presence is required for copper binding. Genomic neighborhood analysis indicates that MbnP proteins, and by extension kynurenine-containing copper sites, are widespread and may play diverse roles in microbial copper homeostasis. 
    more » « less
  4. Abstract

    An important question is how well the models submitted to CASP retain the properties of target structures. We investigate several properties related to binding. First we explore the binding of small molecules as probes, and count the number of interactions between each residue and such probes, resulting in a binding fingerprint. The similarity between two fingerprints, one for the X‐ray structure and the other for a model, is determined by calculating their correlation coefficient. The fingerprint similarity weakly correlates with global measures of accuracy, and GDT_TS higher than 80 is a necessary but not sufficient condition for the conservation of surface binding properties. The advantage of this approach is that it can be carried out without information on potential ligands and their binding sites. The latter information was available for a few targets, and we explored whether the CASP14 models can be used to predict binding sites and to dock small ligands. Finally, we tested the ability of models to reproduce protein–protein interactions by docking both the X‐ray structures and the models to their interaction partners in complexes. The analysis showed that in CASP14 the quality of individual domain models is approaching that offered by X‐ray crystallography, and hence such models can be successfully used for the identification of binding and regulatory sites, as well as for assembling obligatory protein–protein complexes. Success of ligand docking, however, often depends on fine details of the binding interface, and thus may require accounting for conformational changes by simulation methods.

     
    more » « less
  5. Predictions of the structures of stoichiometric, fractional, or nonstoichiometric hydrates of organic molecular crystals are immensely challenging due to the extensive search space of different water contents, host molecular placements throughout the crystal, and internal molecular conformations. However, the dry frameworks of these hydrates, especially for nonstoichiometric or isostructural dehydrates, can often be predicted from a standard anhydrous crystal structure prediction (CSP) protocol. Inspired by developments in the field of drug binding, we introduce an efficient data-driven and topologically aware approach for predicting organic molecular crystal hydrate structures through a mapping of water positions within the crystal structure. The method does not require a priori specification of water content and can, therefore, predict stoichiometric, fractional, and nonstoichiometric hydrate structures. This approach, which we term a mapping approach for crystal hydrates (MACH), establishes a set of rules for systematic determination of favorable positions for water insertion within predicted or experimental crystal structures based on considerations of the chemical features of local environments and void regions. The proposed approach is tested on hydrates of three pharmaceutically relevant compounds that exhibit diverse crystal packing motifs and void environments characteristic of hydrate structures. Overall, we show that our mapping approach introduces an advance in the efficient performance of hydrate CSP through generation of stable hydrate stoichiometries at low cost and should be considered an integral component for CSP workflows. 
    more » « less