Dataset Information

A simple extension to the CMASA method for the prediction of catalytic residues in the presence of single point mutations.

ABSTRACT: The automatic identification of catalytic residues still remains an important challenge in structural bioinformatics. Sequence-based methods are good alternatives when the query shares a high percentage of identity with a well-annotated enzyme. However, when the homology is not apparent, which occurs with many structures from the structural genome initiative, structural information should be exploited. A local structural comparison is preferred to a global structural comparison when predicting functional residues. CMASA is a recently proposed method for predicting catalytic residues based on a local structure comparison. The method achieves high accuracy and a high value for the Matthews correlation coefficient. However, point substitutions or a lack of relevant data strongly affect the performance of the method. In the present study, we propose a simple extension to the CMASA method to overcome this difficulty. Extensive computational experiments are shown as proof of concept instances, as well as for a few real cases. The results show that the extension performs well when the catalytic site contains mutated residues or when some residues are missing. The proposed modification could correctly predict the catalytic residues of a mutant thymidylate synthase, 1EVF. It also successfully predicted the catalytic residues for 3HRC despite the lack of information for a relevant side chain atom in the PDB file.

SUBMITTER: Flores DI

PROVIDER: S-EPMC4182483 | biostudies-literature | 2014

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A simple extension to the CMASA method for the prediction of catalytic residues in the presence of single point mutations.

Flores David I DI Sotelo-Mundo Rogerio R RR Brizuela Carlos A CA

PloS one 20140930 9

The automatic identification of catalytic residues still remains an important challenge in structural bioinformatics. Sequence-based methods are good alternatives when the query shares a high percentage of identity with a well-annotated enzyme. However, when the homology is not apparent, which occurs with many structures from the structural genome initiative, structural information should be exploited. A local structural comparison is preferred to a global structural comparison when predicting f ...[more]

PMID: 25268770

Similar Datasets

Project description:BACKGROUND:Knowledge of catalytic residues can play an essential role in elucidating mechanistic details of an enzyme. However, experimental identification of catalytic residues is a tedious and time-consuming task, which can be expedited by computational predictions. Despite significant development in active-site prediction methods, one of the remaining issues is ranked positions of putative catalytic residues among all ranked residues. In order to improve ranking of catalytic residues and their prediction accuracy, we have developed a meta-approach based method CSmetaPred. In this approach, residues are ranked based on the mean of normalized residue scores derived from four well-known catalytic residue predictors. The mean residue score of CSmetaPred is combined with predicted pocket information to improve prediction performance in meta-predictor, CSmetaPred_poc. RESULTS:Both meta-predictors are evaluated on two comprehensive benchmark datasets and three legacy datasets using Receiver Operating Characteristic (ROC) and Precision Recall (PR) curves. The visual and quantitative analysis of ROC and PR curves shows that meta-predictors outperform their constituent methods and CSmetaPred_poc is the best of evaluated methods. For instance, on CSAMAC dataset CSmetaPred_poc (CSmetaPred) achieves highest Mean Average Specificity (MAS), a scalar measure for ROC curve, of 0.97 (0.96). Importantly, median predicted rank of catalytic residues is the lowest (best) for CSmetaPred_poc. Considering residues ranked ?20 classified as true positive in binary classification, CSmetaPred_poc achieves prediction accuracy of 0.94 on CSAMAC dataset. Moreover, on the same dataset CSmetaPred_poc predicts all catalytic residues within top 20 ranks for ~73% of enzymes. Furthermore, benchmarking of prediction on comparative modelled structures showed that models result in better prediction than only sequence based predictions. These analyses suggest that CSmetaPred_poc is able to rank putative catalytic residues at lower (better) ranked positions, which can facilitate and expedite their experimental characterization. CONCLUSIONS:The benchmarking studies showed that employing meta-approach in combining residue-level scores derived from well-known catalytic residue predictors can improve prediction accuracy as well as provide improved ranked positions of known catalytic residues. Hence, such predictions can assist experimentalist to prioritize residues for mutational studies in their efforts to characterize catalytic residues. Both meta-predictors are available as webserver at: http://14.139.227.206/csmetapred/ .

Project description:Worldwide structural genomics projects continue to release new protein structures at an unprecedented pace, so far nearly 6000, but only about 60% of these proteins have any sort of functional annotation.We explored a range of features that can be used for the prediction of functional residues given a known three-dimensional structure. These features include various centrality measures of nodes in graphs of interacting residues: closeness, betweenness and page-rank centrality. We also analyzed the distance of functional amino acids to the general center of mass (GCM) of the structure, relative solvent accessibility (RSA), and the use of relative entropy as a measure of sequence conservation. From the selected features, neural networks were trained to identify catalytic residues. We found that using distance to the GCM together with amino acid type provide a good discriminant function, when combined independently with sequence conservation. Using an independent test set of 29 annotated protein structures, the method returned 411 of the initial 9262 residues as the most likely to be involved in function. The output 411 residues contain 70 of the annotated 111 catalytic residues. This represents an approximately 14-fold enrichment of catalytic residues on the entire input set (corresponding to a sensitivity of 63% and a precision of 17%), a performance competitive with that of other state-of-the-art methods.We found that several of the graph based measures utilize the same underlying feature of protein structures, which can be simply and more effectively captured with the distance to GCM definition. This also has the added the advantage of simplicity and easy implementation. Meanwhile sequence conservation remains by far the most influential feature in identifying functional residues. We also found that due the rapid changes in size and composition of sequence databases, conservation calculations must be recalibrated for specific reference databases.

Project description:BACKGROUND: A basic question of protein structural studies is to which extent mutations affect the stability. This question may be addressed starting from sequence and/or from structure. In proteomics and genomics studies prediction of protein stability free energy change (DeltaDeltaG) upon single point mutation may also help the annotation process. The experimental DeltaDeltaG values are affected by uncertainty as measured by standard deviations. Most of the DeltaDeltaG values are nearly zero (about 32% of the DeltaDeltaG data set ranges from -0.5 to 0.5 kcal/mole) and both the value and sign of DeltaDeltaG may be either positive or negative for the same mutation blurring the relationship among mutations and expected DeltaDeltaG value. In order to overcome this problem we describe a new predictor that discriminates between 3 mutation classes: destabilizing mutations (DeltaDeltaG<-1.0 kcal/mol), stabilizing mutations (DeltaDeltaG>1.0 kcal/mole) and neutral mutations (-1.0</=DeltaDeltaG</=1.0 kcal/mole). RESULTS: In this paper a support vector machine starting from the protein sequence or structure discriminates between stabilizing, destabilizing and neutral mutations. We rank all the possible substitutions according to a three state classification system and show that the overall accuracy of our predictor is as high as 56% when performed starting from sequence information and 61% when the protein structure is available, with a mean value correlation coefficient of 0.27 and 0.35, respectively. These values are about 20 points per cent higher than those of a random predictor. CONCLUSIONS: Our method improves the quality of the prediction of the free energy change due to single point protein mutations by adopting a hypothesis of thermodynamic reversibility of the existing experimental data. By this we both recast the thermodynamic symmetry of the problem and balance the distribution of the available experimental measurements of free energy changes. This eliminates possible overestimations of the previously described methods trained on an unbalanced data set comprising a number of destabilizing mutations higher than stabilizing ones.

Project description:AimTo develop a fast, low-cost diagnostic strategy to identify single point mutations in highly variable genomes such as hepatitis C virus (HCV).MethodsIn patients with HCV infection, resistance-associated amino acid substitutions within the viral quasispecies prior to therapy can confer decreased susceptibility to direct-acting antiviral agents and lead to treatment failure and virological relapse. One such naturally occurring mutation is the Q80K substitution in the HCV-NS3 protease gene, which confers resistance to PI inhibitors, particularly simeprevir. Low-cost, highly sensitive techniques enabling routine detection of these single point mutations would be useful to identify patients at a risk of treatment failure. LightCycler methods, based on real-time PCR with sequence-specific probe hybridization, have been implemented in most diagnostic laboratories. However, this technique cannot identify single point mutations in highly variable genetic environments, such as the HCV genome. To circumvent this problem, we developed a new method to homogenize all nucleotides present in a region except the point mutation of interest.ResultsUsing nucleotide-specific probes Q, K, and R substitutions at position 80 were clearly identified at a sensitivity of 10% (mutations present at a frequency of at least 10% were detected). The technique was successfully applied to identify the Q80K substitution in 240 HCV G1 serum samples, with performance comparable to that of direct Sanger sequencing, the current standard procedure for this purpose. The new method was then validated in a Catalonian population of 202 HCV G1-infected individuals. Q80K was detected in 14.6% of G1a patients and 0% of G1b in our setting.ConclusionA fast, low-cost diagnostic strategy based on real-time PCR and fluorescence resonance energy transfer probe melting curve analysis has been successfully developed to identify single point mutations in highly variable genomes such as hepatitis C virus. This technique can be adapted to detect any single point mutation in highly variable genomes.

Dataset Information

A simple extension to the CMASA method for the prediction of catalytic residues in the presence of single point mutations.

Publications

A simple extension to the CMASA method for the prediction of catalytic residues in the presence of single point mutations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets