Quantitative Comparison of Catalytic Mechanisms and Overall Reactions in Convergently Evolved Enzymes: Implications for Classification of Enzyme Function

Functionally analogous enzymes are those that catalyze similar reactions on similar substrates but do not share common ancestry, providing a window on the different structural strategies nature has used to evolve required catalysts. Identification and use of this information to improve reaction classification and computational annotation of enzymes newly discovered in the genome projects would benefit from systematic determination of reaction similarities. Here, we quantified similarity in bond changes for overall reactions and catalytic mechanisms for 95 pairs of functionally analogous enzymes (non-homologous enzymes with identical first three numbers of their EC codes) from the MACiE database. Similarity of overall reactions was computed by comparing the sets of bond changes in the transformations from substrates to products. For similarity of mechanisms, sets of bond changes occurring in each mechanistic step were compared; these similarities were then used to guide global and local alignments of mechanistic steps. Using this metric, only 44% of pairs of functionally analogous enzymes in the dataset had significantly similar overall reactions. For these enzymes, convergence to the same mechanism occurred in 33% of cases, with most pairs having at least one identical mechanistic step. Using our metric, overall reaction similarity serves as an upper bound for mechanistic similarity in functional analogs. For example, the four carbon-oxygen lyases acting on phosphates (EC 4.2.3) show neither significant overall reaction similarity nor significant mechanistic similarity. By contrast, the three carboxylic-ester hydrolases (EC 3.1.1) catalyze overall reactions with identical bond changes and have converged to almost identical mechanisms. The large proportion of enzyme pairs that do not show significant overall reaction similarity (56%) suggests that at least for the functionally analogous enzymes studied here, more stringent criteria could be used to refine definitions of EC sub-subclasses for improved discrimination in their classification of enzyme reactions. The results also indicate that mechanistic convergence of reaction steps is widespread, suggesting that quantitative measurement of mechanistic similarity can inform approaches for functional annotation.     

Use of pathway analysis and genome context methods for functional genomics of Mycoplasma pneumoniae nucleotide metabolism

Elementary modes analysis allows one to reveal whether a set of known enzymes is sufficient to sustain functionality of the cell. Moreover, it is helpful in detecting missing reactions and predicting which enzymes could fill these gaps. Here, we perform a comprehensive elementary modes analysis and a genomic context analysis of Mycoplasma pneumoniae nucleotide metabolism, and search for new enzyme activities. The purine and pyrimidine networks are reconstructed by assembling enzymes annotated in the genome or found experimentally. We show that these reaction sets are sufficient for enabling synthesis of DNA and RNA in M. pneumoniae. Special focus is on the key modes for growth. Moreover, we make an educated guess on the nutritional requirements of this micro-organism. For the case that M. pneumoniae does not require adenine as a substrate, we suggest adenylosuccinate synthetase (EC 6.3.4.4), adenylosuccinate lyase (EC 4.3.2.2) and GMP reductase (EC 1.7.1.7) to be operative. GMP reductase activity is putatively assigned to the NRDI_MYCPN gene on the basis of the genomic context analysis. For the pyrimidine network, we suggest CTP synthase (EC 6.3.4.2) to be active. Further experiments on the nutritional requirements are needed to make a decision. Pyrimidine metabolism appears to be more appropriate as a drug target than purine metabolism since it shows lower plasticity.     

ORENZA: a web resource for studying ORphan ENZyme activities

Background  

Despite the current availability of several hundreds of thousands of amino acid sequences, more than 36% of the enzyme activities (EC numbers) defined by the Nomenclature Committee of the International Union of Biochemistry and Molecular Biology (NC-IUBMB) are not associated with any amino acid sequence in major public databases. This wide gap separating knowledge of biochemical function and sequence information is found for nearly all classes of enzymes. Thus, there is an urgent need to explore these sequence-less EC numbers, in order to progressively close this gap.     

Assignment of EC Numbers to Enzymatic Reactions with Reaction Difference Fingerprints

The EC numbers represent enzymes and enzyme genes (genomic information), but they are also utilized as identifiers of enzymatic reactions (chemical information). In the present work (ECAssigner), our newly proposed reaction difference fingerprints (RDF) are applied to assign EC numbers to enzymatic reactions. The fingerprints of reactant molecules minus the fingerprints of product molecules will generate reaction difference fingerprints, which are then used to calculate reaction Euclidean distance, a reaction similarity measurement, of two reactions. The EC number of the most similar training reaction will be assigned to an input reaction. For 5120 balanced enzymatic reactions, the RDF with a fingerprint length at 3 obtained at the sub-subclass, subclass, and main class level with cross-validation accuracies of 83.1%, 86.7%, and 92.6% respectively. Compared with three published methods, ECAssigner is the first fully automatic server for EC number assignment. The EC assignment system (ECAssigner) is freely available via: http://cadd.whu.edu.cn/ecassigner/.     

Energy coupling in type II topoisomerases: why do they hydrolyze ATP?

Type II topoisomerases are essential enzymes in all cells. They help to solve the topological problems of DNA by passing one double helix through a transient break in another, in a reaction coupled to the hydrolysis of ATP. Members of one class of the enzymes, DNA gyrases, are configured to carry out an intramolecular reaction, removing positive supercoiling and introducing negative supercoiling into circular DNA using free energy derived from ATP hydrolysis. The nonsupercoiling class, including bacterial topoisomerase IV and eukaryotic topoisomerase II enzymes, can carry out both intra- and intermolecular reactions, and their primary role is the unlinking (decatenation) of daughter replicons before partition. In these enzymes, ATP hydrolysis is coupled to a reduction in DNA complexity (catenation, supercoiling, and knotting) below the level expected at equilibrium. This review discusses our current understanding of the mechanisms behind the coupling of the energy of ATP hydrolysis to topological changes catalyzed by both of these classes of enzyme.     

Oligonucleotide facilitators enhance the catalytic activity of RNA-cleaving DNA enzymes.

From in vitro selection studies, DNA structures have been found that cleave target RNA sequence specifically and show a certain similarity to the well-investigated hammerhead ribozymes. Such DNA enzymes are more resistant to nuclease-mediated degradation than RNA enzymes. On the other hand, their cleavage activity is lower than the activity of hammerhead ribozymes. In the present study, we improved the activity of DNA enzymes by adding oligonucleotide facilitators complementary to the 5' and the 3' ends of the substrate to the cleavage reaction. DNA enzyme activity in vitro was monitored under multiple turnover conditions using short RNA model substrates. We have shown that oligonucleotide facilitators strongly enhance the multiple turnover activity of the DNA enzyme reaction. In one of our model systems with a suitable facilitator combination, we were able to observe a more than 200-fold enhancement of the k(cat)/Km value. The comparison of two DNA enzyme-substrate systems showed that the principal effects of the facilitators were independent of the substrate sequence. However, the degree of facilitator effect was noticeably dependent on the basic catalytic efficiency of DNA enzymes. Furthermore, the efficiency of the DNA enzyme reaction with facilitator was compared with the reaction of a DNA enzyme with a stem sequence extended by the sequence of the facilitator. The multiple turnover activity of such a "long DNA enzyme" is higher than the activity of the short DNA enzyme without facilitators. However, when compared with the multiple turnover reactions of the short DNA enzyme with facilitator, the reaction with the long DNA enzyme is considerably slower. The results obtained with our model systems demonstrate that oligonucleotide facilitators enable DNA enzymes to act as effective multiple turnover catalysts by cleavage of RNA substrates.     

Automatic Assignment of EC Numbers

A wide range of research areas in molecular biology and medical biochemistry require a reliable enzyme classification system, e.g., drug design, metabolic network reconstruction and system biology. When research scientists in the above mentioned areas wish to unambiguously refer to an enzyme and its function, the EC number introduced by the Nomenclature Committee of the International Union of Biochemistry and Molecular Biology (IUBMB) is used. However, each and every one of these applications is critically dependent upon the consistency and reliability of the underlying data for success. We have developed tools for the validation of the EC number classification scheme. In this paper, we present validated data of 3788 enzymatic reactions including 229 sub-subclasses of the EC classification system. Over 80% agreement was found between our assignment and the EC classification. For 61 (i.e., only 2.5%) reactions we found that their assignment was inconsistent with the rules of the nomenclature committee; they have to be transferred to other sub-subclasses. We demonstrate that our validation results can be used to initiate corrections and improvements to the EC number classification scheme.     

Enzyme-induced inactivation of transminases by acetylenic and vinyl analogues of 4-aminobutyrate.

The reactions of two analogues of 4-aminobutyrate, namely 4-aminohex-5-ynoate and 4-aminohex-5-enoate, with three transaminases were studied. Three pure enzymes were used, aminobutyrate transaminase (EC 2.6.1.19), ornithine transaminase (EC 2.6.1.13) and aspartate transaminase (EC 2.6.1.1), and the course of the reactions was studied by observing changes in the absorption spectrum of the bound coenzyme and by observing loss of activity. All of the enzymes were inactivated by either inhibitor, but amino-hexenoate showed a marked specificity for aminobutyrate transaminase. Aminohexynoate was most potent towards ornithine transaminase, and with this enzyme transamination of the inhibitor is an important factor in protecting the enzyme. Most of the reactions could be analysed as first order, with the observed rate constant showing a hyperbolic dependence on inhibitor concentration.     

Accelerated Search Kinetics Mediated by Redox Reactions of DNA Repair Enzymes

A charge transport (CT) mechanism has been proposed in several articles to explain the localization of base excision repair (BER) enzymes to lesions on DNA. The CT mechanism relies on redox reactions of iron-sulfur cofactors that modify the enzyme's binding affinity. These redox reactions are mediated by the DNA strand and involve the exchange of electrons between BER enzymes along DNA. We propose a mathematical model that incorporates enzyme binding/unbinding, electron transport, and enzyme diffusion along DNA. Analysis of our model within a range of parameter values suggests that the redox reactions can increase desorption of BER enzymes not already bound to lesions, allowing the enzymes to be recycled—thus accelerating the overall search process. This acceleration mechanism is most effective when enzyme copy numbers and enzyme diffusivity along the DNA are small. Under such conditions, we find that CT BER enzymes find their targets more quickly than simple passive enzymes that simply attach to the DNA without desorbing.     

The product-selective blot: a technique for measuring enzyme activities in large numbers of samples and in native electrophoresis gels

A method termed "product-selective" blotting has been developed for screening large numbers of samples for enzyme activity. The technique is particularly well suited to detection of enzymes in native electrophoresis gels. The principle of the method was demonstrated by blotting samples from glutaminase (EC 3.5.1.2) or glutamate synthase (EC 1.4.7.1) reactions into an agarose gel embedded with ion-exchange resin under conditions favoring binding of product (glutamate) over substrates and other substances in the reaction mixture. After washes to remove these unbound substances, the product was measured using either fluorometric staining or radiometric techniques. Glutaminase activity in native electrophoresis gels was visualized by a related procedure in which substrates and products from reactions run in the electrophoresis gel were blotted directly into a resin-containing "image gel." Considering the selective-binding materials available for use in the image gel, along with the possible detection systems, this method has potentially broad application.     

Enzyme degradation during drying

During drying of food materials a multitude of chemical reactions and/or physical changes may occur. In this article attention is focused on one of these, namely, inactivation of enzymes during drying. The prediction of enzyme retention during drying is of interest to the pharmaceutical industry for the production of dry enzyme preparations and to the food processing industry in drying operations of food materials containing enzymes. In this article calculated enzyme retentions are presented for different drying histories and shapes of drying particles. In the numerical calculations it is assumed that enzyme degradation kinetics are first-order reactions, of which reaction constants are known as a function of temperature and water concentration in the drying material. From the calculations, conclusions can be drawn about conditions favorable for high enzyme retentions, or for high enzyme degradations.     

Depolymerization of starch and pectin using superporous matrix supported enzymes

Immobilized enzyme catalyzed biotransformations involving macromolecular substrates and/or products are greatly retarded due to slow diffusion of large substrate molecules in and out of the typical enzyme supports. Slow diffusion of macromolecules into the matrix pores can be speeded up by use of macroporous supports as enzyme carriers. Depolymerization reactions of polysaccharides like starch, pectin, and dextran to their respective low molecular weight products are some of the reactions that can benefit from use of such superporous matrices. In the present work, an indigenously prepared rigid cross-linked cellulose matrix (called CELBEADS) has been used as support for immobilizing alpha amylase (1,4-alpha-D-glucan glucanohydrolase, EC 3.2.1.1.) and pectinase (endo-PG: poly(1,4-alpha-galactouronide) glycanohydrolase, EC 3.2.1.15). The immobilized enzymes were used for starch and pectin hydrolysis respectively, in batch, packed bed and expanded bed modes. The macroporosity of CELBEADS was found to permit through-flow and easy diffusion of substrates pectin and starch to enzyme sites in the porous supports and gave reaction rates comparable to the rates obtained using soluble enzymes.     

Leishmania mexicana: purine-metabolizing enzymes of amastigotes and promastigotes

Cultured promastigote and isolated amastigote forms of Leishmania mexicana mexicana have been surveyed for the presence of enzymes involved in purine metabolism. Quantitative but not qualitative differences between the enzymes of two forms were discovered. There were found to be significant differences between the enzyme content of L. m. mexicana and that reported for L. donovani. Extracts of both parasite forms of L. m. mexicana were found to have higher levels of adenine deaminase (EC 3.5.4.2) and guanine deaminase (EC 3.5.4.3) than adenosine deaminase (EC 3.5.4.4). There appeared to be two distinct nucleosidases (EC 3.2.2.1), one active on nucleosides, the other on deoxynucleosides. Phosphorylase (EC 2.4.2.1) could be detected only in the catabolic direction. Nucleotidases were present, but were more active on 3' (EC 3.1.3.6)- than 5' (EC 3.1.3.5)-nucleotides. Phosphoribosyltransferase (EC 2.4.2.7,.8 and .22) and nucleoside kinase (EC 2.7.1.20) activities were detected in both forms. Nucleotide-interconverting enzymes were found to be present, with IMP dehydrogenase (EC 1.2.1.14) being the most active. Cell fractionation experiments revealed that, in the promastigote, enzyme separation within the parasite may play an important part in regulating cellular purine metabolism.     

Convergent evolution of enzyme active sites is not a rare phenomenon

Since convergent evolution of enzyme active sites was first identified in serine proteases, other individual instances of this phenomenon have been documented. However, a systematic analysis assessing the frequency of this phenomenon across enzyme space is still lacking. This work uses the Query3d structural comparison algorithm to integrate for the first time detailed knowledge about catalytic residues, available through the Catalytic Site Atlas (CSA), with the evolutionary information provided by the Structural Classification of Proteins (SCOP) database. This study considers two modes of convergent evolution: (i) mechanistic analogues which are enzymes that use the same mechanism to perform related, but possibly different, reactions (considered here as sharing the first three digits of the EC number); and (ii) transformational analogues which catalyse exactly the same reaction (identical EC numbers), but may use different mechanisms. Mechanistic analogues were identified in 15% (26 out of 169) of the three-digit EC groups considered, showing that this phenomenon is not rare. Furthermore 11 of these groups also contain transformational analogues. The catalytic triad is the most widespread active site; the results of the structural comparison show that this mechanism, or variations thereof, is present in 23 superfamilies. Transformational analogues were identified for 45 of the 951 four-digit EC numbers present within the CSA and about half of these were also mechanistic analogues exhibiting convergence of their active sites. This analysis has also been extended to the whole Protein Data Bank to provide a complete and manually curated list of the all the transformational analogues whose structure is classified in SCOP. The results of this work show that the phenomenon of convergent evolution is not rare, especially when considering large enzymatic families.     

A comparison of the catalytic properties of cellobiose:quinone oxidoreductase and cellobiose oxidase from Phanerochaete chrysosporium.

Several catalytic properties of the FAD enzyme cellobiose:quinone oxidoreductase (CBQ) and the heme/FAD enzyme, cellobiose oxidase (CBO) have been investigated and compared. Dichlorophenol-indophenol was found to be a very good electron acceptor for cellobiose oxidation by both enzymes. The optimal pH value for this oxidation with dichlorophenol-indophenol as a co-substrate was observed around pH 4 for both enzymes. The turnover numbers of this reaction were also very similar. The Km values for cellobiose oxidation were identical, whereas the Km for CBO with dichlorophenol-indophenol is lower than that of CBQ. Atmospheric oxygen is a very poor electron acceptor for both CBO and CBQ, however, CBO can utilize cytochrome c as an effective electron acceptor, while CBQ cannot. The specific activity of CBO for cytochrome c is thus about 200-times higher than for oxygen. Thus, one way to distinguish the two enzymes is by the cytochrome-c-reducing ability of CBO. Therefore, we propose that the nomenclature for CBO is tentatively changed to cellobiose:cytochrome c oxidoreductase until a rational name can be installed. Both enzymes have radical-reducing activities. The cation radical, derived from 1,2,4,5-tetramethoxybenzene, was reduced by both enzymes at almost the same reaction rate. The phenoxyradical produced by lignin peroxidase, catalyzing the oxidation of acetosyringon, was also reduced by both enzymes. The reduction of phenoxyradicals formed by phenoloxidases (lignin peroxidases, as well as laccases) may be important in preventing repolymerization reactions which we suggest would significantly facilitate lignin degradation.     

Methanol esterification reactions catalyzed by snake venom and bovine intestinal 5'-nucleotide phosphodiesterases. Formation of nucleoside 5'-monophosphate methyl esters from guanosine 5'-triphosphate and other nucleoside 5'-polyphosphates.

It is not known whether the enzymes 5'-nucleotide phosphodiesterase/nucleotide pyrophosphatase (EC 3.1.4.1/EC 3.6.1.9) catalyze the transfer of nucleotides to acceptors other than water. We have investigated the action of snake venom and bovine intestinal mucosa phosphodiesterases on nucleoside 5'-polyphosphates in the presence of methanol. In those conditions, GTP was converted by snake venom phosphodiesterase to a mixture of GMP and another compound with a different retention time in reverse-phase high-performance liquid chromatography. That compound, by ultraviolet, 1H- and 13C-nuclear magnetic resonance spectroscopic analysis, and by enzyme analysis, was characterized as the methyl ester of GMP (GMP-OMe). The molar fraction [GMP-OMe]/[GMP + GMP-OMe] formed was higher than the molar fraction of methanol as a solvent in reaction mixtures. Similar reactions took place at comparable rates with snake venom and bovine intestinal mucosa phosphodiesterases using several nucleoside 5'-polyphosphates as substrates. The ability of 5'-nucleotide phosphodiesterases to catalyze transfer reactions to a non-water acceptor is relevant to the mechanism of the enzymes, to their use as analytical tools, and to their possible use/role in the preparative/in vivo synthesis of nucleotide esters.     

Kinetics of protein-modification reactions. Stoichiometry of modification-produced enzyme inactivation: modification of rhodanese by 2,4,6-trinitrobenzenesulphonic acid.

A mathematical treatment is presented for the dependence of enzyme activity loss on the numbers and reactivities of the groups essential for catalytic function, when enzyme protein modification is carried out by the use of concentrations of protein reactive groups well in excess of that of modifying agent. Experimentally obtained data on the modification of rhodanese (thiosulphate sulphurtransferase, EC 2.8.1.1) by 2,4,6-trinitrobenzenesulphonic acid are presented, and it is shown that, at pH9.00, the fractional concentration of rhodanese groups, or of rhodanese group reactivities, essential for enzyme catalytic function is 0.88; this value is found to decrease with decreasing pH of the reaction medium. The possibility that rhodanese inactivation by 2,4,6-trinitrobenzenesulphonic acid is brought about by modification of groups other than amino groups is ruled out by a comparison of the enzyme-inactivation and protein-modification stoichiometries, for putative reaction models for enzyme and modifying agent.     

Br?nsted-Evans-Polanyi relationships for C–C bond forming and C–C bond breaking reactions in thiamine-catalyzed decarboxylation of 2-keto acids using density functional theory

The concept of generalized enzyme reactions suggests that a wide variety of substrates can undergo enzymatic transformations, including those whose biotransformation has not yet been realized. The use of quantum chemistry to evaluate kinetic feasibility is an attractive approach to identify enzymes for the proposed transformation. However, the sheer number of novel transformations that can be generated makes this impractical as a screening approach. Therefore, it is essential to develop structure/activity relationships based on quantities that are more efficient to calculate. In this work, we propose a structure/activity relationship based on the free energy of binding or reaction of non-native substrates to evaluate the catalysis relative to that of native substrates. While Br?nsted-Evans-Polanyi (BEP) relationships such as that proposed here have found broad application in heterogeneous catalysis, their extension to enzymatic catalysis is limited. We report here on density functional theory (DFT) studies for C–C bond formation and C–C bond cleavage associated with the decarboxylation of six 2-keto acids by a thiamine-containing enzyme (EC 1.2.7.1) and demonstrate a linear relationship between the free energy of reaction and the activation barrier. We then applied this relationship to predict the activation barriers of 17 chemically similar novel reactions. These calculations reveal that there is a clear correlation between the free energy of formation of the transition state and the free energy of the reaction, suggesting that this method can be further extended to predict the kinetics of novel reactions through our computational framework for discovery of novel biochemical transformations.