A QM/MM study of the reaction mechanism of (R)-hydroxynitrile lyases from Arabidopsis thaliana (AtHNL)

Enzyme engineering improves catalytic efficiency and enantioselectivity of hydroxynitrile lyase for promiscuous retro-nitroaldolase activity

Protein engineering to improve promiscuous catalytic activity is important for biocatalytic application of enzymes in green synthesis. We uncovered the significance of binding site residues in Arabidopsis thaliana hydroxynitrile lyase (AtHNL) for promiscuous retro-nitroaldolase activity. Engineering of AtHNL has improved enantioselective retro-nitroaldolase activity, a synthetically important biotransformation, for the production of enantiopure β-nitroalcohols having absolute configuration opposite to that of the stereopreference of the HNL. The variant F179A has shown ∼ 12 fold increased selectivity towards the retro-nitroaldol reaction over cyanogenesis, the natural activity of the parent enzyme. Screening of the two saturation libraries of Phe179 and Tyr14 revealed several variants with higher k_cat, while F179N showed ∼ 2.4-fold k_cat/K_m than the native enzyme towards retro-nitroaldol reaction. Variants F179N, F179M, F179W, F179V, F179I, Y14L, and Y14M have shown > 99% ee in the preparation of (S)-2-nitro-1-phenylethanol (NPE) from the racemic substrate, while F179N has shown the E value of 138 vs. 81 by the wild type. Our molecular docking and dynamics simulations (MDS) studies results provided insights into the molecular basis of higher enantioselectivity by the F179N toward the retro-nitroaldolase activity than the other mutants. Binding energy calculations also showed the higher negative binding free energy in the case of F179N-(R)-NPE compared to other complexes that support our experimental low K_m by the F179N for NPE. A plausible retro-nitroaldol reaction mechanism was proposed based on the MDS study of enzyme-substrate interaction.

Demystifying the catalytic pathway of Mycobacterium tuberculosis isocitrate lyase

Pulmonary tuberculosis, caused by Mycobacterium tuberculosis, is one of the most persistent diseases leading to death in humans. As one of the key targets during the latent/dormant stage of M. tuberculosis, isocitrate lyase (ICL) has been a subject of interest for new tuberculosis therapeutics. In this work, the cleavage of the isocitrate by M. tuberculosis ICL was studied using quantum mechanics/molecular mechanics method at M06-2X/6-31+G(d,p): AMBER level of theory. The electronic embedding approach was applied to provide a better depiction of electrostatic interactions between MM and QM regions. Two possible pathways (pathway I that involves Asp108 and pathway II that involves Glu182) that could lead to the metabolism of isocitrate was studied in this study. The results suggested that the core residues involved in isocitrate catalytic cleavage mechanism are Asp108, Cys191 and Arg228. A water molecule bonded to Mg²⁺ acts as the catalytic base for the deprotonation of isocitrate C(2)–OH group, while Cys191 acts as the catalytic acid. Our observation suggests that the shuttle proton from isocitrate hydroxyl group C(2) atom is favourably transferred to Asp108 instead of Glu182 with a lower activation energy of 6.2 kcal/mol. Natural bond analysis also demonstrated that pathway I involving the transfer of proton to Asp108 has a higher intermolecular interaction and charge transfer that were associated with higher stabilization energy. The QM/MM transition state stepwise catalytic mechanism of ICL agrees with the in vitro enzymatic assay whereby Asp108Ala and Cys191Ser ICL mutants lost their isocitrate cleavage activities.

Plant Defenses Against Pests Driven by a Bidirectional Promoter

The plant defense responses to pests results in the synchronized change of a complex network of interconnected genes and signaling pathways. An essential part of this process is mediated by the binding of transcription factors to the specific responsive cis-elements within in the promoters of phytophagous-responsive genes. In this work, it is reported the identification and characterization of a bidirectional promoter that simultaneously co-regulate two divergent genes, At5g10300 and At5g10290, upon arthropod feeding. Computational analysis identified the presence of cis-elements within the intergenic region between two loci, mainly from the DOF but also from the AP2/ERF, Golden 2-like and bHLH families. The function of the bidirectional promoter was analyzed using two enhanced variants of the GFP and CherryFP reporter genes, in both orientations, in transient tobacco and stably transformed Arabidopsis plants. Promoter activity was tested in response to feeding of Tetranychus urticae and Pieris brassicae, as well as wounding, flagellin and chitin treatments. Using RT-qPCR assays and confocal microscopy, it was shown that all treatments resulted in the induction of both reporter genes. Furthermore, our findings revealed the asymmetric character of the promoter with stronger activity in the forward than in the reverse orientation. This study provides an example of a bidirectional promoter with a strong potential to be used in plant biotechnology in pest control that requires stacking of the defense genes.

Identical active sites in hydroxynitrile lyases show opposite enantioselectivity and reveal possible ancestral mechanism

Evolutionarily related hydroxynitrile lyases from rubber tree (HbHNL) and from Arabidopsis thaliana (AtHNL) follow different catalytic mechanisms with opposite enantioselectivity toward mandelonitrile. We hypothesized that the HbHNL-like mechanism evolved from an enzyme with an AtHNL-like mechanism. We created ancestor-like composite active-sites in each scaffold to elucidate how this transition may have occurred. Surprisingly, a composite active site in HbHNL maintained (S)-selectivity, while the identical set of active site residues in AtHNL maintained (R)-selectivity. Composite active-site mutants that are (S)-selective without the Lys236 and Thr11 that are required for the classical (S)-HNL mechanism suggests a new mechanism. Modeling suggested a possibility for this new mechanism that does not exist in modern enzymes. Thus, the last common ancestor of HbHNL and AtHNL may have used an extinct mechanism, not the AtHNL-like mechanism. Multiple mechanisms are possible with the same catalytic residues and residues outside the active site strongly influence mechanism and enantioselectivity.

Evolution of a Catalytic Mechanism

The means by which superfamilies of specialized enzymes arise by gene duplication and functional divergence are poorly understood. The escape from adaptive conflict hypothesis, which posits multiple copies of a gene encoding a primitive inefficient and highly promiscuous generalist ancestor, receives support from experiments showing that resurrected ancestral enzymes are indeed more substrate-promiscuous than their modern descendants. Here, we provide evidence in support of an alternative model, the innovation-amplification-divergence hypothesis, which posits a single-copied ancestor as efficient and specific as any modern enzyme. We argue that the catalytic mechanisms of plant esterases and descendent acetone cyanohydrin lyases are incompatible with each other (e.g., the reactive substrate carbonyl must bind in opposite orientations in the active site). We then show that resurrected ancestral plant esterases are as catalytically specific as modern esterases, that the ancestor of modern acetone cyanohydrin lyases was itself only very weakly promiscuous, and that improvements in lyase activity came at the expense of esterase activity. These observations support the innovation-amplification-divergence hypothesis, in which an ancestor gains a weak promiscuous activity that is improved by selection at the expense of the ancestral activity, and not the escape from adaptive conflict in which an inefficient generalist ancestral enzyme steadily loses promiscuity throughout the transition to a highly active specialized modern enzyme.

Probing the Catalytic Mechanism Involved in the Isocitrate Lyase Superfamily: Hybrid Quantum Mechanical/Molecular Mechanical Calculations on 2,3-Dimethylmalate Lyase

The isocitrate lyase (ICL) superfamily catalyzes the cleavage of the C(2)-C(3) bond of various α-hydroxy acid substrates. Members of the family are found in bacteria, fungi, and plants and include ICL itself, oxaloacetate hydrolase (OAH), 2-methylisocitrate lyase (MICL), and (2R,3S)-dimethylmalate lyase (DMML) among others. ICL and related targets have been the focus of recent studies to treat bacterial and fungal infections, including tuberculosis. The catalytic process by which this family achieves C(2)-C(3) bond breaking is still not clear. Extensive structural studies have been performed on this family, leading to a number of plausible proposals for the catalytic mechanism. In this paper, we have applied quantum mechanical/molecular mechanical (QM/MM) methods to the most recently reported family member, DMML, to assess whether any of the mechanistic proposals offers a clear energetic advantage over the others. Our results suggest that Arg161 is the general base in the reaction and Cys124 is the general acid, giving rise to a rate-determining barrier of approximately 10 kcal/mol.