Association of putative concave protein-binding sites with the fluctuation behavior of residues

Here, we propose a binding site prediction method based on the high frequency end of the spectrum in the native state of the protein structural dynamics. The spectrum is obtained using an elastic network model (GNM). High frequency vibrating (HFV) residues are determined from the fastest modes dynamics. HFV residue clusters and the associated surface patch residues are tested for their likelihood to locate at the binding interfaces using two different data sets, the Benchmark Set of mainly enzymes and antigen/antibodies and the Cluster Set of more diverse structures. The binding interface is identified to be within 7.5 A of the HFV residue clusters in the Benchmark Set and Cluster Set, for 77% and 70% of the structures, respectively. The success rate increases to 88% and 84%, respectively, by using the surface patches. The results suggest that concave binding interfaces, typically those of enzyme-binding sites, are enriched by the HFV residues. Thus, we expect that the association of HFV residues with the interfaces is mostly for enzymes. If, however, a binding region has invaginations and cavities, as in some of the antigen/antibodies and in cases in the Cluster data set, we expect it would be detected there too. This implies that binding sites possess several (inter-related) properties such as cavities, high packing density, conservation, and disposition for hotspots at binding surfaces. It further suggests that the high frequency vibrating residue-based approach is a potential tool for identification of regions likely to serve as protein-binding sites. The software is available at http://www.prc.boun.edu.tr/PRC/software.html.     

The gaussian network model: precise prediction of residue fluctuations and application to binding problems

The single-parameter Gamma matrix of force constants proposed by the Gaussian Network Model (GNM) is iteratively modified to yield native state fluctuations that agree exactly with experimentally observed values. The resulting optimized Gamma matrix contains residue-specific force constants that may be used for an accurate analysis of ligand binding to single or multiple sites on proteins. Bovine Pancreatic Trypsin Inhibitor (BPTI) is used as an example. The calculated off-diagonal elements of the Gamma matrix, i.e., the optimized spring constants, obey a Lorentzian distribution. The mean value of the spring constants is approximately -0.1, a value much weaker than -1 of the GNM. Few of the spring constants are positive, indicating repulsion between residues. Residue pairs with large number of neighbors have spring constants around the mean, -0.1. Large negative spring constants are between highly correlated pairs of residues. The fluctuations of the distance between anticorrelated pairs of residues are subject to smaller spring constants. The importance of the number of neighbors of residue pairs in determining the elements of the Gamma matrix is pointed out. Allosteric effects of binding on a single or multiple residues of BPTI are illustrated and discussed. Comparison of the predictions of the present model with those of the standard GNM shows that the two models agree at lower modes, i.e., those relating to global motions, but they disagree at higher modes. In the higher modes, the present model points to the important contributions from specific residues whereas the standard GNM fails to do so.     

Basis for half-site ligand binding in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase is an allosterically regulated octameric enzyme composed of four heterodimers of a catalytic IDH2 subunit and a regulatory IDH1 subunit. Despite structural predictions that the enzyme would contain eight isocitrate binding sites, four NAD(+) binding sites, and four AMP binding sites, only half of the sites for each ligand can be measured in binding assays. On the basis of a potential interaction between side chains of Cys-150 residues in IDH2 subunits in each tetramer of the enzyme, ligand binding assays of wild-type (IDH1/IDH2) and IDH1/IDH2(C150S) octameric enzymes were conducted in the presence of dithiothreitol. These assays demonstrated the presence of eight isocitrate and four AMP binding sites for the wild-type enzyme in the presence of dithiothreitol and for the IDH1/IDH2(C150S) enzyme in the absence or presence of this reagent, suggesting that interactions between sulfhydryl side chains of IDH2 Cys-150 residues limit access to these sites. However, only two NAD(+) sites could be measured for either enzyme. A tetrameric form of IDH (an IDH1(G15D)/IDH2 mutant enzyme) demonstrated half-site binding for isocitrate (two sites) in the absence of dithiothreitol and full-site binding (four sites) in the presence of dithiothreitol. Only one NAD(+) site could be measured for the tetramer under both conditions. In the context of the structure of the enzyme, these results suggest that an observed asymmetry between heterotetramers in the holoenzyme contributes to interactions between IDH2 Cys-150 residues and to half-site binding of isocitrate, but that a form of negative cooperativity may limit access to apparently equivalent NAD(+) binding sites.     

A matching algorithm for catalytic residue site selection in computational enzyme design

A loop closure-based sequential algorithm, PRODA_MATCH, was developed to match catalytic residues onto a scaffold for enzyme design in silico. The computational complexity of this algorithm is polynomial with respect to the number of active sites, the number of catalytic residues, and the maximal iteration number of cyclic coordinate descent steps. This matching algorithm is independent of a rotamer library that enables the catalytic residue to take any required conformation during the reaction coordinate. The catalytic geometric parameters defined between functional groups of transition state (TS) and the catalytic residues are continuously optimized to identify the accurate position of the TS. Pseudo-spheres are introduced for surrounding residues, which make the algorithm take binding into account as early as during the matching process. Recapitulation of native catalytic residue sites was used as a benchmark to evaluate the novel algorithm. The calculation results for the test set show that the native catalytic residue sites were successfully identified and ranked within the top 10 designs for 7 of the 10 chemical reactions. This indicates that the matching algorithm has the potential to be used for designing industrial enzymes for desired reactions.     

Residue centrality, functionally important residues, and active site shape: analysis of enzyme and non-enzyme families

The representation of protein structures as small-world networks facilitates the search for topological determinants, which may relate to functionally important residues. Here, we aimed to investigate the performance of residue centrality, viewed as a family fold characteristic, in identifying functionally important residues in protein families. Our study is based on 46 families, including 29 enzyme and 17 non-enzyme families. A total of 80% of these central positions corresponded to active site residues or residues in direct contact with these sites. For enzyme families, this percentage increased to 91%, while for non-enzyme families the percentage decreased substantially to 48%. A total of 70% of these central positions are located in catalytic sites in the enzyme families, 64% are in hetero-atom binding sites in those families binding hetero-atoms, and only 16% belong to protein-protein interfaces in families with protein-protein interaction data. These differences reflect the active site shape: enzyme active sites locate in surface clefts, hetero-atom binding residues are in deep cavities, while protein-protein interactions involve a more planar configuration. On the other hand, not all surface cavities or clefts are comprised of central residues. Thus, closeness centrality identifies functionally important residues in enzymes. While here we focus on binding sites, we expect to identify key residues for the integration and transmission of the information to the rest of the protein, reflecting the relationship between fold and function. Residue centrality is more conserved than the protein sequence, emphasizing the robustness of protein structures.     

A Structure-Based Approach for Detection of Thiol Oxidoreductases and Their Catalytic Redox-Active Cysteine Residues

Cysteine (Cys) residues often play critical roles in proteins, for example, in the formation of structural disulfide bonds, metal binding, targeting proteins to the membranes, and various catalytic functions. However, the structural determinants for various Cys functions are not clear. Thiol oxidoreductases, which are enzymes containing catalytic redox-active Cys residues, have been extensively studied, but even for these proteins there is little understanding of what distinguishes their catalytic redox Cys from other Cys functions. Herein, we characterized thiol oxidoreductases at a structural level and developed an algorithm that can recognize these enzymes by (i) analyzing amino acid and secondary structure composition of the active site and its similarity to known active sites containing redox Cys and (ii) calculating accessibility, active site location, and reactivity of Cys. For proteins with known or modeled structures, this method can identify proteins with catalytic Cys residues and distinguish thiol oxidoreductases from the enzymes containing other catalytic Cys types. Furthermore, by applying this procedure to Saccharomyces cerevisiae proteins containing conserved Cys, we could identify the majority of known yeast thiol oxidoreductases. This study provides insights into the structural properties of catalytic redox-active Cys and should further help to recognize thiol oxidoreductases in protein sequence and structure databases.     

Protein kinase C isozymes and their selectivity towards ruboxistaurin

Protein kinase C (PKC) isozymes are an important class of enzymes in cell signaling and as drug targets. They are involved in specific pathways and have selectivity towards certain ligands, despite their high sequence similarities. Ruboxistaurin is a specific inhibitor of PKC-beta. To understand the molecular determinants for the selectivity of ruboxistaurin, we derived the three-dimensional structures of the kinase domains of PKC-alpha, -betaI, and -zeta using homology modeling. Several binding orientations of ruboxistaurin in the binding sites of these PKC catalytic domains were analyzed, and a putative alternative binding site for PKC-zeta was identified in its kinase domain. The calculated free energy of binding correlates well with the IC(50) of the inhibitor against each PKC isozyme. A residue-based energy decomposition analysis attributed the binding free energy to several key residues in the catalytic sites of these enzymes, revealing potential protein-ligand interactions responsible for ligand binding. The contiguous binding site revealed in the catalytic domain of PKC-zeta provides avenues for selective drug design. The details of structural nuances for specific inhibition of PKC isozymes are presented in the context of the three-dimensional structures of this important class of enzymes.     

Models with energy penalty on interresidue rotation address insufficiencies of conventional elastic network models

In this study, I present a new elastic network model, to our knowledge, that addresses insufficiencies of two conventional models—the Gaussian network model (GNM) and the anisotropic network model (ANM). It has been shown previously that the GNM is not rotation-invariant due to its energy, which penalizes rigid-body rotation (external rotation). As a result, GNM models are found contaminated with rigid-body rotation, especially in the most collective ones. A new model (EPIRM) is proposed to remove such external component in modes. The extracted internal motions result from a potential that penalizes interresidue stretching and rotation in a protein. The new model is shown to pertinently describe crystallographic temperature factors (B-factors) and protein open↔closed transitions. Also, the capability of separating internal and external motions in GNM slow modes permits reexamining important mechanochemical properties in enzyme active sites. The results suggest that catalytic residues stay closer to rigid-body rotation axes than their immediate backbone neighbors. I show that the cumulative density of states for EPIRM and ANM follow different power laws as functions of low-mode frequencies. When using a cutoff distance of 7.5 Å, The cumulative density of states of EPIRM scales faster than that of all-atom normal mode analysis and slower than that of simple lattices.     

A template search reveals mechanistic similarities and differences in beta-ketoacyl synthases (KAS) and related enzymes

A detailed comparison of the active sites in beta-ketoacyl synthases (KAS) and related enzymes has been made. Using three-dimensional templates of the three catalytic residues to scan the protein structural database reveals differences in both the geometry and the catalytic role of equivalent residues in different members of the family. The template based on the catalytic cysteine and two histidines in the KAS I and II is totally specific for this family, with no false hits. However, the role of the histidines in catalysis is different between KAS I/II and thiolase on the one hand and KAS III/chalcone synthase on the other. In contrast, a template comprising only cysteine and one histidine is not specific with many hits including members of the KAS family, metal binding sites, other active sites in nonhomologous proteins, and some "random" nonactive sites.     

Non-catalytic Binding Sites Induce Weaker Long-Range Evolutionary Rate Gradients than Catalytic Sites in Enzymes

Enzymes exhibit a strong long-range evolutionary constraint that extends from their catalytic site and affects even distant sites, where site-specific evolutionary rate increases monotonically with distance. While protein–protein sites in enzymes were previously shown to induce only a weak conservation gradient, a comprehensive relationship between different types of functional sites in proteins and the magnitude of evolutionary rate gradients they induce has yet to be established. Here, we systematically calculate the evolutionary rate (dN/dS) of sites as a function of distance from different types of binding sites in enzymes and other proteins: catalytic sites, non-catalytic ligand binding sites, allosteric binding sites, and protein–protein interaction sites. We show that catalytic sites indeed induce significantly stronger evolutionary rate gradient than all other types of non-catalytic binding sites. In addition, catalytic sites in enzymes with no known allosteric function still induce strong long-range conservation gradients. Notably, the weak long-range conservation gradients induced by non-catalytic binding sites in enzymes is nearly identical in magnitude to those induced by ligand binding sites in non-enzymes. Finally, we show that structural determinants such as local solvent exposure of sites cannot explain the observed difference between catalytic and non-catalytic functional sites. Our results suggest that enzymes and non-enzymes share similar evolutionary constraints only when examined from the perspective of non-catalytic functional sites. Hence, the unique evolutionary rate gradient from catalytic sites in enzymes is likely driven by the optimization of catalysis rather than ligand binding and allosteric functions.     

Looking at enzymes from the inside out: the proximity of catalytic residues to the molecular centroid can be used for detection of active sites and enzyme-ligand interfaces

Analysis of the distances of the exposed residues in 175 enzymes from the centroids of the molecules indicates that catalytic residues are very often found among the 5% of residues closest to the enzyme centroid. This property of catalytic residues is implemented in a new prediction algorithm (named EnSite) for locating the active sites of enzymes and in a new scheme for re-ranking enzyme-ligand docking solutions. EnSite examines only 5% of the molecular surface (represented by surface dots) that is closest to the centroid, identifying continuous surface segments and ranking them by their area size. EnSite ranks the correct prediction 1-4 in 97% of the cases in a dataset of 65 monomeric enzymes (rank 1 for 89% of the cases) and in 86% of the cases in a dataset of 176 monomeric and multimeric enzymes from all six top-level enzyme classifications (rank 1 in 74% of the cases). Importantly, identification of buried or flat active sites is straightforward because EnSite "looks" at the molecular surface from the inside out. Detailed examination of the results indicates that the proximity of the catalytic residues to the centroid is a property of the functional unit, defined as the assembly of domains or chains that form the active site (in most cases the functional unit corresponds to a single whole polypeptide chain). Using the functional unit in the prediction further improves the results. The new property of active sites is also used for re-evaluating enzyme-inhibitor unbound docking results. Sorting the docking solutions by the distance of the interface to the centroid of the enzyme improves remarkably the ranks of nearly correct solutions compared to ranks based on geometric-electrostatic-hydrophobic complementarity scores.     

Enhanced performance in prediction of protein active sites with THEMATICS and support vector machines

Theoretical microscopic titration curves (THEMATICS) is a computational method for the identification of active sites in proteins through deviations in computed titration behavior of ionizable residues. While the sensitivity to catalytic sites is high, the previously reported sensitivity to catalytic residues was not as high, about 50%. Here THEMATICS is combined with support vector machines (SVM) to improve sensitivity for catalytic residue prediction from protein 3D structure alone. For a test set of 64 proteins taken from the Catalytic Site Atlas (CSA), the average recall rate for annotated catalytic residues is 61%; good precision is maintained selecting only 4% of all residues. The average false positive rate, using the CSA annotations is only 3.2%, far lower than other 3D-structure-based methods. THEMATICS-SVM returns higher precision, lower false positive rate, and better overall performance, compared with other 3D-structure-based methods. Comparison is also made with the latest machine learning methods that are based on both sequence alignments and 3D structures. For annotated sets of well-characterized enzymes, THEMATICS-SVM performance compares very favorably with methods that utilize sequence homology. However, since THEMATICS depends only on the 3D structure of the query protein, no decline in performance is expected when applied to novel folds, proteins with few sequence homologues, or even orphan sequences. An extension of the method to predict non-ionizable catalytic residues is also presented. THEMATICS-SVM predicts a local network of ionizable residues with strong interactions between protonation events; this appears to be a special feature of enzyme active sites.     

Structural and functional characterization of binding sites in metallocarboxypeptidases based on Optimal Docking Area analysis

The metallocarboxypeptidases (MCPs) belonging to the clan MC were studied by the Optimal Docking Area (ODA) method to evaluate protein-protein binding sites and to provide a basis for the identification of binding partners for this class of enzymes. The ODA method identifies surface patches with optimal desolvation energy based on the selection of low-energy docking regions, generated from a set of surface points around the protein. With few exceptions, the ODA method identified surface patches with a significant low-energy docking surface for all the MCPs with known three-dimensional structure. Overall, in 14 out of 24 cases, the detected ODA patches were correctly located (i.e. more than 50% of the predicted residues were in known protein-protein binding sites), yielding a global success rate of 58%. More specifically, the success rate increased up to 80% on the ODA patches detected for the catalytic domains of the M14A subfamily, independently on the partner. Interestingly, the ODA residues on the catalytic domain were correctly located in the interface with the N-terminal pro domain in all MCPs. The spatial distribution of the ODA patches for the different members of the family is in relation to the origin and function of the particular MCP, which allowed distinguishing between them. In good agreement with the experimentally characterized protein interfaces, the total average surface area of the theoretically derived ODA patches for the catalytic domain of MCPs is around 1700 A2 and their content in hydrophobic residues is about 40%. As a particular case, the average surface area of the ODA patches in MCPs of crop insect pests is about twice that of the MCPs of vertebrates, which might be related to their particular function. We recognized two binding regions for the catalytic domain of the MCPs, one of them accounting for nearly all the known intermolecular interactions made up by the enzymes. Protein inhibitors seem to have evolved to dock on this subset of ODA patches, evoking the binding mode of the N-terminal pro domains. The second binding region detected, for which no ligands have been identified so far, seems to be related to the acquisition/maintenance of the native structure of the peptidase. Overall, the ODA method has been successful in identifying low-energy docking areas in a set of structurally and functionally related proteins, suggesting that it can be easily extended to other families in the search for protein-protein binding sites and for their functional significance.     

Insertion of new sequences into the catalytic domain of an enzyme

Activities of enzymes can be modified by the replacement of active-site amino acids with residues that strengthen specific interactions with substrates or that alter the specificity. The scope for engineered enzymes would be broadened if additional, new sequences could be inserted into a catalytic domain. Properly designed, these sequences could encode new ligand binding sites, be intermediates in the construction of chimeric enzymes, or alter the internal flexibility and "breathing" modes of the active-site region. As a first step toward this objective, we inserted oligopeptides of up to 14 amino acids into various locations within an 82 amino acid region of the adenylate synthesis domain of Escherichia coli methionyl-tRNA synthetase. These sites include ones that are flanked by sequences that are conserved between the proteins from E. coli and the yeast Saccharomyces cerevisiae and those that are essential for activity and stability. We found that all of the insertional mutants are stable and some have catalytic parameters for adenylate synthesis that are comparable to those of the wild-type enzyme. Thus, such an approach may provide for a variety of novel applications.     

Analysis of glycoside hydrolase family 98: catalytic machinery, mechanism and a novel putative carbohydrate binding module

Glycoside hydrolases (GHs) are diverse enzymes of biotechnological and medical importance. Bioinformatics contributes to our understanding of GH structure and function in various ways, including dissection of their typically modular structures and detection of the distant evolutionary relationships between families that often allow for prediction of catalytic sites. Here these twin strands are applied to the recently described GH98 family, the founder member of which is a blood group glycotope-cleaving endo-beta-galactosidase of potential medical importance from Clostridium perfringens. Three domains can be discerned including a central catalytic TIM barrel domain in which putative catalytic residues can be assigned. Distant homologies and domain contexts suggest that the N-terminal domain is a novel carbohydrate binding module.     

The acceptor specificity of UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferases

The in vitro and in vivo specificity of the family of peptide:N-acetylgalactosaminyltransferases (GalNAcT) is analyzed on the basis of the reactivity and/or inhibitory activity of peptides and protein segments. The transferases appear to be multi-substrate enzymes with extended active sites containing a least nine subsites that interact cooperatively with a linear segment of at least nine amino acid residues on the acceptor polypeptide. Functional acceptor sites are located on the surface of the protein and extended conformations (-strand conformation) are preferred. The acceptor specificity of GalNAc-T can be predicted from the primary structure of the acceptor peptide with an accuracy of 70 to 80%. The same GalNAc-T enzymes catalyze the glycosylation of both serine and threonine residues. The higher in vitro catalytic efficiency toward threonine versus serine is the result of enhanced binding as well as increased reaction velocity, both effects being the result of steric interactions between the active site of the enzyme and the methyl group of threonine. Results from substrate binding studies suggest that GalNAc-T catalyzed transfer proceeds via an ordered sequential mechanism.     

Binding of the transition state analog MgADP-fluoroaluminate to F1-ATPase

Escherichia coli F1-ATPase from mutant betaY331W was potently inhibited by fluoroaluminate plus MgADP but not by MgADP alone. beta-Trp-331 fluorescence was used to measure MgADP binding to catalytic sites. Fluoroaluminate induced a very large increase in MgADP binding affinity at catalytic site one, a smaller increase at site two, and no effect at site three. Mutation of either of the critical catalytic site residues beta-Lys-155 or beta-Glu-181 to Gln abolished the effects of fluoroaluminate on MgADP binding. The results indicate that the MgADP-fluoroaluminate complex is a transition state analog and independently demonstrate that residues beta-Lys-155 and (particularly) beta-Glu-181 are important for generation and stabilization of the catalytic transition state. Dicyclohexylcarbodiimide-inhibited enzyme, with 1% residual steady-state ATPase, showed normal transition state formation as judged by fluoroaluminate-induced MgADP binding affinity changes, consistent with a proposed mechanism by which dicyclohexylcarbodiimide prevents a conformational interaction between catalytic sites but does not affect the catalytic step per se. The fluorescence technique should prove valuable for future transition state studies of F1-ATPase.     

Analysis of the active center of branching enzyme II from maize endosperm

Analysis of the primary structure of mBEII, with those of other branching and amylolytic enzymes as reference, identifies four highly conserved regions which may be involved in substrate binding and in catalysis. When one of the amino acid residues corresponding to the putative catalytic sites of mBEII, i.e., Asp-386, Glu-441, and Asp-509, was replaced, activity disappeared. These putative catalytic residues are located in three different regions (regions 2–4) of the four highly conserved regions (regions 1–4) which exist in the primary structure of most starch hydrolases and related enzymes, including branching enzymes. Region 3, which contains Glu-441 as one of the putative catalytic residues, was located downstream of the carboxyl-terminal position previously reported. The importance of the carboxyl amino acid residues was also demonstrated by chemical modification of the branching enzyme protein using 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide.     

Analysis of catalytic residues of Thermoactinomyces vulgaris R-47 alpha-amylase II (TVA II) by site-directed mutagenesis

To confirm that the catalytic residues (Asp325, Glu354, and Asp421) are necessary for the hydrolysis of starch, pullulan, and cyclodextrins, we constructed TVA II mutated by site-directed mutagenesis. The mutated enzymes (D325N, E354Q, and D421N) had markedly reduced levels of activity, less than 0.006% of the wild type, indicating that these three residues are the catalytic sites for these substrates. Even E354D had reduced levels of activity, less than 0.05% of wild type. These four mutated enzymes retained a trace of activity. From the result of hydrolysis patterns for maltohexaose, in particular, D421N, unlike D325N and E354Q, catalyzed transglycosylation rather than hydrolysis. The results suggest that Asp421 could function to capture water molecules.     

Categoric prediction of metal ion mechanisms in the active sites of 17 select type II restriction endonucleases

Recently determined crystal structures of type II restriction endonucleases have produced a plethora of information on the basis for target site sequence selectivity. The positioning and role of metal ions in DNA recognition sites might reflect important properties of protein-DNA interaction. Although acidic and basic groups in the active sites can be identified, and in some cases divalent-metal binding sites delineated, a convincing picture clarifying the way in which the attacking hydroxide ion is generated, and the leaving group stabilized, has not been elucidated for any of the enzymes. We have examined the interatomic distances between metal ions and proposed key catalytic residues in the binding sites of seventeen type II restriction endonucleases whose crystal structures are documented in literature. The summary and critical evaluation of structural assignments and predictions made earlier have been useful to group these enzymes. All the enzymes used for this study have been categorized on the basis of the number of metal ions identified in their crystal structures. Among 17 experimentally characterized (not putative) type II REases, whose apparently full-length sequences are available in REBASE, we predict 8 (47%) to follow the single metal ion mechanism, 5 to follow the two metal ion mechanism, 2, the three metal ion mechanism, 1, the four metal ion mechanism and 1 the six metal ion mechanism.