Transition States, analogues, and drug development

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.

Review of phage display: A jack-of-all-trades and master of most biomolecule display

Phage display was first described by George P. Smith when it was shown that virus particles were capable of presenting foreign proteins on their surface. The technology has paved the way for the evolution of various biomolecules presentation and diverse selection strategies. This unique feature has been applied as a versatile platform for numerous applications in drug discovery, protein engineering, diagnostics, and vaccine development. Over the decades, the limits of biomolecules displayed on phage particles have expanded from peptides to proteomes and even alternative scaffolds. This has allowed phage display to be viewed as a versatile display platform to accommodate various biomolecules ranging from small peptides to larger proteomes which has significantly impacted advancements in the biomedical industry. This review will explore the vast array of biomolecules that have been successfully employed in phage display technology in biomedical research.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

EzMechanism: an automated tool to propose catalytic mechanisms of enzyme reactions

Over the years, hundreds of enzyme reaction mechanisms have been studied using experimental and simulation methods. This rich literature on biological catalysis is now ripe for use as the foundation of new knowledge-based approaches to investigate enzyme mechanisms. Here, we present a tool able to automatically infer mechanistic paths for a given three-dimensional active site and enzyme reaction, based on a set of catalytic rules compiled from the Mechanism and Catalytic Site Atlas, a database of enzyme mechanisms. EzMechanism (pronounced as 'Easy' Mechanism) is available to everyone through a web user interface. When studying a mechanism, EzMechanism facilitates and improves the generation of hypotheses, by making sure that relevant information is considered, as derived from the literature on both related and unrelated enzymes. We validated EzMechanism on a set of 62 enzymes and have identified paths for further improvement, including the need for additional and more generic catalytic rules.

Design and evaluation of substrate-product analog inhibitors for racemases and epimerases utilizing a 1,1-proton transfer mechanism

Racemases and epimerases catalyze the inversion of stereochemistry at asymmetric carbon atoms to generate stereoisomers that often play important roles in normal and pathological physiology. Consequently, there is interest in developing inhibitors of these enzymes for drug discovery. A strategy for the rational design of substrate-product analog (SPA) inhibitors of racemases and epimerases utilizing a direct 1,1-proton transfer mechanism is elaborated. This strategy assumes that two groups on the asymmetric carbon atom remain fixed at active-site binding determinants, while the hydrogen and third, motile group move during catalysis, with the latter potentially traveling between an R- and S-pocket at the active site. SPAs incorporate structural features of the substrate and product, often with geminal disubstitution on the asymmetric carbon atom to simultaneously present the motile group to both the R- and S-pockets. For racemases operating on substrates bearing three polar groups (glutamate, aspartate, and serine racemases) or with compact, hydrophobic binding pockets (proline racemase), substituent motion is limited and the design strategy furnishes inhibitors with poor or modest binding affinities. The approach is most successful when substrates have a large, motile hydrophobic group that binds at a plastic and/or capacious hydrophobic site. Potent inhibitors were developed for mandelate racemase, isoleucine epimerase, and α-methylacyl-CoA racemase using the SPA inhibitor design strategy, exhibiting binding affinities ranging from substrate-like to exceeding that of the substrate by 100-fold. This rational approach for designing inhibitors of racemases and epimerases having the appropriate active-site architectures is a useful strategy for furnishing compounds for drug development.

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

Phosphoinositide and redox dysregulation by the anticancer methylthioadenosine phosphorylase transition state inhibitor

Methylthio-DADMe-immucillin-A (MTDIA) is an 86 picomolar inhibitor of 5'-methylthioadenosine phosphorylase (MTAP) with potent and specific anti-cancer efficacy. MTAP salvages S-adenosylmethionine (SAM) from 5'-methylthioadenosine (MTA), a toxic metabolite produced during polyamine biosynthesis. Changes in MTAP expression are implicated in cancer growth and development, making MTAP an appealing target for anti-cancer therapeutics. Since SAM is involved in lipid metabolism, we hypothesised that MTDIA alters the lipidomes of MTDIA-treated cells. To identify these effects, we analysed the lipid profiles of MTDIA-treated Saccharomyces cerevisiae using ultra-high resolution accurate mass spectrometry (UHRAMS). MTAP inhibition by MTDIA, and knockout of the Meu1 gene that encodes for MTAP in yeast, caused global lipidomic changes and differential abundance of lipids involved in cell signaling. The phosphoinositide kinase/phosphatase signaling network was specifically impaired upon MTDIA treatment, and was independently validated and further characterised via altered localization of proteins integral to this network. Functional consequences of dysregulated lipid metabolism included a decrease in reactive oxygen species (ROS) levels induced by MTDIA that was contemporaneous with changes in immunological response factors (nitric oxide, tumour necrosis factor-alpha and interleukin-10) in mammalian cells. These results indicate that lipid homeostasis alterations and concomitant downstream effects may be associated with MTDIA mechanistic efficacy.

Distinctive Formation of a DNA-Protein Cross-Link during the Repair of DNA Oxidative Damage: Insights into Human Disease from MD Simulations and QM/MM Calculations

Reactive oxygen species damage DNA and result in health issues. The major damage product, 8-oxo-7,8-dihydroguanine (8oG), is repaired by human adenine DNA glycosylase homologue (MUTYH). Although MUTYH misfunction is associated with a genetic disorder called MUTYH-associated polyposis (MAP) and MUTYH is a potential target for cancer drugs, the catalytic mechanism required to develop disease treatments is debated in the literature. This study uses molecular dynamics simulations and quantum mechanics/molecular mechanics techniques initiated from DNA-protein complexes that represent different stages of the repair pathway to map the catalytic mechanism of the wild-type MUTYH bacterial homologue (MutY). This multipronged computational approach characterizes a DNA-protein cross-linking mechanism that is consistent with all previous experimental data and is a distinct pathway across the broad class of monofunctional glycosylase repair enzymes. In addition to clarifying how the cross-link is formed, accommodated by the enzyme, and hydrolyzed for product release, our calculations rationalize why cross-link formation is favored over immediate glycosidic bond hydrolysis, the accepted mechanism for all other monofunctional DNA glycosylases to date. Calculations on the Y126F mutant MutY highlight critical roles for active site residues throughout the reaction, while investigation of the N146S mutant rationalizes the connection between the analogous N224S MUTYH mutation and MAP. In addition to furthering our knowledge of the chemistry associated with a devastating disorder, the structural information gained about the distinctive MutY mechanism compared to other repair enzymes represents an important step for the development of specific and potent small-molecule inhibitors as cancer therapeutics.

The Origin of Anion-π Autocatalysis

The autocatalysis of epoxide-opening ether cyclizations on the aromatic surface of anion-π catalysts stands out as a leading example of emergent properties expected from the integration of unorthodox interactions into catalysis. A working hypothesis was proposed early on, but the mechanism of anion-π autocatalysis has never been elucidated. Here, we show that anion-π autocatalysis is almost independent of peripheral crowding in substrate and product. Inaccessible asymmetric anion-π autocatalysis and sometimes erratic reproducibility further support that the origin of anion-π autocatalysis is more complex than originally assumed. The apparent long-distance communication without physical contact calls for the inclusion of water between substrate and product on the catalytic aromatic surface. Efficient anion-π autocatalysis around equimolar amounts but poor activity in dry solvents and with excess water indicate that this inclusion of water requires high precision. Computational models suggest that two water molecules transmit dual substrate activation by the product and serve as proton shuttles along antiparallel but decoupled hydrogen-bonded chains to delocalize and stabilize evolving charge density in the transition state by "anion-π double bonds". This new transition-state model of anion-π autocatalysis provides a plausible mechanism that explains experimental results and brings anion-π catalysis to an unprecedented level of sophistication.

Capturing the free energy of transition state stabilization: insights from the inhibition of mandelate racemase

Mandelate racemase (MR) catalyses the Mg²⁺-dependent interconversion of (R)- and (S)-mandelate. To effect catalysis, MR stabilizes the altered substrate in the transition state (TS) by approximately 26 kcal mol^-1 (-ΔG_tx), such that the upper limit of the virtual dissociation constant of the enzyme-TS complex is 2 × 10^-19 M. Designing TS analogue inhibitors that capture a significant amount of ΔG_tx for binding presents a challenge since there are a limited number of protein binding determinants that interact with the substrate and the structural simplicity of mandelate constrains the number of possible isostructural variations. Indeed, current intermediate/TS analogue inhibitors of MR capture less than or equal to 30% of ΔG_tx because they fail to fully capitalize on electrostatic interactions with the metal ion, and the strength and number of all available electrostatic and H-bond interactions with binding determinants present at the TS. Surprisingly, phenylboronic acid (PBA), 2-formyl-PBA, and para-chloro-PBA capture 31-38% of ΔG_tx. The boronic acid group interacts with the Mg²⁺ ion and multiple binding determinants that effect TS stabilization. Inhibitors capable of forming multiple interactions can exploit the cooperative interactions that contribute to optimum binding of the TS. Hence, maximizing interactions with multiple binding determinants is integral to effective TS analogue inhibitor design. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Conformational change of organic cofactor PLP is essential for catalysis in PLP-dependent enzymes

Pyridoxal 5'-phosphate (PLP)-dependent enzymes are ubiquitous, catalyzing various biochemical reactions of approximately 4% of all classified enzymatic activities. They transform amines and amino acids into important metabolites or signaling molecules and are important drug targets in many diseases. In the crystal structures of PLP-dependent enzymes, organic cofactor PLP showed diverse conformations depending on the catalytic step. The conformational change of PLP is essential in the catalytic mechanism. In the study, we review the sophisticated catalytic mechanism of PLP, especially in transaldimination reactions. Most drugs targeting PLP-dependent enzymes make a covalent bond to PLP with the transaldimination reaction. A detailed understanding of organic cofactor PLP will help develop a new drug against PLP-dependent enzymes. [BMB Reports 2022; 55(9): 439-446].

Assessment of Reversibility for Covalent Cysteine Protease Inhibitors Using Quantum Mechanics/Molecular Mechanics Free Energy Surfaces

We have used molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials to investigate the reaction mechanism for covalent inhibition of cathepsin K and assess the reversibility of inhibition. The computed free energy profiles suggest that a nucleophilic attack by the catalytic cysteine on the inhibitor warhead and proton transfer from the catalytic histidine occur in a concerted manner. The results indicate that the reaction is more strongly exergonic for the alkyne-based inhibitors, which bind irreversibly to cathepsin K, than for the nitrile-based inhibitor odanacatib, which binds reversibly. Gas-phase energies were also calculated for the addition of methanethiol to structural prototypes for a number of warheads of interest in cysteine protease inhibitor design in order to assess electrophilicity. The approaches presented in this study are particularly applicable to assessment of novel warheads, and computed transition state geometries can be incorporated into molecular models for covalent docking.

The Dynamic Range of Acidity: Tracking Rules for the Unidirectional Penetration of Cellular Compartments

Labeled ammonium cations with pKa ∼7.4 accumulate in acidic organelles because they can be neutralized transiently to cross the membrane at cytosolic pH 7.2 but not at their internal pH<5.5. Retention in early endosomes with less acidic internal pH was achieved recently using weaker acids of up to pKa 9.8. We report here that primary ammonium cations with higher pKa 10.6, label early endosomes more efficiently. This maximized early endosome tracking coincides with increasing labeling of Golgi networks with similarly weak internal acidity. Guanidinium cations with pKa 13.5 cannot cross the plasma membrane in monomeric form and label the plasma membrane with selectivity for vesicles embarking into endocytosis. Self-assembled into micelles, guanidinium cations enter cells like arginine-rich cell-penetrating peptides and, driven by their membrane potential, penetrate mitochondria unidirectionally despite their high inner pH. The resulting tracking rules with an approximated dynamic range of pKa change ∼3.5 are expected to be generally valid, thus enabling the design of chemistry tools for biology research in the broadest sense. From a practical point of view, most relevant are two complementary fluorescent flipper probes that can be used to image the mechanics at the very beginning of endocytosis.

Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective

The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.

Synthesis and Characterization of Transition-State Analogue Inhibitors against Human DNA Methyltransferase 1

Hypermethylation of CpG regions by human DNA methyltransferase 1 (DNMT1) silences tumor-suppression genes, and inhibition of DNMT1 can reactivate silenced genes. The 5-azacytidines are approved inhibitors of DNMT1, but their mutagenic mechanism limits their utility. A synthon approach from the analogues of S-adenosylhomocysteine, methionine, and deoxycytidine recapitulated the chemical features of the DNMT1 transition state in the synthesis of 16 chemically stable transition-state mimics. Inhibitors causing both full and partial inhibition of purified DNMT1 were characterized. The inhibitors show modest selectivity for DNMT1 versus DNMT3b. Active-site docking predicts inhibitor interactions with S-adenosyl-l-methionine and deoxycytidine regions of the catalytic site, validated by direct binding analysis. Inhibitor action with purified DNMT1 is not reflected in cultured cells. A partial inhibitor activated cellular DNA methylation, and a full inhibitor had no effect on cellular DNA methylation. These compounds provide chemical access to a new family of noncovalent DNMT chemical scaffolds for use in DNA methyltransferases.

Current Status of the Use of Multifunctional Enzymes as Anti-Cancer Drug Targets

Fighting cancer is one of the major challenges of the 21st century. Among recently proposed treatments, molecular-targeted therapies are attracting particular attention. The potential targets of such therapies include a group of enzymes that possess the capability to catalyze at least two different reactions, so-called multifunctional enzymes. The features of such enzymes can be used to good advantage in the development of potent selective inhibitors. This review discusses the potential of multifunctional enzymes as anti-cancer drug targets along with the current status of research into four enzymes which by their inhibition have already demonstrated promising anti-cancer effects in vivo, in vitro, or both. These are PFK-2/FBPase-2 (involved in glucose homeostasis), ATIC (involved in purine biosynthesis), LTA₄H (involved in the inflammation process) and Jmjd6 (involved in histone and non-histone posttranslational modifications). Currently, only LTA₄H and PFK-2/FBPase-2 have inhibitors in active clinical development. However, there are several studies proposing potential inhibitors targeting these four enzymes that, when used alone or in association with other drugs, may provide new alternatives for preventing cancer cell growth and proliferation and increasing the life expectancy of patients.

Design, synthesis, and evaluation of transition-state analogs as inhibitors of the bacterial quorum sensing autoinducer synthase CepI

Quorum sensing is a bacterial signaling system that involves the synthesis, secretion and detection of signal molecules called autoinducers. The main autoinducer in Gram-negative bacteria are acylated homoserine lactones, produced by the LuxI family of autoinducer synthases and detected by the LuxR family of autoinducer receptors. Quorum sensing allows for changes in gene expression and bacterial behaviors in a coordinated, cell density dependent manner. Quorum sensing controls the expression of virulence factors in some human pathogens, making quorum sensing an antibacterial drug target. Here we describe the design and synthesis of transition-state analogs of the autoinducer synthase enzymatic reaction and the evaluation of these compounds as inhibitors of the synthase CepI. One such compound potently inhibits CepI and constitutes a new type of inhibitor against this underdeveloped antibacterial target.

Octahedral Trifluoromagnesate, an Anomalous Metal Fluoride Species, Stabilizes the Transition State in a Biological Motor

Isoelectronic metal fluoride transition state analogue (TSA) complexes, MgF₃^– and AlF₄^–, have proven to be immensely useful in understanding mechanisms of biological motors utilizing phosphoryl transfer. Here we report a previously unobserved octahedral TSA complex, MgF₃(H₂O)⁻, in a 1.5 Å resolution Zika virus NS3 helicase crystal structure. ¹⁹F NMR provided independent validation and also the direct observation of conformational tightening resulting from ssRNA binding in solution. The TSA stabilizes the two conformations of motif V of the helicase that link ATP hydrolysis with mechanical work. DFT analysis further validated the MgF₃(H₂O)⁻ species, indicating the significance of this TSA for studies of biological motors.

A resistant mutant of Plasmodium falciparum purine nucleoside phosphorylase uses wild-type neighbors to maintain parasite survival

Plasmodium falciparum purine nucleoside phosphorylase (PfPNP) catalyzes an essential step in purine salvage for parasite growth. 4′-Deaza-1′-Aza-2′-Deoxy-1′-(9-Methylene)-Immucillin-G (DADMe-ImmG) is a transition state analog inhibitor of this enzyme, and P. falciparum infections in an Aotus primate malaria model can be cleared by oral administration of DADMe-ImmG. P. falciparum cultured under increasing DADMe-ImmG drug pressure exhibited PfPNP gene amplification, increased protein expression, and point mutations involved in DADMe-ImmG binding. However, the weak catalytic properties of the M183L resistance mutation (∼17,000-fold decrease in catalytic efficiency) are inconsistent with the essential function of PfPNP. We hypothesized that M183L subunits may form mixed oligomers of native and mutant PfPNP monomers to give hybrid hexameric enzymes with properties conferring DADMe-ImmG resistance. To test this hypothesis, we designed PfPNP constructs that covalently linked native and the catalytically weak M183L mutant subunits. Engineered hybrid PfPNP yielded trimer-of-dimer hexameric protein with alternating native and catalytically weak M183L subunits. This hybrid PfPNP gave near-native K_m values for substrate, but the affinity for DADMe-ImmG and catalytic efficiency were both reduced approximately ninefold relative to a similar construct of native subunits. Contact between the relatively inactive M183L and native subunits is responsible for altered properties of the hybrid protein. Thus, gene amplification of PfPNP provides adequate catalytic activity while resistance to DADMe-ImmG occurs in the hybrid oligomer to promote parasite survival. Coupled with the slow development of drug resistance, this resistance mechanism highlights the potential for DADMe-ImmG use in antimalarial combination therapies.

Selective Golgi α-mannosidase II inhibitors: N-alkyl substituted pyrrolidines with a basic functional group

Enzymatic assays, molecular modeling and NMR studies of novel 1,4-dideoxy-1,4-imino-l-lyxitols provided new information on the GH38 family enzyme inhibitors and their selectivity.

Tricyclic Nucleobase Analogs and Their Ribosides as Substrates and Inhibitors of Purine-Nucleoside Phosphorylases III. Aminopurine Derivatives

Etheno-derivatives of 2-aminopurine, 2-aminopurine riboside, and 7-deazaadenosine (tubercidine) were prepared and purified using standard methods. 2-Aminopurine reacted with aqueous chloroacetaldehyde to give two products, both exhibiting substrate activity towards bacterial (E. coli) purine-nucleoside phosphorylase (PNP) in the reverse (synthetic) pathway. The major product of the chemical synthesis, identified as 1,N2-etheno-2-aminopurine, reacted slowly, while the second, minor, but highly fluorescent product, reacted rapidly. NMR analysis allowed identification of the minor product as N2,3-etheno-2-aminopurine, and its ribosylation product as N2,3-etheno-2-aminopurine-N2--D-riboside. Ribosylation of 1,N2-etheno-2-aminopurine led to analogous N2--d-riboside of this base. Both enzymatically produced ribosides were readily phosphorolysed by bacterial PNP to the respective bases. The reaction of 2-aminopurine-N9- -D-riboside with chloroacetaldehyde gave one major product, clearly distinct from that obtained from the enzymatic synthesis, which was not a substrate for PNP. A tri-cyclic 7-deazaadenosine (tubercidine) derivative was prepared in an analogous way and shown to be an effective inhibitor of the E. coli, but not of the mammalian enzyme. Fluorescent complexes of amino-purine analogs with E. coli PNP were observed.

Structural Basis for Finding OG Lesions and Avoiding Undamaged G by the DNA Glycosylase MutY

The adenine glycosylase MutY selectively initiates repair of OG:A lesions and, by comparison, avoids G:A mispairs. The ability to distinguish these closely related substrates relies on the C-terminal domain of MutY, which structurally resembles MutT. To understand the mechanism for substrate specificity, we crystallized MutY in complex with DNA containing G across from the high-affinity azaribose transition state analogue. Our structure shows that G is accommodated by the OG site and highlights the role of a serine residue in OG versus G discrimination. The functional significance of Ser308 and its neighboring residues was evaluated by mutational analysis, revealing the critical importance of a β loop in the C-terminal domain for mutation suppression in cells, and biochemical performance in vitro. This loop comprising residues Phe307, Ser308, and His309 (Geobacillus stearothermophilus sequence positions) is conserved in MutY but absent in MutT and other DNA repair enzymes and may therefore serve as a MutY-specific target exploitable by chemical biological probes.

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H₂F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H₂F in driving hydride transfer remains unclear. H₂F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM simulations reveal that hydride transfer is triggered by polarization at the C6 of H₂F. The σ C_β–H bonds contribute to the buildup of the cationic character during the chemical transformation, and hyperconjugation influences the formation of the transition state. Our findings provide key insights into the hydride transfer mechanism of the DHFR-catalyzed reaction, which is a target for antiproliferative drugs and a paradigmatic model in mechanistic enzymology.

Substrate-Differentiated Transition States of SET7/9-Catalyzed Lysine Methylation

Transition state stabilization is essential for rate acceleration of enzymatic reactions. Despite extensive studies on various transition state structures of enzymes, an intriguing puzzle is whether an enzyme can accommodate multiple transition states (TSs) to catalyze a chemical reaction. It is experimentally challenging to study this proposition in terms of the choices of suitable enzymes and the feasibility to distinguish multiple TSs. As a paradigm with the protein lysine methyltransferase (PKMT) SET7/9 paired with its physiological substrates H3 and p53, their TSs were solved with experimental kinetic isotope effects as computational constraints. Remarkably, SET7/9 adopts two structurally distinct TSs, a nearly symmetric S_N2 and an extremely early S_N2, for H3K4 and p53K372 methylation, respectively. The two TSs are also different from those previously revealed for other PKMTs. The setting of multiple TSs is expected to be essential for SET7/9 and likely other PKMTs to act on broad substrates with high efficiency.

Tri-Cyclic Nucleobase Analogs and their Ribosides as Substrates of Purine-Nucleoside Phosphorylases. II Guanine and Isoguanine Derivatives

Etheno-derivatives of guanine, O6-methylguanine, and isoguanine were prepared and purified using standard methods. The title compounds were examined as potential substrates of purine-nucleoside phosphorylases from various sources in the reverse (synthetic) pathway. It was found that 1,N2-etheno-guanine and 1,N6-etheno-isoguanine are excellent substrates for purine-nucleoside phosphorylase (PNP) from E. coli, while O6-methyl-N2,3-etheno-guanine exhibited moderate activity vs. this enzyme. The latter two compounds displayed intense fluorescence in neutral aqueous medium, and so did the corresponding ribosylation products. By contrast, PNP from calf spleens exhibited only modest activity towards 1,N6-etheno-isoguanine; the remaining compounds were not ribosylated by this enzyme. The enzymatic ribosylation of 1,N6-etheno-isoguanine using two forms of calf PNP (wild type and N243D) and E. coli PNP (wild type and D204N) gave three different products, which were identified on the basis of NMR analysis and comparison with the product of the isoguanosine reaction with chloroacetic aldehyde, which gave an essentially single compound, identified unequivocally as N9-riboside. With the wild-type E. coli enzyme as a catalyst, N9--d- and N7--d-ribosides are obtained in proportion ~1:3, while calf PNP produced another riboside, tentatively identified as N6--d-riboside. The potential application of various forms of PNP for synthesis of the tri-cyclic nucleoside analogs is discussed.

Insights into the ubiquinol/dioxygen binding and proton relay pathways of the alternative oxidase

The alternative oxidase (AOX) is a monotopic diiron carboxylate protein which catalyzes the four-electron reduction of dioxygen to water by ubiquinol. Although we have recently determined the crystal structure of Trypanosoma brucei AOX (TAO) in the presence and absence of ascofuranone (AF) derivatives (which are potent mixed type inhibitors) the mechanism by which ubiquinol and dioxygen binds to TAO remain inconclusive. In this article, ferulenol was identified as the first competitive inhibitor of AOX which has been used to probe the binding of ubiquinol. Surface plasmon resonance reveals that AF is a quasi-irreversible inhibitor of TAO whilst ferulenol binding is completely reversible. The structure of the TAO-ferulenol complex, determined at 2.7 Å, provided insights into ubiquinol binding and has also identified a potential dioxygen molecule bound in a side-on conformation to the diiron center for the first time.

Proteases: Pivot Points in Functional Proteomics

Proteases drive the life cycle of all proteins, ensuring the transportation and activation of newly minted, would-be proteins into their functional form while recycling spent or unneeded proteins. Far from their image as engines of protein digestion, proteases play fundamental roles in basic physiology and regulation at multiple levels of systems biology. Proteases are intimately associated with disease and modulation of proteolytic activity is the presumed target for successful therapeutics. "Proteases: Pivot Points in Functional Proteomics" examines the crucial roles of proteolysis across a wide range of physiological processes and diseases. The existing and potential impacts of proteolysis-related activity on drug and biomarker development are presented in detail. All told the decisive roles of proteases in four major categories comprising 23 separate subcategories are addressed. Within this construct, 15 sets of subject-specific, tabulated data are presented that include identification of proteases, protease inhibitors, substrates, and their actions. Said data are derived from and confirmed by over 300 references. Cross comparison of datasets indicates that proteases, their inhibitors/promoters and substrates intersect over a range of physiological processes and diseases, both chronic and pathogenic. Indeed, "Proteases: Pivot Points …" closes by dramatizing this very point through association of (pro)Thrombin and Fibrin(ogen) with: hemostasis, innate immunity, cardiovascular and metabolic disease, cancer, neurodegeneration, and bacterial self-defense.

Studies with Guanidinium- and Amidinium-Based Inhibitors Suggest Minimal Stabilization of Allylic Carbocation Intermediates by Dehydrosqualene and Squalene Synthases

Dehydrosqualene and squalene synthases catalyze the redox neutral and the reductive, head-to-head dimerization of farnesyl diphosphate, respectively. In each case, the reaction is thought to proceed via an initial dissociation of farnesyl diphosphate to form an allylic carbocation-pyrophosphate ion pair. This work describes the synthesis and testing of inhibitors in which a guanidinium or amidinium moiety is flanked by a phosphonylphosphinate group and a hydrocarbon tail. These functional groups bear a planar, delocalized, positive charge and therefore should act as excellent mimics of an allylic carbocation. An inhibitor bearing a neutral urea moiety was also prepared as a control. The positively charged inhibitors acted as competitive inhibitors against Staphylococcus aureus dehydrosqualene synthase with K_i values in the low micromolar range. Surprisingly, the neutral urea inhibitor was the most potent of the three. Similar trends were seen with the first half reaction of human squalene synthase. One interpretation of these results is that the active sites of these enzymes do not directly stabilize the allylic carbocation via electrostatic or π-cation interactions. Instead, it is likely that the enzymes use tight binding to the pyrophosphate and lipid moieties to promote catalysis and that electrostatic stabilization of the carbocation is provided by the bound pyrophosphate product. An alternate possibility is that these inhibitors cannot bind to the "ionization FPP-binding site" of the enzyme and only bind to the "nonionizing FPP-binding site". In either case, all reported attempts to generate potent inhibitors with cationic FPP analogues have been unsuccessful to date.

Protein structure and computational drug discovery

The first protein structures revealed a complex web of weak interactions stabilising the three-dimensional shape of the molecule. Small molecule ligands were then found to exploit these same weak binding events to modulate protein function or act as substrates in enzymatic reactions. As the understanding of ligand–protein binding grew, it became possible to firstly predict how and where a particular small molecule might interact with a protein, and then to identify putative ligands for a specific protein site. Computer-aided drug discovery, based on the structure of target proteins, is now a well-established technique that has produced several marketed drugs. We present here an overview of the various methodologies being used for structure-based computer-aided drug discovery and comment on possible future developments in the field.

Synthesis, Reactivity, Functionalization, and ADMET Properties of Silicon-Containing Nitrogen Heterocycles

Silicon-containing compounds have been largely ignored in drug design and development, despite their potential to improve not only the potency but also the physicochemical and ADMET ( absorption, distribution, metabolism, excretion, toxicity) properties of drug-like candidates because of the unique characteristics of silicon. This deficiency is in large part attributable to a lack of general methods for synthesizing diverse organosilicon structures. Accordingly, a new building block strategy has been developed that diverges from traditional approaches to incorporation of silicon into drug candidates. Flexible, multi-gram-scale syntheses of silicon-containing tetrahydroquinoline and tetrahydroisoquinoline building blocks are described, along with methods by which diversely functionalized silicon-containing nitrogen heterocycles can be rapidly built using common reactions optimized to accommodate the properties of silicon. Furthermore, to better clarify the liabilities and advantages of silicon incorporation, select compounds and their carbon analogues were challenged in ADMET-focused biological studies.

Immucillins in Infectious Diseases

The Immucillins are chemically stable analogues that mimic the ribocation and leaving-group features of N-ribosyltransferase transition states. Infectious disease agents often rely on ribosyltransferase chemistry in pathways involving precursor synthesis for nucleic acids, salvage of nucleic acid precursors, or synthetic pathways with nucleoside intermediates. Here, we review three infectious agents and the use of the Immucillins to taget enzymes essential to the parasites. First, DADMe-Immucillin-G is a purine nucleoside phosphorylase (PNP) inhibitor that blocks purine salvage and shows clinical potential for treatment for the malaria parasite Plasmodium falciparum, a purine auxotroph requiring hypoxanthine for purine nucleotide synthesis. Inhibition of the PNPs in the host and in parasite cells leads to apurinic starvation and death. Second, Helicobacter pylori, a causative agent of human ulcers, synthesizes menaquinone, an essential electron transfer agent, in a pathway requiring aminofutalosine nucleoside hydrolysis. Inhibitors of the H. pylori methylthioadenosine nucleosidase (MTAN) are powerful antibiotics for this organism. Synthesis of menaquinone by the aminofutalosine pathway does not occur in most bacteria populating the human gut microbiome. Thus, MTAN inhibitors provide high-specificity antibiotics for H. pylori and are not expected to disrupt the normal gut bacterial flora. Third, Immucillin-A was designed as a transition state analogue of the atypical PNP from Trichomonas vaginalis. In antiviral screens, Immucillin-A was shown to act as a prodrug. It is active against filoviruses and flaviviruses. In virus-infected cells, Immucillin-A is converted to the triphosphate, is incorporated into the viral transcript, and functions as an atypical chain-terminator for RNA-dependent RNA polymerases. Immucillin-A has entered clinical trials for use as an antiviral. We also summarize other Immucillins that have been characterized in successful clinical trials for T-cell lymphoma and gout. The human trials support the potential development of the Immucillins in infectious diseases.

Potent dialkyl substrate-product analogue inhibitors and inactivators of α-methylacyl-coenzyme A racemase from Mycobacterium tuberculosis by rational design

Rational approaches for the design of enzyme inhibitors furnish powerful strategies for developing pharmaceutical agents and tools for probing biological mechanisms. A new strategy for the development of gem-disubstituted substrate-product analogues as inhibitors of racemases and epimerases is elaborated using α-methylacyl-coenzyme A racemase from Mycobacterium tuberculosis (MtMCR) as a model enzyme. MtMCR catalyzes the epimerization at C2 of acyl-CoA substrates, a key step in the metabolism of branched-chain fatty acids. Moreover, the human enzyme is a potential target for the development of therapeutic agents directed against prostate cancer. We show that rationally designed, N,N-dialkylcarbamoyl-CoA substrate-product analogues inactivate MtMCR. Binding greatly exceeds that of the substrate, (S)-ibuprofenoyl-CoA, up to ∼250-fold and is proportional to the alkyl chain length (4-12 carbons) with the N,N-didecyl and N,N-didodecyl species having competitive inhibition constants with values of 1.9 ± 0.2 μM and 0.42 ± 0.04 μM, respectively. The presence of two decyl chains enhanced binding over a single decyl chain by ∼204-fold. Overall, the results reveal that gem-disubstituted substrate-product analogues can yield extremely potent inhibitors of an epimerase with a capacious active site.

Tricyclic nitrogen base 1,N6-ethenoadenine and its ribosides as substrates for purine-nucleoside phosphorylases: Spectroscopic and kinetic studies

The title compound is an excellent substrate for E. coli PNP, as well as for its D204N mutant. The main product of the synthetic reaction is N9-riboside, but some amount of N7-riboside is also present. Surprisingly, 1,N⁶-ethenoadenine is also ribosylated by both wild-type and mutated (N243D) forms of calf PNP, which catalyze the synthesis of a different riboside, tentatively identified as N6-β-D-ribosyl-1,N⁶-ethenoadenine. All ribosides are susceptible to phosphorolysis by the E. coli PNP (wild type). All the ribosides are fluorescent and can be utilized as analytical probes.

Synthesis of bis-Phosphate Iminoaltritol Enantiomers and Structural Characterization with Adenine Phosphoribosyltransferase

Phosphoribosyl transferases (PRTs) are essential in nucleotide synthesis and salvage, amino acid, and vitamin synthesis. Transition state analysis of several PRTs has demonstrated ribocation-like transition states with a partial positive charge residing on the pentose ring. Core chemistry for synthesis of transition state analogues related to the 5-phospho-α-d-ribosyl 1-pyrophosphate (PRPP) reactant of these enzymes could be developed by stereospecific placement of bis-phosphate groups on an iminoaltritol ring. Cationic character is provided by the imino group and the bis-phosphates anchor both the 1- and 5-phosphate binding sites. We provide a facile synthetic path to these molecules. Cyclic-nitrone redox methodology was applied to the stereocontrolled synthesis of three stereoisomers of a selectively monoprotected diol relevant to the synthesis of transition-state analogue inhibitors. These polyhydroxylated pyrrolidine natural product analogues were bis-phosphorylated to generate analogues of the ribocationic form of 5-phosphoribosyl 1-phosphate. A safe, high yielding synthesis of the key intermediate represents a new route to these transition state mimics. An enantiomeric pair of iminoaltritol bis-phosphates (L-DIAB and D-DIAB) was prepared and shown to display inhibition of Plasmodium falciparum orotate phosphoribosyltransferase and Saccharomyces cerevisiae adenine phosphoribosyltransferase (ScAPRT). Crystallographic inhibitor binding analysis of L- and D-DIAB bound to the catalytic sites of ScAPRT demonstrates accommodation of both enantiomers by altered ring geometry and bis-phosphate catalytic site contacts.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

The mechanism and high-free-energy transition state of lac repressor-lac operator interaction

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)–lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (k_a) and dissociation (k_d) rate constants. RO-destabilizing solutes (urea, KCl) reduce k_a comparably (urea) or more than (KCl) they increase k_d, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce k_d without affecting k_a significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Mechanistic enzymology in drug discovery: a fresh perspective

Given the therapeutic and commercial success of small-molecule enzyme inhibitors, as exemplified by kinase inhibitors in oncology, a major focus of current drug-discovery and development efforts is on enzyme targets. Understanding the course of an enzyme-catalysed reaction can help to conceptualize different types of inhibitor and to inform the design of screens to identify desired mechanisms. Exploiting this information allows the thorough evaluation of diverse compounds, providing the knowledge required to efficiently optimize leads towards differentiated candidate drugs. This review highlights the rationale for conducting high-quality mechanistic enzymology studies and considers the added value in combining such studies with orthogonal biophysical methods.

Optimization of Hydroxyethylamine Transition State Isosteres as Aspartic Protease Inhibitors by Exploiting Conformational Preferences

NMR conformational analysis of a hydroxyethylamine peptide isostere developed as an aspartic protease inhibitor shows that it is a flexible architecture. Cyclization to form pyrrolidines, piperidines, or morpholines results in a preorganization of the whole system in solution. The resulting conformation is similar to the conformation of the inhibitor in the active site of BACE-1. This entropic gain results in increased affinity for the enzyme when compared with the acyclic system. For morpholines 27 and 29, the combination of steric and electronic factors is exploited to orient substituents toward S1, S1', and S2' pockets both in the solution and in the bound states. These highly preorganized molecules proved to be the most potent compounds of the series. Additionally, the morpholines, unlike the pyrrolidine and piperidine analogues, have been found to be brain penetrant BACE-1 inhibitors.