Leonel F. Murga

Selective prediction of interaction sites in protein structures with THEMATICS.

BackgroundMethods are now available for the prediction of interaction sites in protein 3D structures. While many of these methods report high success rates for site prediction, often these predictions are not very selective and have low precision. Precision in site prediction is addressed using Theoretical Microscopic Titration Curves (THEMATICS), a simple computational method for the identification of active sites in enzymes. Recall and precision are measured and compared with other methods for the prediction of catalytic sites.ResultsUsing a test set of 169 enzymes from the original Catalytic Residue Dataset (CatRes) it is shown that THEMATICS can deliver precise, localised site predictions. Furthermore, adjustment of the cut-off criteria can improve the recall rates for catalytic residues with only a small sacrifice in precision. Recall rates for CatRes/CSA annotated catalytic residues are 41.1%, 50.4%, and 54.2% for Z score cut-off values of 1.00, 0.99, and 0.98, respectively. The corresponding precision rates are 19.4%, 17.9%, and 16.4%. The success rate for catalytic sites is higher, with correct or partially correct predictions for 77.5%, 85.8%, and 88.2% of the enzymes in the test set, corresponding to the same respective Z score cut-offs, if only the CatRes annotations are used as the reference set. Incorporation of additional literature annotations into the reference set gives total success rates of 89.9%, 92.9%, and 94.1%, again for corresponding cut-off values of 1.00, 0.99, and 0.98. False positive rates for a 75-protein test set are 1.95%, 2.60%, and 3.12% for Z score cut-offs of 1.00, 0.99, and 0.98, respectively.ConclusionWith a preferred cut-off value of 0.99, THEMATICS achieves a high success rate of interaction site prediction, about 86% correct or partially correct using CatRes/CSA annotations only and about 93% with an expanded reference set. Success rates for catalytic residue prediction are similar to those of other structure-based methods, but with substantially better precision and lower false positive rates. THEMATICS performs well across the spectrum of E.C. classes. The method requires only the structure of the query protein as input. THEMATICS predictions may be obtained via the web from structures in PDB format at: http://pfweb.chem.neu.edu/thematics/submit.html

Partial Order Optimum Likelihood (POOL): Maximum Likelihood Prediction of Protein Active Site Residues Using 3D Structure and Sequence Properties

A new monotonicity-constrained maximum likelihood approach, called Partial Order Optimum Likelihood (POOL), is presented and applied to the problem of functional site prediction in protein 3D structures, an important current challenge in genomics. The input consists of electrostatic and geometric properties derived from the 3D structure of the query protein alone. Sequence-based conservation information, where available, may also be incorporated. Electrostatics features from THEMATICS are combined with multidimensional isotonic regression to form maximum likelihood estimates of probabilities that specific residues belong to an active site. This allows likelihood ranking of all ionizable residues in a given protein based on THEMATICS features. The corresponding ROC curves and statistical significance tests demonstrate that this method outperforms prior THEMATICS-based methods, which in turn have been shown previously to outperform other 3D-structure-based methods for identifying active site residues. Then it is shown that the addition of one simple geometric property, the size rank of the cleft in which a given residue is contained, yields improved performance. Extension of the method to include predictions of non-ionizable residues is achieved through the introduction of environment variables. This extension results in even better performance than THEMATICS alone and constitutes to date the best functional site predictor based on 3D structure only, achieving nearly the same level of performance as methods that use both 3D structure and sequence alignment data. Finally, the method also easily incorporates such sequence alignment data, and when this information is included, the resulting method is shown to outperform the best current methods using any combination of sequence alignments and 3D structures. Included is an analysis demonstrating that when THEMATICS features, cleft size rank, and alignment-based conservation scores are used individually or in combination THEMATICS features represent the single most important component of such classifiers.

Statistical criteria for the identification of protein active sites using theoretical microscopic titration curves

Theoretical Microscopic Titration Curves (THEMATICS) may be used to identify chemically important residues in active sites of enzymes by characteristic deviations from the normal, sigmoidal Henderson–Hasselbalch titration behavior. Clusters of such deviant residues in physical proximity constitute reliable predictors of the location of the active site. Originally the residues with deviant predicted behavior were identified by human observation of the computed titration curves. However, it is preferable to select the unusual residues by mathematically well‐defined criteria, in order to reduce the chance of error, eliminate any possible biases, and substantially speed up the selection process. Here we present some simple statistical tests that constitute such selection criteria. The first derivatives of the predicted titration curves resemble distribution functions and are normalized. The moments of these first derivative functions are computed. It is shown that the third and fourth moments, measures of asymmetry and kurtosis, respectively, are good measures of the deviations from normal behavior. Results are presented for 44 different enzymes. Detailed results are given for 4 enzymes with 4 different types of chemistry: arginine kinase from Limulus polyphemus (horseshoe crab); β‐lactamase from Escherichia coli; glutamate racemase from Aquifex pyrophilus; and 3‐isopropylmalate dehydrogenase from Thiobacillus ferrooxidans. The relationship between the statistical measures of nonsigmoidal behavior in the predicted titration curves and the catalytic activity of the residue is discussed. Proteins 2005.

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.

Prediction of active sites for protein structures from computed chemical properties

MOTIVATION Identification of functional information for a protein from its three-dimensional (3D) structure is a major challenge in genomics. The power of theoretical microscopic titration curves (THEMATICS), when coupled with a statistical analysis, provides a method for high-throughput screening for identification of catalytic sites and binding sites with high accuracy and precision. The method requires only the 3D structure of the query protein as input, but it performs as well as other methods that depend on sequence alignments and structural similarities.

Prediction of interaction sites from apo 3D structures when the holo conformation is different

The predictability of catalytic and binding sites from apo structures is addressed for proteins that undergo significant conformational change upon binding. Theoretical microscopic titration curves (THEMATICS), an electrostatics‐based method for the prediction of functional sites, is performed on a test set of 24 proteins with both apo and holo structures available. For 23 of these 24 proteins (96%), THEMATICS predicts the correct catalytic or binding site for both the apo and holo forms. For only one of the 24 proteins, THEMATICS makes the correct prediction for the holo structure but fails for the apo structure. The metrics used by THEMATICS to identify functional residues generally are larger in absolute value for the functional residues in the holo forms compared to the corresponding residues in the apo forms. However, even in the apo forms, these identifying metrics are still statistically significantly larger for functional residues than for residues not involved in catalysis or binding. This indicates that some of the unusual electrostatic properties of functional residues are preserved in the apo conformation. Evidence is presented that certain residues immediately surrounding the active catalytic and binding residues impart functionally important chemical and electrostatic properties to the active residues. At least parts of these microenvironments exist in the unbound conformations, such that THEMATICS is able to distinguish the functional residues even in the apo structures. Proteins 2008.

Active site prediction for comparative model structures with thematics.

THEMATICS (Theoretical Microscopic Titration Curves) is a simple, reliable computational predictor of the active sites of enzymes from structure. Our method, based on well-established Finite Difference Poisson-Boltzmann techniques, identifies the ionisable residues with anomalous predicted titration behavior. A cluster of two or more such perturbed residues is a very reliable predictor of the active site. The protein does not have to bear any resemblance in sequence or structure to any previously characterized protein, but the method does require the three-dimensional structure. We now present evidence that THEMATICS can also locate the active site in structures built by comparative modeling from similar structures. Results are given for a total of 21 sets of proteins, including 21 templates and 83 comparative model structures. Detailed results are presented for three sets of orthologous proteins (Triosephosphate isomerase, 6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase, and Aspartate aminotransferase) and for one set of human homologues of Aldose reductase with different functions. THEMATICS correctly locates the active site in the model structures. This suggests that the method can be applicable to a much larger set of proteins for which an experimentally determined structure is unavailable. With a few exceptions, the predicted active sites in the comparative model structures are similar to that of the corresponding template structure.

Physicochemical Methods for Prediction of Functional Information for Proteins

Structural genomics initiatives are determining thousands of new protein structures. Many of these structures are of unknown function, and computational methods for the rapid determination of functional information from protein structure are needed. We present details of how functional information is obtained from the structure using THEMATICS (Theoretical Microscopic Titration Curves). THEMATICS is a computational procedure that gives information about chemical reactivity, based on solution of the Poisson-Boltzmann equations for the electrical potential function. We show how anomalies in predicted titration curves are established. We show further that when residues with anomalous predicted titration curves form a cluster in physical space, these residues tend to be very highly conserved across species and such clusters are reliable predictors of the active site. Results are given for ten enzymes; detailed results are shown for the enzymes triosephosphate isomerase (from chicken), 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (from E. coli), and papain (from papaya).

Theory of the Stark Effect in protein systems containing an electron donor–acceptor couple

Abstract Mixed-valency can occur in a variety of biological systems, such as the Cu(I)–Cu(II) pair in hemocyanin, Fe(II)–Fe(III) in many iron–oxo and iron–sulfur proteins, and Mn(II)–Mn(III) or Mn(III)–Mn(IV) in the photosynthetic water oxidase. The characterization of the ground states of such systems often has been controversial. Stark Effect spectroscopy is proving to be a valuable tool for the elucidation of systems of this type. The purpose of the present work is to develop theory for the spectral lineshape for the case where the electron donor and acceptor are coupled directly in a strong electric field. A mixed-valence dimer with an applied electric field aligned along the internuclear axis is studied using a two-site small-polaron model. Potential energy surfaces are calculated in the adiabatic (Born–Oppenheimer) approximation. It is shown that two nuclear coordinates (one totally symmetric and one antisymmetric) are coupled to the electronic motion, whereas only the antisymmetric coordinate is coupled in the absence of an electric field. For a strongly localized system, such as a protein system where electron donor and acceptor sites are separated by large distances, the potential surfaces become highly asymmetrical, but coupling to the totally symmetric mode is not significant. For a localized case corresponding to a valence-trapped two-metal cluster, the displacement along the totally symmetric coordinate is directly proportional to the applied field strength. Along the antisymmetric coordinate, the lowest potential surface is an asymmetric double well. For a delocalized (valence-averaged) two-metal cluster, there is significant displacement along the antisymmetric coordinate, an effect which also vanishes in the absence of an applied field. Contributions to the linewidth are estimated. Localized systems show larger field-induced shift in frequency maximum, whereas delocalized systems show greater field-induced line broadening.

A Hubbard model for the second hyperpolarizability in alternating polymers

Abstract The static second hyperpolarizability in low polymers of the type …ABAB… is studied using a two-band Hubbard Hamiltonian. For fully reduced bridged metallic polymers, there exists a non-zero value for the metal–ligand energy gap which corresponds to a γ of maximum magnitude. For mixed-valence bridged polymers, | γ | may be increased if the metal–ligand energy gap is increased. This means that | γ | may be enhanced significantly by synthetic alteration of the bridging ligand. Large, negative values for γ are predicted for polymers of the von Kameke type for the system in the mixed-valence state where the number of πd electrons equals twice the number of metal atoms minus one.