Navigating the landscape of enzyme design: from molecular simulations to machine learning.

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

ALDELE: All-Purpose Deep Learning Toolkits for Predicting the Biocatalytic Activities of Enzymes.

Rapidly predicting enzyme properties for catalyzing specific substrates is essential for identifying potential enzymes for industrial transformations. The demand for sustainable production of valuable industry chemicals utilizing biological resources raised a pressing need to speed up biocatalyst screening using machine learning techniques. In this research, we developed an all-purpose deep-learning-based multiple-toolkit (ALDELE) workflow for screening enzyme catalysts. ALDELE incorporates both structural and sequence representations of proteins, alongside representations of ligands by subgraphs and overall physicochemical properties. Comprehensive evaluation demonstrated that ALDELE can predict the catalytic activities of enzymes, and particularly, it identifies residue-based hotspots to guide enzyme engineering and generates substrate heat maps to explore the substrate scope for a given biocatalyst. Moreover, our models notably match empirical data, reinforcing the practicality and reliability of our approach through the alignment with confirmed mutation sites. ALDELE offers a facile and comprehensive solution by integrating different toolkits tailored for different purposes at affordable computational cost and therefore would be valuable to speed up the discovery of new functional enzymes for their exploitation by the industry.

Human diet-derived polyphenolic compounds and hepatic diseases: From therapeutic mechanisms to clinical utilization.

This review focuses on the potential ameliorative effects of polyphenolic compounds derived from human diet on hepatic diseases. It discusses the molecular mechanisms and recent advancements in clinical applications. Edible polyphenols have been found to play a therapeutic role, particularly in liver injury, liver fibrosis, NAFLD/NASH, and HCC. In the regulation of liver injury, polyphenols exhibit anti-inflammatory and antioxidant effects, primarily targeting the TGF-ß, NF-κB/TLR4, PI3K/AKT, and Nrf2/HO-1 signaling pathways. In the regulation of liver fibrosis, polyphenolic compounds effectively reverse the fibrotic process by inhibiting the activation of hepatic stellate cells (HSC). Furthermore, polyphenolic compounds show efficacy against NAFLD/NASH by inhibiting lipid oxidation and accumulation, mediated through the AMPK, SIRT, and PPARγ pathways. Moreover, several polyphenolic compounds exhibit anti-HCC activity by suppressing tumor cell proliferation and metastasis. This inhibition primarily involves blocking Akt and Wnt signaling, as well as inhibiting the epithelial-mesenchymal transition (EMT). Additionally, clinical trials and nutritional evidence support the notion that certain polyphenols can improve liver disease and associated metabolic disorders. However, further fundamental research and clinical trials are warranted to validate the efficacy of dietary polyphenols.

Faculty standardized patients versus traditional teaching method to improve clinical competence among traditional Chinese medicine students: a prospective randomized controlled trial.

BACKGROUND: Standardized patients (SPs) simulation training models have been widely used in various fields, the study of using SPs in Traditional Chinese medicine (TCM) is still a new filed. Previous studies have demonstrated the effectiveness of occupational SP for TCM (OSP-TCM), which has an increasingly problem of high time and financial costs. The faculty SPs for TCM (FSP-TCM) simulation training model may provide a better alternative. This study aims to test and determine whether FSP-TCM simulations are more cost-effective than OSP-TCM and traditional educational models to improve the clinical competence of TCM students. METHODS: This study was a single-blind, prospective, randomized controlled trial conducted between February 2023 and October 2023. The participants were randomized into FSP-TCM group, OSP-TCM group and traditionally taught group (TT group) in the ratio of 1:1:1. The duration of this training program was 12 weeks (36 credit hours). Formative and summative assessments were integrated to evaluate the effectiveness of teaching and learning. Three distinct questionnaires were utilized to collect feedback from students, SPs, and teachers at the conclusion of the course. Additionally, analysis of cost comparisons between OSP-TCM and FSP-TCM were performed in the study. RESULTS: The study comprised a total of 90 students, with no dropouts during the research. In the formative evaluation, students assigned to both the FSP-TCM and OSP-TCM groups demonstrated higher overall scores compared to those in the TT group. Notably, their performance in "physical examination" (Pa = 0.01, Pb = 0.04, Pc = 0.93) and "comprehensive ability" (Pa = 0.01, Pb = 0.006, Pc = 0.96) significantly exceeded that of the TT group. In the summary evaluation, both SP-TCM groups students outperforms TT group in the online systematic knowledge test (Pa = 0.019, Pb = 0.04, Pc = 0.97), the application of TCM technology (Pa = 0.01, Pb = 0.03, Pc = 0.93) and real-time assessment (Pa= 0.003, Pb = 0.01, Pc = 0.93). The feedback questionnaire demonstrated that both SP-TCM groups showed higher levels of agreement for this course in "satisfaction with the course" (Pa = 0.03; Pb = 0.02) and "enhanced TCM clinical skills" (Pa = 0.02; Pb = 0.03) than TT group. The SP questionnaire showed that more FSPs than OSPs in "provided professional feedback" (FSPs: strongly agree 30%, agree 50% vs. OSPs: strongly agree 20%, agree 40%. P = 0.69), and in "gave hints" during the course (FSPs: strongly agree 10%, agree 30% vs. OSPs: strongly agree 0%, agree 10%. P = 0.42). It is noteworthy that FSP-TCM was significantly lower than the OSP-TCM in overall expense (FSP-TCM $7590.00 vs. OSP-TCM $17415.60), and teachers have a positive attitude towards the FSP-TCM. CONCLUSION: FSP-TCM training mode showed greater effectiveness than traditional teaching method in improving clinical competence among TCM students. It was feasible, practical, and cost-effective, and may serve as an alternative method to OSP-TCM simulation.

Origin of the enantioselectivity of alcohol dehydrogenase.

Alcohol dehydrogenases (ADH) are a family of enzymes that catalyse the interconversion between ketones/aldehydes and alcohols in the presence of NADPH cofactor. It is challenging to desymmetrise the substituted cyclopentane-1,3-dione by engineering an ADH, while the reaction mechanism of the metal independent ADH remains elusive. Here we measured the conversion of a model substrate 2-benzyl-2-methylcyclopentane-1,3-dione by LbADH and found it predominately gave the (2R,3R) product. Binding mode analysis of the substrate in LbADH from molecular dynamics simulations disclosed the origin of the enantioselectivity of the enzyme; the opening and closing of the loop 191-205 above the substrate are responsible for shaping the binding pocket to orientate the substrate, so as to give different stereoisomer products. Using QM/MM calculations, we elucidated the reaction mechanism of LbADH. Furthermore, we demonstrated the reaction profile corresponding to the production of different stereoisomers, which is in accordance with our experimental observations. This research here will shed a light on the rational engineering of ADH to achieve stereodivergent stereoisomer products.

Taming Tris(bipyridine)ruthenium(II) and Its Reactions in Water by Capture/Release with Shape-Switchable Symmetry-Matched Cyclophanes.

Electron/proton transfers in water proceeding from ground/excited states are the elementary reactions of chemistry. These reactions of an iconic class of molecules─polypyridineRu(II)─are now controlled by capturing or releasing three of them with hosts that are shape-switchable. Reversible erection or collapse of the host walls allows such switchability. Some reaction rates are suppressed by factors of up to 120 by inclusive binding of the metal complexes. This puts nanometric coordination chemistry in a box that can be open or shut as necessary. Such second-sphere complexation can allow considerable control to be exerted on photocatalysis, electrocatalysis, and luminescent sensing involving polypyridineRu(II) compounds. The capturing states of hosts are symmetry-matched to guests for selective binding and display submicromolar affinities. A perching complex, which is an intermediate state between capturing and releasing states, is also demonstrated.

Thebaine is Selectively Demethylated by Thebaine 6-O-Demethylase and Codeine-3-O-demethylase at Distinct Binding Sites: A Computational Study.

Two homologous 2-oxoglutarate-dependent (ODD) nonheme enzymes thebaine 6-O-demethylase (T6ODM) and codeine-3-O-demethylase (CODM), are involved in the morphine biosynthesis pathway from thebaine, catalyzing the O-demethylation reaction with precise regioselectivity at C6 and C3 positions of thebaine respectively. We investigated the origin of the regioselectivity of these enzymes by combined molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations and found that Thebaine binds at the two distinct sites of T6ODM and CODM, which determines the regioselectivity of the enzymes. A remarkable oxo rotation is observed in the decarboxylation process. Starting from the closed pentacoordinate configuration, the C-terminal lid adopts an open conformation in the octahedral Fe(IV) = O complex to facilitate the subsequent demethylation. Phe241 and Phe311 stabilize the substrate in the binding pocket, while Arg219 acts as a gatekeeper residue to stabilize the substrate. Our results unravel the regioselectivity in 2-OG dependent nonheme enzymes and may shed light for exploring the substrate scope of these enzymes and developing novel biotechnology for morphine biosynthesis.

The phosphorylation mechanism of mevalonate diphosphate decarboxylase: a QM/MM study.

Mevalonate diphosphate decarboxylase (MDD) catalyses a crucial step of the mevalonate pathway via Mg2+-ATP-dependent phosphorylation and decarboxylation reactions to ultimately produce isopentenyl diphosphate, the precursor of isoprenoids, which is essential to bacterial functions and provides ideal building blocks for the biosynthesis of isopentenols. However, the metal ion(s) in MDD has not been unambiguously resolved, which limits the understanding of the catalytic mechanism and the exploitation of enzymes for the development of antibacterial therapies or the mevalonate metabolic pathway for the biosynthesis of biofuels. Here by analogizing structurally related kinases and molecular dynamics simulations, we constructed a model of the MDD-substrate-ATP-Mg2+ complex and proposed that MDD requires two Mg2+ ions for maintaining a catalytically active conformation. Subsequent QM/MM studies indicate that MDD catalyses the phosphorylation of its substrate mevalonate diphosphate (MVAPP) via a direct phosphorylation reaction, instead of the previously assumed catalytic base mechanism. The results here would shed light on the active conformation of MDD-related enzymes and their catalytic mechanisms and therefore be useful for developing novel antimicrobial therapies or reconstructing mevalonate metabolic pathways for the biosynthesis of biofuels.

Regio- and stereoselectivity in the CYP450BM3-catalyzed hydroxylation of complex terpenoids: a QM/MM study.

The site-selective C-H oxidation of terpenoids by P450 attracts great attention because of their wide range of biological activities. However, the binding and catalytic mechanism of P450 for the hydroxylation of complex terpenoid substrates remains elusive, which has limited the rational engineering of P450 as a biocatalyst for terpenoid biosynthesis. Here, we studied the origin of the selectivity and reactivity of P450BM3 in the hydroxylation of terpenoids by combining molecular dynamics simulations and QM/MM calculations, using artemisinin as a model compound. We found that the conformational change of the ß1 sheet at the substrate entrance and the displacement of the ß' helix were critical for reshaping the binding pocket to modulate substrate entrance and positioning the C-H to be activated toward the oxidative species of P450 for the subsequent hydrogen abstraction, the rate-determining step of hydroxylation. There is a distinct linear correlation between activation barriers and reaction coordinates, indicating that reaction coordinates can be used as a facile descriptor for predicting the reactivity of P450BM3. These findings would provide valuable guidance for predicting the selectivity and reactivity of P450BM3 for the selective hydroxylation of non-native terpenoid substrates so as to prioritize the rationally designed enzymes for terpenoid biosynthesis.

Determining the Online Measurable Input Variables in Human Joint Moment Intelligent Prediction Based on the Hill Muscle Model.

Introduction: Human joint moment is a critical parameter to rehabilitation assessment and human-robot interaction, which can be predicted using an artificial neural network (ANN) model. However, challenge remains as lack of an effective approach to determining the input variables for the ANN model in joint moment prediction, which determines the number of input sensors and the complexity of prediction. Methods: To address this research gap, this study develops a mathematical model based on the Hill muscle model to determining the online input variables of the ANN for the prediction of joint moments. In this method, the muscle activation, muscle-tendon moment velocity and length in the Hill muscle model and muscle-tendon moment arm are translated to the online measurable variables, i.e. muscle electromyography (EMG), joint angles and angular velocities of the muscle span. To test the predictive ability of these input variables, an ANN model is designed and trained to predict joint moments. The ANN model with the online measurable input variables is tested on the experimental data collected from ten healthy subjects running with the speeds of 2, 3, 4 and 5 m/s on a treadmill. The variance accounted for (VAF) between the predicted and inverse dynamics moment is used to evaluate the prediction accuracy. Results: The results suggested that the method can predict joint moments with a higher accuracy (mean VAF = 89.67±5.56 %) than those obtained by using other joint angles and angular velocities as inputs (mean VAF = 86.27±6.6%) evaluated by jack-knife cross-validation. Conclusions: The proposed method provides us with a powerful tool to predict joint moment based on online measurable variables, which establishes the theoretical basis for optimizing the input sensors and detection complexity of the prediction system. It may facilitate the research on exoskeleton robot control and real-time gait analysis in motor rehabilitation.

Stereodivergent Protein Engineering of a Lipase To Access All Possible Stereoisomers of Chiral Esters with Two Stereocenters.

Enzymatic stereodivergent synthesis to access all possible product stereoisomers bearing multiple stereocenters is relatively undeveloped, although enzymes are being increasingly used in both academic and industrial areas. When two stereocenters and thus four stereoisomeric products are involved, obtaining stereodivergent enzyme mutants for individually accessing all four stereoisomers would be ideal. Although significant success has been achieved in directed evolution of enzymes in general, stereodivergent engineering of one enzyme into four highly stereocomplementary variants for obtaining the full complement of stereoisomers bearing multiple stereocenters remains a challenge. Using Candida antarctica lipase B (CALB) as a model, we report the protein engineering of this enzyme into four highly stereocomplementary variants needed for obtaining all four stereoisomers in transesterification reactions between racemic acids and racemic alcohols in organic solvents. By generating and screening less than 25 variants of each isomer, we achieved >90% selectivity for all of the four possible stereoisomers in the model reaction. This difficult feat was accomplished by developing a strategy dubbed "focused rational iterative site-specific mutagenesis" (FRISM) at sites lining the enzyme's binding pocket. The accumulation of single mutations by iterative site-specific mutagenesis using a restricted set of rationally chosen amino acids allows the formation of ultrasmall mutant libraries requiring minimal screening for stereoselectivity. The crystal structure of all stereodivergent CALB variants, flanked by MD simulations, uncovered the source of selectivity.

Catalytic mechanism of mevalonate kinase revisited, a QM/MM study.

Mevalonate Kinase (MVK) catalyses the ATP-Mg2+ mediated phosphate transfer of mevalonate to produce mevalonate 5-phosphate and is a key kinase in the mevalonate pathway in the biosynthesis of isopentenyl diphosphate, the precursor of isoprenoid-based biofuels. However, the crystal structure in complex with the native substrate mevalonate, ATP and Mg2+ has not been resolved, which has limited the understanding of its reaction mechanism and therefore its application in the production of isoprenoid-based biofuels. Here using molecular docking, molecular dynamics (MD) simulations and a hybrid QM/MM study, we revisited the location of Mg2+ resolved in the crystal structure of MVK and determined a catalytically competent MVK structure in complex with the native substrate mevalonate and ATP. We demonstrated that significant conformational change on a flexible loop connecting the α6 and α7 helix is induced by the substrate binding. Further, we found that Asp204 is coordinated to the Mg2+ ion. Arg241 plays a crucial role in organizing the triphosphoryl tail of ATP for in-line phosphate transfer and stabilizing the negative charge that accumulates at the ß,γ-bridging oxygen of ATP upon bond cleavage. Remarkably, we revealed that the phosphorylation of mevalonate catalyzed by MVK occurs via a direct phosphorylation mechanism, instead of the conventionally postulated catalytic base mechanism. The catalytically competent complex structure of MVK as well as the mechanism of reaction will pave the way for the rational engineering of MVK to exploit its applications in the production of biofuels.

Are the Apo Proteins Suitable for the Rational Discovery of Allosteric Drugs?

Allosteric modulators, by targeting the less-conserved allosteric sites, represent an innovative strategy in drug discovery. These modulators have a distinctive advantage over orthosteric ligands that attach to the conserved, functional orthosteric sites. However, in structure-based drug design, it remains unclear whether allosteric protein structures determined without orthosteric ligand binding are suitable for allosteric drug screening. In this study, we performed large-scale conformational samplings of six representative allosteric proteins uncomplexed ( apo) and complexed ( holo) with orthosteric ligands to explore the effect of orthosteric site binding on the conformational dynamics of allosteric sites. The results, coupled with the redocking evaluation of allosteric modulators to their apo and holo proteins using their MD trajectories, indicated that orthosteric site binding had an effect on the dynamics of the allosteric sites and allosteric modulators preferentially bound to their holo proteins. According to the analysis data, we constructed a new correlation model for quantifying the allosteric site change driven by substrate binding to the orthosteric site. These results highlight the strong demand to select holo allosteric proteins as initial inputs in structure-based allosteric drug screening when the distance between orthosteric and allosteric sites in the protein is below 5 Å, which is expected to contribute to allosteric drug discovery.

Improving the Activity and Stability of Human Galactokinase for Therapeutic and Biotechnological Applications.

Galactokinase catalyses the site- and stereospecific phosphorylation of α-d-galactose. As such it has attracted interest as a biocatalyst for the introduction of phosphate groups into monosaccharides. However, attempts to broaden the substrate range of human galactokinase have generally resulted in substantially reduced activity. The enzyme also has biotechnological potential in enzyme replacement therapy (ERT) for type II galactosaemia. The return-to-consensus approach can be used to identify residues that can be altered to increase protein stability and enzyme activity. This approach identified six residues of potential interest in human galactokinase. Some of the single consensus variants (M60V, D268E, A334S and G373S) increased the catalytic turnover of the enzyme, but none resulted in improved stability. When all six changes were introduced into the protein (M60V/M180V/D268E/A334S/R366Q/G373S), thermal stability was increased. Molecular dynamics simulations suggested that these changes altered the protein's conformation at key sites. The number of salt bridges and hydrogen bonds was also increased. Combining the six consensus variations with Y379W (a variant with greater substrate promiscuity) increased the stability of this variant and its turnover towards some substrates. Thus, the six consensus variants can be used to stabilise catalytically interesting variants of human galactokinase and might also be useful if the protein were to be used in ERT.

Water-mediated network in the resistance mechanism of fosfomycin.

Fosfomycin Resistance Kinase A (FomA) catalyzes the phosphorylation of fosfomycin, which is an effective antibiotic for treating urinary tract infections. Understanding the chemical reaction mechanism is essential for developing strategies to counter the resistance of fosfomycin in clinical settings. Here the catalytic mechanism of FomA was investigated using molecular dynamic simulations in conjunction with quantum mechanics/molecular mechanics calculations (B97d/AMBER99). Our QM/MM study disclosed that the phosphorylation reaction catalyzed by FomA follows a dissociative mechanism, in contrast to the previously proposed associative mechanism. In addition, we found that His58, a characteristic residue in the AAK family, plays a key role in positioning the phosphate group of fosfomycin in the transition state. Molecular dynamic simulations revealed the important roles of Lys9 and Lys18 in arranging the nucleotide for phosphate transfer. Furthermore, we identified a four-membered water network mediated by Asp171 and Ser13 that is critical in ordering ATP for phosphate transfer. The active structure and reaction mechanism of FomA will provide valuable insights for developing new strategies to tackle the resistance to Fosfomycin-based antibiotic therapies.

Spatial requirement for PAMO for transformation of non-native linear substrates.

Phenylacetone monooxygenase is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenases family, and therefore it is an ideal candidate for the synthesis of industrially relevant ester or lactone compounds. However, its limited substrate scope has largely limited its industrial applications. Linear substrates are interesting from an industrial point of view, it is thus necessary to identify the essential spatial requirement for achieving high conversions for non-native linear substrates. Here using molecular dynamics simulations, we compared the conversion of a non-native linear substrate 2-octanone and the native substrate phenylacetone, catalyzed by the WT enzyme and a quadruple variant P253F/G254A/R258M/L443F that exhibits significantly improved activity towards 2-octanone. We uncovered that a remarkable movement of L289 is crucial for a reshaping of the active site of the quadruple variant so as to prevent the aliphatic substrate from moving away from the C4a-peroxyflavin, thus enabling it to keep a catalytically relevant pose during the oxygenation process. By performing steady-state kinetic analysis of two single-mutation variants at position 258, we further validated that the L289 reposition is attributed to the combined effect of quadruple mutations. In order to further explore the substrate scope of PAMO we also studied the binding of cyclopentanone and 2-phenylcyclohexanone, which are the typical substrates of CPMO in group I and CHMO in group III, respectively. Our study provides fundamental atomic-level insights in rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Modulation of the mobility of a key region in human galactokinase: Impacts on catalysis and stability.

Galactokinase catalyses the phosphorylation of α-d-galactose and some structurally related monosaccharides. The enzyme is of interest due to its potential as a biocatalyst for the production of sugar 1-phosphates and due to its involvement in the inherited metabolic disease type II galactosemia. It has been previously shown that a region (residues 231-245) in human galactokinase often has altered mobility when active site residues are varied. We hypothesised that the reverse may be true and that designing changes to this region might affect the functioning of the active site of the enzyme. Focussing on four residues (Leu-231, Gln-242, Glu-244 and Glu-245) we conducted molecular dynamics simulations to explore the effects of changing these residues to glycine or serine. In most cases the variations resulted in local changes to the 231-245 region and global changes to the root mean squared fluctuation (RMSF) of the protein. The four serine variants were expressed as recombinant proteins. All had altered steady state enzyme kinetic parameters with α-d-galactose as a substrate. However, these changes were generally less than ten-fold in magnitude. Changes were also observed with 2-deoxy-α-d-galactose, α-d-galactosamine and α-d-talose as substrates, including (in some cases) loss of detectable activity, suggesting that these variations can tune the specificity of the enzyme. This study demonstrates that activity and specificity of human galactokinase can be modulated by variations designed to affect active site flexibility. It is likely that this principle can be generalised to other enzymes.

Insight into the mechanism of galactokinase: Role of a critical glutamate residue and helix/coil transitions.

Galactokinase, the enzyme which catalyses the first committed step in the Leloir pathway, has attracted interest due to its potential as a biocatalyst and as a possible drug target in the treatment of type I galactosemia. The mechanism of the enzyme is not fully elucidated. Molecular dynamics (MD) simulations of galactokinase with the active site residues Arg-37 and Asp-186 altered predicted that two regions (residues 174-179 and 231-240) had different dynamics as a consequence. Interestingly, the same two regions were also affected by alterations in Arg-105, Glu-174 and Arg-228. These three residues were identified as important in catalysis in previous computational studies on human galactokinase. Alteration of Arg-105 to methionine resulted in a modest reduction in activity with little change in stability. When Arg-228 was changed to methionine, the enzyme's interaction with both ATP and galactose was affected. This variant was significantly less stable than the wild-type protein. Changing Glu-174 to glutamine (but not to aspartate) resulted in no detectable activity and a less stable enzyme. Overall, these combined in silico and in vitro studies demonstrate the importance of a negative charge at position 174 and highlight the critical role of the dynamics in to key regions of the protein. We postulate that these regions may be critical for mediating the enzyme's structure and function.