Analysis of catalytic residues in enzyme active sites

We present an analysis of the residues directly involved in catalysis in 178 enzyme active sites. Specific criteria were derived to define a catalytic residue, and used to create a catalytic residue dataset, which was then analysed in terms of properties including secondary structure, solvent accessibility, flexibility, conservation, quaternary structure and function. The results indicate the dominance of a small set of amino acid residues in catalysis and give a picture of a general active site environment. It is hoped that this information will provide a better understanding of the molecular mechanisms involved in catalysis and a heuristic basis for predicting catalytic residues in enzymes of unknown function.     

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included α-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein. Figure Collective motions in C^α atoms of the active site of cold-active xylanase Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

The human TREX2 3' -> 5'-exonuclease structure suggests a mechanism for efficient nonprocessive DNA catalysis

The 3' --> 5'-exonucleases process DNA ends in many DNA repair pathways of human cells. Determination of the human TREX2 structure is the first of a dimeric 3'-deoxyribonuclease and indicates how this highly efficient nonprocessive enzyme removes nucleotides at DNA 3' termini. Symmetry in the TREX2 dimer positions the active sites at opposite outer edges providing open access for the DNA. Adjacent to each active site is a flexible region containing three arginines positioned appropriately to bind DNA and to control its entry into the active site. Mutation of these three arginines to alanines reduces the DNA binding capacity by approximately 100-fold with no effect on catalysis. The human TREX2 catalytic residues overlay with the bacterial DnaQ family of 3'-exonucleases confirming the structural conservation of the catalytic sites despite limited sequence identity, and mutations of these residues decrease the still measurable activity by approximately 10(5)-fold, confirming their catalytic role.     

Role of active site rigidity in activity: MD simulation and fluorescence study on a lipase mutant

Relationship between stability and activity of enzymes is maintained by underlying conformational flexibility. In thermophilic enzymes, a decrease in flexibility causes low enzyme activity while in less stable proteins such as mesophiles and psychrophiles, an increase in flexibility is associated with enhanced enzyme activity. Recently, we identified a mutant of a lipase whose stability and activity were enhanced simultaneously. In this work, we probed the conformational dynamics of the mutant and the wild type lipase, particularly flexibility of their active site using molecular dynamic simulations and time-resolved fluorescence techniques. In contrast to the earlier observations, our data show that active site of the mutant is more rigid than wild type enzyme. Further investigation suggests that this lipase needs minimal reorganization/flexibility of active site residues during its catalytic cycle. Molecular dynamic simulations suggest that catalytically competent active site geometry of the mutant is relatively more preserved than wild type lipase, which might have led to its higher enzyme activity. Our study implies that widely accepted positive correlation between conformation flexibility and enzyme activity need not be stringent and draws attention to the possibility that high enzyme activity can still be accomplished in a rigid active site and stable protein structures. This finding has a significant implication towards better understanding of involvement of dynamic motions in enzyme catalysis and enzyme engineering through mutations in active site.     

Protein flexibility and rigidity predicted from sequence

Structural flexibility has been associated with various biological processes such as molecular recognition and catalytic activity. In silico studies of protein flexibility have attempted to characterize and predict flexible regions based on simple principles. B-values derived from experimental data are widely used to measure residue flexibility. Here, we present the most comprehensive large-scale analysis of B-values. We used this analysis to develop a neural network-based method that predicts flexible-rigid residues from amino acid sequence. The system uses both global and local information (i.e., features from the entire protein such as secondary structure composition, protein length, and fraction of surface residues, and features from a local window of sequence-consecutive residues). The most important local feature was the evolutionary exchange profile reflecting sequence conservation in a family of related proteins. To illustrate its potential, we applied our method to 4 different case studies, each of which related our predictions to aspects of function. The first 2 were the prediction of regions that undergo conformational switches upon environmental changes (switch II region in Ras) and the prediction of surface regions, the rigidity of which is crucial for their function (tunnel in propeller folds). Both were correctly captured by our method. The third study established that residues in active sites of enzymes are predicted by our method to have unexpectedly low B-values. The final study demonstrated how well our predictions correlated with NMR order parameters to reflect motion. Our method had not been set up to address any of the tasks in those 4 case studies. Therefore, we expect that this method will assist in many attempts at inferring aspects of function.     

Role of protein conformational dynamics in the catalysis by 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

Enzymatic catalysis has conflicting structural requirements of the enzyme. In order for the enzyme to form a Michaelis complex, the enzyme must be in an open conformation so that the substrate can get into its active center. On the other hand, in order to maximize the stabilization of the transition state of the enzymatic reaction, the enzyme must be in a closed conformation to maximize its interactions with the transition state. The conflicting structural requirements can be resolved by a flexible active center that can sample both open and closed conformational states. For a bisubstrate enzyme, the Michaelis complex consists of two substrates in addition to the enzyme. The enzyme must remain flexible upon the binding of the first substrate so that the second substrate can get into the active center. The active center is fully assembled and stabilized only when both substrates bind to the enzyme. However, the side-chain positions of the catalytic residues in the Michaelis complex are still not optimally aligned for the stabilization of the transition state, which lasts only approximately 10(-13) s. The instantaneous and optimal alignment of catalytic groups for the transition state stabilization requires a dynamic enzyme, not an enzyme which undergoes a large scale of movements but an enzyme which permits at least a small scale of adjustment of catalytic group positions. This review will summarize the structure, catalytic mechanism, and dynamic properties of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase and examine the role of protein conformational dynamics in the catalysis of a bisubstrate enzymatic reaction.     

Effect of a new non-cleavable substrate analog on wild-type and serine mutants in the signature sequence of adenylosuccinate lyase of Bacillus subtilis and Homo sapiens

Adenylosuccinate lyase (ASL) catalyzes two beta-elimination reactions in purine biosynthesis, leading to the question of whether the two substrates occupy the same or different active sites. Kinetic studies of Bacillus subtilis and human ASL with a new substrate analog, adenosine phosphonobutyric acid, 2'(3'), 5'-diphosphate (APBADP), show that it acts as a competitive inhibitor with respect to either substrate (K(I) approximately 0.1 microM), indicating that the two substrates occupy the same active site. Binding studies show that both the B. subtilis and human ASLs bind up to 4 mol of APBADP per mole of enzyme tetramer and that both enzymes exhibit cooperativity: negative for B. subtilis ASL and positive for human ASL. Mutant B. subtilis ASLs, with replacements for residues previously identified as critical for catalysis, bind the substrate analog similarly to wild-type ASL. Two serines in a flexible loop of ASL have been proposed to play roles in catalysis because they are close to the substrate in the crystal structure of Escherichia coli ASL. We have now mutated the corresponding serines to alanines in B. subtilis and human ASL to evaluate their involvement in enzyme function. Kinetic data reveal that human Ser(289) and B. subtilis Ser(262) and Ser(263) are essential for catalysis, while the ability of these Ser mutants to bind APBADP suggests that they do not contribute to substrate affinity. Although these serines are not visible in the crystal structure of human adenylosuccinate lyase complexed with substrate or products (PDB #2VD6), they may be interacting with the active sites.     

Sequence and structural differences between enzyme and nonenzyme homologs

To improve our understanding of the evolution of novel functions, we performed a sequence, structural, and functional analysis of homologous enzymes and nonenzymes of known three-dimensional structure. In most examples identified, the nonenzyme is derived from an ancestral catalytic precursor (as opposed to the reverse evolutionary scenario, nonenzyme to enzyme), and the active site pocket has been disrupted in some way, owing to the substitution of critical catalytic residues and/or steric interactions that impede substrate binding and catalysis. Pairwise sequence identity is typically insignificant, and almost one-half of the enzyme and nonenzyme pairs do not share any similarity in function. Heterooligomeric enzymes comprising homologous subunits in which one chain is catalytically inactive and enzyme polypeptides that contain internal catalytic and noncatalytic duplications of an ancient enzyme domain are also discussed.     

A tale of two isomerases: compact versus extended active sites in ketosteroid isomerase and phosphoglucose isomerase

Understanding the catalytic efficiency and specificity of enzymes is a fundamental question of major practical and conceptual importance in biochemistry. Although progress in biochemical and structural studies has enriched our knowledge of enzymes, the role in enzyme catalysis of residues that are not nearest neighbors of the reacting substrate molecule is largely unexplored experimentally. Here computational active site predictors, THEMATICS and POOL, were employed to identify functionally important residues that are not in direct contact with the reacting substrate molecule. These predictions then guided experiments to explore the active sites of two isomerases, Pseudomonas putida ketosteroid isomerase (KSI) and human phosphoglucose isomerase (PGI), as prototypes for very different types of predicted active sites. Both KSI and PGI are members of EC 5.3 and catalyze similar reactions, but they represent significantly different degrees of remote residue participation, as predicted by THEMATICS and POOL. For KSI, a compact active site of mostly first-shell residues is predicted, but for PGI, an extended active site in which residues in the first, second, and third layers around the reacting substrate are predicted. Predicted residues that have not been previously tested experimentally were investigated by site-directed mutagenesis and kinetic analysis. In human PGI, single-point mutations of the predicted second- and third-shell residues K362, H100, E495, D511, H396, and Q388 show significant decreases in catalytic activity relative to that of the wild type. The results of these experiments demonstrate that, as predicted, remote residues are very important in PGI catalysis but make only small contributions to catalysis in KSI.     

活性部位的柔性

比较酶在变性过程中构象和活力变化,发现在活性完全丧失时尚无可察 觉的整体构象变化。排除变性剂抑制和寡聚酶解聚等可能性之后,提出了酶活性部位柔性假说。随后用多种实验方法直接证实了活性部位的构象变化先于分子整体构象变化,并与活性丧失同步,根据催化过程中活性部位构旬变化,以及限制活性部位构象变化对酶活性的影响,提出了酶活性部位柔性为酶充分表现其催化活性所必需的设想。     

Analysis and prediction of the location of catalytic residues in enzymes

The catalytic residues of an enzyme are defined as the amino acids directly involved in chemical catalysis. They mainly act as a general acid--base, electrophilic or nucleophilic catalyst or they polarize and stabilize the transition state. An analysis of the structural features of 36 catalytic residues in 17 enzymes of known structure and with defined mechanism is reported. Residues that bind metal ions (Zn2+ and Cu2+) are considered separately. The features examined are: residue type, location in secondary structure, separation between the residues, accessibility to solvent, intra-protein electrostatic interactions, mobility as evaluated from crystallographic temperature factors, polarity of the environment and the sequence conservation between homologous enzymes of residues that were sequentially or spatially close to the catalytic residue. In general the environment of catalytic residues is similar to that of polar side chains that have low accessibility to solvent. Two algorithms have been developed to identify probable catalytic residues. Scanning an alignment of homologous enzyme sequences for peaks of sequence conservation identifies 13 out of the 16 catalytic residues with 50 residues overpredicted. When the conservation of the spatially close residues is used instead, a different set of 13 residues are identified with 47 residues overpredicted. A combination of the two algorithms identifies 11 residues with 36 residues overpredicted.     

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.     

Conformational mobility in the active site of a heme peroxidase

Conformational mobility of the distal histidine residue has been implicated for several different heme peroxidase enzymes, but unambiguous structural evidence is not available. In this work, we present mechanistic, spectroscopic, and structural evidence for peroxide- and ligand-induced conformational mobility of the distal histidine residue (His-42) in a site-directed variant of ascorbate peroxidase (W41A). In this variant, His-42 binds "on" to the heme in the oxidized form, duplicating the active site structure of the cytochromes b but, in contrast to the cytochromes b, is able to swing "off" the iron during catalysis. This conformational flexibility between the on and off forms is fully reversible and is used as a means to overcome the inherently unreactive nature of the on form toward peroxide, so that essentially complete catalytic activity is maintained. Contrary to the widely adopted view of heme enzyme catalysis, these data indicate that strong coordination of the distal histidine to the heme iron does not automatically undermine catalytic activity. The data add a new dimension to our wider appreciation of structure/activity correlations in other heme enzymes.     

Conformational dynamics of free and catalytically active thermolysin are indistinguishable by hydrogen/deuterium exchange mass spectrometry

Conformational dynamics are thought to be a prerequisite for the catalytic activity of enzymes. However, the exact relationship between structural fluctuations and function is not well understood. In this work hydrogen/deuterium exchange (HDX) and electrospray ionization mass spectrometry (ESI-MS) are used for exploring the conformational dynamics of thermolysin. Amide HDX reflects the internal mobility of proteins; regions that undergo frequent unfolding-refolding show faster exchange than segments that are highly stable. Thermolysin is a zinc protease with an active site that is located between two lobes. Substrate turnover is associated with hinge bending that leads to a closed conformation. Product release regenerates the open form, such that steady-state catalysis involves a continuous closing/opening cycle. HDX/ESI-MS with proteolytic peptide mapping in the absence of substrate shows that elements in the periphery of the two lobes are most mobile. A comparison with previous X-ray data suggests that these peripheral regions undergo quite pronounced structural changes during the catalytic cycle. In contrast, active site residues exhibit only a moderate degree of backbone flexibility, and the central zinc appears to be in a fairly rigid environment. The presence of both rigid and moderately flexible elements in the active site may reflect a carefully tuned balance that is required for function. Interestingly, the HDX behavior of catalytically active thermolysin is indistinguishable from that of the free enzyme. This result is consistent with the view that catalytically relevant motions preexist in the resting state and that enzyme function can only be performed within the limitations given by the intrinsic dynamics of the protein. The data presented in this work indicate the prevalence of stochastic elements in the function of thermolysin, rather than supporting a deterministic mechanism.     

Flexibility analysis of enzyme active sites by crystallographic temperature factors

Protein flexibility is inherent to protein structural behavior. Experimental evidence for protein flexibility is extensive both in solution and in the solid state. A major question is whether the flexibility observed in enzymes is simply an inherent property of proteins that must always be borne in mind or is essential for catalysis or substrate binding. The temperature factors or B-values, as determined crystallographically, are linearly related to the mean square displacement of an atom and give an indication of atomic flexibility in the crystalline state. In this paper, we describe the frequency distributions of the normalized B-factor (B'-factor) for the active site and non-active site residues in the selected 69 apo-enzymes. This analysis was performed over the entire sequences and for different structural subsets defined by the three-dimensional structure of proteins, as alpha-helices, beta-structures and coil conformation and buried and non-buried residues. The results show that in all cases, the active site residues predominantly occur in region of low B'-factor and the non-active site residues have a tendency to exist in the high B'-factor region. This observation suggests that the active site residues, in general, are less flexible than the non-active site residues and therefore the vibrational and the fast collective motions of the C(alpha) atoms of proteins appear not to have clear biological significance.     

The enzyme activity allosteric regulation model based on the composite nature of catalytic and regulatory sites concept

A new kinetic model of enzymatic catalysis is proposed, which postulates that enzyme solutions are equilibrium systems of oligomers differing in the number of subunits and in the mode of their assembly. It is suggested that the catalytic and regulatory sites of allosteric enzymes are of composite nature and appear as a result of subunits joining. Two possible joining modes are postulated at each oligomerization step. Catalytic site may arise on oligomer formed only by one of these modes. Effector acts by fastening together components of certain oligomeric form and increases the life time of this form. It leads to a shift of oligomer equilibrium and increases a proportion of effector-binding oligomers. Effectors-activators bind the oligomers carrying composite catalytic sites and effectors-inhibitors bind the oligomers, which do not carry active catalytic sites. Thus, catalytic activity control in such system is explained by effector-induced changes of a catalytic sites number, but not of a catalytic site activity caused by changes of subunit's tertiary structure. The postulates of the model do not contradict available experimental data and lead to a new type of general rate equation, which allows to describe and understand the specific kinetic behavior of allosteric enzymes as well as Michaelis type enzymes. All known rate equations of allosteric The equation was tested by modeling the kinetics of human erythrocyte phosphofructokinase. It enabled to reproduce quantitatively the 66 kinetic curves experimentally obtained for this enzyme under different reaction conditions.     

Conservation of Flexible Residue Clusters among Structural and Functional Enzyme Homologues

Conformational flexibility between structural ensembles is an essential component of enzyme function. Although the broad dynamical landscape of proteins is known to promote a number of functional events on multiple time scales, it is yet unknown whether structural and functional enzyme homologues rely on the same concerted residue motions to perform their catalytic function. It is hypothesized that networks of contiguous and flexible residue motions occurring on the biologically relevant millisecond time scale evolved to promote and/or preserve optimal enzyme catalysis. In this study, we use a combination of NMR relaxation dispersion, model-free analysis, and ligand titration experiments to successfully capture and compare the role of conformational flexibility between two structural homologues of the pancreatic ribonuclease family: RNase A and eosinophil cationic protein (or RNase 3). In addition to conserving the same catalytic residues and structural fold, both homologues show similar yet functionally distinct clusters of millisecond dynamics, suggesting that conformational flexibility can be conserved among analogous protein folds displaying low sequence identity. Our work shows that the reduced conformational flexibility of eosinophil cationic protein can be dynamically and functionally reproduced in the RNase A scaffold upon creation of a chimeric hybrid between the two proteins. These results support the hypothesis that conformational flexibility is partly required for catalytic function in homologous enzyme folds, further highlighting the importance of dynamic residue sectors in the structural organization of proteins.     

On the relationship between thermal stability and catalytic power of enzymes

The possible relationship between the thermal stability and the catalytic power of enzymes is of great current interest. In particular, it has been suggested that thermophilic or hyperthermophilic (Tm) enzymes have lower catalytic power at a given temperature than the corresponding mesophilic (Ms) enzymes, because the thermophilic enzymes are less flexible (assuming that flexibility and catalysis are directly correlated). These suggestions presume that the reduced dynamics of the thermophilic enzymes is the reason for their reduced catalytic power. The present paper takes the specific case of dihydrofolate reductase (DHFR) and explores the validity of the above argument by simulation approaches. It is found that the Tm enzymes have restricted motions in the direction of the folding coordinate, but this is not relevant to the chemical process, since the motions along the reaction coordinate are perpendicular to the folding motions. Moreover, it is shown that the rate of the chemical reaction is determined by the activation barrier and the corresponding reorganization energy, rather than by dynamics or flexibility in the ground state. In fact, as far as flexibility is concerned, we conclude that the displacement along the reaction coordinate is larger in the Tm enzyme than in the Ms enzyme and that the general trend in enzyme catalysis is that the best catalyst involves less motion during the reaction than the less optimal catalyst. The relationship between thermal stability and catalysis appears to reflect the fact that to obtain small electrostatic reorganization energy it is necessary to invest some folding energy in the overall preorganization process. Thus, the optimized catalysts are less stable. This trend is clearly observed in the DHFR case.     

Understanding nature's catalytic toolkit

Enzymes catalyse numerous reactions in nature, often causing spectacular accelerations in the catalysis rate. One aspect of understanding how enzymes achieve these feats is to explore how they use the limited set of residue side chains that form their 'catalytic toolkit'. Combinations of different residues form 'catalytic units' that are found repeatedly in different unrelated enzymes. Most catalytic units facilitate rapid catalysis in the enzyme active site either by providing charged groups to polarize substrates and to stabilize transition states, or by modifying the pKa values of other residues to provide more effective acids and bases. Given recent efforts to design novel enzymes, the rise of structural genomics and subsequent efforts to predict the function of enzymes from their structure, these units provide a simple framework to describe how nature uses the tools at her disposal, and might help to improve techniques for designing and predicting enzyme function.     

Intermediates in serine recombinase-mediated site-specific recombination

Site-specific recombinases are enzymes that promote precise rearrangements of DNA sequences. They do this by cutting and rejoining the DNA strands at specific positions within a pair of target sites recognized and bound by the recombinase. One group of these enzymes, the serine recombinases, initiates strand exchange by making double-strand breaks in the DNA of the two sites, in an intermediate built around a catalytic tetramer of recombinase subunits. However, these catalytic steps are only the culmination of a complex pathway that begins when recombinase subunits recognize and bind to their target sites as dimers. To form the tetramer-containing reaction intermediate, two dimer-bound sites are brought together by protein dimer-dimer interactions. During or after this initial synapsis step, the recombinase subunit and tetramer conformations change dramatically by repositioning of component subdomains, bringing about a transformation of the enzyme from an inactive to an active configuration. In natural serine recombinase systems, these steps are subject to elaborate regulatory mechanisms in order to ensure that cleavage and rejoining of DNA strands only happen when and where they should, but we and others have identified recombinase mutants that have lost dependence on this regulation, thus facilitating the study of the basic steps leading to catalysis. We describe how our studies on activated mutants of two serine recombinases, Tn3 resolvase and Sin, are providing us with insights into the structural changes that occur before catalysis of strand exchange, and how these steps in the reaction pathway are regulated.