Computational Design of Enantiocomplementary Epoxide Hydrolases for Asymmetric Synthesis of Aliphatic and Aromatic Diols.

The use of enzymes in preparative biocatalysis often requires tailoring enzyme selectivity by protein engineering. Herein we explore the use of computational library design and molecular dynamics simulations to create variants of limonene epoxide hydrolase that produce enantiomeric diols from meso-epoxides. Three substrates of different sizes were targeted: cis-2,3-butene oxide, cyclopentene oxide, and cis-stilbene oxide. Most of the 28 designs tested were active and showed the predicted enantioselectivity. Excellent enantioselectivities were obtained for the bulky substrate cis-stilbene oxide, and enantiocomplementary mutants produced (S,S)- and (R,R)-stilbene diol with >97 % enantiomeric excess. An (R,R)-selective mutant was used to prepare (R,R)-stilbene diol with high enantiopurity (98 % conversion into diol, >99 % ee). Some variants displayed higher catalytic rates (kcat ) than the original enzyme, but in most cases KM values increased as well. The results demonstrate the feasibility of computational design and screening to engineer enantioselective epoxide hydrolase variants with very limited laboratory screening.

A robust cosolvent-compatible halohydrin dehalogenase by computational library design.

To improve the applicability of halohydrin dehalogenase as a catalyst for reactions in the presence of organic cosolvents, we explored a computational library design strategy (Framework for Rapid Enzyme Stabilization by Computational libraries) that involves discovery and in silico evaluation of stabilizing mutations. Energy calculations, disulfide bond predictions and molecular dynamics simulations identified 218 point mutations and 35 disulfide bonds with predicted stabilizing effects. Experiments confirmed 29 stabilizing point mutations, most of which were located in two distinct regions, whereas introduction of disulfide bonds was not effective. Combining the best mutations resulted in a 12-fold mutant (HheC-H12) with a 28°C higher apparent melting temperature and a remarkable increase in resistance to cosolvents. This variant also showed a higher optimum temperature for catalysis while activity at low temperature was preserved. Mutant H12 was used as a template for the introduction of mutations that enhance enantioselectivity or activity. Crystal structures showed that the structural changes in the H12 mutant mostly agreed with the computational predictions and that the enhanced stability was mainly due to mutations that redistributed surface charges and improved interactions between subunits, the latter including better interactions of water molecules at the subunit interfaces.

X-ray crystallographic validation of structure predictions used in computational design for protein stabilization.

Protein engineering aimed at enhancing enzyme stability is increasingly supported by computational methods for calculation of mutant folding energies and for the design of disulfide bonds. To examine the accuracy of mutant structure predictions underlying these computational methods, crystal structures of thermostable limonene epoxide hydrolase variants obtained by computational library design were determined. Four different predicted effects indeed contributed to the obtained stabilization: (i) enhanced interactions between a flexible loop close to the N-terminus and the rest of the protein; (ii) improved interactions at the dimer interface; (iii) removal of unsatisfied hydrogen bonding groups; and (iv) introduction of additional positively charged groups at the surface. The structures of an eightfold and an elevenfold mutant showed that most mutations introduced the intended stabilizing interactions, and side-chain conformations were correctly predicted for 72-88% of the point mutations. However, mutations that introduced a disulfide bond in a flexible region had a larger influence on the backbone conformation than predicted. The enzyme active sites were unaltered, in agreement with the observed preservation of catalytic activities. The structures also revealed how a c-Myc tag, which was introduced for facile detection and purification, can reduce access to the active site and thereby lower the catalytic activity. Finally, sequence analysis showed that comprehensive mutant energy calculations discovered stabilizing mutations that are not proposed by the consensus or B-FIT methods.

Computational library design for increasing haloalkane dehalogenase stability.

We explored the use of a computational design framework for the stabilization of the haloalkane dehalogenase LinB. Energy calculations, disulfide bond design, molecular dynamics simulations, and rational inspection of mutant structures predicted many stabilizing mutations. Screening of these in small mutant libraries led to the discovery of seventeen point mutations and one disulfide bond that enhanced thermostability. Mutations located in or contacting flexible regions of the protein had a larger stabilizing effect than mutations outside such regions. The combined introduction of twelve stabilizing mutations resulted in a LinB mutant with a 23 °C increase in apparent melting temperature (Tm,app , 72.5 °C) and an over 200-fold longer half-life at 60 °C. The most stable LinB variants also displayed increased compatibility with co-solvents, thus allowing substrate conversion and kinetic resolution at much higher concentrations than with the wild-type enzyme.

Computationally designed libraries for rapid enzyme stabilization.

The ability to engineer enzymes and other proteins to any desired stability would have wide-ranging applications. Here, we demonstrate that computational design of a library with chemically diverse stabilizing mutations allows the engineering of drastically stabilized and fully functional variants of the mesostable enzyme limonene epoxide hydrolase. First, point mutations were selected if they significantly improved the predicted free energy of protein folding. Disulfide bonds were designed using sampling of backbone conformational space, which tripled the number of experimentally stabilizing disulfide bridges. Next, orthogonal in silico screening steps were used to remove chemically unreasonable mutations and mutations that are predicted to increase protein flexibility. The resulting library of 64 variants was experimentally screened, which revealed 21 (pairs of) stabilizing mutations located both in relatively rigid and in flexible areas of the enzyme. Finally, combining 10-12 of these confirmed mutations resulted in multi-site mutants with an increase in apparent melting temperature from 50 to 85°C, enhanced catalytic activity, preserved regioselectivity and a >250-fold longer half-life. The developed Framework for Rapid Enzyme Stabilization by Computational libraries (FRESCO) requires far less screening than conventional directed evolution.

Biocatalytic and structural properties of a highly engineered halohydrin dehalogenase.

Two highly engineered halohydrin dehalogenase variants were characterized in terms of their performance in dehalogenation and epoxide cyanolysis reactions. Both enzyme variants outperformed the wild-type enzyme in the cyanolysis of ethyl (S)-3,4-epoxybutyrate, a conversion yielding ethyl (R)-4-cyano-3-hydroxybutyrate, an important chiral building block for statin synthesis. One of the enzyme variants, HheC2360, displayed catalytic rates for this cyanolysis reaction enhanced up to tenfold. Furthermore, the enantioselectivity of this variant was the opposite of that of the wild-type enzyme, both for dehalogenation and for cyanolysis reactions. The 37-fold mutant HheC2360 showed an increase in thermal stability of 8 °C relative to the wild-type enzyme. Crystal structures of this enzyme were elucidated with chloride and ethyl (S)-3,4-epoxybutyrate or with ethyl (R)-4-cyano-3-hydroxybutyrate bound in the active site. The observed increase in temperature stability was explained in terms of a substantial increase in buried surface area relative to the wild-type HheC, together with enhanced interfacial interactions between the subunits that form the tetramer. The structures also revealed that the substrate binding pocket was modified both by substitutions and by backbone movements in loops surrounding the active site. The observed changes in the mutant structures are partly governed by coupled mutations, some of which are necessary to remove steric clashes or to allow backbone movements to occur. The importance of interactions between substitutions suggests that efficient directed evolution strategies should allow for compensating and synergistic mutations during library design.

Aminoacyl-coenzyme A synthesis catalyzed by a CoA ligase from Penicillium chrysogenum.

Coenzyme A ligases play an important role in metabolism by catalyzing the activation of carboxylic acids. In this study we describe the synthesis of aminoacyl-coenzyme As (CoAs) catalyzed by a CoA ligase from Penicillium chrysogenum. The enzyme accepted medium-chain length fatty acids as the best substrates, but the proteinogenic amino acids L-phenylalanine and L-tyrosine, as well as the non-proteinogenic amino acids D-phenylalanine, D-tyrosine and (R)- and (S)-ß-phenylalanine were also accepted. Of these amino acids, the highest activity was found for (R)-ß-phenylalanine, forming (R)-ß-phenylalanyl-CoA. Homology modeling suggested that alanine 312 is part of the active site cavity, and mutagenesis (A312G) yielded a variant that has an enhanced catalytic efficiency with ß-phenylalanines and D-α-phenylalanine.

The Penicillium chrysogenum aclA gene encodes a broad-substrate-specificity acyl-coenzyme A ligase involved in activation of adipic acid, a side-chain precursor for cephem antibiotics.

Activation of the cephalosporin side-chain precursor to the corresponding CoA-thioester is an essential step for its incorporation into the beta-lactam backbone. To identify an acyl-CoA ligase involved in activation of adipate, we searched in the genome database of Penicillium chrysogenum for putative structural genes encoding acyl-CoA ligases. Chemostat-based transcriptome analysis was used to identify the one presenting the highest expression level when cells were grown in the presence of adipate. Deletion of the gene renamed aclA, led to a 32% decreased specific rate of adipate consumption and a threefold reduction of adipoyl-6-aminopenicillanic acid levels, but did not affect penicillin V production. After overexpression in Escherichia coli, the purified protein was shown to have a broad substrate range including adipate. Finally, protein-fusion with cyan-fluorescent protein showed co-localization with microbody-borne acyl-transferase. Identification and functional characterization of aclA may aid in developing future metabolic engineering strategies for improving the production of different cephalosporins.

Characterization of a phenylacetate-CoA ligase from Penicillium chrysogenum.

Enzymatic activation of PAA (phenylacetic acid) to phenylacetyl-CoA is an important step in the biosynthesis of the beta-lactam antibiotic penicillin G by the fungus Penicillium chrysogenum. CoA esters of PAA and POA (phenoxyacetic acid) act as acyl donors in the exchange of the aminoadipyl side chain of isopenicillin N to produce penicillin G or penicillin V. The phl gene, encoding a PCL (phenylacetate-CoA ligase), was cloned in Escherichia coli as a maltose-binding protein fusion and the biochemical properties of the enzyme were characterized. The recombinant fusion protein converted PAA into phenylacetyl-CoA in an ATP- and magnesium-dependent reaction. PCL could also activate POA, but the catalytic efficiency of the enzyme was rather low with k(cat)/K(m) values of 0.23+/-0.06 and 7.8+/-1.2 mM(-1).s(-1) for PAA and POA respectively. Surprisingly, PCL was very efficient in catalysing the conversion of trans-cinnamic acids to the corresponding CoA thioesters [k(cat)/K(m)=(3.1+/-0.4)x10(2) mM(-1).s(-1) for trans-cinnamic acid]. Of all the substrates screened, medium-chain fatty acids, which also occur as the side chains of the natural penicillins F, DF, H and K, were the best substrates for PCL. The high preference for fatty acids could be explained by a homology model of PCL that was constructed on the basis of sequence similarity with the Japanese firefly luciferase. The results suggest that PCL has evolved from a fatty-acid-activating ancestral enzyme that may have been involved in the beta-oxidation of fatty acids.

Saturation mutagenesis reveals the importance of residues alphaR145 and alphaF146 of penicillin acylase in the synthesis of beta-lactam antibiotics.

Penicillin acylase (PA) from Escherichia coli can catalyze the coupling of an acyl group to penicillin- and cephalosporin-derived beta-lactam nuclei, a conversion that can be used for the industrial synthesis of beta-lactam antibiotics. The modest synthetic properties of the wild-type enzyme make it desirable to engineer improved mutants. Analysis of the crystal structure of PA has shown that residues alphaR145 and alphaF146 undergo extensive repositioning upon binding of large ligands to the active site, suggesting that these residues may be good targets for mutagenesis aimed at improving the catalytic performance of PA. Therefore, site-saturation mutagenesis was performed on both positions and a complete set of all 38 variants was subjected to rapid HPLC screening for improved ampicillin synthesis. Not less than 33 mutants showed improved synthesis, indicating the importance of the mutated residues in PA-catalyzed acyl transfer kinetics. In several mutants at low substrate concentrations, the maximum level of ampicillin production was increased up to 1.5-fold, and the ratio of the synthetic rate over the hydrolytic rate was increased 5-15-fold. Moreover, due to increased tendency of the acyl-enzyme intermediate to react with beta-lactam nucleophile instead of water, mutants alphaR145G, alphaR145S and alphaR145L demonstrated an enhanced synthetic yield over wild-type PA at high substrate concentrations. This was accompanied by an increased conversion of 6-APA to ampicillin as well as a decreased undesirable hydrolysis of the acyl donor. Therefore, these mutants are interesting candidates for the enzymatic production of semi-synthetic beta-lactam antibiotics.

Acetobacter turbidans alpha-amino acid ester hydrolase: how a single mutation improves an antibiotic-producing enzyme.

The alpha-amino acid ester hydrolase (AEH) from Acetobacter turbidans is a bacterial enzyme catalyzing the hydrolysis and synthesis of beta-lactam antibiotics. The crystal structures of the native enzyme, both unliganded and in complex with the hydrolysis product D-phenylglycine are reported, as well as the structures of an inactive mutant (S205A) complexed with the substrate ampicillin, and an active site mutant (Y206A) with an increased tendency to catalyze antibiotic production rather than hydrolysis. The structure of the native enzyme shows an acyl binding pocket, in which D-phenylglycine binds, and an additional space that is large enough to accommodate the beta-lactam moiety of an antibiotic. In the S205A mutant, ampicillin binds in this pocket in a non-productive manner, making extensive contacts with the side chain of Tyr(112), which also participates in oxyanion hole formation. In the Y206A mutant, the Tyr(112) side chain has moved with its hydroxyl group toward the catalytic serine. Because this changes the properties of the beta-lactam binding site, this could explain the increased beta-lactam transferase activity of this mutant.

The patatin-like protein from the latex of Hevea brasiliensis (Hev b 7) is not a vacuolar protein.

Upon centrifugation, rubber latex is divided into a layer of rubber particles, the cytosol, and the lutoid-body fraction, which is of vacuolar origin. One of the proteins isolated from the lutoid-body fraction is a protein with a molecular mass of 43 kDa, which has esterase activity on p-nitrophenylpalmitate and which shows significant sequence similarity with patatin, a vacuolar protein with esterase activity from potato (Solanum tuberosum). This protein is a major allergen in rubber latex products (Hev b 7) and can also be isolated from the cytosol fraction of rubber latex. The mature protein isolated from lutoid-bodies has no structural features expected for a vacuolar protein: the N-terminal methionine in the cDNA-derived sequence is cleaved off, the second residue is N-acetylated, and the C-terminal sequence is identical to that in the cDNA-derived sequence. Thus the patatin-like protein in Hevea brasiliensis is not a vacuolar protein, but may be associated with not yet characterized particles in the cytoplasm, which either sediment with lutoid-bodies or remain in the cytosol fraction, depending on the centrifugation conditions.

The sequence and crystal structure of the alpha-amino acid ester hydrolase from Xanthomonas citri define a new family of beta-lactam antibiotic acylases.

alpha-Amino acid ester hydrolases (AEHs) catalyze the hydrolysis and synthesis of esters and amides with an alpha-amino group. As such, they can synthesize beta-lactam antibiotics from acyl compounds and beta-lactam nuclei obtained from the hydrolysis of natural antibiotics. This article describes the gene sequence and the 1.9-A resolution crystal structure of the AEH from Xanthomonas citri. The enzyme consists of an alpha/beta-hydrolase fold domain, a helical cap domain, and a jellyroll beta-domain. Structural homology was observed to the Rhodococcus cocaine esterase, indicating that both enzymes belong to the same class of bacterial hydrolases. Docking of a beta-lactam antibiotic in the active site explains the substrate specificity, specifically the necessity of an alpha-amino group on the substrate, and explains the low specificity toward the beta-lactam nucleus.

X-ray analysis of two antibiotic-synthesizing bacterial ester hydrolases: preliminary results.

Alpha-amino-acid ester hydrolases are multimeric enzymes of potential use in antibiotic production. Knowledge of their structure could help to engineer these enzymes into economically viable biocatalysts. The alpha-amino-acid ester hydrolases from Xanthomonas citri and Acetobacter turbidans have been crystallized. The X. citri enzyme crystallizes in a primitive monoclinic space group (unit-cell parameters a = 90.1, b = 125.8, c = 132.1 A, beta = 90.9 degrees ). The A. turbidans enzyme crystallizes in both a primitive orthorhombic (a = 99.1, b = 104.9, c = 284.9 A) and a body-centred cubic space group with a = b = c = 342.2 A. From both enzymes, diffraction-quality crystals (resolution 3.0 A or better) were obtained. Data-collection statistics are reported for data sets from both enzymes.

Pancreatic-type ribonuclease 1 gene duplications in rat species.

Mammalian pancreatic-type ribonucleases (RNases) 1 represent single-copy genes in the genome of most investigated mammalian species, including Mus musculus and other murid rodents. However, in six species belonging to the genus Rattus and closely related taxa, several paralogous gene products were identified by Southern blotting and PCR amplifications of genomic sequences. Phylogenies of nucleotide and derived amino acid sequences were reconstructed by several procedures, with three Mus species as outgroup. Duplications of the RNase 1 occurred after the divergence of Niviventer cremoriventer and Leopoldamys edwardsi from the other investigated species. Four groups of paralogous genes could be identified from specific amino acid sequence features in each of them. Low ratios of nonsynonymous-to-synonymous substitutions and the paucity of pseudogene features suggest functional gene products. One of the RNase 1 genes of R. norvegicus is expressed in the pancreas. RNases 1 were isolated from pancreatic tissues of R. rattus and R. exulans and submitted to N-terminal amino acid sequence analysis. In R. rattus, the orthologue of the expressed gene of R. norvegicus was identified, but in R. exulans, two paralogous gene products were found. The gene encoding for one of these had not yet been found by PCR amplification of genomic DNA. A well-defined group of orthologous sequences found in five investigated species codes for very basic RNases. Northern blot analysis showed expression of messenger RNA for this RNase in the spleen of R. norvegicus, but the protein product could not be identified. Evolutionary rates of RNase 1, expressed as nucleotide substitutions per site per 10(3) million years (Myr), vary between 5 and 9 in the lines leading to Mus, Niviventer, and Lepoldamys (on the basis of an ancestral date of mouse/rat divergence of 12.2 Myr) and between 20 and 50 in the lines to the other sequences after divergence from Niviventer and Leopoldamys (5.5 Myr).

Identification of the catalytic residues of alpha-amino acid ester hydrolase from Acetobacter turbidans by labeling and site-directed mutagenesis.

The alpha-amino acid ester hydrolase from Acetobacter turbidans ATCC 9325 is capable of hydrolyzing and synthesizing the side chain peptide bond in beta-lactam antibiotics. Data base searches revealed that the enzyme contains an active site serine consensus sequence Gly-X-Ser-Tyr-X-Gly that is also found in X-prolyl dipeptidyl aminopeptidase. The serine hydrolase inhibitor p-nitrophenyl-p'-guanidino-benzoate appeared to be an active site titrant and was used to label the alpha-amino acid ester hydrolase. Electrospray mass spectrometry and tandem mass spectrometry analysis of peptides from a CNBr digest of the labeled protein showed that Ser(205), situated in the consensus sequence, becomes covalently modified by reaction with the inhibitor. Extended sequence analysis showed alignment of this Ser(205) with the catalytic nucleophile of some alpha/beta-hydrolase fold enzymes, which posses a catalytic triad composed of a nucleophile, an acid, and a base. Based on the alignments, 10 amino acids were selected for site-directed mutagenesis (Arg(85), Asp(86), Tyr(143), Ser(156), Ser(205), Tyr(206), Asp(338), His(370), Asp(509), and His(610)). Mutation of Ser(205), Asp(338,) or His(370) to an alanine almost fully inactivated the enzyme, whereas mutation of the other residues did not seriously affect the enzyme activity. Circular dichroism measurements showed that the inactivation was not caused by drastic changes in the tertiary structure. Therefore, we conclude that the catalytic domain of the alpha-amino acid ester hydrolase has an alpha/beta-hydrolase fold structure with a catalytic triad of Ser(205), Asp(338), and His(370). This distinguishes the alpha-amino acid ester hydrolase from the Ntn-hydrolase family of beta-lactam antibiotic acylases.

Cloning, sequence analysis, and expression in Escherichia coli of the gene encoding an alpha-amino acid ester hydrolase from Acetobacter turbidans.

The alpha-amino acid ester hydrolase from Acetobacter turbidans ATCC 9325 is capable of hydrolyzing and synthesizing beta-lactam antibiotics, such as cephalexin and ampicillin. N-terminal amino acid sequencing of the purified alpha-amino acid ester hydrolase allowed cloning and genetic characterization of the corresponding gene from an A. turbidans genomic library. The gene, designated aehA, encodes a polypeptide with a molecular weight of 72,000. Comparison of the determined N-terminal sequence and the deduced amino acid sequence indicated the presence of an N-terminal leader sequence of 40 amino acids. The aehA gene was subcloned in the pET9 expression plasmid and expressed in Escherichia coli. The recombinant protein was purified and found to be dimeric with subunits of 70 kDa. A sequence similarity search revealed 26% identity with a glutaryl 7-ACA acylase precursor from Bacillus laterosporus, but no homology was found with other known penicillin or cephalosporin acylases. There was some similarity to serine proteases, including the conservation of the active site motif, GXSYXG. Together with database searches, this suggested that the alpha-amino acid ester hydrolase is a beta-lactam antibiotic acylase that belongs to a class of hydrolases that is different from the Ntn hydrolase superfamily to which the well-characterized penicillin acylase from E. coli belongs. The alpha-amino acid ester hydrolase of A. turbidans represents a subclass of this new class of beta-lactam antibiotic acylases.