Enzyme/non-enzyme discrimination and prediction of enzyme active site location using charge-based methods

Calculations of charge interactions complement analysis of a characterised active site, rationalising pH-dependence of activity and transition state stabilisation. Prediction of active site location through large DeltapK(a)s or electrostatic strain is relevant for structural genomics. We report a study of ionisable groups in a set of 20 enzymes, finding that false positives obscure predictive potential. In a larger set of 156 enzymes, peaks in solvent-space electrostatic properties are calculated. Both electric field and potential match well to active site location. The best correlation is found with electrostatic potential calculated from uniform charge density over enzyme volume, rather than from assignment of a standard atom-specific charge set. Studying a shell around each molecule, for 77% of enzymes the potential peak is within that 5% of the shell closest to the active site centre, and 86% within 10%. Active site identification by largest cleft, also with projection onto a shell, gives 58% of enzymes for which the centre of the largest cleft lies within 5% of the active site, and 70% within 10%. Dielectric boundary conditions emphasise clefts in the uniform charge density method, which is suited to recognition of binding pockets embedded within larger clefts. The variation of peak potential with distance from active site, and comparison between enzyme and non-enzyme sets, gives an optimal threshold distinguishing enzyme from non-enzyme. We find that 87% of the enzyme set exceeds the threshold as compared to 29% of the non-enzyme set. Enzyme/non-enzyme homologues, "structural genomics" annotated proteins and catalytic/non-catalytic RNAs are studied in this context.     

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large‐scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post‐translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus‐responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.     

Dbmodeling

Genome sequencing efforts are providing us with complete genetic blueprints for hundreds of organisms. We are now faced with assigning, understanding, and modifying the functions of proteins encoded by these genomes. DBMODELING is a relational database of annotated comparative protein structure models and their metabolic pathway characterization, when identified. This procedure was applied to complete genomes such as Mycobacterium tuberculosis and Xylella fastidiosa. The main interest in the study of metabolic pathways is that some of these pathways are not present in humans, which makes them selective targets for drug design, decreasing the impact of drugs in humans. In the database, there are currently 1116 proteins from two genomes. It can be accessed by any researcher at http://www. biocristalografia.df.ibilce. unesp.br/tools/. This project confirms that homology modeling is a useful tool in structural bioinformatics and that it can be very valuable in annotating genome sequence information, contributing to structural and functional genomics, and analyzing protein-ligand docking.     

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.     

Computational active site analysis of molecular pathways to improve functional classification of enzymes

This study describes a method to computationally assess the function of homologous enzymes through small molecule binding interaction energy. Three experimentally determined X-ray structures and four enzyme models from ornithine cyclo-deaminase, alanine dehydrogenase, and mu-crystallin were used in combination with nine small molecules to derive a function score (FS) for each enzyme-model combination. While energy values varied for a single molecule-enzyme combination due to differences in the active sites, we observe that the binding energies for the entire pathway were proportional for each set of small molecules investigated. This proportionality of energies for a reaction pathway appears to be dependent on the amino acids in the active site and their direct interactions with the small molecules, which allows a function score (FS) to be calculated to assess the specificity of each enzyme. Potential of mean force (PMF) calculations were used to obtain the energies, and the resulting FS values demonstrate that a measurement of function may be obtained using differences between these PMF values. Additionally, limitations of this method are discussed based on: (a) larger substrates with significant conformational flexibility; (b) low homology enzymes; and (c) open active sites. This method should be useful in accurately predicting specificity for single enzymes that have multiple steps in their reactions and in high throughput computational methods to accurately annotate uncharacterized proteins based on active site interaction analysis.     

Structural genomics of proteins from conserved biochemical pathways and processes

During the past year, X-ray crystallographers and solution NMR spectroscopists have made significant progress towards the complete structural characterization of conserved biochemical pathways and processes. Some of these advances were made in the context of nascent structural genomics programs, which promise to accelerate structural studies of biologically and medically important proteins. The results of high-throughput protein production, crystallization, structure determination, homology modeling and functional annotation published by two such programs have provided insight into the evolution and function of enzymes in the isoprenoid biosynthesis and ribulose monophosphate pathways.     

Structural and functional analysis of tetracenomycin F2 cyclase from Streptomyces glaucescens. A type II polyketide cyclase

Tetracenomycin F2 cyclase (tcmI gene product), catalyzes an aromatic rearrangement in the biosynthetic pathway for tetracenomycin C in Streptomyces glaucescens. The x-ray structure of this small enzyme has been determined to 1.9-A resolution together with an analysis of site-directed mutants of potential catalytic residues. The protein exhibits a dimeric betaalphabeta ferredoxin-like fold that utilizes strand swapping between subunits in its assembly. The fold is dominated by four strands of antiparallel sheet and a layer of alpha-helices, which creates a cavity that is proposed to be the active site. This type of secondary structural arrangement has been previously observed in polyketide monooxygenases and suggests an evolutionary relationship between enzymes that catalyze adjacent steps in these biosynthetic pathways. Mutational analysis of all of the obvious catalytic bases within the active site suggests that the enzyme functions to steer the chemical outcome of the cyclization rather than providing a specific catalytic group. Together, the structure and functional analysis provide insight into the structural framework necessary to perform the complex rearrangements catalyzed by this class of polyketide cyclases.     

The emerging structural understanding of transglutaminase 3

Transglutaminases (TGase; protein-glutamine: amine gamma-glutamyl-transferase) are a family of calcium-dependent acyl-transfer enzymes ubiquitously expressed in mammalian cells and responsible for catalyzing covalent cross-links between proteins or peptides. A series of recent crystal structures have revealed the overall architecture of TGase enzymes, and provided a deep look at their active site, calcium and magnesium ions, and the manner by which guanine nucleotides interact with this enzyme. These structures, backed with extensive biochemical studies, are providing new insights as to how access to the enzyme's active site may be gated through the coordinated changes in cellular calcium and magnesium concentrations and GTP/GDP. Calcium-activated TGase 3 can bind, hydrolyze, and is inhibited by GTP, despite lacking structural homology with other GTP binding proteins. A structure based sequence homology among the TGase enzyme family shows that these essential structural features are shared among other members of the TGase family.     

Structural basis of kinetic variations in Sucrose Phosphate Phosphatase (SPP) of rice and Anabaena unveiled through computational analysis

Sucrose is an important storage form of assimilated carbon in many plant species. Unlike other sucrose biosynthetic enzymes, Sucrose Phosphate Phosphatase (SPP), the terminal enzyme in sucrose biosynthetic pathway, is the least understood. SPPs from different organisms have different kinetic properties. The current study focuses on the structural differences among SPP homologues and unveils the probable structural basis of kinetic variations. We have employed computational methods of molecular modeling and structure comparisons and identified structural variations in some of the substrate binding residues, amino acid substitutions in regions that are lining the active site and minute structural differences that can enhance the nucleophilicity of a catalytic nucleophile (Asp 9 ). We report a structurally and hence functionally important amino acid substitution (Asp 159 by Alanine) in one of the rice SPP isoforms, which can result in the disruption of a H-bond that helps in binding of sucrose at the active site of the enzyme. In this paper we discuss the structural basis of enhanced catalytic efficiency of rice SPP in comparison with a cyanobacterium (Anabaena variabilis). The natural mutations identified in our analysis of the SPP catalytic domain would be useful in re-designing the enzyme for enhanced catalytic efficiency and higher sucrose production.     

In silico quest for putative drug targets in <Emphasis Type="Italic">Helicobacter pylori</Emphasis> HPAG1: molecular modeling of candidate enzymes from lipopolysaccharide biosynthesis pathway

Aimed at identification and structural characterization of novel putative therapeutic targets in H. pylori, the etiological agent of numerous gastrointestinal diseases including peptic ulcer and gastric cancer, the present study comprised of three phases. First, through subtractive analysis of metabolic pathways of Helicobacter pylori HPAG1 and human, as documented in the KEGG database, 11 pathogen-specific pathways were identified. Next, all proteins involved in these pathogen-specific pathways were scrutinized in search of promising targets and the study yielded 25 candidate target proteins that are likely to be essential for the pathogen viability, but have no homolog in human. The lipopolysaccharide (LPS) biosynthesis pathway was found to be the largest contributor (nine proteins) to this list of candidate proteins. Considering the importance of LPS in H. pylori virulence, 3D structural models of three predicted target enzymes of this pathway, namely 2-dehydro-3-deoxy-phosphooctonate aldolase, UDP-3-O-[3-hydroxymyristoyl] N-acetylglucosamine deacetylase and Phosphoheptose isomerase, were then built up using the homology modeling approaches. Binding site analysis and docking of the known biological substrate PEP to 2-dehydro-3-deoxyphosphooctonate aldolase revealed the potential binding pocket present in the single monomeric form of the enzyme and identified 11 amino acid residues that might play the key roles in this protein-ligand interaction.     

Broad‐substrate screen as a tool to identify substrates for bacterial Gcn5‐related N‐acetyltransferases with unknown substrate specificity

Due to a combination of efforts from individual laboratories and structural genomics centers, there has been a surge in the number of members of the Gcn5‐related acetyltransferasesuperfamily that have been structurally determined within the past decade. Although the number of three‐dimensional structures is increasing steadily, we know little about the individual functions of these enzymes. Part of the difficulty in assigning functions for members of this superfamily is the lack of information regarding how substrates bind to the active site of the protein. The majority of the structures do not show ligand bound in the active site, and since the substrate‐binding domain is not strictly conserved, it is difficult to predict the function based on structure alone. Additionally, the enzymes are capable of acetylating a wide variety of metabolites and many may exhibit promiscuity regarding their ability to acetylate multiple classes of substrates, possibly having multiple functions for the same enzyme. Herein, we present an approach to identify potential substrates for previously uncharacterized members of the Gcn5‐related acetyltransferase superfamily using a variety of metabolites including polyamines, amino acids, antibiotics, peptides, vitamins, catecholamines, and other metabolites. We have identified potential substrates for eight bacterial enzymes of this superfamily. This information will be used to further structurally and functionally characterize them.     

Non-homologous isofunctional enzymes: A systematic analysis of alternative solutions in enzyme evolution

Background  

Evolutionarily unrelated proteins that catalyze the same biochemical reactions are often referred to as analogous - as opposed to homologous - enzymes. The existence of numerous alternative, non-homologous enzyme isoforms presents an interesting evolutionary problem; it also complicates genome-based reconstruction of the metabolic pathways in a variety of organisms. In 1998, a systematic search for analogous enzymes resulted in the identification of 105 Enzyme Commission (EC) numbers that included two or more proteins without detectable sequence similarity to each other, including 34 EC nodes where proteins were known (or predicted) to have distinct structural folds, indicating independent evolutionary origins. In the past 12 years, many putative non-homologous isofunctional enzymes were identified in newly sequenced genomes. In addition, efforts in structural genomics resulted in a vastly improved structural coverage of proteomes, providing for definitive assessment of (non)homologous relationships between proteins.     

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.     

Rat submaxillary gland serine protease, tonin. Structure solution and refinement at 1.8 A resolution

Tonin is a mammalian serine protease that is capable of generating the vasoconstrictive agent, angiotensin II, directly from its precursor protein, angiotensinogen, a process that normally requires two enzymes, renin and angiotensin-converting enzyme. The X-ray crystallographic structure determination and refinement of tonin at 1.8 A resolution and the analysis of the resulting model are reported. The initial phases were obtained by the method of molecular replacement using as the search model the structure of bovine trypsin. The refined model of tonin consists of 227 amino acid residues out of the 235 in the complete molecule, 149 water molecules, and one zinc ion. The R-factor (R = sigma Fo - Fc/sigma Fo) is 0.196 for the 14,997 measured data between 8 and 1.8 A resolution with I greater than or equal to sigma (I). It is estimated that the overall root-mean-square error in the coordinates is about 0.3 A. The structure of tonin that has been determined is not in its active conformation, but one that has been perturbed by the binding of Zn2+ in the active site. Zn2+ was included in the buffer to aid the crystallization. Nevertheless, the structure of tonin that is described is for the most part similar to its native form as indicated by the close tertiary structural homology with kallikrein. The differences in the structures of the two enzymes are concentrated in several loop regions; these structural differences are probably responsible for the differences in their reactivities and specificities.     

High-throughput screening for potent and selective inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase

Plasmodium falciparum is the causative agent of the most serious and fatal malarial infections, and it has developed resistance to commonly employed chemotherapeutics. The de novo pyrimidine biosynthesis enzymes offer potential as targets for drug design, because, unlike the host, the parasite does not have pyrimidine salvage pathways. Dihydroorotate dehydrogenase (DHODH) is a flavin-dependent mitochondrial enzyme that catalyzes the fourth reaction in this essential pathway. Coenzyme Q (CoQ) is utilized as the oxidant. Potent and species-selective inhibitors of malarial DHODH were identified by high-throughput screening of a chemical library, which contained 220,000 drug-like molecules. These novel inhibitors represent a diverse range of chemical scaffolds, including a series of halogenated phenyl benzamide/naphthamides and urea-based compounds containing napthyl or quinolinyl substituents. Inhibitors in these classes with IC(50) values below 600 nm were purified by high pressure liquid chromatography, characterized by mass spectroscopy, and subjected to kinetic analysis against the parasite and human enzymes. The most active compound is a competitive inhibitor of CoQ with an IC(50) against malarial DHODH of 16 nm, and it is 12,500-fold less active against the human enzyme. Site-directed mutagenesis of residues in the CoQ-binding site significantly reduced inhibitor potency. The structural basis for the species selective enzyme inhibition is explained by the variable amino acid sequence in this binding site, making DHODH a particularly strong candidate for the development of new anti-malarial compounds.     

Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: defining the entry point into the novel benzoate oxidation (box) pathway

Xenobiotic aromatic compounds represent one of the most significant classes of environmental pollutants. A novel benzoate oxidation (box) pathway has been identified recently in Burkholderia xenovorans LB400 (referred to simply as LB400) that is capable of assimilating benzoate and intimately tied to the degradation of polychlorinated biphenyls (PCBs). The box pathway in LB400 is present in two paralogous copies (boxM and boxC) and encodes eight enzymes with the first committed step catalyzed by benzoate CoA ligase (BCL). As a first step towards delineating the biochemical role of the box pathway in LB400, we have carried out functional studies of the paralogous BCL enzymes (BCLM and BCLC) with 20 different putative substrates. We have established a structural rationale for the observed substrate specificities on the basis of a 1.84 A crystal structure of BCLM in complex with benzoate. These data show that, while BCLM and BCLC display similar overall substrate specificities, BCLM is significantly more active towards benzoate and 2-aminobenzoate with tighter binding (Km) and a faster reaction rate (Vmax). Despite these clear functional differences, the residues that define the substrate-binding site in BCLM are completely conserved in BCLC, suggesting that second shell residues may play a significant role in substrate recognition and catalysis. Furthermore, comparison of the active site of BCLM with the recently solved structures of 4-chlorobenzoate CoA ligase and 2, 3-dihydroxybenzoate CoA ligase offers additional insight into the molecular features that mediate substrate binding in adenylate-forming enzymes. This study provides the first biochemical characterization of a Box enzyme from LB400 and the first structural characterization of a Box enzyme from any organism, and further substantiates the concept of distinct roles for the two paralogous box pathways in LB400.     

Structural analysis of two enzymes catalysing reverse metabolic reactions implies common ancestry

The crystal structure of the dimeric anthranilate phosphoribosyltransferase (AnPRT) reveals a new category of phosphoribosyltransferases, designated as class III. The active site of this enzyme is located within the flexible hinge region of its two-domain structure. The pyrophosphate moiety of phosphoribosylpyrophosphate is co-ordinated by a metal ion and is bound by two conserved loop regions within this hinge region. With the structure of AnPRT available, structural analysis of all enzymatic activities of the tryptophan biosynthesis pathway is complete, thereby connecting the evolution of its enzyme members to the general development of metabolic processes. Its structure reveals it to have the same fold, topology, active site location and type of association as class II nucleoside phosphorylases. At the level of sequences, this relationship is mirrored by 13 structurally invariant residues common to both enzyme families. Taken together, these data imply common ancestry of enzymes catalysing reverse biological processes--the ribosylation and deribosylation of metabolic pathway intermediates. These relationships establish new links for enzymes involved in nucleotide and amino acid metabolism.     

Complete amino acid sequence of the medium-chain S-acyl fatty acid synthetase thio ester hydrolase from rat mammary gland

The complete amino acid sequence of the medium-chain S-acyl fatty acid synthetase thio ester hydrolase (thioesterase II) from rat mammary gland is presented. Most of the sequence was derived by analysis of peptide fragments produced by cleavage at methionyl, glutamyl, lysyl, arginyl, and tryptophanyl residues. A small section of the sequence was deduced from a previously analyzed cDNA clone. The protein consists of 260 residues and has a blocked amino-terminal methionine and calculated Mr of 29,212. The carboxy-terminal sequence, verified by Edman degradation of the carboxy-terminal cyanogen bromide fragment and carboxypeptidase Y digestion of the intact thioesterase II, terminates with a serine residue and lacks three additional residues predicted by the cDNA sequence. The native enzyme contains three cysteine residues but no disulfide bridges. The active site serine residue is located at position 101. The rat mammary gland thioesterase II exhibits approximately 40% homology with a thioesterase from mallard uropygial gland, the sequence of which was recently determined by cDNA analysis [Poulose, A.J., Rogers, L., Cheesbrough, T. M., & Kolattukudy, P. E. (1985) J. Biol. Chem. 260, 15953-15958]. Thus the two enzymes may share similar structural features and a common evolutionary origin. The location of the active site in these thioesterases differs from that of other serine active site esterases; indeed, the enzymes do not exhibit any significant homology with other serine esterases, suggesting that they may constitute a separate new family of serine active site enzymes.     

Molecular cloning and amino acid sequence of rat enkephalinase

cDNA clones encoding rat enkephalinase (neutral endopeptidase, EC 3.4.24.11) have been isolated in lambda gt10 libraries from both brain and kidney mRNAs and the complete 742 amino acid sequence of rat enkephalinase is presented. The enzyme possesses a single transmembrane spanning domain near the N-terminal of the molecule but lacks a signal sequence. Because enkephalinase has it active site located extracellularly and is thus an ectopeptidase, we suggest that the N-terminal transmembrane region of the enzyme anchors the protein in membranes and that the majority of the protein, including the carboxy terminus, is extracellular. Enkephalinase, a zinc-containing metallo enzyme, displays homology with other zinc metallo enzymes such as carboxypeptidase A, B and E, suggesting enzymatic similarities in these enzymes.