Recent Trends in Quantum Chemical Modeling of Enzymatic Reactions

Asymmetric photoenzymatic incorporation of fluorinated motifs into olefins

Enzymes capable of assimilating fluorinated feedstocks are scarce. This situation poses a challenge for the biosynthesis of fluorinated compounds used in pharmaceuticals, agrochemicals, and materials. We developed a photoenzymatic hydrofluoroalkylation that integrates fluorinated motifs into olefins. The photoinduced promiscuity of flavin-dependent ene-reductases enables the generation of carbon-centered radicals from iodinated fluoroalkanes, which are directed by the photoenzyme to engage enantioselectively with olefins. This approach facilitates stereocontrol through interaction between a singular fluorinated unit and the enzyme, securing high enantioselectivity at β, γ, or δ positions of fluorinated groups through enzymatic hydrogen atom transfer-a process that is notably challenging with conventional chemocatalysis. This work advances enzymatic strategies for integrating fluorinated chemical feedstocks and opens avenues for asymmetric synthesis of fluorinated compounds.

Efficient and easible biocatalysts: Strategies for enzyme improvement. A review

Owing to the environmental friendliness and vast advantages that enzymes offer in the biotechnology and industry fields, biocatalysts are a prolific investigation field. However, the low catalytic activity, stability, and specific selectivity of the enzyme limit the range of the reaction enzymes involved in. A comprehensive understanding of the protein structure and dynamics in terms of molecular details enables us to tackle these limitations effectively and enhance the catalytic activity by enzyme engineering or modifying the supports and solvents. Along with different strategies including computational, enzyme engineering based on DNA recombination, enzyme immobilization, additives, chemical modification, and physicochemical modification approaches can be promising for the wide spread of industrial enzyme usage. This is attributed to the successful application of biocatalysts in industrial and synthetic processes requires a system that exhibits stability, activity, and reusability in a continuous flow process, thereby reducing the production cost. The main goal of this review is to display relevant approaches for improving enzyme characteristics to overcome their industrial application.

Computational Study of the Fries Rearrangement Catalyzed by Acyltransferase from Pseudomonas protegens

The acyltransferase from Pseudomonas protegens (PpATase) catalyzes in nature the reversible transformation of monoacetylphloroglucinol to diacetylphloroglucinol and phloroglucinol. Interestingly, this enzyme has been shown to catalyze the promiscuous transformation of 3-hydroxyphenyl acetate to 2',4'-dihydroxyacetophenone, representing a biological version of the Fries rearrangement. In the present study, we report a mechanistic investigation of this activity of PpATase using quantum chemical calculations. A detailed mechanism is proposed, and the energy profile for the reaction is presented. The calculations show that the acylation of the enzyme is highly exothermic, while the acetyl transfer back to the substrate is only slightly exothermic. The deprotonation of the C6-H of the substrate is rate-limiting, and a remote aspartate residue (Asp137) is proposed to be the general base group in this step. Analysis of the binding energies of various acetyl acceptors shows that PpATase can promote both intramolecular and intermolecular Fries rearrangement towards diverse compounds.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Quantum chemical modeling of enantioselective sulfoxidation and epoxidation reactions by indole monooxygenase VpIndA1

Indole monooxygenases (IMOs) are enzymes from the family of Group E monooxygenases, requiring flavin adenine dinucleotide (FAD) for their activities. IMOs play important roles in both sulfoxidation and epoxidation reactions. The broad substrate range and high selectivity of IMOs make them promising biocatalytic tools for synthesizing chiral compounds. In the present study, quantum chemical calculations using the cluster approach were performed to investigate the reaction mechanism and the enantioselectivity of the IMO from Variovorax paradoxus EPS (VpIndA1). The sulfoxidation of methyl phenyl sulfide (MPS) and the epoxidation of indene were chosen as the representative reactions. The calculations confirmed that the FADOOH intermediate is the catalytic species in the VpIndA1 reactions. The oxidation of MPS adopts a one-step mechanism involving the direct oxygen-transfer from FADOOH to the substrate and the proton transfer from the -OH group back to FAD, while the oxidation of indene follows a stepwise mechanism involving a carbocation intermediate. It was computationally predicted that VpIndA1 prefers the formation of (S)-product for the MPS sulfoxidation and (1S,2R)-product for the indene epoxidation, consistent with the experimental observations. Importantly, the factors controlling the stereo-preference of the two reactions are identified. The findings in the present study provide valuable insights into the VpIndA1-catalyzed reactions, which are essential for the rational design of this enzyme and other IMOs for industrial applications. It is also worth emphasizing that the quantum chemical cluster approach is again demonstrated to be powerful in studying the enantioselectivity of enzymatic reactions.

Origin of Chemoselectivity of Halohydrin Dehalogenase-Catalyzed Epoxide Ring-Opening Reactions

By performing molecular dynamics (MD), quantum mechanical/molecular mechanical (QM/MM) calculations, and QM cluster calculations, the origin of chemoselectivity of halohydrin dehalogenase (HHDH)-catalyzed ring-opening reactions of epoxide with the nucleophilic reagent NO2- has been explored. Four possible chemoselective pathways were considered, and the computed results indicate that the pathway associated with the nucleophilic attack on the Cα position of epoxide by NO2- is most energetically favorable and has an energy barrier of 12.9 kcal/mol, which is close to 14.1 kcal/mol derived from experimental kinetic data. A hydrogen bonding network formed by residues Ser140, Tyr153, and Arg157 can strengthen the electrophilicity of the active site of the epoxide substrate to affect chemoselectivity. To predict the energy barrier trends of the chemoselective transition states, multiple analyses including distortion analysis and electrophilic Parr function (Pk+) analysis were carried out with or without an enzyme environment. The obtained insights should be valuable for the rational design of enzyme-catalyzed and biomimetic organocatalytic epoxide ring-opening reactions with special chemoselectivity.

Reduction of molecular oxygen in flavodiiron proteins - Catalytic mechanism and comparison to heme-copper oxidases

The family of flavodiiron proteins (FDPs) plays an important role in the scavenging and detoxification of both molecular oxygen and nitric oxide. Using electrons from a flavin mononucleotide cofactor molecular oxygen is reduced to water and nitric oxide is reduced to nitrous oxide and water. While the mechanism for NO reduction in FDPs has been studied extensively, there is very little information available about O2 reduction. Here we use hybrid density functional theory (DFT) to study the mechanism for O2 reduction in FDPs. An important finding is that a proton coupled reduction is needed after the O2 molecule has bound to the diferrous diiron active site and before the OO bond can be cleaved. This is in contrast to the mechanism for NO reduction, where both NN bond formation and NO bond cleavage occurs from the same starting structure without any further reduction, according to both experimental and computational results. This computational result for the O2 reduction mechanism should be possible to evaluate experimentally. Another difference between the two substrates is that the actual OO bond cleavage barrier is low, and not involved in rate-limiting the reduction process, while the barrier connected with bond cleavage/formation in the NO reduction process is of similar height as the rate-limiting steps. We suggest that these results may be part of the explanation for the generally higher activity for O2 reduction as compared to NO reduction in most FDPs. Comparisons are also made to the O2 reduction reaction in the family of heme‑copper oxidases.

Computation-Driven Rational Design of Self-Assembled Short Peptides for Catalytic Hydrogen Production

Self-assembling peptides represent a captivating area of study in nanotechnology and biomaterials. This interest is largely driven by their unique properties and the vast application potential across various fields such as catalytic functions. However, design complexities, including high-dimensional sequence space and structural diversity, pose significant challenges in the study of such systems. In this work, we explored the possibility of self-assembled peptides to catalyze the hydrolysis of hydrosilane for hydrogen production using ab initio calculations and carried out wet-lab experiments to confirm the feasibility of these catalytic reactions under ambient conditions. Further, we delved into the nuanced interplay between sequence, structural conformation, and catalytic activity by combining modeling with experimental techniques such as transmission electron microscopy and nuclear magnetic resonance and proposed a dual mode of the microstructure of the catalytic center. Our results reveal that although research in this area is still at an early stage, the development of self-assembled peptide catalysts for hydrogen production has the potential to provide a more sustainable and efficient alternative to conventional hydrogen production methods. In addition, this work also demonstrates that a computation-driven rational design supplemented by experimental validation is an effective protocol for conducting research on functional self-assembled peptides.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Biocatalytic strategy for the construction of sp3-rich polycyclic compounds from directed evolution and computational modelling

Catalysis with engineered enzymes has provided more efficient routes for the production of active pharmaceutical agents. However, the potential of biocatalysis to assist in early-stage drug discovery campaigns remains largely untapped. In this study, we have developed a biocatalytic strategy for the construction of sp3-rich polycyclic compounds via the intramolecular cyclopropanation of benzothiophenes and related heterocycles. Two carbene transferases with complementary regioisomer selectivity were evolved to catalyse the stereoselective cyclization of benzothiophene substrates bearing diazo ester groups at the C2 or C3 position of the heterocycle. The detailed mechanisms of these reactions were elucidated by a combination of crystallographic and computational analyses. Leveraging these insights, the substrate scope of one of the biocatalysts could be expanded to include previously unreactive substrates, highlighting the value of integrating evolutionary and rational strategies to develop enzymes for new-to-nature transformations. The molecular scaffolds accessed here feature a combination of three-dimensional and stereochemical complexity with 'rule-of-three' properties, which should make them highly valuable for fragment-based drug discovery campaigns.

Insight into chemically reactive metabolites of aliphatic amine pollutants: A de novo prediction strategy and case study of sertraline

The uncommon metabolic pathways of organic pollutants are easily overlooked, potentially leading to idiosyncratic toxicity. Prediction of their biotransformation associated with the toxic effects is the very purpose that this work focuses, to develop a de novo method to mechanistically predict the reactive toxicity pathways of uncommon metabolites from start aliphatic amine molecules, which employed sertraline triggered by CYP450 enzymes as a model system, as there are growing concerns about the effects on human health posed by antidepressants in the aquatic environment. This de novo prediction strategy combines computational and experimental methods, involving DFT calculations upon sequential growth, in vitro and in vivo assays, dissecting chemically reactive mechanism relevant to toxicity, and rationalizing the fundamental factors. Significantly, desaturation and debenzylation-aromatization as the emerging metabolic pathways of sertraline have been elucidated, with the detection of DNA adducts of oxaziridine metabolite in mice, highlighting the potential reactive toxicity. Molecular orbital analysis supports the reactivity preference for toxicological-relevant C-N desaturation over N-hydroxylation of sertraline, possibly extended to several other aliphatic amines based on the Bell-Evans-Polanyi principle. It was further validated toward some other wide-concerned aliphatic amine pollutants involving atrazine, ε-caprolactam, 6PPD via in silico and in vitro assays, thereby constituting a complete path for de novo prediction from case study to general applications.

Role of Explicit Solvation in Computational Modeling of Chemical Reactions: Mechanism of Cu(I) Transfer Between Thiolate-Based Chelators in Water

Solvation plays important roles in controlling the thermodynamic and kinetic aspects of chemical reactions. The conventional approaches to treat solvation via electronic structure methods are likely to become inadequate, when the reacting solutes have strong electrostatic and hydrogen bonding interactions with the solvent and undergo significant structural changes during the course of the reaction. In this article, we present evidence of such solvent and structural effects in the computational study of the Cu(I) transfer reaction between thiolate-based chelators dithiobutylamine (DTBA) and dithiotheritol (DTT) in water, inspired from biological copper trafficking phenomena. We propose a general solution to the problem by combining classical molecular dynamics (MD) simulations of the bulk system and static quantum chemistry calculations. The fluctuating solvation shell was estimated from MD, and energetics was assessed by averaging QM energies of a series of molecular clusters constructed from the MD snapshots. Applying this approach, we propose a reaction pathway with estimates of relative intermediate stabilities and barriers, which suggest the overall reaction to be reversible in nature and likely to go through both two and three coordinated intermediates, confirming previous studies of similar protein analogues. An interesting fact that emerged from our study was the strong indication that the rate-determining step is the deprotonation of initial thiol bound Cu(I) complex, without involving any Cu(I)-S bonds. The proposed method will lead to a better treatment of solvations, and these mechanistic insights will aid our understanding of biological copper(I) trafficking.

Simulating chemical reactions promoted by self-assembled peptides with catalytic properties

Peptides that self-assemble exhibit distinct three-dimensional structures and attributes, positioning them as promising candidates for biocatalysts. Exploring their catalytic processes enhances our comprehension of the catalytic actions inherent to self-assembling peptides, laying a theoretical foundation for creating novel biocatalysts. The investigation into the intricate reaction mechanisms of these entities is rendered challenging due to the vast variability in peptide sequences, their aggregated formations, supportive elements, structures of active sites, types of catalytic reactions, and the interplay between these variables. This complexity hampers the elucidation of the linkage between sequence, structure, and catalytic efficiency in self-assembling peptide catalysts. This chapter delves into the latest progress in understanding the mechanisms behind peptide self-assembly, serving as a catalyst in hydrolysis and oxidation reactions, and employing computational analyses. It discusses the establishment of models, selection of computational strategies, and analysis of computational procedures, emphasizing the application of modeling techniques in probing the catalytic mechanisms of peptide self-assemblies. It also looks ahead to the potential future trajectories within this research domain. Despite facing numerous obstacles, a thorough investigation into the structural and catalytic mechanisms of peptide self-assemblies, combined with the ongoing advancement in computational simulations and experimental methodologies, is set to offer valuable theoretical insights for the development of new biocatalysts, thereby significantly advancing the biocatalysis field.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

Proxy-approach in understanding the bisubstrate activity of strictosidine synthases

Besides tryptamine (1) and secologanin (2), non-cognate substrates also undergo a Pictet-Spengler reaction (PSR) catalyzed by strictosidine synthases (STR) with differing catalytic properties. We characterized the bisubstrate binding aspect of catalysis - order, affinity, and cooperativity - with STR orthologs from Rauvolfia serpentina (RsSTR) and Ophiorrhiza pumila (OpSTR) by an isothermal titration calorimetry (ITC) based 'proxy approach' that employed a non-reactive tryptamine analog (m1) to capture its inert ternary complexes with STRs and (2). ITC studies with OpSTR and (2) revealed 'tryptamine-first' cooperative binding with (1) and a simultaneous cooperative binding with (m1). Binding cooperativity among (m1) and (2) towards OpSTR was higher than RsSTR. Crystallographic study of RsSTR-(m1) complex helped to understand the unreactive binding of (m1) in terms of orientation and interactions in the RsSTR pocket. PSR with (m1) was revealed to be energetically unfeasible by the density functional theory (DFT) scans of the first hydrogen abstraction by RsSTR. The effect of pH on the bisubstrate binding to OpSTR was deciphered by molecular dynamics simulations (MDS), which also provided a molecular basis for the stability of complex of OpSTR with (m1) and (2). Therefore, we investigated STRs from a substrate binding perspective to inform drug-design and rational enzyme engineering efforts.

Molecular insight into the enhanced performance of CALB toward PBDF degradation

Poly(butylene diglycolate-co-furandicarboxylate) (PBDF) is a newly developed biodegradable copolyester. Candida antarctica lipase B (CALB) has been identified as an effective catalyst for PBDF degradation. The mechanism is elucidated using a combination of molecular dynamics simulations and quantum chemistry approaches. The findings unveil a four-step catalytic reaction pathway. Furthermore, bond analysis, charge and interaction analysis are conducted to gain a more comprehensive understanding of the PBDF degradation process. Additionally, through the introduction of single-point mutations to crucial residues in CALB's active sites, two mutants, T138I and D134I, are discovered to exhibit improved catalytic efficiency. These significant findings contribute to the advancement of our comprehension concerning the molecular mechanism of underlying copolyesters degradation, while also presenting a novel approach for expediting the degradation rate by the CALB enzyme modification.

Uncovering the Importance of Ligand Mobility on Cellular Uptake of Nanoparticles: Insights from Experimental, Computational, and Theoretical Investigations

The cellular uptake of nanoparticles (NPs) by biological cells is an important and fundamental process in drug delivery. Previous studies reveal that the physicochemical properties of nanoparticles as well as those of functionalized ligands can both critically affect the uptake behaviors. However, the effect of the conjugation strategy (i.e., the "bond" between the ligand and the NP) on the cellular uptake is overlooked and remains largely elusive. Here, by taking the broadly employed gold nanoparticle as an example, we comprehensively assessed the relationship between the conjugation strategy and uptake behaviors by introducing three ligands with the same functional terminal but different anchoring sites. As revealed by in vitro cell experiments and multiscale molecular simulations, the uptake efficiency of gold NPs was positively correlated with the strength of the "bond" and more specifically the ligand mobility on the NP surface. Moreover, we validated the results presented above by proposing a thermodynamic theory for the wrapping of NPs with mobile ligands. Further, we also showed that the endocytic pathway of NPs was highly dependent on ligand mobility. Overall, this study uncovered a vital role of conjugation strategy in the cellular uptake and may provide useful guidelines for tailoring the biobehaviors of nanoparticles.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

Enzymatic Synthesis of Noncanonical α-Amino Acids Containing γ-Tertiary Alcohols

Noncanonical amino acids (ncAAs) containing tertiary alcohols are valuable as precursors of natural products and active pharmaceutical ingredients. However, the assembly of such ncAA scaffolds from simple material by C-C bond formation remains a challenging task due to the presence of multiple stereocenters and large steric hindrance. In this study, we present a novel solution to this problem through highly selective enzymatic decarboxylative aldol addition. This method allows for the streamlined assembly of multifunctionalized ncAAs with γ-tertiary alcohols from readily available materials, such as L -aspartatic acid and isatins, vicinal diones and keto esters. The modularity of electrophiles furnished four classes of ncAAs with decent efficiency as well as excellent site and stereocontrol. Computational modeling was employed to gain detailed insight into the catalytic mechanism and to provide a rationale for the observed selectivities. The method offers a single-step approach to producing multifunctionalized ncAAs, which can be directly utilized in peptide synthesis and bioactivity assessment.

Remote stereocontrol with azaarenes via enzymatic hydrogen atom transfer

Strategies for achieving asymmetric catalysis with azaarenes have traditionally fallen short of accomplishing remote stereocontrol, which would greatly enhance accessibility to distinct azaarenes with remote chiral centres. The primary obstacle to achieving superior enantioselectivity for remote stereocontrol has been the inherent rigidity of the azaarene ring structure. Here we introduce an ene-reductase system capable of modulating the enantioselectivity of remote carbon-centred radicals on azaarenes through a mechanism of chiral hydrogen atom transfer. This photoenzymatic process effectively directs prochiral radical centres located more than six chemical bonds, or over 6 Å, from the nitrogen atom in azaarenes, thereby enabling the production of a broad array of azaarenes possessing a remote γ-stereocentre. Results from our integrated computational and experimental investigations underscore that the hydrogen bonding and steric effects of key amino acid residues are important for achieving such high stereoselectivities.

Enzyme Engineering Renders Chlorinase the Activity of Fluorinase

Organofluorine compounds have attracted substantial attention owing to their wide application in agrochemistry. Fluorinase (FlA) is a unique enzyme in nature that can incorporate fluorine into an organic molecule. Chlorinase (SalL) has a similar mechanism as fluorinase and can use chloride but not fluoride as a substrate to generate 5'-chloro-deoxyadenosine (5'-ClDA) from S-adenosyl-l-methionine (SAM). Therefore, identifying the features that lead to this selectivity for halide ions is highly important. Here, we engineered SalL to gain the function of FlA. We found that residue Tyr70 plays a key role in this conversion through alanine scanning. Site-saturation mutagenesis experiments demonstrated that Y70A/C/S/T/G all exhibited obvious fluorinase activity. The G131S mutant of SalL, in which the previously thought crucial residue Ser158 for fluoride binding in FlA was introduced, did not exhibit fluorination activity. Compared with the Y70T single mutant, the double mutant Y70T/W129F increased 5'-fluoro-5-deoxyadenosine (5'-FDA) production by 76%. The quantum mechanics (QM)/molecular mechanics (MM) calculations suggested that the lower energy barriers and shorter nucleophilic distance from F- to SAM in the mutants than in the SalL wild-type may contribute to the activity. Therefore, our study not only renders SalL the activity of FlA but also sheds light on the enzyme selectivity between fluoride versus chloride.

Catalytic Effects of Active Site Conformational Change in the Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Imidazole glycerol phosphate synthase (IGPS) is a class-I glutamine amidotransferase (GAT) that hydrolyzes glutamine. Ammonia is produced and transferred to a second active site, where it reacts with N1-(5'-phosphoribosyl)-formimino-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) to form precursors to purine and histidine biosynthesis. Binding of PrFAR over 25 Å away from the active site increases glutaminase efficiency by ∼4500-fold, primarily altering the glutamine turnover number. IGPS has been the focus of many studies on allosteric communication; however, atomic details for how the glutamine hydrolysis rate increases in the presence of PrFAR are lacking. We present a density functional theory study on 237-atom active site cluster models of IGPS based on crystallized structures representing the inactive and allosterically active conformations and investigate the multistep reaction leading to thioester formation and ammonia production. The proposed mechanism is supported by similar, well-studied enzyme mechanisms, and the corresponding energy profile is consistent with steady-state kinetic studies of PrFAR + IGPS. Additional active site models are constructed to examine the relationship between active site structural change and transition-state stabilization via energy decomposition schemes. The results reveal that the inactive IGPS conformation does not provide an adequately formed oxyanion hole structure and that repositioning of the oxyanion strand relative to the substrate is vital for a catalysis-competent oxyanion hole, with or without the hVal51 dihedral flip. These findings are valuable for future endeavors in modeling the IGPS allosteric mechanism by providing insight into the atomistic changes required for rate enhancement that can inform suitable reaction coordinates for subsequent investigations.

Binding of Per- and Polyfluoroalkyl Substances (PFAS) to the PPARγ/RXRα-DNA Complex

Nuclear receptors are the fundamental building blocks of gene expression regulation and the focus of many drug targets. While binding to DNA, nuclear receptors act as transcription factors, governing a multitude of functions in the human body. Peroxisome proliferator-activator receptor γ (PPARγ) and the retinoid X receptor α (RXRα) form heterodimers with unique properties and have a primordial role in insulin sensitization. This PPARγ/RXRα heterodimer has been shown to be impacted by per- and polyfluoroalkyl substances (PFAS) and linked to a variety of significant health conditions in humans. Herein, a selection of the most common PFAS (legacy and emerging) was studied utilizing molecular dynamics simulations for PPARγ/RXRα. The local and global structural effects of PFAS binding on the known ligand binding pockets of PPARγ and RXRα as well as the DNA binding domain (DBD) of RXRα were inspected. The binding free energies were predicted computationally and were compared between the different binding pockets. In addition, two electronic structure approaches were utilized to model the interaction of PFAS within the DNA binding domain, density functional theory (DFT) and domain-based pair natural orbital coupled cluster with perturbative triples (DLPNO-CCSD(T)) approaches, with implicit solvation. Residue decomposition and hydrogen-bonding analysis were also performed, detailing the role of prominent residues in molecular recognition. The role of l-carnitine is explored as a potential in vivo remediation strategy for PFAS interaction with the PPARγ/RXRα heterodimer. In this work, it was found that PFAS can bind and act as agonists for all of the investigated pockets. For the first time in the literature, PFAS are postulated to bind to the DNA binding domain in a nonspecific manner. In addition, for the PPARγ ligand binding domain, l-carnitine shows promise in replacing smaller PFAS from the pocket.

Electronic origins of the stereochemistry in β-lactam formed through the Staudinger reaction catalyzed by a nucleophile

This paper evaluates the electronic effects of molecular substituents on the stereoselectivity of the umpolung Staudinger catalytic reaction. This is especially important because experimental studies on constructing the β-lactam ring, a core structure of most antibiotics, through catalyzed Staudinger reactions have been massively progressing over the last century. Yet, there is a necessity for an in-depth understanding of the reaction mechanisms to help chemists, working on the well-established discoveries, improve these to optimize the stereoselectivity and yield of synthetic methods. Access to practical and effective advancements in forming optically pure β-lactam is paramount in the field of medical chemistry. This paper specifically investigates how changing the N-protecting group in the imine fragment can switch the stereoselectivity of the PPY-catalyzed Staudinger reaction. To do so, we employed the density functional theory (DFT) for geometry optimization and electronic analysis at the B3LYP/6-31G(d) level of theory to examine and compare the role of N-tosyl (N-Ts) and N-triflyl (N-Tf) imine on the mechanism pathways, i.e., imine-first or ketene-first, and stereochemistry of the reaction, i.e., cis or trans β-lactam. Our results show that the reaction mechanism pathway cannot be simply switched from ketene-first to imine-first by changing the substituent on the imine nitrogen atom, which is contrary to the reported experimental results, and both imines go through the ketene-first mechanism with different stereochemistries, which is cis selective for imine-Ts and trans selective for imine-Tf. Based on electronic analyses, the reversal in diastereoselectivity in the N-triflyl imine system could be attributed to the charge transfers and electron-density distribution over the transition states. Therefore, the cis/trans selectivity of the PPY-catalyzed Staudinger reaction could be effectively controlled by the electronic characteristics of the molecular substituents in the reactants. A N-protecting group in imine with a more electron-withdrawing nature seems to accelerate the stereo-determining step, ring closure, and increase the stabilization charge transfers in the transition state, leading to a preference for trans β-lactam formation. It seems that using a N-substituent with a higher electron-withdrawing nature can initially activate the imine by the nucleophilic catalyst in competition with ketene, i.e., imine-first versus ketene-first. These results can provide an insight into select proper substituents for the fragments to synthesis β-lactam with a desired stereochemistry. Also, a comprehensive comparison was performed between calculations with and without dispersion.

Stimuli-responsive probes for amplification-based imaging of miRNAs in living cells

MicroRNAs (miRNAs) have emerged as important biomarkers in biomedicine and bioimaging due to their roles in various physiological and pathological processes. Real-time and in situ monitoring of dynamic fluctuation of miRNAs in living cells is crucial for understanding these processes. However, current miRNA imaging probes still have some limitations, including the lack of effective amplification methods for low abundance miRNAs bioanalysis and uncontrollable activation, leading to background signals and potential false-positive results. Therefore, researchers have been integrating activatable devices with miRNA amplification techniques to design stimuli-responsive nanoprobes for "on-demand" and precise imaging of miRNAs in living cells. In this review, we summarize recent advances of stimuli-responsive probes for the amplification-based imaging of miRNAs in living cells and discuss the future challenges and opportunities in this field, aiming to provide valuable insights for accurate disease diagnosis and monitoring.

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Role of mutations in a chemoenzymatic enantiodivergent C(sp3)-H insertion: exploring the mechanism and origin of stereoselectivity

New-to-nature enzymes have emerged as powerful catalysts in recent years for streamlining various stereoselective organic transformations. While synthetic strategies employing engineered enzymes have witnessed proliferating success, there is limited clarity on the mechanistic front and more so when considering molecular-level insights into the role of selected mutations, dramatically escalating catalytic competency and selectivity. We have investigated the mechanism and correlation between mutations and exquisite stereoselectivity of a lactone carbene insertion into the C(sp3)-H bond of substituted aniline, catalyzed by two mutants of a cytochrome P450 variant, "P411" (engineered through directed evolution) in which the axial cysteine has been mutated to serine, utilizing various computational tools. The pivotal role of S264 and L/R328 mutations in the active site has been delineated computationally using two cluster models, thus rationalizing the enantiodivergence. This report provides much-needed insights into the origin of enantiodivergence, furnishing a mechanistic framework for understanding the anchoring effects of H-bond donor residues with the lactone ring. This study is expected to have important implications in the rational design of stereodivergent enzymes and toward successful in silico enzyme designing.

Prediction of Enzyme Catalysis by Computing Reaction Energy Barriers via Steered QM/MM Molecular Dynamics Simulations and Machine Learning

The prediction of enzyme activity is one of the main challenges in catalysis. With computer-aided methods, it is possible to simulate the reaction mechanism at the atomic level. However, these methods are usually expensive if they are to be used on a large scale, as they are needed for protein engineering campaigns. To alleviate this situation, machine learning methods can help in the generation of predictive-decision models. Herein, we test different regression algorithms for the prediction of the reaction energy barrier of the rate-limiting step of the hydrolysis of mono-(2-hydroxyethyl)terephthalic acid by the MHETase ofIdeonella sakaiensis. As a training data set, we use steered quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulation snapshots and their corresponding pulling work values. We have explored three algorithms together with three chemical representations. As an outcome, our trained models are able to predict pulling works along the steered QM/MM MD simulations with a mean absolute error below 3 kcal mol-1 and a score value above 0.90. More challenging is the prediction of the energy maximum with a single geometry. Whereas the use of the initial snapshot of the QM/MM MD trajectory as input geometry yields a very poor prediction of the reaction energy barrier, the use of an intermediate snapshot of the former trajectory brings the score value above 0.40 with a low mean absolute error (ca. 3 kcal mol-1). Altogether, we have faced in this work some initial challenges of the final goal of getting an efficient workflow for the semiautomatic prediction of enzyme-catalyzed energy barriers and catalytic efficiencies.

Nucleolar Essential Protein 1 (Nep1): Elucidation of enzymatic catalysis mechanism by molecular dynamics simulation and quantum mechanics study

The Nep1 protein is essential for the formation of eukaryotic and archaeal small ribosomal subunits, and it catalyzes the site-directed SAM-dependent methylation of pseudouridine (Ψ) during pre-rRNA processing. It possesses a non-trivial topology, namely, a 31 knot in the active site. Here, we address the issue of seemingly unfeasible deprotonation of Ψ in Nep1 active site by a distant aspartate residue (D101 in S. cerevisiae), using a combination of bioinformatics, computational, and experimental methods. We identified a conserved hydroxyl-containing amino acid (S233 in S. cerevisiae, T198 in A. fulgidus) that may act as a proton-transfer mediator. Molecular dynamics simulations, based on the crystal structure of S. cerevisiae, and on a complex generated by molecular docking in A. fulgidus, confirmed that this amino acid can shuttle protons, however, a water molecule in the active site may also serve this role. Quantum-chemical calculations based on density functional theory and the cluster approach showed that the water-mediated pathway is the most favorable for catalysis. Experimental kinetic and mutational studies reinforce the requirement for the aspartate D101, but not S233. These findings provide insight into the catalytic mechanisms underlying proton transfer over extended distances and comprehensively elucidate the mode of action of Nep1.

Applications and Potential of In Silico Approaches for Psychedelic Chemistry

Molecular-level investigations of the Central Nervous System have been revolutionized by the development of computational methods, computing power, and capacity advances. These techniques have enabled researchers to analyze large amounts of data from various sources, including genomics, in vivo, and in vitro drug tests. In this review, we explore how computational methods and informatics have contributed to our understanding of mental health disorders and the development of novel drugs for neurological diseases, with a special focus on the emerging field of psychedelics. In addition, the use of state-of-the-art computational methods to predict the potential of drug compounds and bioinformatic tools to integrate disparate data sources to create predictive models is also discussed. Furthermore, the challenges associated with these methods, such as the need for large datasets and the diversity of in vitro data, are explored. Overall, this review highlights the immense potential of computational methods and informatics in Central Nervous System research and underscores the need for continued development and refinement of these techniques and more inclusion of Quantitative Structure-Activity Relationships (QSARs).

How to Compute Atomistic Insight in DFT Clusters: The REG-IQA Approach

The relative energy gradient (REG) method is paired with the topological energy partitioning method interacting quantum atoms (IQA), as REG-IQA, to provide detailed and unbiased knowledge on the intra- and interatomic interactions. REG operates on a sequence of geometries representing a dynamical change of a system. Its recent application to peptide hydrolysis of the human immunodeficiency virus-1 (HIV-1) protease (PDB code: 4HVP) has demonstrated its full potential in recovering reaction mechanisms and through-space electrostatic and exchange-correlation effects, making it a compelling tool for analyzing enzymatic reactions. In this study, the computational efficiency of the REG-IQA method for the 133-atom HIV-1 protease quantum mechanical system is analyzed in every detail and substantially improved by means of three different approaches. The first approach of smaller integration grids for IQA integrations reduces the computational overhead by about a factor of 3. The second approach uses the line-simplification Ramer-Douglas-Peucker (RDP) algorithm, which outputs the minimal number of geometries necessary for the REG-IQA analysis for a predetermined root mean squared error (RMSE) tolerance. This cuts the computational time of the whole REG analysis by a factor of 2 if an RMSE of 0.5 kJ/mol is considered. The third approach consists of a "biased" or "unbiased" selection of a specific subset of atoms of the whole initial quantum mechanical model wave-function, which results in more than a 10-fold speed-up per geometry for the IQA calculation, without deterioration of the outcome of the REG-IQA analysis. Finally, to show the capability of these approaches, the findings gathered from the HIV-1 protease system are also applied to a different system named haloalcohol dehalogenase (HheC). In summary, this study takes the REG-IQA method to a computationally feasible and highly accurate level, making it viable for the analysis of a multitude of enzymatic systems.

Uncovering the Mechanism of Azepino-Indole Skeleton Formation via Pictet-Spengler Reaction by Strictosidine Synthase: A Quantum Chemical Investigation

Strictosidine synthase (STR) catalyzes the Pictet-Spengler (PS) reaction of tryptamine and secologanin to produce strictosidine. Recent studies demonstrated that the enzyme can also catalyze the reaction of non-natural substrates to form new alkaloid skeletons. For example, the PS condensation of 1H-indole-4-ethanamine with secologanin could be promoted by the STR from Rauvolfia serpentina (RsSTR) to generate a rare class of skeletons with a seven-membered ring, namely azepino-[3,4,5-cd]-indoles, which are precursors for the synthesis of new compounds displaying antimalarial activity. In the present study, the detailed reaction mechanism of RsSTR-catalyzed formation of the rare seven-membered azepino-indole skeleton through the PS reaction was revealed at the atomic level by quantum chemical calculations. The structures of the transition states and intermediates involved in the reaction pathway were optimized, and the energetics of the complete reaction were analyzed. Based on our calculation results, the most likely pathway of the enzyme-catalyzed reaction was determined, and the rate-determining step of the reaction was clarified. The mechanistic details obtained in the present study are important in understanding the promiscuous activity of RsSTR in the formation of the rare azepino-indole skeleton molecule and are also helpful in designing STR enzymes for the synthesis of other new alkaloid skeleton molecules.

Pre-oxidation coupled with charged covalent organic framework membranes for highly efficient removal of organic chloramines precursors in algae-containing water treatment

Organic chloramines in water would pose both chemical and microbiological risks. It is essential to remove the precursors of organic chloramine (amino acids and decomposed peptides/proteins) to limit its formation in disinfection. In our work, nanofiltration was chosen to remove organic chloramines precursors. To solve the "trade-off" effect and low rejection of small molecules in algae organic matter, we synthesized a thin film composite (TFC) nanofiltration (NF) membrane with a crumpled polyamide (PA) layer via interfacial polymerization on polyacrylonitrile (PAN) composite support loaded with covalent organic framework (COF) nanoparticles (TpPa-SO₃H). The obtained NF membrane (PA-TpPa-SO₃H/PAN) increased the permeance from 10.2 to 28.2 L m^-2 h^-1 bar^-1 and the amino acid rejection from 24% to 69% compared to the control NF membrane. The addition of TpPa-SO₃H nanoparticles decreased the thickness of PA layers, increased the hydrophilicity of the membrane, and increased the transition energy barrier for amino acids transferring through the membrane, which was identified by scanning electron microscope, contact angle test, and density functional theory computations, respectively. Finally, pre-oxidation coupled with PA-TpPa-SO₃H/PAN membrane nanofiltration on the limitation of organic chloramines formation was evaluated. We found that the combined application of KMnO₄ pre-oxidation and PA-TpPa-SO₃H/PAN membranes nanofiltration in algae-containing water treatment could minimize the formation of organic chloramines in subsequent chlorination and maintain a high flux during filtration. Our work provides an effective way for algae-containing water treatment and organic chloramines control.

Theoretical Study of the Catalytic Mechanism of the Cu-Only Superoxide Dismutase

The catalytic mechanisms for the wild-type and the mutated Cu-only superoxide dismutase were studied using the hybrid density functional B3LYP and a quantum chemical cluster approach. Optimal protonation states of the active site were examined for each stage of the catalytic cycle. For both the reductive and the oxidative half-reactions, the arrival of the substrate O2•- was found to be accompanied by a charge-compensating H⁺ with exergonicities of -15.4 kcal·mol and -4.7 kcal·mol, respectively. The second-sphere Glu-110 and first-sphere His-93 were suggested to be the transient protonation site for the reductive and the oxidative half-reactions, respectively, which collaborates with the hydrogen bonding water chain to position the substrate near the redox-active copper center. For the reductive half-reaction, the rate-limiting step was found to be the inner-sphere electron transfer from the partially coordinated O2•- to Cu^II with a barrier of 8.1 kcal·mol. The formed O₂ is released from the active site with an exergonicity of -14.9 kcal·mol. For the oxidative half-reaction, the inner-sphere electron transfer from Cu^I to the partially coordinated O2•- was found to be accompanied by the proton transfer from the protonated His-93 and barrierless. The rate-limiting step was found to be the second proton transfer from the protonated Glu-110 to HO2- with a barrier of 7.3 kcal·mol. The barriers are reasonably consistent with experimental activities, and a proton-transfer rate-limiting step in the oxidative half-reaction could explain the experimentally observed pH-dependence. For the E110Q CuSOD, Asp-113 was suggested to be likely to serve as the transient protonation site in the reductive half-reaction. The rate-limiting barriers were found to be 8.0 and 8.6 kcal·mol, respectively, which could explain the slightly lower performance of E110X mutants. The results were found to be stable, with respect to the percentage of exact exchange in B3LYP.

Tautomer-Specific Deacylation and Ω-Loop Flexibility Explain the Carbapenem-Hydrolyzing Broad-Spectrum Activity of the KPC-2 β-Lactamase

KPC-2 (Klebsiella pneumoniae carbapenemase-2) is a globally disseminated serine-β-lactamase (SBL) responsible for extensive β-lactam antibiotic resistance in Gram-negative pathogens. SBLs inactivate β-lactams via a mechanism involving a hydrolytically labile covalent acyl-enzyme intermediate. Carbapenems, the most potent β-lactams, evade the activity of many SBLs by forming long-lived inhibitory acyl-enzymes; however, carbapenemases such as KPC-2 efficiently deacylate carbapenem acyl-enzymes. We present high-resolution (1.25-1.4 Å) crystal structures of KPC-2 acyl-enzymes with representative penicillins (ampicillin), cephalosporins (cefalothin), and carbapenems (imipenem, meropenem, and ertapenem) obtained utilizing an isosteric deacylation-deficient mutant (E166Q). The mobility of the Ω-loop (residues 165-170) negatively correlates with antibiotic turnover rates (kcat), highlighting the role of this region in positioning catalytic residues for efficient hydrolysis of different β-lactams. Carbapenem-derived acyl-enzyme structures reveal the predominance of the Δ1-(2R) imine rather than the Δ2 enamine tautomer. Quantum mechanics/molecular mechanics molecular dynamics simulations of KPC-2:meropenem acyl-enzyme deacylation used an adaptive string method to differentiate the reactivity of the two isomers. These identify the Δ1-(2R) isomer as having a significantly (7 kcal/mol) higher barrier than the Δ2 tautomer for the (rate-determining) formation of the tetrahedral deacylation intermediate. Deacylation is therefore likely to proceed predominantly from the Δ2, rather than the Δ1-(2R) acyl-enzyme, facilitated by tautomer-specific differences in hydrogen-bonding networks involving the carbapenem C-3 carboxylate and the deacylating water and stabilization by protonated N-4, accumulating a negative charge on the Δ2 enamine-derived oxyanion. Taken together, our data show how the flexible Ω-loop helps confer broad-spectrum activity upon KPC-2, while carbapenemase activity stems from efficient deacylation of the Δ2-enamine acyl-enzyme tautomer.

Emerging Metabolic Profiles of Sulfonamide Antibiotics by Cytochromes P450: A Computational-Experimental Synergy Study on Emerging Pollutants

Metabolism, especially by CYP450 enzymes, is the main reason for mediating the toxification and detoxification of xenobiotics in humans, while some uncommon metabolic pathways, especially for emerging pollutants, probably causing idiosyncratic toxicity are easily overlooked. The pollution of sulfonamide antibiotics in aqueous system has attracted increasing public attention. Hydroxylation of the central amine group can trigger a series of metabolic processes of sulfonamide antibiotics in humans; however, this work parallelly reported the coupling and fragmenting initiated by amino H-abstraction of sulfamethoxazole (SMX) catalyzed by human CYP450 enzymes. Elucidation of the emerging metabolic profiles was mapped via a multistep synergy between computations and experiments, involving preliminary DFT computations and in vitro and in vivo assays, profiling adverse effects, and rationalizing the fundamental factors via targeted computations. Especially, the confirmed SMX dimer was shown to potentially act as a metabolism disruptor in humans, while spin aromatic delocalization resulting in the low electron donor ability of amino radicals was revealed as the fundamental factor to enable coupling of sulfonamide antibiotics by CYP450 through the nonconventional nonrebound pathway. This work may further strengthen the synergistic use of computations prior to experiments to avoid wasteful experimental screening efforts in environmental chemistry and toxicology.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (ε_bulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (ε_bulk^') from intrinsic KIEs on the rate-determining steps and the observed ε_bulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and S_N2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase

Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of β-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-β-amino acids. This enzyme has the potential to be used as a tool in the production of β-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol^-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.

QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase

The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b₅₉₅. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.

Alleviating the unwanted effects of oxidative stress on Aβ clearance: a review of related concepts and strategies for the development of computational modelling

Treatment for Alzheimer's disease (AD) can be more effective in the early stages. Although we do not completely understand the aetiology of the early stages of AD, potential pathological factors (amyloid beta [Aβ] and tau) and other co-factors have been identified as causes of AD, which may indicate some of the mechanism at work in the early stages of AD. Today, one of the primary techniques used to help delay or prevent AD in the early stages involves alleviating the unwanted effects of oxidative stress on Aβ clearance. 4-Hydroxynonenal (HNE), a product of lipid peroxidation caused by oxidative stress, plays a key role in the adduction of the degrading proteases. This HNE employs a mechanism which decreases catalytic activity. This process ultimately impairs Aβ clearance. The degradation of HNE-modified proteins helps to alleviate the unwanted effects of oxidative stress. Having a clear understanding of the mechanisms associated with the degradation of the HNE-modified proteins is essential for the development of strategies and for alleviating the unwanted effects of oxidative stress. The strategies which could be employed to decrease the effects of oxidative stress include enhancing antioxidant activity, as well as the use of nanozymes and/or specific inhibitors. One area which shows promise in reducing oxidative stress is protein design. However, more research is needed to improve the effectiveness and accuracy of this technique. This paper discusses the interplay of potential pathological factors and AD. In particular, it focuses on the effect of oxidative stress on the expression of the Aβ-degrading proteases through adduction of the degrading proteases caused by HNE. The paper also elucidates other strategies that can be used to alleviate the unwanted effects of oxidative stress on Aβ clearance. To improve the effectiveness and accuracy of protein design, we explain the application of quantum mechanical/molecular mechanical approach.

Sequential C-H Methylation Catalyzed by the B12 -Dependent SAM Enzyme TokK: Comprehensive Theoretical Study of Selectivities

TokK is a B₁₂ -dependent radical SAM enzyme involved in the biosynthesis of the β-lactam antibiotic asparenomycin A. It can catalyze three methylations on different sp³ -hybridized carbon positions to introduce an isopropyl side chain at the β-lactam ring of pantetheinylated carbapenem. Herein, we report a quantum chemical study of the reaction mechanism of TokK. A stepwise ''push-pull'' radical relay mechanism is proposed for each methylation: a 5'-deoxyadenosine radical first abstracts a hydrogen atom from the substrate in the active site, then methylcobalamin directionally donates a methyl group to the substrate. More importantly, calculations were able to uncover the origin of observed chemoselectivity and stereoselectivity for the first methylation and regioselectivity for the following two methylations. Further detailed distortion/interaction analysis can help to unravel the main factors controlling the selectivities. Our findings of sequential methylations by TokK could have profound implications for studying other B₁₂ -dependent radical SAM enzymes.

Revisiting the Burden Borne by Fumarase: Enzymatic Hydration of an Olefin

Fumarate hydratase (FH) is a remarkable catalyst that decreases the free energy of the catalyzed reaction by 30 kcal mol^-1, much larger than most exceptional enzymes with extraordinary catalytic rates. Two classes of FH are observed in nature: class-I and class-II, which have different folds, yet catalyze the same reversible hydration/dehydration reaction of the dicarboxylic acids fumarate/malate, with equal efficiencies. Using class-I FH from the hyperthermophilic archaeon Methanocaldococcus jannaschii (Mj) as a model along with comparative analysis with the only other available class-I FH structure from Leishmania major (Lm), we provide insights into the molecular mechanism of catalysis in this class of enzymes. The structure of MjFH apo-protein has been determined, revealing that large intersubunit rearrangements occur across apo- and holo-protein forms, with a largely preorganized active site for substrate binding. Site-directed mutagenesis of active site residues, kinetic analysis, and computational studies, including density functional theory (DFT) and natural population analysis, together show that residues interacting with the carboxylate group of the substrate play a pivotal role in catalysis. Our study establishes that an electrostatic network at the active site of class-I FH polarizes the substrate fumarate through interactions with its carboxylate groups, thereby permitting an easier addition of a water molecule across the olefinic bond. We propose a mechanism of catalysis in FH that occurs through transition-state stabilization involving the distortion of the electronic structure of the substrate olefinic bond mediated by the charge polarization of the bound substrate at the enzyme active site.

Mechanistic Investigation of Lysine-Targeted Covalent Inhibition of PI3Kδ via ONIOM QM:QM Computations

Phosphoinositide 3-kinase (PI3K) enzymes are important drug targets, especially in oncology, and several inhibitors are currently under investigation in clinical trials for the treatment of lymphocytic leukemia, follicular lymphoma, breast, thyroid, colorectal, and lung cancer. Targeted covalent inhibitors hold significant promise for drug discovery research especially for kinases. Targeting the lysine residues attracts attention as a new strategy in designing targeted covalent inhibitors, since the lysine residue provides several advantages over the traditional cysteine residue. Recently, new highly selective covalent inhibitors of PI3Kδ with activated ester warheads, targeting the conserved Lys779 residue, were reported. Based on the observed kinetics, a covalent inhibition mechanism was proposed, but the atomistic details of the reaction are still not understood. Therefore, in the present work, we have conducted quantum chemical ONIOM M06-2X/6-31+G(d,p):PM6 calculations on the active site cluster structure of PI3Kδ to elucidate the microscopic details of the mechanism of the aminolysis reaction between Lys779 and the ester inhibitors. Our calculations clearly discriminate the noncovalent methyl ester inhibitor and the covalent inhibitors with activated phenolic esters. For the representative p-NO₂, p-F, p-H, and p-OCH₃ phenolic esters, the Gibbs free energy profiles of alternative mechanistic paths through either Asp782 or Asp911 demonstrate the modulatory role of active site aspartate residues. The most plausible path alters depending on the electron-withdrawing/donating nature of the p-substituted phenolate leaving group. Inhibitors with sufficiently strong electron-withdrawing group prefer direct dissociation of the leaving group from the tetrahedral zwitterion intermediate, while the ones with electron-donating group favor the formation of a neutral tetrahedral intermediate prior to the dissociation. The relative Gibbs free energy barriers of p-NO₂ < p-F < p-H < p-OCH₃ substituted phenyl esters display the same qualitative trend as the experimentally measured k_inact/K_I values. Our results provide in depth insight into the mechanism, which can pave the way for optimizing the inhibitor efficiency.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Perspective on the Current State-of-the-Art of Quantum Computing for Drug Discovery Applications

Computational chemistry is an essential tool in the pharmaceutical industry. Quantum computing is a fast evolving technology that promises to completely shift the computational capabilities in many areas of chemical research by bringing into reach currently impossible calculations. This perspective illustrates the near-future applicability of quantum computation of molecules to pharmaceutical problems. We briefly summarize and compare the scaling properties of state-of-the-art quantum algorithms and provide novel estimates of the quantum computational cost of simulating progressively larger embedding regions of a pharmaceutically relevant covalent protein-drug complex involving the drug Ibrutinib. Carrying out these calculations requires an error-corrected quantum architecture that we describe. Our estimates showcase that recent developments on quantum phase estimation algorithms have dramatically reduced the quantum resources needed to run fully quantum calculations in active spaces of around 50 orbitals and electrons, from estimated over 1000 years using the Trotterization approach to just a few days with sparse qubitization, painting a picture of fast and exciting progress in this nascent field.

A Combined Computational–Experimental Study on the Substrate Binding and Reaction Mechanism of Salicylic Acid Decarboxylase

Salicylic acid decarboxylase (SDC) from the amidohydrolase superfamily (AHS) catalyzes the reversible decarboxylation of salicylic acid to form phenol. In this study, the substrate binding mode and reaction mechanism of SDC were investigated using computational and crystallographic methods. Quantum chemical calculations show that the enzyme follows the general mechanism of AHS decarboxylases. Namely, the reaction begins with proton transfer from a metal-coordinated aspartic acid residue (Asp298 in SDC) to the C1 of salicylic acid, which is followed by the C–C bond cleavage, to generate the phenol product and release CO2. Interestingly, the calculations show that SDC is a Mg-dependent enzyme rather than the previously proposed Zn-dependent, and the substrate is shown to be bidentately coordinated to the metal center in the catalysis, which is also different from the previous proposal. These predictions are corroborated by the crystal structure of SDC solved in complex with the substrate analogue 2-nitrophenol. The mechanistic insights into SDC in the present study provide important information for the rational design of the enzyme.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

The Mechanism of Biochemical NO‐Sensing: Insights from Computational Chemistry

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO‐receptor in humans is soluble guanylyl cyclase (sGC) – a histidine‐ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO‐signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high‐level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H‐NOX protein from Shewanella oneidensis (So H‐NOX) homologous to the NO‐binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six‐coordinate intermediate heme‐NO complex, followed by relatively facile His decoordination yielding a five‐coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.