Mutations in distant residues moderately increase the enantioselectivity of Pseudomonas fluorescens esterase towards methyl 3bromo-2-methylpropanoate and ethyl 3phenylbutyrate

Directed evolution combined with saturation mutagenesis identified six different point mutations that each moderately increases the enantioselectivity of an esterase from Pseudomonas fluorescens (PFE) towards either of two chiral synthons. Directed evolution identified a Thr230Ile mutation that increased the enantioselectivity from 12 to 19 towards methyl (S)-3-bromo-2-methylpropanoate. Saturation mutagenesis at Thr230 identified another mutant, Thr230Pro, with higher-than-wild-type enantioselectivity (E=17). Previous directed evolution identified mutants Asp158Asn and Leu181Gln that increased the enantioselectivity from 3.5 to 5.8 and 6.6, respectively, towards ethyl (R)-3-phenylbutyrate. In this work, saturation mutagenesis identified other mutations that further increase the enantioselectivity to 12 (Asp158Leu) and 10 (Leu181Ser). A homology model of PFE indicates that all mutations lie outside the active site, 12-14 A from the substrate and suggests how the distant mutations might indirectly change the substrate-binding site. Since proteins contain many more residues far from the active site than close to the active site, random mutagenesis is strongly biased in favor of distant mutations. Directed evolution rarely screens all mutations, so it usually finds the distant mutations because they are more common, but probably not the most effective.     

Computational design of protein–ligand binding: Modifying the specificity of asparaginyl‐tRNA synthetase

A method for computational design of protein–ligand interactions is implemented and tested on the asparaginyl‐ and aspartyl‐tRNA synthetase enzymes (AsnRS, AspRS). The substrate specificity of these enzymes is crucial for the accurate translation of the genetic code. The method relies on a molecular mechanics energy function and a simple, continuum electrostatic, implicit solvent model. As test calculations, we first compute AspRS‐substrate binding free energy changes due to nine point mutations, for which experimental data are available; we also perform large‐scale redesign of the entire active site of each enzyme (40 amino acids) and compare to experimental sequences. We then apply the method to engineer an increased binding of aspartyl‐adenylate (AspAMP) into AsnRS. Mutants are obtained using several directed evolution protocols, where four or five amino acid positions in the active site are randomized. Promising mutants are subjected to molecular dynamics simulations; Poisson‐Boltzmann calculations provide an estimate of the corresponding, AspAMP, binding free energy changes, relative to the native AsnRS. Several of the mutants are predicted to have an inverted binding specificity, preferring to bind AspAMP rather than the natural substrate, AsnAMP. The computed binding affinities are significantly weaker than the native, AsnRS:AsnAMP affinity, and in most cases, the active site structure is significantly changed, compared to the native complex. This almost certainly precludes catalytic activity. One of the designed sequences has a higher affinity and more native‐like structure and may represent a valid candidate for Asp activity. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Expanding the substrate scope of enzymes: combining mutations obtained by CASTing

In a previous paper, the combinatorial active-site saturation test (CAST) was introduced as an effective strategy for the directed evolution of enzymes toward broader substrate acceptance. CASTing comprises the systematic design and screening of focused libraries around the complete binding pocket, but it is only the first step of an evolutionary process because only the initial libraries of mutants are considered. In the present study, a simple method is presented for further optimization of initial hits by combining the mutational changes obtained from two different libraries. Combined lipase mutants were screened for hydrolytic activity against six notoriously difficult substrates (bulky carboxylic acid esters) and improved mutants showing significantly higher activity were identified. The enantioselectivity of the mutants in the hydrolytic kinetic resolution of two substrates was also studied, with the best mutant-substrate combination resulting in a selectivity factor of E=49. Finally, the catalytic profile of the evolved mutants in the hydrolysis of simple nonbranched carboxylic acid esters, ranging from acetate to palmitate, was studied for theoretical reasons.     

Directed evolution of epoxide hydrolase from A. radiobacter toward higher enantioselectivity by error-prone PCR and DNA shuffling

The enantioselectivity of epoxide hydrolase from Agrobacterium radiobacter (EchA) was improved using error-prone PCR and DNA shuffling. An agar plate assay was used to screen the mutant libraries for activity. Screening for improved enantioselectivity was subsequently done by spectrophotometric progress curve analysis of the conversion of para-nitrophenyl glycidyl ether (pNPGE). Kinetic resolutions showed that eight mutants were obtained with up to 13-fold improved enantioselectivity toward pNPGE and at least three other epoxides. The large enhancements in enantioselectivity toward epichlorohydrin and 1,2-epoxyhexane indicated that pNPGE acts as an epoxyalkane mimic. Active site mutations were found in all shuffled mutants, which can be explained by an interaction of the affected amino acid with the epoxide oxygen or the hydrophobic moiety of the substrate. Several mutations in the shuffled mutants had additive effects.     

Directed evolution of an enantioselective lipase

BACKGROUND: The biocatalytic production of enantiopure compounds is of steadily increasing importance to the chemical and biotechnological industry. In most cases, however, it is impossible to identify an enzyme that possesses the desired enantioselectivity. Therefore, there is a strong need to create by molecular biological methods novel enzymes which display high enantioselectivity. RESULTS: A bacterial lipase from Pseudomonas aeruginosa (PAL) was evolved to catalyze with high enantioselectivity the hydrolysis of the chiral model substrate 2-methyldecanoic acid p-nitrophenyl ester. Successive rounds of random mutagenesis by ep-PCR and saturation mutagenesis resulted in an increase in enantioselectivity from E=1.1 for the wild-type enzyme to E=25.8 for the best variant which carried five amino acid substitutions. The recently solved three-dimensional structure of PAL allowed us to analyze the structural consequences of these substitutions. CONCLUSIONS: A highly enantioselective lipase was created by increasing the flexibility of distinct loops of the enzyme. Our results demonstrate that enantioselective enzymes can be created by directed evolution, thereby opening up a large area of novel applications in biotechnology.     

Enzyme redesign: two mutations cooperate to convert cycloartenol synthase into an accurate lanosterol synthase

Efforts to modify the catalytic specificity of enzymes consistently show that it is easier to broaden the substrate or product specificity of an accurate enzyme than to restrict the selectivity of one that is promiscuous. Described herein are experiments in which cycloartenol synthase was redesigned to become a highly accurate lanosterol synthase. Several single mutants have been described that modify the catalytic specificity of cycloartenol to form some lanosterol. Modeling studies were undertaken to identify combinations of mutations that cooperate to decrease the formation of products other than lanosterol. A double mutant was constructed and characterized and was shown to cyclize oxidosqualene accurately to lanosterol (99%). This catalytic change entailed both relocating polarity with a His477Asn mutation and modifying steric constraints with an Ile481Val mutation.     

Enantioselective enzymes for organic synthesis created by directed evolution

A new concept for the creation of enzymes displaying improved enantioselectivity in a given reaction is described; it is based on "evolution in the test tube". Accordingly, proper molecular biological methods for random mutagenesis, gene expression, and high-throughput screening systems for the rapid assay of enantioselectivity are combined. Several rounds of mutagenesis and screening are generally necessary in order to create mutant enzymes that show high degrees of enantioselectivity, as in the case of the lipase-catalyzed hydrolytic kinetic resolution of a chiral ester in which the original enantioselectivity of 2 % ee (E = 1) increases to > 90% ee (E = 25).     

Systematic characterization of adenosine triphosphate response to lung cancer epidermal growth factor receptor missense mutations: A molecular insight into “generic” drug resistance mutations

Many missense mutations in human epidermal growth factor receptor (EGFR) are clinically involved in lung cancer and may cause acquired resistance to tyrosine kinase inhibitors. Traditionally, the resistance is considered to be established by impairing inhibitor affinity due to the mutations. However, it was found that, instead of blocking inhibitor binding, the gatekeeper mutation T790M can improve the kinase affinity for its natural substrate adenosine triphosphate (ATP), which is thus regarded as a “generic” resistance mutation that will reduce the potency of any ATP-competitive reversible kinase inhibitor. In this study, we attempt to systematically investigate the binding behavior of ATP to clinically observed EGFR missense mutants in nonsmall-cell lung cancer to identify those substantial mutations that may significantly increase (or decrease) ATP affinity. Several substantial mutations are excluded because they are also involved in kinase's catalytic activity or directly influence inhibitor binding, thus largely complicating the multiple dependent relationships of kinase, ATP, and inhibitor. Two new “generic” resistance mutations, A839T and E758G, are identified, which can improve ATP affinity by forming a favorable hydrogen bond and by eliminating unfavorable electrostatic effect between the kinase and ATP, respectively.     

Computational Study of Mechanism and Enantioselectivity of Imine Reductase from Amycolatopsis orientalis

Imine reductases (IREDs) are NADPH-dependent enzymes (NADPH=nicotinamide adenine dinucleotide phosphate) that catalyze the reduction of imines to amines. They exhibit high enantioselectivity for a broad range of substrates, making them of interest for biocatalytic applications. In this work, we have employed density functional theory (DFT) calculations to elucidate the reaction mechanism and the origins of enantioselectivity of IRED from Amycolatopsis orientalis. Two substrates are considered, namely 1-methyl-3,4-dihydroisoquinoline and 2-propyl-piperideine. A model of the active site is built on the basis of the available crystal structure. For both substrates, different binding modes are first evaluated, followed by calculation of the hydride transfer transition states from each complex. We have also investigated the effect of mutations of certain important active site residues (Tyr179Ala and Asn241Ala) on the enantioselectivity. The calculated energies are consistent with the experimental observations and the analysis of transition states geometries provides insights into the origins of enantioselectivity of this enzyme.     

Engineering temperature-sensitive SH3 domains.

BACKGROUND: The ability to control specific protein-protein interactions conditionally in vivo would be extremely helpful for analyzing protein-protein interaction networks. SH3 (Src homology 3) modular protein binding domains are found in many signaling proteins and they play a crucial role in signal transduction by binding to proline-rich sequences. RESULTS: Random in vitro mutagenesis coupled with yeast two-hybrid screening was used to identify mutations in the second SH3 domain of Nck that render interaction with its ligand temperature sensitive. Four of the mutants were functionally temperature sensitive in mammalian cells, where temperature sensitivity was correlated with a pronounced instability of the mutant domains at the nonpermissive temperature. Two of the mutations affect conserved residues in the hydrophobic core (Val133 and Val160), suggesting a general strategy for engineering temperature-sensitive SH3-containing proteins. Indeed mutagenesis of the corresponding positions in another SH3 domain, that of Crk-1, rendered the full-length Crk-1 protein temperature sensitive for function and stability in mammalian cells. CONCLUSIONS: Construction of temperature-sensitive SH3 domains is a novel approach to regulating the function of SH3 domains in vivo. Such mutants will be valuable in dissecting SH3-mediated signaling pathways. Furthermore, the methodology described here to isolate temperature-sensitive domains should be widely applicable to any domain involved in protein-protein interactions.     

Multiscale simulation reveals multiple pathways for H2 and O2 transport in a [NiFe]-hydrogenase

Hydrogenases are enzymes that catalyze the reversible conversion of hydrogen molecules to protons and electrons. The mechanism by which the gas molecules reach the active site is important for understanding the function of the enzyme and may play a role in the selectivity for hydrogen over inhibitor molecules. Here, we develop a general multiscale molecular simulation approach for the calculation of diffusion rates and determination of pathways by which substrate or inhibitor gases can reach the protein active site. Combining kinetic data from both equilibrium simulations and enhanced sampling, we construct a master equation describing the movement of gas molecules within the enzyme. We find that the time-dependent gas population of the active site can be fit to the same phenomenological rate law used to interpret experiments, with corresponding diffusion rates in very good agreement with experimental data. However, in contrast to the conventional picture, in which the gases follow a well-defined hydrophobic tunnel, we find that there is a diverse network of accessible pathways by which the gas molecules can reach the active site. The previously identified tunnel accounts for only about 60% of the total flux. Our results suggest that the dramatic decrease in the diffusion rate for mutations involving the residue Val74 could be in part due to the narrowing of the passage Val74-Arg476, immediately adjacent to the binding site, explaining why mutations of Leu122 had only a negligible effect in experiment. Our method is not specific to the [NiFe]-hydrogenase and should be generally applicable to the transport of small molecules in proteins.     

New Structural Motif for Carboxylic Acid Perhydrolases

Some serine hydrolases also catalyze a promiscuous reaction— reversible perhydrolysis of carboxylic acids to make peroxycarboxylic acids. Five X‐ray crystal structures of these carboxylic acid perhydrolases show a proline in the oxyanion loop. Here, we test whether this proline is essential for high perhydrolysis activity using Pseudomonas fluorescens esterase (PFE). The L29P variant of this esterase catalyzes perhydrolysis 43‐fold faster (k_cat comparison) than the wild type. Surprisingly, saturation mutagenesis at the 29 position of PFE identified six other amino acid substitutions that increase perhydrolysis of acetic acid at least fourfold over the wild type. The best variant, L29I PFE, catalyzed perhydrolysis 83‐times faster (k_cat comparison) than wild‐type PFE and twice as fast as L29P PFE. Despite the different amino acid in the oxyanion loop, L29I PFE shows a similar selectivity for hydrogen peroxide over water as L29P PFE (β₀=170 vs. 160 M ^?1), and a similar fast formation of acetyl‐enzyme (140 vs. 62 U mg^?1). X‐ray crystal structures of L29I PFE with and without bound acetate show an unusual mixture of two different oxyanion loop conformations. The type II β‐turn conformation resembles the wild‐type structure and is unlikely to increase perhydrolysis, but the type I β‐turn conformation creates a binding site for a second acetate. Modeling suggests that a previously proposed mechanism for L29P PFE can be extended to include L29I PFE, so that an acetate accepts a hydrogen bond to promote faster formation of the acetyl‐enzyme.     

Inverting the enantioselectivity of a carbonyl reductase via substrate-enzyme docking-guided point mutation

Substrate-enzyme docking-guided point mutation of a carbonyl reductase from Sporobolomyces salmonicolor led to mutant enzymes, which reversed the enantiopreference and enhanced the enantioselectivity toward the reduction of para-substituted acetophenones. Such a dramatic change in the enantioselectivity indicates that the 245 residue in the catalytic site plays a critical role in determining the enantioselectivity of these ketone reductions, providing valuable insight into our understanding of how residues involved in substrate binding affect the orientation of bound substrate and thus control the reduction stereoselectivity.     

Creation of a broad-range and highly stereoselective D-amino acid dehydrogenase for the one-step synthesis of D-amino acids

Using both rational and random mutagenesis, we have created the first known broad substrate range, nicotinamide cofactor dependent, and highly stereoselective d-amino acid dehydrogenase. This new enzyme is capable of producing d-amino acids via the reductive amination of the corresponding 2-keto acid with ammonia. This biocatalyst was the result of three rounds of mutagenesis and screening performed on the enzyme meso-diaminopimelate d-dehydrogenase. The first round targeted the active site of the wild-type enzyme and produced mutants that were no longer strictly dependent on the native substrate. The second and third rounds produced mutants that had an increased substrate range including straight- and branched-aliphatic amino acids and aromatic amino acids. The very high selectivity toward the d-enantiomer (95 to >99% ee) was shown to be preserved even after the addition of the five mutations found in the three rounds of mutagenesis and screening. This new enzyme could complement and improve upon current methods for d-amino acid synthesis.     

Increasing the enantioselectivity of cyclopentanone monooxygenase (CPMO): profile of new CPMO mutants

A series of cyclohexanones substituted at the 4-position with a selection of hydrophobic and hydrophilic groups were used as substrates in the evaluation of six new cyclopentanone monooxygenase (CPMO) mutants. These mutants were obtained through evolutionary modifications in two specific regions of the CPMO’s putative active site. Several mutant enzymes with improved enantioselectivity were identified. Analysis of the results, in terms of a diamond model, illustrates how a family of cyclohexanone substrates may be used to explore putative active sites of Baeyer–Villiger monooxygenases (BVMOs) and to design productive mutations for specific substrates.     

A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.     

羰基还原酶辅酶结合域位点突变对其不对称催化性能的影响

基于同源模型的比较和分析,发现羰基还原酶SCR1辅酶结合域P124和W125位点对辅酶NADPH的结合形成了一定的空间位阻效应.通过对该位点进行小侧基氨基酸的取代突变,该酶的底物专一性和立体选择性均发生了不同程度的改变,表明该位点是酶与辅酶有效结合的关键位点,而且它与辅酶结合的空间效应进一步影响了底物结合域活性中心对不同构型的底物及其对映体产物的亲和作用.在底物专一性方面,野生型酶对2-羟基苯乙酮和2-溴苯乙酮及其衍生物等底物表现出较高的催化活性,而突变株W125A,W125G,P124A/W125A和P124G/W125G对苯乙酮及其部分衍生物和2-辛酮等底物的催化活性均有所提高.对于酶的立体选择性,部分突变株发生了转化产物对映体构型反转的现象,突变株P124A/W125A和P124G/W125G催化还原2-羟基苯乙酮和4-氯乙酰乙酸乙酯均生成了(R)-型产物.     

Insight into catalytically relevant correlated motions in human purine nucleoside phosphorylase

The catalytic site of the homotrimeric enzyme human purine nucleoside phosphorylase enzyme (hPNP) features residue F200 and the 241-265 loop directly skirting the purine base and a residue belonging to the adjacent monomer, F159, immediately conterminous to the ribosyl moiety. Crystallographic B-factors of apo human purine nucleoside phosphorylase, and hPNP complexed with substrate/transition state (TS) analogues, show that residue E250 is the centroid of a highly mobile loop region. Furthermore, superimposition of apo hPNP and hPNP complexed with TS analogue Immucillin-H shows a tightening of the active site, caused by the ligand-dependent 241-265 loop rearrangement taking place upon substrate/inhibitor binding, suggesting a putative dynamic role of the loop in binding/catalysis. However, crystallographic structures reveal only average atomic positions, and more detailed information is needed to discern the dynamic behavior of hPNP. The Essential Dynamics (ED) method is used here to investigate the existence of correlated motions in hPNP and consequently proposes mutagenesis assays to estimate the relative importance of these motions in the phosphorolytic efficiency of the reaction catalyzed by hPNP. We compare the concerted motions obtained from multiple molecular dynamics simulations of apo and Michaelis complex of hPNP both in vacuo and in solution. The results of the principal component analysis for the apo hPNP indicate the existence of strong correlations predominantly in the vicinity of residue F159. However, for the Michaelis complex, concerted motions are seen mostly around both active site residue F200 and loop residue E250. Additionally, for a simulation depicting the relaxation of tight complexed hPNP with a TS analogue, toward its relaxed apo form (after removal of the TS analog), a combination of the apo hPNP and Michaelis complex motions is found, with prominent concerted modes centered around neighboring residues F159, F200, and E250. Finally, we probed the extent to which these concerted motions bear an intrinsic catalytic role by performing experimental site-directed mutagenesis on some residues, followed by kinetic analysis. The F159G and F200G mutants displayed a strong increase in K(M) and modest decrease in k(cat), suggesting that these concerted motions may provide dynamical roles in substrate binding and/or catalysis. However, further structural data for the hPNP mutants are needed to confirm our hypothesis.     

Mapping the substrate selectivity of new hydrolases using colorimetric screening: lipases from Bacillus thermocatenulatus and Ophiostoma piliferum,esterases from Pseudomonas fluorescens and Streptomyces diastatochromogenes

Recent advances in biochemistry and molecular biology have simplified the discovery and preparation of new hydrolases. Although these hydrolases might solve problems in organic synthesis, measuring their selectivity, especially enantioselectivity, remains tedious and time consuming. Recently, we developed a colorimetric screening method to measure the enantioselectivity of hydrolases. Here we apply this rapid screening method to map the substrate selectivity of four new hydrolases: lipases from the thermophilic Bacillus thermocatenulatus (DSM 730, BTL2) and a filamentous fungus Ophiostoma piliferum (NRRL 18917, OPL) and esterases from two bacteria, Pseudomonas fluorescens (SIK-W1, esterase I, PFE) and Streptomyces diastatochromogenes (Tü 20, SDE). We screened a general library of 29 substrates and a chiral library of 23 pairs of enantiomers. All four hydrolases catalysed the hydrolysis of unnatural substrates, but the two lipases accepted a broader range of substrates than the two esterases. As expected, the two lipases favoured more hydrophobic substrates, while the two esterases showed a preference for smaller substrates. Several moderately enantioselective reactions were identified for the solketal esters: BTL2, butyrate, E=7.9 (R); octanoate, E=4.9 (R) and 3-bromo-2-methyl propionate methyl esters, PFE, E=12 (S); SDE, E=5.6 (S). OPL showed low enantioselectivity toward all substrates tested. The current colorimetric screen could not measure the selectivity for several slow-reacting substrates. Traditional screening identified high enantioselectivity of BTL2 and PFE toward one of these slow substrates, 1-phenylethyl acetate (E>50).     

Ferrous binding to the multicopper oxidases Saccharomyces cerevisiae Fet3p and human ceruloplasmin: contributions to ferroxidase activity

The multicopper oxidases are a family of enzymes that couple the reduction of O(2) to H(2)O with the oxidation of a range of substrates. Saccharomyces cerevisiae Fet3p and human ceruloplasmin (hCp) are members of this family that exhibit ferroxidase activity. Their high specificity for Fe(II) has been attributed to the existence of a binding site for iron. In this study, mutations at the E185 and Y354 residues, which are putative ligands for iron in Fet3p, have been generated and characterized. The effects of these mutations on the electronic structure of the T1 Cu site have been assessed, and the reactivities of this site toward 1,4-hydroquinone (a weak binding substrate) and Fe(II) have been evaluated and interpreted in terms of the semiclassical Marcus theory for electron transfer. The electronic and geometric structure of the Fe(II) substrate bound to Fet3p and hCp has been studied for the first time, using variable-temperature variable field magnetic circular dichroism (VTVH MCD) spectroscopy. The iron binding sites in Fet3p and hCp appear to be very similar in nature, and their contributions to the ferroxidase activity of these proteins have been analyzed. It is found that these iron binding sites play a major role in tuning the reduction potential of iron to provide a large driving force for the ferroxidase reaction, while still supporting the delivery of the Fe(III) product to the acceptor protein. Finally, the analysis of possible electron-transfer (ET) pathways from the protein-bound Fe(II) to the T1 Cu site indicates that the E185 residue not only plays a role in iron binding, but also provides the dominant ET pathway to the T1 Cu site.