Conformational Variation in Enzyme Catalysis: A Structural Study on Catalytic Residues

Conformational variation in catalytic residues can be captured as alternative snapshots in enzyme crystal structures. Addressing the question of whether active site flexibility is an intrinsic and essential property of enzymes for catalysis, we present a comprehensive study on the 3D variation of active sites of 925 enzyme families, using explicit catalytic residue annotations from the Mechanism and Catalytic Site Atlas and structural data from the Protein Data Bank. Through weighted pairwise superposition of the functional atoms of active sites, we captured structural variability at single-residue level and examined the geometrical changes as ligands bind or as mutations occur. We demonstrate that catalytic centres of enzymes can be inherently rigid or flexible to various degrees according to the function they perform, and structural variability most often involves a subset of the catalytic residues, usually those not directly involved in the formation or cleavage of bonds. Moreover, data suggest that 2/3 of active sites are flexible, and in half of those, flexibility is only observed in the side chain. The goal of this work is to characterise our current knowledge of the extent of flexibility at the heart of catalysis and ultimately place our findings in the context of the evolution of catalysis as enzymes evolve new functions and bind different substrates.     

Non-empirical study of the phosphorylation reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase: relevance of the theory of intermolecular interactions

The subject of this study was an analysis of the role of active site residues in the phosphoryl transfer reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase (ThiK). The ThiK-catalyzed reaction is of special interest due to the lack of a highly conserved aspartate residue serving as a catalytic base. ONIOM(B3LYP:PM3) models of stationary points along the reaction pathway consisted of reactants, two magnesium ions and several highly conserved ThiK active site residues. The results indicate that an S_N2-like mechanism of ThiK, with γ-phosphate acting as an alcohol-activating base is reasonable. Geometries of substrates, transition state and products were utilized in the non-empirical analysis of the physical nature of catalytic interactions taking place in the ThiK active site. The role of particular residues was investigated in terms of their ability to preferentially stabilize the transition state relative to substrates (differential transition state stabilization, DTSS) or products (differential product stabilization, DPS). It seems that Mg2, Glu126 and Cys198 play a major catalytic role, whereas Mg1 and the same Cys198 are responsible for product release. It is remarkable that no dominant role of an electrostatic term in the interactions involved in catalytic activity is observed for product release. Determination of catalytic fields expressing differential electrostatic potential of the transition state with respect to substrates revealed the optimal electrostatic features of an ideal catalyst for the studied reaction. The predicted catalytic environment is in agreement with experimental data showing increased catalytic activity of ThiK upon mutation of Cys198 to aspartate. Figure Catalytic fields for ThiK-catalyzed reaction juxtaposed with the positions of active site residues of a model system. Magnesium ions are considered part of the transition state/reactants. The surface of constant electronic density is colored according to differential electrostatic potential of transition state with respect to reactants. The sign of the differential potential reflects the electrostatic properties of a complementary molecular environment. Red (green) color denotes regions where a negative (positive) charge would be optimal for catalytic activity     

Uniquely Localized Intra-Molecular Amino Acid Concentrations at the Glycolytic Enzyme Catalytic/Active Centers of Archaea, Bacteria and Eukaryota are Associated with Their Proposed Temporal Appearances on Earth

The distributions of amino acids at most-conserved sites nearest catalytic/active centers (C/AC) in 4,645 sequences of ten enzymes of the glycolytic Embden-Meyerhof-Parnas pathway in Archaea, Bacteria and Eukaryota are similar to the proposed temporal order of their appearance on Earth. Glycine, isoleucine, leucine, valine, glutamic acid and possibly lysine often described as prebiotic, i.e., existing or occurring before the emergence of life, were localized in positional and conservational defined aggregations in all enzymes of all Domains. The distributions of all 20 biologic amino acids in most-conserved sites nearest their C/ACs were quite different either from distributions in sites less-conserved and further from their C/ACs or from all amino acids regardless of their position or conservation. The major concentrations of glycine, e.g., perhaps the earliest prebiotic amino acid, occupies ≈16 % of all the most-conserved sites within a volume of ≈7–8 Å radius from their C/ACs and decreases linearly towards the molecule’s peripheries. Spatially localized major concentrations of isoleucine, leucine and valine are in the mid-conserved and mid-distant sites from their C/ACs in protein interiors. Lysine and glutamic acid comprise ≈25–30 % of all amino acids within an irregular volume bounded by ≈24–28 Å radii from their C/ACs at the most-distant least-conserved sites. The unreported characteristics of these amino acids: their spatially and conservationally identified concentrations in Archaea, Bacteria and Eukaryota, suggest some common structural organization of glycolytic enzymes that may be relevant to their evolution and that of other proteins. We discuss our data in relation to enzyme evolution, their reported prebiotic putative temporal appearances on Earth, abundances, biological “cost”, neighbor-sequence preferences or “ordering” and some thermodynamic parameters.     

The role of group bulkiness in the catalytic activity of psychrophile cold-active protein tyrosine phosphatase

The cold-active protein tyrosine phosphatase found in psychrophilic Shewanella species exhibits high catalytic efficiency at low temperatures as well as low thermostability, both of which are characteristics shared by many cold-active enzymes. The structure of cold-active protein tyrosine phosphatase is notable for the presence of three hydrophobic sites (termed the CA, Zn-1 and Zn-2 sites) behind the loop structures comprising the catalytic region. To identify the structural components responsible for specific enzyme characteristics, we determined the structure of wild-type cold-active protein tyrosine phosphatase at high resolution (1.1 A) and measured the catalytic efficiencies of enzymes containing mutations in the three hydrophobic sites. The bulkiness of the amino acid side chains in the core region of the Zn-1 site strongly affects the thermostability and the catalytic efficiency at low temperatures. The mutant enzyme I115M possessed a higher kcat at low temperatures. Elucidation of the crystal structure of I115M at a resolution of 1.5 A revealed that the loop structures involved in retaining the nucleophilic group and the acid catalyst are more flexible than in the wild-type enzyme.     

A computational method for design of connected catalytic networks in proteins

Computational design of new active sites has generally proceeded by geometrically defining interactions between the reaction transition state(s) and surrounding side‐chain functional groups which maximize transition‐state stabilization, and then searching for sites in protein scaffolds where the specified side‐chain–transition‐state interactions can be realized. A limitation of this approach is that the interactions between the side chains themselves are not constrained. An extensive connected hydrogen bond network involving the catalytic residues was observed in a designed retroaldolase following directed evolution. Such connected networks could increase catalytic activity by preorganizing active site residues in catalytically competent orientations, and enabling concerted interactions between side chains during catalysis, for example, proton shuffling. We developed a method for designing active sites in which the catalytic side chains, in addition to making interactions with the transition state, are also involved in extensive hydrogen bond networks. Because of the added constraint of hydrogen‐bond connectivity between the catalytic side chains, to find solutions, a wider range of interactions between these side chains and the transition state must be considered. Our new method starts from a ChemDraw‐like two‐dimensional representation of the transition state with hydrogen‐bond donors, acceptors, and covalent interaction sites indicated, and all placements of side‐chain functional groups that make the indicated interactions with the transition state, and are fully connected in a single hydrogen‐bond network are systematically enumerated. The RosettaMatch method can then be used to identify realizations of these fully‐connected active sites in protein scaffolds. The method generates many fully‐connected active site solutions for a set of model reactions that are promising starting points for the design of fully‐preorganized enzyme catalysts.     

Conserved Structural Chemistry for Incision Activity in Structurally Non-homologous Apurinic/Apyrimidinic Endonuclease APE1 and Endonuclease IV DNA Repair Enzymes

Non-coding apurinic/apyrimidinic (AP) sites in DNA form spontaneously and as DNA base excision repair intermediates are the most common toxic and mutagenic in vivo DNA lesion. For repair, AP sites must be processed by 5′ AP endonucleases in initial stages of base repair. Human APE1 and bacterial Nfo represent the two conserved 5′ AP endonuclease families in the biosphere; they both recognize AP sites and incise the phosphodiester backbone 5′ to the lesion, yet they lack similar structures and metal ion requirements. Here, we determined and analyzed crystal structures of a 2.4 Å resolution APE1-DNA product complex with Mg²⁺ and a 0.92 Å Nfo with three metal ions. Structural and biochemical comparisons of these two evolutionarily distinct enzymes characterize key APE1 catalytic residues that are potentially functionally similar to Nfo active site components, as further tested and supported by computational analyses. We observe a magnesium-water cluster in the APE1 active site, with only Glu-96 forming the direct protein coordination to the Mg²⁺. Despite differences in structure and metal requirements of APE1 and Nfo, comparison of their active site structures surprisingly reveals strong geometric conservation of the catalytic reaction, with APE1 catalytic side chains positioned analogously to Nfo metal positions, suggesting surprising functional equivalence between Nfo metal ions and APE1 residues. The finding that APE1 residues are positioned to substitute for Nfo metal ions is supported by the impact of mutations on activity. Collectively, the results illuminate the activities of residues, metal ions, and active site features for abasic site endonucleases.     

Active site identification through geometry-based and sequence profile-based calculations: burial of catalytic clefts

Electrostatics calculations with proteins that are uniformly charged over volume can aid enzyme/non-enzyme discrimination. For known enzymes, such methods locate active sites to within 5% on the enzyme surface, in 77% of a test set. We now report that removing the dielectric boundary improves active site location to 80%, with optimal discrimination between enzymes and non-enzymes of around 80% specificity and 80% sensitivity. This calculation quantifies burial of solvent-accessible regions. Many of the true enzymes incorrectly assigned as non-enzymes have active sites at subunit boundaries. These are missed in monomer-based calculations. Catalytic and non-catalytic antibodies are studied in this context of active/binding site burial. Whilst catalytic antibodies, on average, have marginally higher active site burial than non-catalytic antibodies, these values are generally smaller than for non-antibody enzymes, possibly contributing to their relatively low turnover. Prediction of active site location improves further when sequence profile-based weights replace the uniform charge distribution, so that a combination of burial and amino acid conservation is assessed. Accuracy rises to 93% of active sites to within 5%, in the test set, for the optimal profile weights scheme. The equivalent value in a separate validation set is 89% to within 5%. Enzyme/non-enzyme and enzyme functional site predictions are made for structural genomics proteins, suggesting that a substantial majority of these are non-enzymes.     

Mechanism of post-translational quinone formation in copper amine oxidases and its relationship to the catalytic turnover

Copper amine oxidases (CAOs) post-translationally construct a redox-active quinone from an amino acid side chain in their polypeptide chain. As such, these enzymes illustrate how nature is able to expand upon naturally-occurring side chains to create new, catalytically powerful functionalities. The active sites of the CAOs are highly unusual in their ability to catalyze two very different reactions: single-turnover, oxygen-dependent quinone formation, followed by catalytic oxidation (formally dehydrogenation) of amines. This review summarizes our current understanding of the pathway whereby the 2,4,5-trihydroxyphenylalanyl quinone (TPQ) cofactor is generated from the phenolic side chain of tyrosine. This reaction occurs spontaneously intermediates in the presence of O(2) and active site bound Cu(II), without the assistance of other proteins or cofactors. Ongoing work has focused on uncovering the details of the TPQ formation mechanism. A larger goal is to understand how a single active site is capable of supporting both quinone formation and subsequent catalytic turnover.     

Distinctive structural motifs co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole in the alpha/beta‐hydrolase fold enzymes

The alpha/beta‐hydrolases (ABH) are among the largest structural families of proteins that are found in nature. Although they vary in their sequence and function, the ABH enzymes use a similar acid–base‐nucleophile catalytic mechanism to catalyze reactions on different substrates. Because ABH enzymes are biocatalysts with a wide range of potential applications, protein engineering has taken advantage of their catalytic versatility to develop enzymes with industrial applications. This study is a comprehensive analysis of 40 ABH enzyme families focusing on two identified substructures: the nucleophile zone and the oxyanion zone, which co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole, and independently reported as critical for the enzymatic activity. We also frequently observed an aromatic cluster near the nucleophile and oxyanion zones, and opposite the ligand‐binding site. The nucleophile zone, the oxyanion zone and the residue cluster enriched in aromatic side chains comprise a three‐dimensional structural organization that shapes the active site of ABH enzymes and plays an important role in the enzymatic function by structurally stabilizing the catalytic nucleophile and the residues of the oxyanion hole. The structural data support the notion that the aromatic cluster can participate in co‐ordination of the catalytic histidine loop, and properly place the catalytic histidine next to the catalytic nucleophile.     

Enzyme catalysis: removing chemically 'essential' residues by site-directed mutagenesis

Enzymatic catalysis relies on the action of the amino acid side chains arrayed in the enzyme active sites. Usually, only two or three 'essential' residues are directly involved in the bond making and breaking steps leading to product formation. For the past 20 years, enzymologists have been addressing the role of such residues by changing them into chemically inert side chains. Removal of an 'essential' group often does not abolish activity, but can significantly alter the catalytic mechanism. Such results underscore the sophistication of enzyme catalysis and the functional plasticity of enzyme active sites.     

Bioinformatics and molecular modeling in chemical enzymology. Active sites of hydrolases

Comparison and multiple alignments of amino acid sequences of a representative number of related enzymes demonstrate the existence of certain positions of amino acid residues which are permanently reproducible in all members of the whole family. The use of the bioinformatic approach revealed conservative residues in each of the related enzymes and ranked amino acid conservatism for the overall enzymatic catalysis. Glycine and aspartic acid residues were shown to be the most essential for structure and catalytic activity of enzymes. Amino acid residues forming catalytic subsite of the active site of enzymes are always highly conservative. Analysis revealed that aspartic acid carboxyl group is the most frequently employed nucleophilic (in deprotonated form) and electrophilic (in protonated form) agent involved in activation of molecules by the mechanism of general base and acidic catalyses in the catalytic sites of enzymes. Glycine is a unique amino acid possessing the highest possibilities for rotation along C–C and C–N bonds of the polypeptide chain. The conservative fixation of the glycine residue in polypeptide chains of related enzymes provides a possibility for directed assembly of amino acid residues into the catalytic subsite structure. It is possible that the conservative glycines provide known conformational mobility of the protein and the active site. Methods of molecular modeling were used for analysis of structural substitutions of conservative and non-conservative glycines and their effects on geometry of catalytic site of typical hydrolases. The substitution of glycine(s) for alanine significantly altered the catalytic site structures.     

The Enzyme Activity Allosteric Regulation Model Based on the Composite Nature of Catalytic and Regulatory Sites Concept

Abstract

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.     

Structure of the complex between trypanosomal triosephosphate isomerase and N-hydroxy-4-phosphono-butanamide: binding at the active site despite an "open" flexible loop conformation.

The structure of triosephosphate isomerase from Trypanosoma brucei complexed with the competitive inhibitor N-hydroxy-4-phosphono-butanamide was determined by X-ray crystallography to a resolution of 2.84 A. Full occupancy binding of the inhibitor is observed only at one of the active sites of the homodimeric enzyme where the flexible loop is locked in a completely open conformation by crystal contacts. There is evidence that the inhibitor also binds to the second active site of the enzyme, but with low occupancy. The hydroxamyl group of the inhibitor forms hydrogen bonds to the side chains of Asn 11, Lys 13, and His 95, whereas each of its three methylene units is involved in nonpolar interactions with the side chain of the flexible loop residue Ile 172. Interactions between the hydroxamyl and the catalytic base Glu 167 are absent. The binding of this phosphonate inhibitor exhibits three unusual features: (1) the flexible loop is open, in contrast with the binding mode observed in eight other complexes between triosephosphate isomerase and various phosphate and phosphonate compounds; (2) compared with these complexes the present structure reveals a 1.5-A shift of the anion-binding site; (3) this is the first phosphonate inhibitor that is not forced by the enzyme into an eclipsed conformation about the P-CH2 bond. The results are discussed with respect to an ongoing drug design project aimed at the selective inhibition of glycolytic enzymes of T. brucei.     

M42 aminopeptidase catalytic site: the structural and functional role of a strictly conserved aspartate residue

The M42 aminopeptidases are a family of dinuclear aminopeptidases widely distributed in Prokaryotes. They are potentially associated to the proteasome, achieving complete peptide destruction. Their most peculiar characteristic is their quaternary structure, a tetrahedron-shaped particle made of twelve subunits. The catalytic site of M42 aminopeptidases is defined by seven conserved residues. Five of them are involved in metal ion binding which is important to maintain both the activity and the oligomeric state. The sixth conserved residue, a glutamate, is the catalytic base deprotonating the water molecule during peptide bond hydrolysis. The seventh residue is an aspartate whose function remains poorly understood. This aspartate residue, however, must have a critical role as it is strictly conserved in all MH clan enzymes. It forms some kind of catalytic triad with the histidine residue and the metal ion of the M2 binding site. We assess its role in TmPep1050, an M42 aminopeptidase of Thermotoga maritima, through a mutational approach. Asp-62 was substituted with alanine, asparagine, or glutamate residue. The Asp-62 substitutions completely abolished TmPep1050 activity and impeded dodecamer formation. They also interfered with metal ion binding as only one cobalt ion is bound per subunit instead of two. The structure of Asp62Ala variant was solved at 1.5 Å showing how the substitution has an impact on the active site fold. We propose a structural role for Asp-62, helping to stabilize a crucial loop in the active site and to position correctly the catalytic base and a metal ion ligand of the M1 site.     

Metal complexes as artificial proteases: toward catalytic drugs

Catalytic cleavage of the backbone of a protein related to a disease may cure the disease. Owing to the catalytic nature of the protein inactivation, the drug dosage and the side effects can be reduced with the catalytic drugs. Catalytic drugs can be designed even for proteins lacking active sites. Effective artificial proteases have been designed for proteins or oligomers of oligopeptides. The Co(III) complex of cyclen has been used as the catalytic center for peptide hydrolysis. Binding sites of the catalysts that recognize the targets have been searched by using various kinds of chemical libraries. Some of the artificial metalloproteases reported till date offer a new therapeutic option for amyloidoses (e.g. Alzheimer's disease, type 2 diabetes mellitus, Parkinson's disease, mad cow's disease, etc.).     

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.     

The PD...(D/E)XK motif in restriction enzymes: a link between function and conformation

The active sites of Mg(II)-dependent nucleases feature a cluster of conserved charged residues which includes both acidic (Asp and Glu) and basic (Lys) side chains. In restriction enzymes, these side chains are part of the conserved PD...(D/E)XK functional sequence motif which has been implicated as being important in metal ion binding and catalytic steps. Recent work revealing the unusual behavior of the active site variant D58A of the representative PvuII endonuclease prompted speculation that the array of charged groups in the nuclease active site may also be linked to conformational behavior [Dupureur, C. M., and Conlan, L. H. (2000) Biochemistry 39, 10921-10927]. To address this issue, we analyzed the conformational behavior of active site variants of PvuII endonuclease using both NMR spectroscopic and thermodynamic methods. NMR spectroscopic analysis via (19)F and (1)H-(15)N HSQC experiments indicates that a number of side chain and backbone amide groups are perturbed upon Ala substitution at conserved active site residues Asp58, Glu68, and Lys70. Spectral changes are particularly pronounced for the lowest-activity mutants (D58A and K70A). These changes are accompanied by perturbations in conformational stability. Ala substitution at each of these positions results in 2-5 kcal/mol of stabilization over the wild-type enzyme at pH 7.7, changes which constitute increases in DeltaG(d)(H2O) of 20-50%. The pH dependencies of mutant enzyme stabilities are distinct from those of the wild type, results which confirm that these ionizable groups strongly influence stability. Wild-type enzyme stability is correlated with the ionization of groups shown to be important to metal ion binding and orientation. Correlations between spectral changes and conformational stability indicate that the latter measurements may prove useful in the evaluation of site-directed mutant restriction enzymes. More importantly, these results indicate that structure-function relationships in restriction enzyme active sites can be complex, and that the ensemble of conserved charged residues which mediate DNA hydrolysis in Mg(II)-dependent nucleases constitutes a critical link between function and conformation.     

The unswapped chain of bovine seminal ribonuclease: Crystal structure of the free and liganded monomeric derivative

Bovine seminal ribonuclease, a homodimeric enzyme joined covalently by two interchain disulphide bonds, is an equilibrium mixture of two conformational isomers, MxM and M=M. The major form, MxM, whose crystal structure has been previously determined at 1.9 A resolution, presents the swapping of the N-terminal segments (residues 1-15) and composite active sites formed by residues of different chains. The three-dimensional domain swapping does not occur in the M=M form. The different fold of each N-terminal tail is directed by the hinge loop (residue 16-22) connecting the swapping domain to the body of the protein. Reduction and alkylation of interchain disulphide bridges produce a monomeric derivative and a noncovalent swapped dimer, which are both active. The free and nucleotide-bound forms of the monomer have been crystallized at an alkaline pH and refined at 1.45 and 1.65 A resolution, respectively. In both cases, the N-terminal fragment is folded on the main body of the protein to produce an intact active site and a chain architecture very similar to that of bovine pancreatic ribonuclease. In this new fold of the seminal chain, the hinge loop is disordered. Despite the difference between the tertiary structure of the monomer and that of the chains in the MxM form, the active sites of the two enzymes are virtually indistinguishable. Furthermore, the structure of the liganded enzyme represents the first example of a ribonuclease complex studied at an alkaline pH and provides new information on the binding of a nucleotide when the catalytic histidines are deprotonated.     

Comparison of designed and randomly generated catalysts for simple chemical reactions

There has been recent success in designing enzymes for simple chemical reactions using a two-step protocol. In the first step, a geometric matching algorithm is used to identify naturally occurring protein scaffolds at which predefined idealized active sites can be realized. In the second step, the residues surrounding the transition state model are optimized to increase transition state binding affinity and to bolster the primary catalytic side chains. To improve the design methodology, we investigated how the set of solutions identified by the design calculations relate to the overall set of solutions for two different chemical reactions. Using a TIM barrel scaffold in which catalytically active Kemp eliminase and retroaldolase designs were obtained previously, we carried out activity screens of random libraries made to be compositionally similar to active designs. A small number of active catalysts were found in screens of 10(3) variants for each of the two reactions, which differ from the computational designs in that they reuse charged residues already present in the native scaffold. The results suggest that computational design considerably increases the frequency of catalyst generation for active sites involving newly introduced catalytic residues, highlighting the importance of interaction cooperativity in enzyme active sites.     

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.