A Dual-Focus Workflow for Simultaneously Engineering High Activity and Thermal Stability in Methyl Parathion Hydrolase

Industrial fermentation applications typically require enzymes that exhibit high stability and activity at high temperatures. However, efforts to simultaneously improve these properties are usually limited by a trade-off between stability and activity. This report describes a computational strategy to enhance both activity and thermal stability of the mesophilic organophosphate-degrading enzyme, methyl parathion hydrolase (MPH). To predict hotspot mutation sites, we assembled a library of features associated with the target properties for each residue and then prioritized candidate sites by hierarchical clustering. Subsequent in silico screening with multiple algorithms to simulate selective pressures yielded a subset of 23 candidate mutations. Iterative parallel screening of mutations that improved thermal stability and activity yielded, MPHase-m5b, which exhibited 13.3 °C higher Tm and 4.2 times higher catalytic activity than wild-type (WT) MPH over a wide temperature range. Systematic analysis of crystal structures, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations revealed a wider entrance to the active site that increased substrate access with an extensive network of interactions outside the active site that reinforced αβ/βα sandwich architecture to improve thermal stability. This study thus provides an advanced, rational design framework to improve efficiency in engineering highly active, thermostable biocatalysts for industrial applications.

Evolving ω-amine transaminase AtATA guided by substrate-enzyme binding free energy for enhancing activity and stability against non-natural substrates

In the field of chiral amine synthesis, ω-amine transaminase (ω-ATA) is one of the most established enzymes capable of asymmetric amination under optimal conditions. However, the applicability of ω-ATA toward more non-natural complex molecules remains limited due to its low transamination activity, thermostability, and narrow substrate scope. Here, by employing a combined approach of computational virtual screening strategy and combinatorial active-site saturation test/iterative saturation mutagenesis strategy, we have constructed the best variant M14C3-V5 (M14C3-V62A-V116S-E117I-L118I-V147F) with improved ω-ATA from Aspergillus terreus (AtATA) activity and thermostability toward non-natural substrate 1-acetylnaphthalene, which is the ketone precursor for producing the intermediate (R)-(+)-1-(1-naphthyl)ethylamine [(R)-NEA] of cinacalcet hydrochloride, showing activity enhancement of up to 3.4-fold compared to parent enzyme M14C3 (AtATA-F115L-M150C-H210N-M280C-V149A-L182F-L187F). The computational tools YASARA, Discovery Studio, Amber, and FoldX were applied for predicting mutation hotspots based on substrate-enzyme binding free energies and to show the possible mechanism with features related to AtATA structure, catalytic activity, and stability in silico analyses. M14C3-V5 achieved 71.8% conversion toward 50 mM 1-acetylnaphthalene in a 50 mL preparative-scale reaction for preparing (R)-NEA. Moreover, M14C3-V5 expanded the substrate scope toward aromatic ketone compounds. The generated virtual screening strategy based on the changes in binding free energies has successfully predicted the AtATA activity toward 1-acetylnaphthalene and related substrates. Together with experimental data, these approaches can serve as a gateway to explore desirable performances, expand enzyme-substrate scope, and accelerate biocatalysis.IMPORTANCEChiral amine is a crucial compound with many valuable applications. Their asymmetric synthesis employing ω-amine transaminases (ω-ATAs) is considered an attractive method. However, most ω-ATAs exhibit low activity and stability toward various non-natural substrates, which limits their industrial application. In this work, protein engineering strategy and computer-aided design are performed to evolve the activity and stability of ω-ATA from Aspergillus terreus toward non-natural substrates. After five rounds of mutations, the best variant, M14C3-V5, is obtained, showing better catalytic efficiency toward 1-acetylnaphthalene and higher thermostability than the original enzyme, M14C3. The robust combinational variant acquired displayed significant application value for pushing the asymmetric synthesis of aromatic chiral amines to a higher level.

QM/MM and MM MD Simulations on Enzymatic Degradation of the Nerve Agent VR by Phosphotriesterase

V-type nerve agents are hardly degraded by phosphotriesterase (PTE). Interestingly, the PTE variant of BHR-73MNW can effectively improve the hydrolytic efficiency of VR, especially for its Sp-enantiomer. Here, the whole enzymatic degradation of both Sp and Rp enantiomers of VR by the wild-type PTE and its variant BHR-73MNW was investigated by quantum mechanics/molecular mechanics (QM/MM) calculations and MM molecular dynamics simulations. Present results indicate that the degradation of VR can be initiated by the nucleophilic attack of the bridging OH- and the zinc-bound water molecule. The QM/MM-predicted energy barriers for the hydrolytic process of Sp-VR are 19.8 kcal mol-1 by the variant with water as a nucleophile and 22.0 kcal mol-1 by the wild-type PTE with OH- as a nucleophile, and corresponding degraded products are bound to the dinuclear metal site in monodentate and bidentate coordination modes, respectively. The variant effectively increases the volume of the large pocket, allowing more water molecules to enter the active pocket and resulting in the improvement of the degradation efficiency of Sp-VR. The hydrolysis of Rp-VR is triggered only by the hydroxide with an energy span of 20.6 kcal mol-1 for the wild-type PTE and 20.7 kcal mol-1 for the variant BHR-73-MNW PTE. Such mechanistic insights into the stereoselective degradation of VR by PTE and the role of water may inspire further studies to improve the catalytic efficiency of PTE toward the detoxification of nerve agents.

Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1

Metalloproteinase-1 (MMP-1) catalyzed collagen degradation is essential for a wide variety of normal physiological processes, while at the same time contributing to several diseases in humans. Therefore, a comprehensive understanding of this process is of great importance. Although crystallographic and spectroscopic studies provided fundamental information about the structure and function of MMP-1, the precise mechanism of collagen degradation especially considering the complex and flexible structure of the substrate, remains poorly understood. In addition, how the protein environment dynamically reorganizes at the atomic scale into a catalytically active state capable of collagen hydrolysis remains unknown. In this study, we applied experimentally-guided multiscale molecular modeling methods including classical molecular dynamics (MD), well-tempered (WT) classical metadynamics (MetD), combined quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM MetD simulations to explore and characterize the early catalytic events of MMP-1 collagenolysis. Importantly the study provided a complete atomic and dynamic description of the transition from the open to the closed form of the MMP-1•THP complex. Notably, the formation of catalytically active Michaelis complex competent for collagen cleavage was characterized. The study identified the changes in the coordination state of the catalytic zinc(II) associated with the conformational transformation and the formation of catalytically productive ES complex. Our results confirm the essential role of the MMP-1 catalytic domain's α-helices (hA, hB and hC) and the linker region in the transition to the catalytically competent ES complex. Overall, the results provide unique mechanistic insight into the conformational transformations and associated changes in the coordination state of the catalytic zinc(II) that would be important for the design of effective MMP-1 inhibitors.

Water-Regulated Mechanisms for Degradation of Pesticides Paraoxon and Parathion by Phosphotriesterase: Insight from QM/MM and MD Simulations

The enzymatic degradation of pesticides paraoxon (PON) and parathion (PIN) by phosphotriesterase (PTE) has been investigated by QM/MM calculations and MD simulations. In the PTE-PON complex, Zn α and Zn β in the active site are five- and six-coordinated, respectively, while both zinc ions are six coordinated with the Zn α -bound water molecule (WT1) for the PTE-PIN system. The hydrolytic reactions for PON and PIN are respectively driven by the nucleophilic attack of the bridging-OH - and the Zn α -bound water molecule on the phosphorus center of substrate, and the two-step hydrolytic process is predicted to be the rate-limiting step with the energy spans of 13.8 and 14.4 kcal/mol for PON and PIN, respectively. The computational studies reveal that the presence of the Zn α -bound water molecule depends on the structural feature of substrates characterized by P=O and P=S, which determines the hydrolytic mechanism and efficiency for the degradation of organophosphorus pesticides by PTE.