Enzymatic Transition State Theory and Transition State Analogue Design

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Snapshots of the Reaction Coordinate of a Thermophilic 2'-Deoxyribonucleoside/ribonucleoside Transferase

Nucleosides are ubiquitous to life and are required for the synthesis of DNA, RNA, and other molecules crucial for cell survival. Despite the notoriously difficult organic synthesis of nucleosides, 2'-deoxynucleoside analogues can interfere with natural DNA replication and repair and are successfully employed as anticancer, antiviral, and antimicrobial compounds. Nucleoside 2'-deoxyribosyltransferase (dNDT) enzymes catalyze transglycosylation via a covalent 2'-deoxyribosylated enzyme intermediate with retention of configuration, having applications in the biocatalytic synthesis of 2'-deoxynucleoside analogues in a single step. Here, we characterize the structure and function of a thermophilic dNDT, the protein from Chroococcidiopsis thermalis (CtNDT). We combined enzyme kinetics with structural and biophysical studies to dissect mechanistic features in the reaction coordinate, leading to product formation. Bell-shaped pH-rate profiles demonstrate activity in a broad pH range of 5.5-9.5, with two very distinct pKa values. A pronounced viscosity effect on the turnover rate indicates a diffusional step, likely product (nucleobase1) release, to be rate-limiting. Temperature studies revealed an extremely curved profile, suggesting a large negative activation heat capacity. We trapped a 2'-fluoro-2'-deoxyarabinosyl-enzyme intermediate by mass spectrometry and determined high-resolution structures of the protein in its unliganded, substrate-bound, ribosylated, 2'-difluoro-2'-deoxyribosylated, and in complex with probable transition-state analogues. We reveal key features underlying (2'-deoxy)ribonucleoside selection, as CtNDT can also use ribonucleosides as substrates, albeit with a lower efficiency. Ribonucleosides are the building blocks of RNA and other key intracellular metabolites participating in energy and metabolism, expanding the scope of use of CtNDT in biocatalysis.

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

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.

Association of polymicrobial interactions with dental caries development and prevention

Dental caries is a common oral disease. In many cases, disruption of the ecological balance of the oral cavity can result in the occurrence of dental caries. There are many cariogenic microbiota and factors, and their identification allows us to take corresponding prevention and control measures. With the development of microbiology, the caries-causing bacteria have evolved from the traditional single Streptococcus mutans to the discovery of oral symbiotic bacteria. Thus it is necessary to systematically organized the association of polymicrobial interactions with dental caries development. In terms of ecology, caries occurs due to an ecological imbalance of the microbiota, caused by the growth and reproduction of cariogenic microbiota due to external factors or the disruption of homeostasis by one's own factors. To reduce the occurrence of dental caries effectively, and considering the latest scientific viewpoints, caries may be viewed from the perspective of ecology, and preventive measures can be taken; hence, this article systematically summarizes the prevention and treatment of dental caries from the aspects of ecological perspectives, in particular the ecological biofilm formation, bacterial quorum sensing, the main cariogenic microbiota, and preventive measures.

Enzymatic Synthesis of 2-Chloropurine Arabinonucleosides with Chiral Amino Acid Amides at the C6 Position and an Evaluation of Antiproliferative Activity In Vitro

A number of purine arabinosides containing chiral amino acid amides at the C6 position of the purine were synthesized using a transglycosylation reaction with recombinant E. coli nucleoside phosphorylases. Arsenolysis of 2-chloropurine ribosides with chiral amino acid amides at C6 was used for the enzymatic synthesis, and the reaction equilibrium shifted towards the synthesis of arabinonucleosides. The synthesized nucleosides were shown to be resistant to the action of E. coli adenosine deaminase. The antiproliferative activity of the synthesized nucleosides was studied on human acute myeloid leukemia cell line U937. Among all the compounds, the serine derivative exhibited an activity level (IC50 = 16 μM) close to that of Nelarabine (IC50 = 3 μM) and was evaluated as active.

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'.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation (TDHDX-MS) has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend TDHDX from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine or DAA) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different HDX-MS properties for the two ligands, as expected from their 10⁶-fold differences in strength of binding. By contrast, analogous patterns for TDHDX-MS emerge in the presence of both DAA and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended TDHDX to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility vs. site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

Comparative QM/MM Molecular Dynamics and Umbrella Sampling Simulations: Interaction of the Zinc-Bound Intermediate Gem-Diolate Trapoxin A Inhibitor and Acetyl-l-lysine Substrate with Histone Deacetylase 8

Targeting the genetic material without destruction is a priority to develop safe anticancer drugs. Histone deacetylase 8 (HDAC8), which is proved to be involved in carcinogenesis, is an enzyme associated with the chromatin for post-translational deacetylation of acetylated lysine. In this study, HDAC8 co-crystallized with the intermediate state tetrapeptide Trapoxin A (TA) inhibitor and the holoenzyme are utilized to find their conformational ensembles. Furthermore, the co-crystallized intermediate gem-diolate TA was used to find optimum interactions with the active site residues by conventional molecular dynamics (MD) simulation and QM/MM umbrella sampling. Finally, the intermediate state of the acetyl-l-lysine substrate was explored by QM/MM steered MD and compared to the binding of the intermediate state of the inhibitor. This research showed that HDAC8 is flexible and exists in conformational ensembles in its holoenzyme state. Binding of the intermediate state TA stabilizes its conformation. The optimum binding to the active site of HDAC8 for structures of gem-diolate TA (intermediate state) and acetyl-l-lysine (intermediate state) was determined according to the corresponding energy profiles. The use of these models will aid in the design of potentially reversible, potent, and selective inhibitors of HDAC8 for cancer treatment.

Inhibiting HTLV-1 Protease: A Viable Antiviral Target

Human T-cell lymphotropic virus type 1 (HTLV-1) is a retrovirus that can cause severe paralytic neurologic disease and immune disorders as well as cancer. An estimated 20 million people worldwide are infected with HTLV-1, with prevalence reaching 30% in some parts of the world. In stark contrast to HIV-1, no direct acting antivirals (DAAs) exist against HTLV-1. The aspartyl protease of HTLV-1 is a dimer similar to that of HIV-1 and processes the viral polyprotein to permit viral maturation. We report that the FDA-approved HIV-1 protease inhibitor darunavir (DRV) inhibits the enzyme with 0.8 μM potency and provides a scaffold for drug design against HTLV-1. Analogs of DRV that we designed and synthesized achieved submicromolar inhibition against HTLV-1 protease and inhibited Gag processing in viral maturation assays and in a chronically HTLV-1 infected cell line. Cocrystal structures of these inhibitors with HTLV-1 protease highlight opportunities for future inhibitor design. Our results show promise toward developing highly potent HTLV-1 protease inhibitors as therapeutic agents against HTLV-1 infections.

LuxS quorum sensing system and biofilm formation of oral microflora: A short review article

The LuxS quorum sensing system is considered as the main system that most of the oral bacteria use to communicate in order to create biofilms. Here we identified 11 of the most important biofilm formers that utilize the LuxS system and presented current and recent information regarding this system. Though different bacterial species are able to communicate thorough the LuxS system, it was also found that cross kingdom communication can occur between bacteria and fungi and bacteria and epithelial cells. Immune response also plays and important role in mitigating the effects of biofilms. Here we identified 6 of the most important molecules that are involved in the immune response to biofilms. These immune molecules maintain the stability in the oral cavity by preventing bacteria from overwhelming the space and simultaneously minimizing the immune response in order not to cause tissue damage. Here we also discuss current research being done in order to maintain the balance in the oral cavity via inhibiting biofilm formation without eradicating oral bacteria in order to prevent the overgrowth of other organisms such as Candida albicans. One approach being used is inhibiting AI-2 intermediates which leads to lack of quorum sensing communication between bacteria through the use of intermediate analogues. Another approach that found success is the utilization of D forms of sugars where D-ribose and D-galactose have been proven to inhibit the LuxS system and subsequently preventing the process of quorum sensing leading to the reduction in biofilm formation.

Hydrazides Are Potent Transition-State Analogues for Glutaminyl Cyclase Implicated in the Pathogenesis of Alzheimer's Disease

Amyloidogenic plaques are hallmarks of Alzheimer's disease (AD) and typically consist of high percentages of modified Aβ peptides bearing N-terminally cyclized glutamate residues. The human zinc(II) enzyme glutaminyl cyclase (QC) was shown in vivo to catalyze the cyclization of N-terminal glutamates of Aβ peptides in a pathophysiological side reaction establishing QC as a druggable target for therapeutic treatment of AD. Here, we report crystallographic snapshots of human QC catalysis acting on the neurohormone neurotensin that delineate the stereochemical course of catalysis and suggest that hydrazides could mimic the transition state of peptide cyclization and deamidation. This hypothesis is validated by a sparse-matrix inhibitor screening campaign that identifies hydrazides as the most potent metal-binding group compared to classic Zn binders. The structural basis of hydrazide inhibition is illuminated by X-ray structure analysis of human QC in complex with a hydrazide-bearing peptide inhibitor and reveals a pentacoordinated Zn complex. Our findings inform novel strategies in the design of potent and highly selective QC inhibitors by employing hydrazides as the metal-binding warhead.

Dissecting the Dynamics during Enzyme Catalysis: A Case Study of Pin1 Peptidyl-Prolyl Isomerase

Free energy surfaces have played a central role in studying protein conformational changes and enzymatic reactions over decades. Yet, free energy barriers and kinetics are highly dependent on the coordinates chosen to define the surface, and furthermore, the dynamics during the reactions are often overlooked. Our recent study on the Pin1-catalyzed isomerization reaction has indicated that the isomerization transition events remarkably deviate from the free energy path, highlighting the need to understand the reaction dynamics in more detail. To this end, here we investigate the reaction coordinates that describe the transition states of the free energy and transition pathways by minimizing the cross-entropy function. We show that the isomerization transition events can be expressed by the concerted changes in the improper torsion angle ζ and nearby backbone torsional angles of the ligand, whereas the transition state of the free energy surface involves changes in a broad range of coordinates including multiple protein-ligand interactions. The current result supports the previous finding that the isomerization transitions occur quickly from the conformational excited states, which is in sharp contrast to the slow and collective changes suggested from the free energy path. Our results further indicate that the coordinates derived from the transition trajectories are not sufficient for finding the transition states on the free energy surfaces due to the lack of information from conformational excited states.

Droplet-based optofluidic systems for measuring enzyme kinetics

The study of enzyme kinetics is of high significance in understanding metabolic networks in living cells and using enzymes in industrial applications. To gain insight into the catalytic mechanisms of enzymes, it is necessary to screen an enormous number of reaction conditions, a process that is typically laborious, time-consuming, and costly when using conventional measurement techniques. In recent times, droplet-based microfluidic systems have proved themselves to be of great utility in large-scale biological experimentation, since they consume a minimal sample, operate at high analytical throughput, are characterized by efficient mass and heat transfer, and offer high levels of integration and automation. The primary goal of this review is the introduction of novel microfluidic tools and detection methods for use in high-throughput and sensitive analysis of enzyme kinetics. The first part of this review focuses on introducing basic concepts of enzyme kinetics and describing most common microfluidic approaches, with a particular focus on segmented flow. Herein, the key advantages include accurate control over the flow behavior, efficient mass and heat transfer, multiplexing, and high-level integration with detection modalities. The second part describes the current state-of-the-art platforms for high-throughput and sensitive analysis of enzyme kinetics. In addition to our categorization of recent advances in measuring enzyme kinetics, we have endeavored to critically assess the limitations of each of these detection approaches and propose strategies to improve measurements in droplet-based microfluidics. Graphical abstract.

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.

Plasmodium Purine Metabolism and Its Inhibition by Nucleoside and Nucleotide Analogues

Malaria still affects around 200 million people and is responsible for more than 400,000 deaths per year, mostly children in subequatorial areas. This disease is caused by parasites of the Plasmodium genus. Only a few WHO-recommended treatments are available to prevent or cure plasmodial infections, but genetic mutations in the causal parasites have led to onset of resistance against all commercial antimalarial drugs. New drugs and targets are being investigated to cope with this emerging problem, including enzymes belonging to the main metabolic pathways, while nucleoside and nucleotide analogues are also a promising class of potential drugs. This review highlights the main metabolic pathways targeted for the development of potential antiplasmodial therapies based on nucleos(t)ide analogues, as well as the different series of purine-containing nucleoside and nucleotide derivatives designed to inhibit Plasmodium falciparum purine metabolism.

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.

Conformational Excitation and Nonequilibrium Transition Facilitate Enzymatic Reactions: Application to Pin1 Peptidyl-Prolyl Isomerase

Conformational flexibility of protein is essential for enzyme catalysis. Yet, how protein's conformational rearrangements and dynamics contribute to catalysis remains highly controversial. To unravel protein's role in catalysis, it is inevitable to understand the static and dynamic mechanisms simultaneously. To this end, here the Pin1-catalyzed isomerization reaction is studied from the two perspectives. The static view indicates that the hydrogen bonds involving Pin1 rearrange in a tightly coupled manner with isomerization. In sharp contrast, the isomerization dynamics are found to be very rapid; protein's slow conformational rearrangements thus cannot occur simultaneously with isomerization, and the reaction proceeds in a nonequilibrium manner. The distinctive protein conformations necessary to stabilize the transition state are prepared a priori, i.e., as conformational excited states. The present result suggests that enzymatic reaction is not a simple thermal activation from equilibrium directly to the transition state, thus adding a novel perspective to Pauling's view.

Chemical and Biochemical Perspectives of Protein Lysine Methylation

Protein lysine methylation is a distinct posttranslational modification that causes minimal changes in the size and electrostatic status of lysine residues. Lysine methylation plays essential roles in regulating fates and functions of target proteins in an epigenetic manner. As a result, substrates and degrees (free versus mono/di/tri) of protein lysine methylation are orchestrated within cells by balanced activities of protein lysine methyltransferases (PKMTs) and demethylases (KDMs). Their dysregulation is often associated with neurological disorders, developmental abnormalities, or cancer. Methyllysine-containing proteins can be recognized by downstream effector proteins, which contain methyllysine reader domains, to relay their biological functions. While numerous efforts have been made to annotate biological roles of protein lysine methylation, limited work has been done to uncover mechanisms associated with this modification at a molecular or atomic level. Given distinct biophysical and biochemical properties of methyllysine, this review will focus on chemical and biochemical aspects in addition, recognition, and removal of this posttranslational mark. Chemical and biophysical methods to profile PKMT substrates will be discussed along with classification of PKMT inhibitors for accurate perturbation of methyltransferase activities. Semisynthesis of methyllysine-containing proteins will also be covered given the critical need for these reagents to unambiguously define functional roles of protein lysine methylation.

Transition State Analysis of Adenosine Triphosphate Phosphoribosyltransferase

Adenosine triphosphate phosphoribosyltransferase (ATP-PRT) catalyzes the first step in histidine biosynthesis, a pathway essential to microorganisms and a validated target for antimicrobial drug design. The ATP-PRT enzyme catalyzes the reversible substitution reaction between phosphoribosyl pyrophosphate and ATP. The enzyme exists in two structurally distinct forms, a short- and a long-form enzyme. These forms share a catalytic core dimer but bear completely different allosteric domains and thus distinct quaternary assemblies. Understanding enzymatic transition states can provide essential information on the reaction mechanisms and insight into how differences in domain structure influence the reaction chemistry, as well as providing a template for inhibitor design. In this study, the transition state structures for ATP-PRT enzymes from Campylobacter jejuni and Mycobacterium tuberculosis (long-form enzymes) and from Lactococcus lactis (short-form) were determined and compared. Intrinsic kinetic isotope effects (KIEs) were obtained at reaction sensitive positions for the reverse reaction using phosphonoacetic acid, an alternative substrate to the natural substrate pyrophosphate. The experimental KIEs demonstrated mechanistic similarities between the three enzymes and provided experimental boundaries for quantum chemical calculations to characterize the transition states. Predicted transition state structures support a dissociative reaction mechanism with a D_N*A_N^‡ transition state. Weak interactions from the incoming nucleophile and a fully dissociated ATP adenine are predicted regardless of the difference in overall structure and quaternary assembly. These studies establish that despite significant differences in the quaternary assembly and regulatory machinery between ATP-PRT enzymes from different sources, the reaction chemistry and catalytic mechanism are conserved.

The Transition-State Structure for Human MAT2A from Isotope Effects

Human methionine S-adenosyltransferase (MAT2A) catalyzes the formation of S-adenosylmethionine (SAM) from ATP and methionine. Synthetic lethal genetic analysis has identified MAT2A as an anticancer target in tumor cells lacking expression of 5'-methylthioadenosine phosphorylase (MTAP). Approximately 15% of human cancers are MTAP^-/-. The remainder can be rendered MTAP^- through MTAP inhibitors. We used kinetic isotope effect (KIE), commitment factor (C_f), and binding isotope effect (BIE) measurements combined with quantum mechanical (QM) calculations to solve the transition state structure of human MAT2A. The reaction is characterized by an advanced S_N2 transition state. The bond forming from the nucleophilic methionine sulfur to the 5'-C of ATP is 2.03 Å at the transition state (bond order of 0.67). Departure of the leaving group triphosphate of ATP is well advanced and forms a 2.32 Å bond between the 5'-C of ATP and the oxygen of the triphosphate (bond order of 0.23). Interaction of MAT2A with its MAT2B regulatory subunit causes no change in the intrinsic KIEs, indicating the same transition state structure. The transition state for MAT2A is more advanced along the reaction coordinate (more product-like) than that from the near-symmetrical transition state of methionine adenosyltransferase from E. coli.

Transition State Analogue Inhibitors of 5'-Deoxyadenosine/5'-Methylthioadenosine Nucleosidase from Mycobacterium tuberculosis

Mycobacterium tuberculosis 5'-deoxyadenosine/5'-methylthioadenosine nucleosidase (Rv0091) catalyzes the N-riboside hydrolysis of its substrates 5'-methylthioadenosine (MTA) and 5'-deoxyadenosine (5'-dAdo). 5'-dAdo is the preferred substrate, a product of radical S-adenosylmethionine-dependent enzyme reactions. Rv0091 is characterized by a ribocation-like transition state, with low N-ribosidic bond order, an N7-protonated adenine leaving group, and an activated but weakly bonded water nucleophile. DADMe-Immucillins incorporating 5'-substituents of the substrates 5'-dAdo and MTA were synthesized and characterized as inhibitors of Rv0091. 5'-Deoxy-DADMe-Immucillin-A was the most potent among the 5'-dAdo transition state analogues with a dissociation constant of 640 pM. Among the 5'-thio substituents, hexylthio-DADMe-Immucillin-A was the best inhibitor at 87 pM. The specificity of Rv0091 for the Immucillin transition state analogues differs from those of other bacterial homologues because of an altered hydrophobic tunnel accepting the 5'-substituents. Inhibitors of Rv0091 had weak cell growth effects on M. tuberculosis or Mycobacterium smegmatis but were lethal toward Helicobacter pylori, where the 5'-methylthioadenosine nucleosidase is essential in menaquinone biosynthesis. We propose that Rv0091 plays a role in 5'-deoxyadenosine recycling but is not essential for growth in these Mycobacteria.

Heat Capacity Changes for Transition-State Analogue Binding and Catalysis with Human 5'-Methylthioadenosine Phosphorylase

Human 5'-methylthioadenosine phosphorylase (MTAP) catalyzes the phosphorolysis of 5'-methylthioadenosine (MTA). Its action regulates cellular MTA and links polyamine synthesis to S-adenosylmethionine (AdoMet) salvage. Transition state analogues with picomolar dissociation constants bind to MTAP in an entropically driven process at physiological temperatures, suggesting increased hydrophobic character or dynamic structure for the complexes. Inhibitor binding exhibits a negative heat capacity change (-ΔC_p), and thus the changes in enthalpy and entropy upon binding are strongly temperature-dependent. The ΔC_p of inhibitor binding by isothermal titration calorimetry does not follow conventional trends and is contrary to that expected from the hydrophobic effect. Thus, ligands of increasing hydrophobicity bind with increasing values of ΔC_p. Crystal structures of MTAP complexed to transition-state analogues MT-DADMe-ImmA, BT-DADMe-ImmA, PrT-ImmA, and a substrate analogue, MT-tubercidin, reveal similar active site contacts and overall protein structural parameters, despite large differences in ΔC_p for binding. In addition, ΔC_p values are not correlated with K_d values. Temperature dependence of presteady state kinetics revealed the chemical step for the MTAP reaction to have a negative heat capacity for transition state formation (-ΔC_p^‡). A comparison of the ΔC_p^‡ for MTAP presteady state chemistry and ΔC_p for inhibitor binding revealed those transition-state analogues most structurally and thermodynamically similar to the transition state. Molecular dynamics simulations of MTAP apoenzyme and complexes with MT-DADMe-ImmA and MT-tubercidin show small, but increased dynamic motion in the inhibited complexes. Variable temperature CD spectroscopy studies for MTAP-inhibitor complexes indicate remarkable protein thermal stability (to T_m = 99 °C) in complexes with transition-state analogues.

Structural Snapshots for Mechanism-Based Inactivation of a Glycoside Hydrolase by Cyclopropyl Carbasugars

Glycoside hydrolases (GHs) have attracted considerable attention as targets for therapeutic agents, and thus mechanism-based inhibitors are of great interest. We report the first structural analysis of a carbocyclic mechanism-based GH inactivator, the results of which show that the two Michaelis complexes are in 2 H3 conformations. We also report the synthesis and reactivity of a fluorinated analogue and the structure of its covalently linked intermediate (flattened 2 H3 half-chair). We conclude that these inactivator reactions mainly involve motion of the pseudo-anomeric carbon atom, knowledge that should stimulate the design of new transition-state analogues for use as chemical biology tools.

Conformational and electrostatic analysis of SN1 donor analogue glycomimetic inhibitors of ST3Gal-I mammalian sialyltransferase

Mammalian sialyltransferases play a role in the metastasis of various cancers in humans. Inhibitors of these enzymes will in principle be able to directly inhibit aberrant sialylation in cancer. Inhibitors of ST3Gal-I resembling the donor component of S_N1 Transition State structures were previously evaluated as part of a kinetics study. Here, using classical dynamics simulations and free energy perturbation calculations, we rationalize the performance of three of these donor analogue ST3Gal-I enzyme inhibitors. We find to inhibit the mammalian ST3Gal-I enzyme a donor analogue requires configurationally limited functionality. This is mediated by the binding of the inhibitor to the enzyme. The inhibitor's ability to interact with Y194 and T272 through a charged group such as a carboxylate is especially important. Furthermore, a conformational rigid form approximating the donor substrate is central. Here this is achieved by an intramolecular hydrogen bond formed between the carboxylate group and one of the ribose hydroxyl groups of the cytidine monophosphate (CMP) leaving group. This intramolecular interaction results in the donor substrate conformer complimenting the form of the catalytic binding site. Finally the carboxylate charge is essential for electrostatic pairing with the binding site. Substituting this group for an alcohol or amide results in severe weakening of the ligand binding. The carboxylate thus proves an to be an irreplaceable functional group and an essential pharmacophore.

Transition State Structure and Inhibition of Rv0091, a 5'-Deoxyadenosine/5'-methylthioadenosine Nucleosidase from Mycobacterium tuberculosis

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) is a bacterial enzyme that catalyzes the hydrolysis of the N-ribosidic bond in 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (SAH). MTAN activity has been linked to quorum sensing pathways, polyamine biosynthesis, and adenine salvage. Previously, the coding sequence of Rv0091 was annotated as a putative MTAN in Mycobacterium tuberculosis. Rv0091 was expressed in Escherichia coli, purified to homogeneity, and shown to be a homodimer, consistent with MTANs from other microorganisms. Substrate specificity for Rv0091 gave a preference for 5'-deoxyadenosine relative to MTA or SAH. Intrinsic kinetic isotope effects (KIEs) for the hydrolysis of [1'-(3)H], [1'-(14)C], [5'-(3)H2], [9-(15)N], and [7-(15)N]MTA were determined to be 1.207, 1.038, 0.998, 1.021, and 0.998, respectively. A model for the transition state structure of Rv0091 was determined by matching KIE values predicted via quantum chemical calculations to the intrinsic KIEs. The transition state shows a substantial loss of C1'-N9 bond order, well-developed oxocarbenium character of the ribosyl ring, and weak participation of the water nucleophile. Electrostatic potential surface maps for the Rv0091 transition state structure show similarity to DADMe-immucillin transition state analogues. DADMe-immucillin transition state analogues showed strong inhibition of Rv0091, with the most potent inhibitor (5'-hexylthio-DADMe-immucillinA) displaying a Ki value of 87 pM.

Transition state for the NSD2-catalyzed methylation of histone H3 lysine 36

Nuclear receptor SET domain containing protein 2 (NSD2) catalyzes the methylation of histone H3 lysine 36 (H3K36). It is a determinant in Wolf-Hirschhorn syndrome and is overexpressed in human multiple myeloma. Despite the relevance of NSD2 to cancer, there are no potent, selective inhibitors of this enzyme reported. Here, a combination of kinetic isotope effect measurements and quantum chemical modeling was used to provide subangstrom details of the transition state structure for NSD2 enzymatic activity. Kinetic isotope effects were measured for the methylation of isolated HeLa cell nucleosomes by NSD2. NSD2 preferentially catalyzes the dimethylation of H3K36 along with a reduced preference for H3K36 monomethylation. Primary Me-(14)C and (36)S and secondary Me-(3)H3, Me-(2)H3, 5'-(14)C, and 5'-(3)H2 kinetic isotope effects were measured for the methylation of H3K36 using specifically labeled S-adenosyl-l-methionine. The intrinsic kinetic isotope effects were used as boundary constraints for quantum mechanical calculations for the NSD2 transition state. The experimental and calculated kinetic isotope effects are consistent with an SN2 chemical mechanism with methyl transfer as the first irreversible chemical step in the reaction mechanism. The transition state is a late, asymmetric nucleophilic displacement with bond separation from the leaving group at (2.53 Å) and bond making to the attacking nucleophile (2.10 Å) advanced at the transition state. The transition state structure can be represented in a molecular electrostatic potential map to guide the design of inhibitors that mimic the transition state geometry and charge.

Tight binding enantiomers of pre-clinical drug candidates

MTDIA is a picomolar transition state analogue inhibitor of human methylthioadenosine phosphorylase and a femtomolar inhibitor of Escherichia coli methylthioadenosine nucleosidase. MTDIA has proven to be a non-toxic, orally available pre-clinical drug candidate with remarkable anti-tumour activity against a variety of human cancers in mouse xenografts. The structurally similar compound MTDIH is a potent inhibitor of human and malarial purine nucleoside phosphorylase (PNP) as well as the newly discovered enzyme, methylthioinosine phosphorylase, isolated from Pseudomonas aeruginosa. Since the enantiomers of some pharmaceuticals have revealed surprising biological activities, the enantiomers of MTDIH and MTDIA, compounds 1 and 2, respectively, were prepared and their enzyme binding properties studied. Despite binding less tightly to their target enzymes than their enantiomers compounds 1 and 2 are nanomolar inhibitors.

Chemoenzymatic Synthesis of (36)S Isotopologues of Methionine and S-Adenosyl-L-methionine

Substrates containing isotope labels at specific atoms are required for transition-state analysis based on the measurement of multiple kinetic isotope effects.(36)S-labeled l-methionine and S-adenosyl-l-methionine were synthesized from elemental sulfur using a chemoenzymatic approach with >98% (36)S enrichment. This method provides access to previously inaccessible sulfur isotope-labeled substrates for sulfur kinetic isotope effect studies.

Targeting metallo-carbapenemases via modulation of electronic properties of cephalosporins

The global proliferation of metallo-carbapenemase-producing Enterobacteriaceae has created an unmet need for inhibitors of these enzymes. The rational design of metallo-carbapenemase inhibitors requires detailed knowledge of their catalytic mechanisms. Nine cephalosporins, structurally identical except for the systematic substitution of electron-donating and withdrawing groups in the para position of the styrylbenzene ring, were synthesized and utilized to probe the catalytic mechanism of New Delhi metallo-β-lactamase (NDM-1). Under steady-state conditions, K(m) values were all in the micromolar range (1.5-8.1 μM), whereas k(cat) values varied widely (17-220 s(-1)). There were large solvent deuterium isotope effects for all substrates under saturating conditions, suggesting a proton transfer is involved in the rate-limiting step. Pre-steady-state UV-visible scans demonstrated the formation of short-lived intermediates for all compounds. Hammett plots yielded reaction constants (ρ) of -0.34 ± 0.02 and -1.15 ± 0.08 for intermediate formation and breakdown, respectively. Temperature-dependence experiments yielded ΔG(‡) values that were consistent with the Hammett results. These results establish the commonality of the formation of an azanide intermediate in the NDM-1-catalysed hydrolysis of a range cephalosporins with differing electronic properties. This intermediate is a promising target for judiciously designed β-lactam antibiotics that are poor NDM-1 substrates and inhibitors with enhanced active-site residence times.

Theoretical study on the dehydrogenation reaction of dihydrogen bonded phenol–borane-trimethylamine in the excited state

No dehydrogenation reaction occurs in the ground state of dihydrogen bonded phenol–BTMA. TS-S₁₀ points to the formation of a hydrogen molecule, while TS-S₁₁ points to the B atom. The dehydrogenation reaction along TS-S₁₀ is energy favorable, unlike that along TS-S₁₁.

How quickly can a β-hairpin fold from its transition state?

Understanding the structural nature of the free energy bottleneck(s) encountered in protein folding is essential to elucidating the underlying dynamics and mechanism. For this reason, several techniques, including Φ-value analysis, have previously been developed to infer the structural characteristics of such high free-energy or transition states. Herein we propose that one (or few) appropriately placed backbone and/or side chain cross-linkers, such as disulfides, could be used to populate a thermodynamically accessible conformational state that mimics the folding transition state. Specifically, we test this hypothesis on a model β-hairpin, Trpzip4, as its folding mechanism has been extensively studied and is well understood. Our results show that cross-linking the two β-strands near the turn region increases the folding rate by an order of magnitude, to about (500 ns)⁻¹, whereas cross-linking the termini results in a hyperstable β-hairpin that has essentially the same folding rate as the uncross-linked peptide. Taken together, these findings suggest that cross-linking is not only a useful strategy to manipulate folding free energy barriers, as shown in other studies, but also, in some cases, it can be used to stabilize a folding transition state analogue and allow for direct assessment of the folding process on the downhill side of the free energy barrier. The calculated free energy landscape of the cross-linked Trpzip4 also supports this picture. An empirical analysis further suggests, when folding of β-hairpins does not involve a significant free energy barrier, the folding time (τ) follows a power law dependence on the number of hydrogen bonds to be formed (n_H), namely, τ = τ₀n_H^α, with τ₀ = 20 ns and α = 2.3.

Transition-state inhibitors of purine salvage and other prospective enzyme targets in malaria

Malaria is a leading cause of human death within the tropics. The gradual generation of drug resistance imposes an urgent need for the development of new and selective antimalarial agents. Kinetic isotope effects coupled to computational chemistry have provided the relevant details on geometry and charge of enzymatic transition states to facilitate the design of transition-state analogs. These features have been reproduced into chemically stable mimics through synthetic chemistry, generating inhibitors with dissociation constants in the pico- to femto-molar range. Transition-state analogs are expected to contribute to the control of malaria.

Transition state analogues of Plasmodium falciparum and human orotate phosphoribosyltransferases

The survival and proliferation of Plasmodium falciparum parasites and human cancer cells require de novo pyrimidine synthesis to supply RNA and DNA precursors. Orotate phosphoribosyltransferase (OPRT) is an indispensible component in this metabolic pathway and is a target for antimalarials and antitumor drugs. P. falciparum (Pf) and Homo sapiens (Hs) OPRTs are characterized by highly dissociative transition states with ribocation character. On the basis of the geometrical and electrostatic features of the PfOPRT and HsOPRT transition states, analogues were designed, synthesized, and tested as inhibitors. Iminoribitol mimics of the ribocation transition state in linkage to pyrimidine mimics using methylene or ethylene linkers gave dissociation constants (Kd) as low as 80 nM. Inhibitors with pyrrolidine groups as ribocation mimics displayed slightly weaker binding affinities for OPRTs. Interestingly, p-nitrophenyl riboside 5'-phosphate bound to OPRTs with Kd values near 40 nM. Analogues designed with a C5-pyrimidine carbon-carbon bond to ribocation mimics gave Kd values in the range of 80-500 nM. Acyclic inhibitors with achiral serinol groups as the ribocation mimics also displayed nanomolar inhibition against OPRTs. In comparison with the nucleoside derivatives, inhibition constants of their corresponding 5'-phosphorylated transition state analogues are largely unchanged, an unusual property for a nucleotide-binding site. In silico docking of the best inhibitor into the HsOPRT active site supported an extensive hydrogen bond network associated with the tight binding affinity. These OPRT transition state analogues identify crucial components of potent inhibitors targeting OPRT enzymes. Despite their tight binding to the targets, the inhibitors did not kill cultured P. falciparum.