Entropy-driven binding of picomolar transition state analogue inhibitors to human 5'-methylthioadenosine phosphorylase

A Phenylboronic Acid-Based Transition State Analogue Yields Nanomolar Inhibition of Mandelate Racemase

Mandelate racemase (MR) catalyzes the Mg²⁺-dependent interconversion of (R)- and (S)-mandelate by stabilizing the altered substrate in the transition state (TS) by ∼26 kcal/mol. The enzyme has been employed as a model to explore the limits to which the free energy of TS stabilization may be captured by TS analogues to effect strong binding. Herein, we determined the thermodynamic parameters accompanying binding of a series of bromo-, chloro-, and fluoro-substituted phenylboronic acids (PBAs) by MR and found that binding was predominately driven by favorable entropy changes. 3,4-Dichloro-PBA was discovered to be the most potent inhibitor yet identified for MR, binding with a K_d^app value of 11 ± 2 nM and exceeding the binding of the substrate by ∼72,000-fold. The ΔC_p value accompanying binding (-488 ± 18 cal·mol^-1 K^-1) suggested that dispersion forces contribute significantly to the binding. The pH-dependence of the inhibition revealed that MR preferentially binds the anionic, tetrahedral form of 3,4-dichloro-PBA with a pH-independent K_i value of 5.7 ± 0.5 nM, which was consistent with the observed upfield shift of the ¹¹B NMR signal. The linear free energy relationship between log(k_cat/K_m) and log(1/K_i) for wild-type and 11 MR variants binding 3,4-dichloro-PBA had a slope of 0.8 ± 0.2, indicating that MR recognizes the inhibitor as an analogue of the TS. Hence, halogen substitution may be utilized to capture additional free energy of TS stabilization arising from dispersion forces to enhance the binding of boronic acid inhibitors by MR.

MYCN upregulates the transsulfuration pathway to suppress the ferroptotic vulnerability in MYCN-amplified neuroblastoma

Ferroptosis is an iron-dependent, oxidative form of cell death that is countered mainly by glutathione peroxidase 4 (GPX4) and the production of glutathione (GSH), which is formed from cysteine. The identification of the cancers that may benefit from pharmacological ferroptotic induction is just emerging. We recently demonstrated that inducing ferroptosis genetically or pharmacologically in MYCN-amplified neuroblastoma (NB) is a novel and effective way to kill these cells. MYCN increases iron metabolism and subsequent hydroxyl radicals through increased expression of the transferrin receptor 1 (TfR1) and low levels of the ferroportin receptor. To counter increased hydroxyl radicals, MYCN binds to the promoter of SLC3A2 (solute carrier family 3 member 2). SLC3A2 is a subunit of system Xc-, which is the cysteine-glutamate antiporter that exports glutamate and imports cystine. Cystine is converted to cysteine intracellularly. Here, we investigated other ways MYCN may increase cysteine levels. By performing metabolomics in a syngeneic NB cell line either expressing MYCN or GFP, we demonstrate that the transsulfuration pathway is activated by MYCN. Furthermore, we demonstrate that MYCN-amplified NB cell lines and tumors have higher levels of cystathionine beta-synthase (CBS), the rate-limiting enzyme in transsulfuration, which leads to higher levels of the thioether cystathionine (R-S-(2-amino-2-carboxyethyl)-l-homocysteine). In addition, MYCN-amplified NB tumors have high levels of methylthioadenosine phosphorylase (MTAP), an enzyme that helps salvage methionine following polyamine metabolism. MYCN directly binds to the promoter of MTAP. We propose that MYCN orchestrates both enhanced cystine uptake and enhanced activity of the transsulfuration pathway to counteract increased reactive oxygen species (ROS) from iron-induced Fenton reactions, ultimately contributing to a ferroptosis vulnerability in MYCN-amplified neuroblastoma.

Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations

Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target in the treatment of T-cell cancers and autoimmune disorders. Experimental studies on the inhibition of the hPNP observe that the first ligand bound to one of three subunits effectively inhibits the protein, while the binding of more ligands to the subsequent sites shows negative cooperativities. In this work, we performed extensive end-point and alchemical free energy calculations to determine the binding thermodynamics of the trimeric protein-ligand system. 13 Immucillin inhibitors with experimental results are under calculation. Two widely accepted charge schemes for small molecules including AM1-BCC and RESP are adopted for ligands. The results of RESP are in better agreement with the experimental reference. Further investigations of the interaction networks in the protein-ligand complexes reveal that several residues play significant roles in stabilizing the complex structure. The most commonly observed ones include PHE200, GLU201, MET219, and ASN243. The conformations of the protein in different protein-ligand complexes are observed to be similar. We expect these insights to aid the development of potent drugs targeting hPNP.

Entropy-driven binding of gut bacterial β-glucuronidase inhibitors ameliorates irinotecan-induced toxicity

Irinotecan inhibits cell proliferation and thus is used for the primary treatment of colorectal cancer. Metabolism of irinotecan involves incorporation of β-glucuronic acid to facilitate excretion. During transit of the glucuronidated product through the gastrointestinal tract, an induced upregulation of gut microbial β-glucuronidase (GUS) activity may cause severe diarrhea and thus force many patients to stop treatment. We herein report the development of uronic isofagomine (UIFG) derivatives that act as general, potent inhibitors of bacterial GUSs, especially those of Escherichia coli and Clostridium perfringens. The best inhibitor, C6-nonyl UIFG, is 23,300-fold more selective for E. coli GUS than for human GUS (Ki = 0.0045 and 105 μM, respectively). Structural evidence indicated that the loss of coordinated water molecules, with the consequent increase in entropy, contributes to the high affinity and selectivity for bacterial GUSs. The inhibitors also effectively reduced irinotecan-induced diarrhea in mice without damaging intestinal epithelial cells.

Selective Inhibitors of Helicobacter pylori Methylthioadenosine Nucleosidase and Human Methylthioadenosine Phosphorylase

Bacterial 5'-methylthioadenosine/ S-adenosylhomocysteine nucleosidase (MTAN) hydrolyzes adenine from its substrates to form S-methyl-5-thioribose and S-ribosyl-l-homocysteine. MTANs are involved in quorum sensing, menaquinone synthesis, and 5'-methylthioadenosine recycling to S-adenosylmethionine. Helicobacter pylori uses MTAN in its unusual menaquinone pathway, making H. pylori MTAN a target for antibiotic development. Human 5'-methylthioadenosine phosphorylase (MTAP), a reported anticancer target, catalyzes phosphorolysis of 5'-methylthioadenosine to salvage S-adenosylmethionine. Transition-state analogues designed for HpMTAN and MTAP show significant overlap in specificity. Fifteen unique transition-state analogues are described here and are used to explore inhibitor specificity. Several analogues of HpMTAN bind in the picomolar range while inhibiting human MTAP with orders of magnitude weaker affinity. Structural analysis of HpMTAN shows inhibitors extending through a hydrophobic channel to the protein surface. The more enclosed catalytic sites of human MTAP require the inhibitors to adopt a folded structure, displacing the phosphate nucleophile from the catalytic site.

Enzymatic Transition States and Drug Design

Transition state theory teaches that chemically stable mimics of enzymatic transition states will bind tightly to their cognate enzymes. Kinetic isotope effects combined with computational quantum chemistry provides enzymatic transition state information with sufficient fidelity to design transition state analogues. Examples are selected from various stages of drug development to demonstrate the application of transition state theory, inhibitor design, physicochemical characterization of transition state analogues, and their progress in drug development.

Glutamine Hydrolysis by Imidazole Glycerol Phosphate Synthase Displays Temperature Dependent Allosteric Activation

The enzyme imidazole glycerol phosphate synthase (IGPS) is a model for studies of long-range allosteric regulation in enzymes. Binding of the allosteric effector ligand N'-[5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) stimulates millisecond (ms) timescale motions in IGPS that enhance its catalytic function. We studied the effect of temperature on these critical conformational motions and the catalytic mechanism of IGPS from the hyperthermophile Thermatoga maritima in an effort to understand temperature-dependent allostery. Enzyme kinetic and NMR dynamics measurements show that apo and PRFAR-activated IGPS respond differently to changes in temperature. Multiple-quantum Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments performed at 303, 323, and 343 K (30, 50, and 70°C) reveal that millisecond flexibility is enhanced to a higher degree in apo IGPS than in the PRFAR-bound enzyme as the sample temperature is raised. We find that the flexibility of the apo enzyme is nearly identical to that of its PRFAR activated state at 343 K, whereas conformational motions are considerably different between these two forms of the enzyme at room temperature. Arrhenius analyses of these flexible sites show a varied range of activation energies that loosely correlate to allosteric communities identified by computational methods and reflect local changes in dynamics that may facilitate conformational sampling of the active conformation. In addition, kinetic assays indicate that allosteric activation by PRFAR decreases to 65-fold at 343 K, compared to 4,200-fold at 303 K, which mirrors the decreased effect of PRFAR on ms motions relative to the unactivated enzyme. These studies indicate that at the growth temperature of T. maritima, PFRAR is a weaker allosteric activator than it is at room temperature and illustrate that the allosteric mechanism of IGPS is temperature dependent.

A Paradigm for CH Bond Cleavage: Structural and Functional Aspects of Transition State Stabilization by Mandelate Racemase

Mandelate racemase (MR) from Pseudomonas putida catalyzes the Mg²⁺-dependent, 1,1-proton transfer reaction that racemizes (R)- and (S)-mandelate. MR shares a partial reaction (i.e., the metal ion-assisted, Brønsted base-catalyzed proton abstraction of the α-proton of carboxylic acid substrates) and structural features ((β/α)₇β-barrel and N-terminal α + β capping domains) with a vast group of homologous, yet functionally diverse, enzymes in the enolase superfamily. Mechanistic and structural studies have developed this enzyme into a paradigm for understanding how enzymes such as those of the enolase superfamily overcome kinetic and thermodynamic barriers to catalyze the abstraction of an α-proton from a carbon acid substrate with a relatively high pK_a value. Structural studies on MR bound to intermediate/transition state analogues have delineated those structural features that MR uses to stabilize transition states and enhance reaction rates of proton abstraction. Kinetic, site-directed mutagenesis, and structural studies have also revealed that the phenyl ring of the substrate migrates through the hydrophobic cavity within the active site during catalysis and that the Brønsted acid-base catalysts (Lys 166 and His 297) may be utilized as binding determinants for inhibitor recognition. In addition, structural studies on the adduct formed from the irreversible inhibition of MR by 3-hydroxypyruvate revealed that MR can form and deprotonate a Schiff-base with 3-hydroxypyruvate to yield an enol(ate)-aldehyde adduct, suggesting a possible evolutionary link between MR and the Schiff-base forming aldolases. As the archetype of the enolase superfamily, mechanistic and structural studies on MR will continue to enhance our understanding of enzyme catalysis and furnish insights into the evolution of enzyme function.

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.

Continuous Fluorescence Assays for Reactions Involving Adenine

5'-Methylthioadenosine phosphorylase (MTAP) and 5'-methylthioadenosine nucleosidase (MTAN) catalyze the phosphorolysis and hydrolysis of 5'-methylthioadenosine (MTA), respectively. Both enzymes have low KM values for their substrates. Kinetic assays for these enzymes are challenging, as the ultraviolet absorbance spectra for reactant MTA and product adenine are similar. We report a new assay using 2-amino-5'-methylthioadenosine (2AMTA) as an alternative substrate for MTAP and MTAN enzymes. Hydrolysis or phosphorolysis of 2AMTA forms 2,6-diaminopurine, a fluorescent and easily quantitated product. We kinetically characterize 2AMTA with human MTAP, bacterial MTANs and use 2,6-diaminopurine as a fluorescent substrate for yeast adenine phosphoribosyltransferase. 2AMTA was used as the substrate to kinetically characterize the dissociation constants for three-transition-state analogue inhibitors of MTAP and MTAN. Kinetic values obtained from continuous fluorescent assays with MTA were in good agreement with previously measured literature values, but gave smaller experimental errors. Chemical synthesis from ribose and 2,6-dichloropurine provided crystalline 2AMTA as the oxalate salt. Chemo-enzymatic synthesis from ribose and 2,6-diaminopurine produced 2-amino-S-adenosylmethionine for hydrolytic conversion to 2AMTA. Interaction of 2AMTA with human MTAP was also characterized by pre-steady-state kinetics and by analysis of the crystal structure in a complex with sulfate as a catalytically inert analogue of phosphate. This assay is suitable for inhibitor screening by detection of fluorescent product, for quantitative analysis of hits by rapid and accurate measurement of inhibition constants in continuous assays, and pre-steady-state kinetic analysis of the target enzymes.

An additional role for the Brønsted acid-base catalysts of mandelate racemase in transition state stabilization

Mandelate racemase (MR) catalyzes the interconversion of the enantiomers of mandelate and serves as a paradigm for understanding the enzyme-catalyzed abstraction of an α-proton from a carbon acid substrate with a high pKa. The enzyme utilizes a two-base mechanism with Lys 166 and His 297 acting as Brønsted acid and base catalysts, respectively, in the R → S reaction direction. In the S → R reaction direction, their roles are reversed. Using isothermal titration calorimetry (ITC), MR is shown to bind the intermediate/transition state (TS) analogue inhibitor benzohydroxamate (BzH) in an entropy-driven process with a value of ΔCp equal to -358 ± 3 cal mol(-1) K(-1), consistent with an increased number of hydrophobic interactions. However, MR binds BzH with an affinity that is ∼2 orders of magnitude greater than that predicted solely on the basis of hydrophobic interactions [St. Maurice, M., and Bearne, S. L. (2004) Biochemistry 43, 2524], suggesting that additional specific interactions contribute to binding. To test the hypothesis that cation-π/NH-π interactions between the side chains of Lys 166 and His 297 and the aromatic ring and/or the hydroxamate/hydroximate moiety of BzH contribute to the binding of BzH, site-directed mutagenesis was used to generate the MR variants K166M, K166C, H297N, and K166M/H297N and their binding affinity for various ligands determined using ITC. Comparison of the binding affinities of these MR variants with the intermediate/TS analogues BzH and cyclohexanecarbohydroxamate revealed that cation-π/NH-π interactions between His 297 and the hydroxamate/hydroximate moiety and the phenyl ring of BzH contribute approximately 0.26 and 0.91 kcal/mol to binding, respectively, while interactions with Lys 166 contribute approximately 1.74 and 1.74 kcal/mol, respectively. Similarly, comparison of the binding affinities of these mutants with substrate analogues revealed that Lys 166 contributes >2.93 kcal/mol to the binding of (R)-atrolactate, and His 297 contributes 2.46 kcal/mol to the binding of (S)-atrolactate. These results are consistent with Lys 166 and His 297 playing dual roles in catalysis: they act as Brønsted acid-base catalysts, and they stabilize both the enolate moiety and phenyl ring of the altered substrate in the TS.

A structural and energetic model for the slow-onset inhibition of the Mycobacterium tuberculosis enoyl-ACP reductase InhA

Slow-onset enzyme inhibitors are of great interest for drug discovery programs since the slow dissociation of the inhibitor from the drug-target complex results in sustained target occupancy leading to improved pharmacodynamics. However, the structural basis for slow-onset inhibition is often not fully understood, hindering the development of structure-kinetic relationships and the rational optimization of drug-target residence time. Previously we demonstrated that slow-onset inhibition of the Mycobacterium tuberculosis enoyl-ACP reductase InhA correlated with motions of a substrate-binding loop (SBL) near the active site. In the present work, X-ray crystallography and molecular dynamics simulations have been used to map the structural and energetic changes of the SBL that occur upon enzyme inhibition. Helix-6 within the SBL adopts an open conformation when the inhibitor structure or binding kinetics is substrate-like. In contrast, slow-onset inhibition results in large-scale local refolding in which helix-6 adopts a closed conformation not normally populated during substrate turnover. The open and closed conformations of helix-6 are hypothesized to represent the EI and EI* states on the two-step induced-fit reaction coordinate for enzyme inhibition. These two states were used as the end points for nudged elastic band molecular dynamics simulations resulting in two-dimensional potential energy profiles that reveal the barrier between EI and EI*, thus rationalizing the binding kinetics observed with different inhibitors. Our findings indicate that the structural basis for slow-onset kinetics can be understood once the structures of both EI and EI* have been identified, thus providing a starting point for the rational control of enzyme-inhibitor binding kinetics.

Catalytic site cooperativity in dimeric methylthioadenosine nucleosidase

5′-Methylthioadenosine/S-adenosylhomocysteine nucleosidases (MTANs) are bacterial enzymes that catalyze hydrolysis of the N-ribosidic bonds of 5′-methylthioadenosine (MTA) and S-adenosylhomocysteine (SAH) to form adenine and 5-thioribosyl groups. MTANs are involved in AI-1 and AI-2 bacterial quorum sensing and the unusual futalosine-based menaquinone synthetic pathway in Streptomyces,Helicobacter, and Campylobacter species. Crystal structures show MTANs to be homodimers with two catalytic sites near the dimer interface. Here, we explore the cooperative ligand interactions in the homodimer of Staphylococcus aureus MTAN (SaMTAN). Kinetic analysis indicated negative catalytic cooperativity. Titration of SaMTAN with the transition-state analogue MT-DADMe-ImmA gave unequal catalytic site binding, consistent with negative binding cooperativity. Thermodynamics of MT-DADMe-ImmA binding also gave negative cooperativity, where the first site had different enthalpic and entropic properties than the second site. Cysteine reactivity in a single-cysteine catalytic site loop construct of SaMTAN is reactive in native enzyme, less reactive when inhibitor is bound to one subunit, and nonreactive upon saturation with inhibitor. A fusion peptide heterodimer construct with one inactive subunit (E173Q) and one native subunit gave 25% of native SaMTAN activity, similar to native SaMTAN with MT-DADMe-ImmA at one catalytic site. Pre-steady-state kinetics showed fast chemistry at one catalytic site, consistent with slow adenine release before catalysis occurs at the second catalytic site. The results support the two catalytic sites acting sequentially, with negative cooperativity and product release being linked to motion of a catalytic site loop contributed by the neighboring subunit.

Change in heat capacity for enzyme catalysis determines temperature dependence of enzyme catalyzed rates

The increase in enzymatic rates with temperature up to an optimum temperature (Topt) is widely attributed to classical Arrhenius behavior, with the decrease in enzymatic rates above Topt ascribed to protein denaturation and/or aggregation. This account persists despite many investigators noting that denaturation is insufficient to explain the decline in enzymatic rates above Topt. Here we show that it is the change in heat capacity associated with enzyme catalysis (ΔC(‡)p) and its effect on the temperature dependence of ΔG(‡) that determines the temperature dependence of enzyme activity. Through mutagenesis, we demonstrate that the Topt of an enzyme is correlated with ΔC(‡)p and that changes to ΔC(‡)p are sufficient to change Topt without affecting the catalytic rate. Furthermore, using X-ray crystallography and molecular dynamics simulations we reveal the molecular details underpinning these changes in ΔC(‡)p. The influence of ΔC(‡)p on enzymatic rates has implications for the temperature dependence of biological rates from enzymes to ecosystems.

Thermodynamic analysis of transition-state features in picomolar inhibitors of human 5'-methylthioadenosine phosphorylase

Human 5'-methylthioadenosine phosphorylase (MTAP) is solely responsible for 5'-methylthioadenosine (MTA) metabolism to permit S-adenosylmethionine salvage. Transition-state (TS) analogues of MTAP are in development as anticancer candidates. TS analogues of MTAP incorporate a cationic nitrogen and a protonated 9-deazaadenine leaving group, which are mimics of the ribocation transition state. MT-ImmA and MT-DADMe-ImmA are two examples of these TS analogues. Thermodynamic analysis of MTA, inhibitor, and phosphate binding reveals the cationic nitrogen to provide -2.6 and -3.6 kcal/mol binding free energy for MT-ImmA and MT-DADMe-ImmA, respectively. The protonated deazaadenine provides an additional -1.3 (MT-ImmA) to -1.7 kcal/mol (MT-DADMe-ImmA). MT-DADMe-ImmA is a better match in TS geometry than MT-ImmA and is thermodynamically favored. Binding of TS analogues to the MTAP/phosphate complex is fully entropic, in contrast to TS analogue binding to the related human purine nucleoside phosphorylase/phosphate complex, which is fully enthalpic (Guan, R., Ho, M. C., Brenowitz, M., Tyler, P. C., Evans, G. B., Almo, S. C., and Schramm, V. L. (2011) Biochemistry 50, 10408-10417). The binding thermodynamics of phosphate or TS analogues alone to MTAP are fully dominated by enthalpy. Phosphate anchored in the catalytic site forms an ion pair with the cationic TS analogue to cause stabilization of the enzyme structure in the ternary complex. The ternary-induced conformational changes convert the individual enthalpic binding energies to entropy, resulting in a presumed shift of the protein architecture toward the transition state. Formation of the ternary TS analogue complex with MTAP induces a remarkable increase in thermal stability (ΔTm 28 °C). The enthalpic, entropic, and protein-stability features of TS analogue binding to human MTAP are resolved in these studies.

Crystal structures of the Helicobacter pylori MTAN enzyme reveal specific interactions between S-adenosylhomocysteine and the 5'-alkylthio binding subsite

The bacterial 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) enzyme is a multifunctional enzyme that catalyzes the hydrolysis of the N-ribosidic bond of at least four different adenosine-based metabolites: S-adenosylhomocysteine (SAH), 5'-methylthioadenosine (MTA), 5'-deoxyadenosine (5'-DOA), and 6-amino-6-deoxyfutalosine. These activities place the enzyme at the hub of seven fundamental bacterial metabolic pathways: S-adenosylmethionine (SAM) utilization, polyamine biosynthesis, the purine salvage pathway, the methionine salvage pathway, the SAM radical pathways, autoinducer-2 biosynthesis, and menaquinone biosynthesis. The last pathway makes MTAN essential for Helicobacter pylori viability. Although structures of various bacterial and plant MTANs have been described, the interactions between the homocysteine moiety of SAH and the 5'-alkylthiol binding site of MTAN have never been resolved. We have determined crystal structures of an inactive mutant form of H. pylori MTAN bound to MTA and SAH to 1.63 and 1.20 Å, respectively. The active form of MTAN was also crystallized in the presence of SAH, allowing the determination of the structure of a ternary enzyme-product complex resolved at 1.50 Å. These structures identify interactions between the homocysteine moiety and the 5'-alkylthiol binding site of the enzyme. This information can be leveraged for the development of species-specific MTAN inhibitors that prevent the growth of H. pylori.

Extreme entropy-enthalpy compensation in a drug-resistant variant of HIV-1 protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

Femtomolar inhibitors bind to 5'-methylthioadenosine nucleosidases with favorable enthalpy and entropy

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) catalyzes the hydrolytic cleavage of adenine from methylthioadenosine (MTA). Inhibitor design and synthesis informed by transition state analysis have developed femtomolar inhibitors for MTANs, among the most powerful known noncovalent enzyme inhibitors. Thermodynamic analyses of the inhibitor binding reveals a combination of highly favorable contributions from enthalpic (-24.7 to -4.0 kcal mol(-1)) and entropic (-10.0 to 6.4 kcal mol(-1)) interactions. Inhibitor binding to similar MTANs from different bacterial species gave distinct energetic contributions from similar catalytic sites. Thus, binding of four transition state analogues to EcMTAN and SeMTAN is driven primarily by enthalpy, while binding to VcMTAN is driven primarily by entropy. Human MTA phosphorylase (hMTAP) has a transition state structure closely related to that of the bacterial MTANs, and it binds tightly to some of the same transition state analogues. However, the thermodynamic signature of binding of an inhibitor to hMTAP differs completely from that with MTANs. We conclude that factors other than first-sphere catalytic residue contacts contribute to binding of inhibitors because the thermodynamic signature differs between bacterial species of the same enzyme.