Redesigning the PheA domain of gramicidin synthetase leads to a new understanding of the enzyme's mechanism and selectivity

The PheA domain of gramicidin synthetase A, a non-ribosomal peptide synthetase, selectively binds phenylalanine along with ATP and Mg2+ and catalyzes the formation of an aminoacyl adenylate. In this study, we have used a novel protein redesign algorithm, K*, to predict mutations in PheA that should exhibit improved binding for tyrosine. Interestingly, the introduction of two predicted mutations to PheA did not significantly improve KD, as measured by equilibrium fluorescence quenching. However, the mutations improved the specificity of the enzyme for tyrosine (as measured by kcat/KM), primarily driven by a 56-fold improvement in KM, although the improvement did not make tyrosine the preferred substrate over phenylalanine. Using stopped-flow fluorometry, we examined binding of different amino acid substrates to the wild-type and mutant enzymes in the pre-steady state in order to understand the improvement in KM. Through these investigations, it became evident that substrate binding to the wild-type enzyme is more complex than previously described. These experiments show that the wild-type enzyme binds phenylalanine in a kinetically selective manner; no other amino acids tested appeared to bind the enzyme in the early time frame examined (500 ms). Furthermore, experiments with PheA, phenylalanine, and ATP reveal a two-step binding process, suggesting that the PheA-ATP-phenylalanine complex may undergo a conformational change toward a catalytically relevant intermediate on the pathway to adenylation; experiments with PheA, phenylalanine, and other nucleotides exhibit only a one-step binding process. The improvement in KM for the mutant enzyme toward tyrosine, as predicted by K*, may indicate that redesigning the side-chain binding pocket allows the substrate backbone to adopt productive conformations for catalysis but that further improvements may be afforded by modeling an enzyme:ATP:substrate complex, which is capable of undergoing conformational change.     

Progress in computational protein design

Current progress in computational structure-based protein design is reviewed in the areas of methodology and applications. Foundational advances include new potential functions, more efficient ways of computing energetics, flexible treatments of solvent, and useful energy function approximations, as well as ensemble-based approaches to scoring designs for inclusion of entropic effects, improvements to guaranteed and to stochastic search techniques, and methods to design combinatorial libraries for screening and selection. Applications include new approaches and successes in the design of specificity for protein folding, binding, and catalysis, in the redesign of proteins for enhanced binding affinity, and in the application of design technology to study and alter enzyme catalysis. Computational protein design continues to mature and advance.     

A viable amino acid editing activity in the leucyl-tRNA synthetase CP1-splicing domain is not required in the yeast mitochondria

Aminoacyl-tRNA synthetases are a family of enzymes that are responsible for translating the genetic code in the first step of protein synthesis. Some aminoacyl-tRNA synthetases have editing activities to clear their mistakes and enhance fidelity. Leucyl-tRNA synthetases have a hydrolytic active site that resides in a discrete amino acid editing domain called CP1. Mutational analysis within yeast mitochondrial leucyl-tRNA synthetase showed that the enzyme has maintained an editing active site that is competent for post-transfer editing of mischarged tRNA similar to other leucyl-tRNA synthetases. These mutations that altered or abolished leucyl-tRNA synthetase editing were introduced into complementation assays. Cell viability and mitochondrial function were largely unaffected in the presence of high levels of non-leucine amino acids. In contrast, these editing-defective mutations limited cell viability in Escherichia coli. It is possible that the yeast mitochondria have evolved to tolerate lower levels of fidelity in protein synthesis or have developed alternate mechanisms to enhance discrimination of leucine from non-cognate amino acids that can be misactivated by leucyl-tRNA synthetase.     

Computational design of Candida boidinii xylose reductase for altered cofactor specificity

In this study we introduce a computationally‐driven enzyme redesign workflow for altering cofactor specificity from NADPH to NADH. By compiling and comparing data from previous studies involving cofactor switching mutations, we show that their effect cannot be explained as straightforward changes in volume, hydrophobicity, charge, or BLOSUM62 scores of the residues populating the cofactor binding site. Instead, we find that the use of a detailed cofactor binding energy approximation is needed to adequately capture the relative affinity towards different cofactors. The implicit solvation models Generalized Born with molecular volume integration and Generalized Born with simple switching were integrated in the iterative protein redesign and optimization (IPRO) framework to drive the redesign of Candida boidinii xylose reductase (CbXR) to function using the non‐native cofactor NADH. We identified 10 variants, out of the 8,000 possible combinations of mutations, that improve the computationally assessed binding affinity for NADH by introducing mutations in the CbXR binding pocket. Experimental testing revealed that seven out of ten possessed significant xylose reductase activity utilizing NADH, with the best experimental design (CbXR‐GGD) being 27‐fold more active on NADH. The NADPH‐dependent activity for eight out of ten predicted designs was either completely abolished or significantly diminished by at least 90%, yielding a greater than 10⁴‐fold change in specificity to NADH (CbXR‐REG). The remaining two variants (CbXR‐RTT and CBXR‐EQR) had dual cofactor specificity for both nicotinamide cofactors.     

Steric and thermodynamic limits of design for the incorporation of large unnatural amino acids in aminoacyl‐tRNA synthetase enzymes

Orthogonal aminoacyl‐tRNA synthetase/tRNA pairs from archaea have been evolved to facilitate site specific in vivo incorporation of unnatural amino acids into proteins in Escherichia coli. Using this approach, unnatural amino acids have been successfully incorporated with high translational efficiency and fidelity. In this study, CHARMM‐based molecular docking and free energy calculations were used to evaluate rational design of specific protein–ligand interactions for aminoacyl‐tRNA synthetases. A series of novel unnatural amino acid ligands were docked into the p‐benzoyl‐L ‐phenylalanine tRNA synthetase, which revealed that the binding pocket of the enzyme does not provide sufficient space for significantly larger ligands. Specific binding site residues were mutated to alanine to create additional space to accommodate larger target ligands, and then mutations were introduced to improve binding free energy. This approach was used to redesign binding sites for several different target ligands, which were then tested against the standard 20 amino acids to verify target specificity. Only the synthetase designed to bind Man‐α‐O‐Tyr was predicted to be sufficiently selective for the target ligand and also thermodynamically stable. Our study suggests that extensive redesign of the tRNA synthatase binding pocket for large bulky ligands may be quite thermodynamically unfavorable. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

The 2 A crystal structure of leucyl-tRNA synthetase and its complex with a leucyl-adenylate analogue

Leucyl-, isoleucyl- and valyl-tRNA synthetases are closely related large monomeric class I synthetases. Each contains a homologous insertion domain of approximately 200 residues, which is thought to permit them to hydrolyse ('edit') cognate tRNA that has been mischarged with a chemically similar but non-cognate amino acid. We describe the first crystal structure of a leucyl-tRNA synthetase, from the hyperthermophile Thermus thermophilus, at 2.0 A resolution. The overall architecture is similar to that of isoleucyl-tRNA synthetase, except that the putative editing domain is inserted at a different position in the primary structure. This feature is unique to prokaryote-like leucyl-tRNA synthetases, as is the presence of a novel additional flexibly inserted domain. Comparison of native enzyme and complexes with leucine and a leucyl- adenylate analogue shows that binding of the adenosine moiety of leucyl-adenylate causes significant conformational changes in the active site required for amino acid activation and tight binding of the adenylate. These changes are propagated to more distant regions of the enzyme, leading to a significantly more ordered structure ready for the subsequent aminoacylation and/or editing steps.     

C1-Tetrahydrofolate synthase from rabbit liver. Structural and kinetic properties of the enzyme and its two domains

C1-Tetrahydrofolate synthase is a multifunctional enzyme which catalyzes three reactions in 1-carbon metabolism: 10-formyltetrahydrofolate synthetase; 5,10-methenyltetrahydrofolate cyclohydrolase; 5,10-methylenetetrahydrofolate dehydrogenase. A rapid 1-day purification procedure has been developed which gives 40 mg of pure enzyme from 10 rabbit livers. The 10-formyltetrahydrofolate synthetase activity of this trifunctional enzyme has a specific activity that is 4-fold higher than the enzyme previously purified from rabbit liver. Conditions have been developed for the rapid isolation of a tryptic fragment of the enzyme which contains the methylenetetrahydrofolate dehydrogenase and methenyltetrahydrofolate cyclohydrolase activities. This fragment is a monomer exhibiting a subunit and native molecular weight of 36,000 in most buffers. However, in phosphate buffers the native molecular weight suggests that the fragment is a dimer. Conditions are also given whereby chymotryptic digestion allows the simultaneous isolation from the native enzyme of a large fragment containing the 10-formyltetrahydrofolate synthetase activity and a smaller fragment containing the dehydrogenase and cyclohydrolase activities. The large fragment is a dimer with a subunit molecular weight of 66,000. The small fragment retains all of the dehydrogenase and cyclohydrolase activities of the native enzyme. The large fragment is unstable but retains most of the 10-formyltetrahydrofolate synthetase activity. Km values of substrates for the two fragments are the same as the values for the native enzyme. The 10-formyltetrahydrofolate synthetase activity of the native enzyme requires ammonium or potassium ions for expression of full catalytic activity. The effect of these two ions on the catalytic activity of the large chymotryptic fragment is the same as with the native enzyme. We have shown by differential scanning calorimetry that the native enzyme contains two protein domains which show thermal transitions at 47 and 60 degrees C. Evidence is presented that the two domains are related to the two protein fragments generated by proteolysis of the native enzyme. The larger of the two domains contains the active site for the 10-formyltetrahydrofolate synthetase activity while the smaller domain contains the active site which catalyzes the dehydrogenase and cyclohydrolase reactions. Replacement of sodium ion buffers with either ammonium or potassium ions results in an increase in stability of the large domain of the native enzyme. This change in stability is not accompanied by a change in the quaternary structure of the enzyme.(ABSTRACT TRUNCATED AT 400 WORDS)     

Systematic optimization model and algorithm for binding sequence selection in computational enzyme design

A systematic optimization model for binding sequence selection in computational enzyme design was developed based on the transition state theory of enzyme catalysis and graph‐theoretical modeling. The saddle point on the free energy surface of the reaction system was represented by catalytic geometrical constraints, and the binding energy between the active site and transition state was minimized to reduce the activation energy barrier. The resulting hyperscale combinatorial optimization problem was tackled using a novel heuristic global optimization algorithm, which was inspired and tested by the protein core sequence selection problem. The sequence recapitulation tests on native active sites for two enzyme catalyzed hydrolytic reactions were applied to evaluate the predictive power of the design methodology. The results of the calculation show that most of the native binding sites can be successfully identified if the catalytic geometrical constraints and the structural motifs of the substrate are taken into account. Reliably predicting active site sequences may have significant implications for the creation of novel enzymes that are capable of catalyzing targeted chemical reactions.     

Probing electrostatic interactions and ligand binding in aspartyl-tRNA synthetase through site-directed mutagenesis and computer simulations

Faithful genetic code translation requires that each aminoacyl-tRNA synthetase recognise its cognate amino acid ligand specifically. Aspartyl-tRNA synthetase (AspRS) distinguishes between its negatively-charged Asp substrate and two competitors, neutral Asn and di-negative succinate, using a complex network of electrostatic interactions. Here, we used molecular dynamics simulations and site-directed mutagenesis experiments to probe these interactions further. We attempt to decrease the Asp/Asn binding free energy difference via single, double and triple mutations that reduce the net positive charge in the active site of Escherichia coli AspRS. Earlier, Glutamine 199 was changed to a negatively-charged glutamate, giving a computed reduction in Asp affinity in good agreement with experiment. Here, Lysine 198 was changed to a neutral leucine; then, Lys198 and Gln199 were mutated simultaneously. Both mutants are predicted to have reduced Asp binding and improved Asn binding, but the changes are insufficient to overcome the initial, high specificity of the native enzyme, which retains a preference for Asp. Probing the aminoacyl-adenylation reaction through pyrophosphate exchange experiments, we found no detectable activity for the mutant enzymes, indicating weaker Asp binding and/or poorer transition state stabilization. The simulations show that the mutations' effect is partly offset by proton uptake by a nearby histidine. Therefore, we performed additional simulations where the nearby Histidines 448 and 449 were mutated to neutral or negative residues: (Lys198Leu, His448Gln, His449Gln), and (Lys198Leu, His448Glu, His449Gln). This led to unexpected conformational changes and loss of active site preorganization, suggesting that the AspRS active site has a limited structural tolerance for electrostatic modifications. The data give insights into the complex electrostatic network in the AspRS active site and illustrate the difficulty in engineering charged-to-neutral changes of the preferred ligand.     

Design,activity, and structure of a highly specific artificial endonuclease

We have generated an artificial highly specific endonuclease by fusing domains of homing endonucleases I-DmoI and I-CreI and creating a new 1400 A(2) protein interface between these domains. Protein engineering was accomplished by combining computational redesign and an in vivo protein-folding screen. The resulting enzyme, E-DreI (Engineered I-DmoI/I-CreI), binds a long chimeric DNA target site with nanomolar affinity, cleaving it precisely at a rate equivalent to its natural parents. The structure of an E-DreI/DNA complex demonstrates the accuracy of the protein interface redesign algorithm and reveals how catalytic function is maintained during the creation of the new endonuclease. These results indicate that it may be possible to generate novel highly specific DNA binding proteins from homing endonucleases.     

Cysteinyl-tRNA synthetase from Bacillus stearothermophilus. A structural and functional monomer.

A procedure is described for the purification of cysteinyl-tRNA synthetase as a side product of a multi-enzyme isolation from Bacillus stearothermophilus. The native and denatured enzyme are both shown to have a molecular weight of 54000 by gel filtration and sodium dodecyl sulphate/polyacrylamide gel electrophoresis respectively. Fingerprinting and peptide counting indicate that the polypeptide chain has a nonrepeating primary structure. The enzyme has only one binding site for each of its substrates (cysteine, ATP and tRNACys) as judged by equilibrium dialysis, active-site titration and fluorescence quenching. No evidence for the dimerisation of the enzyme in the presence of these substrates could be found. We conclude that cysteinyl-tRNA synthetase, which is the smallest aminoacyl-tRNA synthetase yet described, is both structurally and functionally monomeric.     

Effect of alanine-293 replacement on the activity, ATP binding, and editing of Escherichia coli leucyl-tRNA synthetase

Leucyl-tRNA synthetase (LeuRS) is a class I aminoacyl-tRNA synthetase that catalyzes leucylation of tRNA(Leu). Several mutants in the CP1 domain of Escherichia coli LeuRS were obtained by introduction of restriction endonuclease sites into its gene, leuS. Of these mutants, only LeuRS-A293F had decreased activity (46%) compared to the native enzyme. To investigate the effect of A293 on enzyme function, A293 was mutated to Y, G, I, R, or D. The mutants were impaired in activity and editing function to varying extents. The decrease in K(m) values for three substrates showed that the binding of ATP to these mutants became much stronger. The inhibition of ATP binding to most of the mutants was also stronger. In particular, LeuRS-A293D had the lowest activity, the strongest ATP binding, and the most impaired editing function. A red shift of the fluorescence emission maximum of LeuRS-A293D indicated a less hydrophobic chromophore environment and a relatively more flexible dynamic conformation. The change in T(m) of LeuRS-A293D was higher than that of all other substitutions. Evidence from sequence alignment and crystal structure of LeuRS from Thermus thermophilus shows that A293 was conserved as R (K) or A and is located at a small helix in the editing domain of the enzyme facing the active site. Hence, any amino acid substitution of A293 may affect the stability of the helix, which may lead to impaired editing function and aminoacylation activity and may be indirectly involved in ATP binding.     

Structural basis for the function of Bacillus subtilis phosphoribosyl-pyrophosphate synthetase

Here we report the first three-dimensional structure of a phosphoribosylpyrophosphate (PRPP) synthetase. PRPP is an essential intermediate in several biosynthetic pathways. Structures of the Bacillus subtilis PRPP synthetase in complex with analogs of the activator phosphate and the allosteric inhibitor ADP show that the functional form of the enzyme is a hexamer. The individual subunits fold into two domains, both of which resemble the type I phosphoribosyltransfereases. The active site is located between the two domains and includes residues from two subunits. Phosphate and ADP bind to the same regulatory site consisting of residues from three subunits of the hexamer. In addition to identifying residues important for binding substrates and effectors, the structures suggest a novel mode of allosteric regulation.     

Isolation and binding properties of leucyl-tRNA synthetase from Escherichia coli MRE 600

Summary A procedure for the large-scale isolation of leucyl-tRNA synthetase from E. coli MRE 600 is described: The enzyme was purified about 320-fold to homogeneity by precipitation with cetyl-trimethyl-ammonium bromide, two consecutive chromatographies on DEAE-cellulose and three on hydroxyapatite with an over-all yield of 4%.The molecular weight of leucyl-tRNA synthetase from E. coli MRE 600 was found to be 99 000 daltons. Binding studies by ultracentrifugation and equilibrium partition showed that the enzyme binds leucine, leucyl-adenylate and tRNA^Leu, each in a 1 : 1 stoichiometry. For ATP only a very weak binding to the enzyme could be observed, which did not allow the evaluation of the complex stoichiometry. The presence of ATP was not required for the binding of leucine or tRNA to leucyl-tRNA synthetase from E. coli MRE 600.     

Expression in Escherichia coli of N- and C-terminally deleted human holocarboxylase synthetase. Influence of the N-terminus on biotinylation and identification of a minimum functional protein

Biotin functions as a covalently bound cofactor of biotindependent carboxylases. Biotin attachment is catalyzed by biotin protein ligases, called holocarboxylase synthetase in mammals and BirA in prokaryotes. These enzymes show a high degree of sequence similarity in their biotinylation domains but differ markedly in the length and sequence of their N terminus. BirA is also the repressor of the biotin operon, and its DNA attachment site is located in its N terminus. The function of the eukaryotic N terminus is unknown. Holocarboxylase synthetase with N- and C-terminal deletions were evaluated for the ability to catalyze biotinylation after expression in Escherichia coli using bacterial and human acceptor substrates. We showed that the minimum functional protein is comprised of the last 349 of the 726-residue protein, which includes the biotinylation domain. Significantly, enzyme containing intermediate length, N-terminal deletions interfered with biotin transfer and interaction with different peptide acceptor substrates. We propose that the N terminus of holocarboxylase synthetase contributes to biotinylation through N- and C-terminal interactions and may affect acceptor substrate recognition. Our findings provide a rationale for the biotin responsiveness of patients with point mutations in the N-terminal sequence of holocarboxylase synthetase. Such mutant enzyme may respond to biotin-mediated stabilization of the substrate-bound complex.     

Vanadium containing bromoperoxidase - Insights into the enzymatic mechanism using X-ray crystallography

The X-ray crystal structure of the vanadium bromoperoxidase from the red algae Corallina pilulifera has been solved in the presence of the known substrates, phenol red and phloroglucinol. A putative substrate binding site has been observed in the active site channel of the enzyme. In addition bromide has been soaked into the crystals and it has been shown to bind unambiguously within the enzyme active site by using the technique of single anomalous dispersion. A specific leucine amino acid is seen to move towards the bromide ion in the wild-type enzyme to produce a hydrophobic environment within the active site. A mutant of the enzyme where arginine 397 has been changed to tryptophan, shows a different behaviour on bromide binding. These results have increased our understanding of the mechanism of the vanadium bromoperoxidases and have demonstrated that the substrate and bromide are specifically bound to the enzyme active site.     

Evidence for the importance of histidine at the active site of argininosuccinate synthetase from soybean

Studies to determine the role of histidine in catalysis by L-argininosuccinate synthetase (EC 6. 3. 4. 5) were carried out with the enzyme isolated from soybean cell suspension cultures. These experiments utilized analogues of the substrates citrulline and aspartate to investigate substrate binding, and to determine which portion of the molecule were required for binding at the active site of the enzyme. Photooxidation studies using rose bengal were carried out to define the importance of histidine residues for catalysis. These studies suggest that an active site histidine residue has an important role to play in the formation of argininosuccinate by this enzyme.     

A Pareto-Optimal Refinement Method for Protein Design Scaffolds

Computational design of protein function involves a search for amino acids with the lowest energy subject to a set of constraints specifying function. In many cases a set of natural protein backbone structures, or “scaffolds”, are searched to find regions where functional sites (an enzyme active site, ligand binding pocket, protein – protein interaction region, etc.) can be placed, and the identities of the surrounding amino acids are optimized to satisfy functional constraints. Input native protein structures almost invariably have regions that score very poorly with the design force field, and any design based on these unmodified structures may result in mutations away from the native sequence solely as a result of the energetic strain. Because the input structure is already a stable protein, it is desirable to keep the total number of mutations to a minimum and to avoid mutations resulting from poorly-scoring input structures. Here we describe a protocol using cycles of minimization with combined backbone/sidechain restraints that is Pareto-optimal with respect to RMSD to the native structure and energetic strain reduction. The protocol should be broadly useful in the preparation of scaffold libraries for functional site design.     

Crystal structures of a ligand-free and malonate-bound human caspase-1: implications for the mechanism of substrate binding

Caspase-1, a mediator of the posttranslational processing of IL-1beta and IL-18, requires an aspartic acid in the P1 position of its substrates. The mechanisms of caspase-1 activation remain poorly understood despite numerous structures of the enzyme complexed with aspartate-based inhibitors. Here we report a crystal structure of ligand-free caspase-1 that displays dramatic rearrangements of loops defining the active site to generate a closed conformation that is incompatible with substrate binding. A structure of the enzyme complexed with malonate shows the protein in its open (active-site ligand-bound) conformation in which malonate reproduces the hydrogen bonding network observed in structures with covalent inhibitors. These results illustrate the essential function of the obligatory aspartate recognition element that opens the active site of caspase-1 to substrates and may be the determinant responsible for the conformational changes between ligand-free and -bound forms of the enzyme, and suggest a new approach for identifying novel aspartic acid mimetics.     

Directed evolution of the promiscuous esterase activity of carbonic anhydrase II

A promiscuous activity of an existing enzyme can confer an evolutionary advantage by providing an immediate response to a new selection pressure and a starting point for the divergence of a new enzyme. This work seeks to examine how this process might take place. Human carbonic anhydrase II (hCAII) is an enzyme that evolved to catalyze the reversible hydration of CO(2) and performs this task at a remarkable rate (k(cat) approximately 10(6) s(-)(1)). hCAII also exhibits promiscuous activity toward highly activated esters such as 4-nitrophenyl acetate. We describe a much weaker esterase activity of hCAII toward the bulkier and much less activated ester substrate 2-naphthyl acetate (2NA). Directed evolution of hCAII produced a variant with 40-fold higher rates toward 2NA, owing to two mutations: one within the active site (Ala65Val) and one at its mouth (Thr200Ala). Structure-activity studies suggest that these mutations led to adaptation of the active site for bulkier substrates and for the catalysis of nonactivated esters. The mutations did not, however, significantly alter the native activity of hCAII. Our results support the notion that the evolution of a new function can be driven by mutations that increase a promiscuous function (which serves as the starting point for the evolutionary process) but do not harm the native function.