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.

Molecular Mechanism of Substrate Oxidation in Lytic Polysaccharide Monooxygenases: Insight from Theoretical Investigations

Lytic polysaccharide monooxygenases (LPMOs) are copper enzymes that today comprise a large enzyme superfamily, grouped into the distinct members AA9-AA17 (with AA12 exempted). The LPMOs have the potential to facilitate the upcycling of biomass waste products by boosting the breakdown of cellulose and other recalcitrant polysaccharides. The cellulose biopolymer is the main component of biomass waste and thus comprises a large, unexploited resource. The LPMOs work through a catalytic, oxidative reaction whose mechanism is still controversial. For instance, the nature of the intermediate performing the oxidative reaction is an open question, and the same holds for the employed co-substrate. Here we review theoretical investigations addressing these questions. The applied theoretical methods are usually based on quantum mechanics (QM), often combined with molecular mechanics (QM/MM). We discuss advantages and disadvantages of the employed theoretical methods and comment on the interplay between theoretical and experimental results.

Pharmacophore Modeling of Targets Infested with Activity Cliffs via Molecular Dynamics Simulation Coupled with QSAR and Comparison with other Pharmacophore Generation Methods: KDR as Case Study

Activity cliffs (ACs) are defined as pairs of structurally similar compounds with large difference in their potencies against certain biotarget. We recently proposed that potent AC members induce significant entropically-driven conformational modifications of the target that unveil additional binding interactions, while their weakly-potent counterparts are enthalpically-driven binders with little influence on the protein target. We herein propose to extract pharmacophores for ACs-infested target(s) from molecular dynamics (MD) frames of purely "enthalpic" potent binder(s) complexed within the particular target. Genetic function algorithm/machine learning (GFA/ML) can then be employed to search for the best possible combination of MD pharmacophore(s) capable of explaining bioactivity variations within a list of inhibitors. We compared the performance of this approach with established ligand-based and structure-based methods. Kinase inserts domain receptor (KDR) was used as a case study. KDR plays a crucial role in angiogenic signaling and its inhibitors have been approved in cancer treatment. Interestingly, GFA/ML selected, MD-based, pharmacophores were of comparable performances to ligand-based and structure-based pharmacophores. The resulting pharmacophores and QSAR models were used to capture hits from the national cancer institute list of compounds. The most active hit showed anti-KDR IC50 of 2.76 µM.

Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force  field, acting on the CRO. We  find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO-environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism.  We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO-environment polarization in quantitative calculations of absorption spectra in FPs.

Functional Characterization and Protein Engineering of a Triterpene 3-/6-/2'-O-Glycosyltransferase Reveal a Conserved Residue Critical for the Regiospecificity

Engineering the function of triterpene glucosyltransferases (GTs) is challenging due to the large size of the sugar acceptors. In this work, we identified a multifunctional glycosyltransferase AmGT8 catalyzing triterpene 3-/6-/2'-O-glycosylation from the medicinal plant Astragalus membranaceus. To engineer its regiospecificity, a small mutant library was built based on semi-rational design. Variants A394F, A394D, and T131V were found to catalyze specific 6-O, 3-O, and 2'-O glycosylation, respectively. The origin of regioselectivity of AmGT8 and its A394F variant was studied by molecular dynamics and hydrogen deuterium exchange mass spectrometry. Residue 394 is highly conserved as A/G and is critical for the regiospecificity of the C- and O-GTs TcCGT1 and GuGT10/14. Finally, astragalosides III and IV were synthesized by mutants A394F, T131V and P192E. This work reports biocatalysts for saponin synthesis and gives new insights into protein engineering of regioselectivity in plant GTs.

Exploring Dihydroflavonol-4-Reductase Reactivity and Selectivity by QM/MM-MD Simulations

Flavonoids are secondary metabolites ubiquitously found in plants. Their antioxidant properties make them highly interesting natural compounds for use in pharmacology. Therefore, unravelling the mechanisms of flavonoid biosynthesis is an important challenge. Among all the enzymes involved in this biosynthetic pathway, dihydroflavonol-4-reductase (DFR) plays a key role in the production of anthocyanins and proanthocyanidins. Here, we provide new information on the mechanism of action of this enzyme by using QM/MM-MD simulations applied to both dihydroquercetin (DHQ) and dihydrokaempferol (DHK) substrates. The consideration of these very similar compounds shed light on the major role played by the enzyme on the stabilization of the transition state but also on the activation of the substrate before the reaction through near-attack conformer effects.

Spin-Regulated Electron Transfer and Exchange-Enhanced Reactivity in Fe4 S4 -Mediated Redox Reaction of the Dph2 Enzyme During the Biosynthesis of Diphthamide

The [4Fe-4S]-dependent radical S-adenosylmethionine (SAM) proteins is one of large families of redox enzymes that are able to carry a panoply of challenging transformations. Despite the extensive studies of structure-function relationships of radical SAM (RS) enzymes, the electronic state-dependent reactivity of the [4Fe-4S] cluster in these enzymes remains elusive. Using combined MD simulations and QM/MM calculations, we deciphered the electronic state-dependent reactivity of the [4Fe-4S] cluster in Dph2, a key enzyme involved in the biosynthesis of diphthamide. Our calculations show that the reductive cleavage of the S-C_(γ) bond is highly dependent on the electronic structure of [4Fe-4S]. Interestingly, the six electronic states can be classified into a low-energy and a high-energy groups, which are correlated with the net spin of Fe4 atom ligated to SAM. Due to the driving force of Fe4-C_(γ) bonding, the net spin on the Fe4 moiety dictate the shift of the opposite spin electron from the Fe1-Fe2-Fe3 block to SAM. Such spin-regulated electron transfer results in the exchange-enhanced reactivity in the lower-energy group compared with those in the higher-energy group. This reactivity principle provides fundamental mechanistic insights into reactivities of [4Fe-4S] cluster in RS enzymes.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

Lewis Acid Transition-Metal-Catalyzed Hydrogen Activation: Structures, Mechanisms, and Reactivities

As a new type of bifunctional catalyst, the Lewis acid transition-metal (LA-TM) catalysts have been widely applied for hydrogen activation. This study presents a mechanistic framework to understand the LA-TM-catalyzed H₂ activation through DFT studies. The mer(trans)-homolytic cleavage, the fac(cis)-homolytic cleavage, the synergetic heterolytic cleavage, and the dissociative heterolytic cleavage should be taken as general mechanisms for the field of LA-TM catalysis. Four typical LA-TM catalysts, the Z-type κ⁴ -L₃ B-Rh complex tri(azaindolyl)borane-Rh, the X-type κ³ -L₂ B-Co complex bis-phosphino-boryl (PBP)-Co, the η² -BC-type κ³ -L₂ B-Pd complex diphosphine-borane (DPB)-Pd, and the Z-type κ² -LB-Pt complex (boryl)iminomethane (BIM)-Pt are selected as representative models to systematically illustrate their mechanistic features and explore the influencing factors on mechanistic variations. Our results indicate that the tri(azaindolyl)borane-Rh catalyst favors the synergetic heterolytic mechanism; the PBP-Co catalyst prefers the mer(trans)-homolytic mechanism; the DPB-Pd catalyst operates through the fac(cis)-homolytic mechanism, whereas the BIM-Pt catalyst tends to undergo the dissociative heterolytic mechanism. The mechanistic variations are determined by the coordination geometry, the LA-TM bonding nature, the electronic structure of the TM center, and the flexibility or steric effect of the LA ligands. The presented mechanistic framework should provide helpful guidelines for LA-TM catalyst design and reaction developments.

Docking of Covalent Ligands: Challenges and Approaches

Covalent ligands have recently regained considerable attention in drug discovery. The rational design of such ligands, however, is still faced with particular challenges, mostly related to the fact that covalent bond formation is a quantum mechanical phenomenon which cannot adequately be handled by the force fields or empirical approaches typically used for noncovalent protein-ligand interactions. Although the necessity for quantum chemical approaches is clear, they cannot yet routinely be applied on large data sets of ligands or for a broader exploration of binding modes in docking calculations. On the other hand, technical solutions for performing docking calculations with covalent ligands are available, but their scope is normally quite limited. Scoring functions typically neglect the contribution from covalent bond formation completely. In this situation, the question arises how to approach covalent ligands and which methods to choose for their docking and design.

Mutagenesis of Sequence Determinants of Truncated Porcine ALOX15 Induces Changes in the Reaction Specificity by Altering the Catalytic Mechanism of Initial Hydrogen Abstraction

The reaction specificity of lipoxygenases is of physiological relevance since the various oxygenation products exhibit different biological activities. Among mammalian ALOX15 orthologs there are arachidonic acid 12- and 15-lipoxygenating enzymes and recent studies suggested an evolutionary switch in that reaction specificity during late primate development. Previous reports showed that 12-lipoxygenating ALOX15 orthologs can be converted to 15-lipoxygenating enzymes by site-directed mutagenesis of some sequence determinants. Unfortunately, the molecular basis for those alterations are not well understood. Here, the arachidonic acid 12-lipoxygenating N-terminal truncation variant of pig ALOX15, for which a crystal structure is available, was used to explore the catalytic mechanism of the specificity switch induced by mutagenesis of Val418 and Val419 sequence determinants. We found that Val418Ile+Val419Met double mutant is dominantly 15-lipoxygenating. Docking and MD simulations, and quantum mechanics/molecular mechanics calculations indicated that the wildtype energy barrier for arachidonic acid 15-lipoxygenation is 3.4 kcal mol^-1 higher than for 12-lipoxygenation. In contrast, for the Val418Ile+Val419Met double mutant the energy barrier for 12-lipoxygenation is 6.0 kcal mol^-1 higher than for 15-lipoxygenation. Our data suggest that enzyme-substrate complex geometries determine the value of these energy barriers and, as a consequence, the reaction specificity of ALOX15 orthologs.

Computational Chemistry: The Fate of Current Methods and Future Challenges

"Where do we go from here?" is the underlying question regarding the future (perhaps foreseeable) developments in computational chemistry. Although this young discipline has already permeated practically all of chemistry, it is likely to become even more powerful with the rapid development of computational hard- and software.

Metal Fluorides as Analogues for Studies on Phosphoryl Transfer Enzymes

The 1994 structure of a transition-state analogue with AlF₄^- and GDP complexed to G1α, a small G protein, heralded a new field of research into the structure and mechanism of enzymes that manipulate the transfer of phosphoryl (PO₃^- ) groups. The number of enzyme structures in the PDB containing metal fluorides (MF_x ) as ligands that imitate either a phosphoryl or a phosphate group was 357 at the end of 2016. They fall into three distinct geometrical classes: 1) Tetrahedral complexes based on BeF₃^- that mimic ground-state phosphates; 2) octahedral complexes, primarily based on AlF₄^- , which mimic "in-line" anionic transition states for phosphoryl transfer; and 3) trigonal bipyramidal complexes, represented by MgF₃^- and putative AlF₃⁰ moieties, which mimic the geometry of the transition state. The interpretation of these structures provides a deeper mechanistic understanding into the behavior and manipulation of phosphate monoesters in molecular biology. This Review provides a comprehensive overview of these structures, their uses, and their computational development.

On the Influence of the Protonation States of Active Site Residues on AChE Reactivation: A QM/MM Approach

Acetylcholinesterase (AChE), an enzyme of the serine hydrolase superfamily, is a mediator of signal transmission at cholinergic synapses by catalyzing acetylcholine cleavage into acetate and choline. This enzyme is vulnerable to covalent inhibition by organophosphate compounds (like VX). Covalent inhibition of AChE does not revert spontaneously. Known reactivator compounds have limited action in restoring catalytic activity. QM/MM simulations of VX-inhibited AChE reactivation by pralidoxime (2-PAM), a classical reactivator, were performed. These afforded a broad view of the effect of protonation states of active-site residues, and provide evidence for the role of Glu202, which needs to be protonated for reactivation to occur. In situ deprotonation of 2-PAM for both protonation states of Glu202 showed that His447 is able to deprotonate 2-PAM with the assistance of Glu202. Because the active site of serine hydrolases is highly conserved, this work provides new insights on the interplay between the catalytic triad residues and this glutamate, newly identified as protonatable.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

Fighting Cancer with Transition Metal Complexes: From Naked DNA to Protein and Chromatin Targeting Strategies

Many transition metal complexes have unique physicochemical properties that can be efficiently exploited in medicinal chemistry for cancer treatment. Traditionally, double-stranded DNA has been assumed to be the main binding target; however, recent studies have shown that nucleosomal DNA as well as proteins can act as dominant molecular binding partners. This has raised new questions about the molecular determinants that govern DNA versus protein binding selectivity, and has offered new ways to rationalize their biological activity and possible side effects. To address these questions, molecular simulations at an atomistic level of detail have been used to complement, support, and rationalize experimental data. Herein we review some relevant studies-focused on platinum and ruthenium compounds-to illustrate the power of state-of-the-art molecular simulation techniques and to demonstrate how the interplay between molecular simulations and experiments can make important contributions to elucidating the target preferences of some promising transition metal anticancer agents. This contribution aims at providing relevant information that may help in the rational design of novel drug-discovery strategies.

Kinetic and Theoretical Studies of Beta-Lactone Reactivity-A Quantitative Scale for Biological Application

Natural products comprise a rich source for bioactive molecules with medicinal relevance. Many of these contain electrophilic scaffolds that bind conserved enzyme active sites covalently. Prominent examples include beta-lactams and beta-lactones, which specifically acylate serine residues in diverse peptidases. Although these scaffolds appear similar, their bioactivities and corresponding protein targets vary. To quantify and dissect these differences in bioactivities, the kinetics of the reactions of beta-butyrolactone with a set of reference amines in buffered aqueous solution at 37 °C have been analyzed. Different product ratios of C1 versus C3 attack on the beta-butyrolactone have been observed, depending on the aliphatic or aromatic nature of the standard amine used. Quantum mechanics/molecular mechanics (QM/MM) calculations reveal that a H₃ O⁺ molecule has a crucial role in stabilizing C3 attack by aniline, through coordination of the lactone ring oxygen. In agreement with their weak proteome reactivity, monocyclic beta-lactams did not react with the set of standard nucleophiles studied herein. Bicyclic beta-lactams, however, exhibited a lower activation barrier, and thus, reacted with standard nucleophiles. This study represents a starting point for semiquantitative reactivity scales for natural products, which, in analogy to chemical reactivity scales, will provide predictions for electrophilic modifications in biological systems.

Isoprenoid Biosynthesis in Pathogenic Bacteria: Nuclear Resonance Vibrational Spectroscopy Provides Insight into the Unusual [4Fe-4S] Cluster of the E. coli LytB/IspH Protein

The LytB/IspH protein catalyzes the last step of the methylerythritol phosphate (MEP) pathway which is used for the biosynthesis of essential terpenoids in most pathogenic bacteria. Therefore, the MEP pathway is a target for the development of new antimicrobial agents as it is essential for microorganisms, yet absent in humans. Substrate-free LytB has a special [4Fe-4S](2+) cluster with a yet unsolved structure. This motivated us to use synchrotron-based nuclear resonance vibrational spectroscopy (NRVS) in combination with quantum chemical-molecular mechanical (QM/MM) calculations to gain more insight into the structure of substrate-free LytB. The apical iron atom of the [4Fe-4S](2+) is clearly linked to three water molecules. We additionally present NRVS data of LytB bound to its natural substrate, (E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate (HMBPP) and to the inhibitors (E)-4-amino-3-methylbut-2-en-1-yl diphosphate and (E)-4-mercapto-3-methylbut-2-en-1-yl diphosphate.

Valence bond and enzyme catalysis: a time to break down and a time to build up

Understanding enzyme catalysis and developing ability to control of it are two great challenges in biochemistry. A few successful examples of computational-based enzyme design have proved the fantastic potential of computational approaches in this field, however, relatively modest rate enhancements have been reported and the further development of complementary methods is still required. Herein we propose a conceptually simple scheme to identify the specific role that each residue plays in catalysis. The scheme is based on a breakdown of the total catalytic effect into contributions of individual protein residues, which are further decomposed into chemically interpretable components by using valence bond theory. The scheme is shown to shed light on the origin of catalysis in wild-type haloalkane dehalogenase (wt-DhlA) and its mutants. Furthermore, the understanding gained through our scheme is shown to have great potential in facilitating the selection of non-optimal sites for catalysis and suggesting effective mutations to enhance the enzymatic rate.

The extraordinary catalytic ability of peroxiredoxins: a combined experimental and QM/MM study on the fast thiol oxidation step

Peroxiredoxins (Prxs) catalyze the reduction of peroxides, a process of key relevance in a variety of cellular processes. The first step in the catalytic cycle of all Prxs is the oxidation of a cysteine residue to sulfenic acid, which occurs 10(3)-10(7) times faster than in free cysteine. We present an experimental kinetics and hybrid QM/MM investigation to explore the reaction of Prxs with H2O2 using alkyl hydroperoxide reductase E from Mycobacterium tuberculosis as a Prx model. We report for the first time the thermodynamic activation parameters of H2O2 reduction using Prx, which show that protein significantly lowers the activation enthalpy, with an unfavourable entropic effect, compared to the uncatalyzed reaction. The QM/MM simulations show that the remarkable catalytic effects responsible for the fast H2O2 reduction in Prxs are mainly due to an active-site arrangement, which establishes a complex hydrogen bond network activating both reactive species.

The Fenna-Matthews-Olson protein revisited: a fully polarizable (TD)DFT/MM description

We report a combined molecular dynamics and quantum mechanics (QM)/molecular mechanics (MM) analysis of the excitonic properties of the Fenna-Matthews-Olson (FMO) protein by using a polarizable MM model combined with a time-dependent density functional theory description. Overall, our results indicate that structural fluctuations, electrostatic interactions, and short-range quantum effects can significantly modulate the model Hamiltonian parameters (site energies and couplings). We find that the specific interactions with the axial ligand and the hydrogen-bonded residues are responsible for the energy ladder, with their effects being mainly due to electrostatic interactions, but with short-range quantum contributions that are not negligible. In addition, a striking modulation of the screening effects experienced by the BChl pairs, due to the heterogeneous polarizability of the FMO and solvent environment, is observed. Finally, we find that the exciton model gives a reliable description of the delocalized excited states in the complex.

Computational reference data for the photochemistry of cyclobutane pyrimidine dimers

The cis-syn cyclobutane pyrimidine dimer is one of the major classes of carcinogenic UV-induced DNA photoproducts. In this work, diverse high-level quantum-chemical methods were used to determine the spectroscopic properties of neutral (singlet and triplet) and charged (cation and anion) species of thymine dimers. Maps of potential energy, charge distribution, electron affinity, and ionization potential of the thymidine dimers were computed along the two dimerization coordinates for neutral and charged species, as well as for the singlet excited state. This set of data aims at providing consistent results computed with the same methods as for photodamage and repair. Based on these results, several different photo-, heat-, and charge-induced mechanisms of dimerization and repair are characterized and discussed. Additionally, a new stable dimer with methylmethylidene-hexahydropyrimidine structure was found in the S0 state.