Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant unusually preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employ static and dynamic structural methods and observe that, compared to R132H, the R132Q active site adopts a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling reveals a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Catalytically distinct IDH1 mutants tune phenotype severity in tumor models

Mutations in isocitrate dehydrogenase 1 (IDH1) impart a neomorphic reaction that produces the oncometabolite D-2-hydroxyglutarate (D2HG), which can inhibit DNA and histone demethylases to drive tumorigenesis via epigenetic changes. Though heterozygous point mutations in patients primarily affect residue R132, there are myriad D2HG-producing mutants that display unique catalytic efficiency of D2HG production. Here, we show that catalytic efficiency of D2HG production is greater in IDH1 R132Q than R132H mutants, and expression of IDH1 R132Q in cellular and mouse xenograft models leads to higher D2HG concentrations in cells, tumors, and sera compared to R132H-expressing models. Reduced representation bisulfite sequencing (RRBS) analysis of xenograft tumors shows expression of IDH1 R132Q relative to R132H leads to hypermethylation patterns in pathways associated with DNA damage. Transcriptome analysis indicates that the IDH1 R132Q mutation has a more aggressive pro-tumor phenotype, with members of EGFR, Wnt, and PI3K signaling pathways differentially expressed, perhaps through non-epigenetic routes. Together, these data suggest that the catalytic efficiency of IDH1 mutants modulate D2HG levels in cellular and in vivo models, resulting in unique epigenetic and transcriptomic consequences where higher D2HG levels appear to be associated with more aggressive tumors.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant uniquely preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employed static and dynamic structural methods and found that, compared to R132H, the R132Q active site adopted a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling revealed a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant uniquely preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employed static and dynamic structural methods and found that, compared to R132H, the R132Q active site adopted a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling revealed a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Insights into Nitrosoalkane Binding to Myoglobin Provided by Crystallography of Wild-Type and Distal Pocket Mutant Derivatives

Nitrosoalkanes (R-N═O; R = alkyl) are biological intermediates that form from the oxidative metabolism of various amine (RNH₂) drugs or from the reduction of nitroorganics (RNO₂). RNO compounds bind to and inhibit various heme proteins. However, structural information on the resulting Fe-RNO moieties remains limited. We report the preparation of ferrous wild-type and H64A sw Mb^II-RNO derivatives (λ_max 424 nm; R = Me, Et, Pr, iPr) from the reactions of Mb^III-H₂O with dithionite and nitroalkanes. The apparent extent of formation of the wt Mb derivatives followed the order MeNO > EtNO > PrNO > iPrNO, whereas the order was the opposite for the H64A derivatives. Ferricyanide oxidation of the Mb^II-RNO derivatives resulted in the formation of the ferric Mb^III-H₂O precursors with loss of the RNO ligands. X-ray crystal structures of the wt Mb^II-RNO derivatives at 1.76-2.0 Å resoln. revealed N-binding of RNO to Fe and the presence of H-bonding interactions between the nitroso O-atoms and distal pocket His64. The nitroso O-atoms pointed in the general direction of the protein exterior, and the hydrophobic R groups pointed toward the protein interior. X-ray crystal structures for the H64A mutant derivatives were determined at 1.74-1.80 Å resoln. An analysis of the distal pocket amino acid surface landscape provided an explanation for the differences in ligand orientations adopted by the EtNO and PrNO ligands in their wt and H64A structures. Our results provide a good baseline for the structural analysis of RNO binding to heme proteins possessing small distal pockets.

Capturing the Dynamic Conformational Changes of Human Isocitrate Dehydrogenase 1 (IDH1) upon Ligand and Metal Binding Using Hydrogen-Deuterium Exchange Mass Spectrometry

Human isocitrate dehydrogenase 1 (IDH1) is a highly conserved metabolic enzyme that catalyzes the interconversion of isocitrate and α-ketoglutarate. Kinetic and structural studies with IDH1 have revealed evidence of striking conformational changes that occur upon binding of its substrates, isocitrate and NADP⁺, and its catalytic metal cation. Here, we used hydrogen-deuterium exchange mass spectrometry (HDX-MS) to build a comprehensive map of the dynamic conformational changes experienced by IDH1 upon ligand binding. IDH1 proved well-suited for HDX-MS analysis, allowing us to capture profound changes in solvent accessibility at substrate binding sites and at a known regulatory region, as well as at more distant local subdomains that appear to support closure of this protein into its active conformation. HDX-MS analysis suggested that IDH1 is primarily purified with NADP(H) bound in the absence of its metal cation. Subsequent metal cation binding, even in the absence of isocitrate, was critical for driving large conformational changes. WT IDH1 folded into its fully closed conformation only when the full complement of substrates and metal was present. Finally, we show evidence supporting a previously hypothesized partially open conformation that forms prior to the catalytically active state, and we propose this conformation is driven by isocitrate binding in the absence of metal.

Readying students for careers in industry: A guided inquiry activity to prepare students for success in biotechnology and pharmaceutical industry positions

While science students are well prepared for careers in biotechnology and pharmaceutical sciences in terms of technical expertise and critical thinking, they rarely have an opportunity to practice the due diligence required for success in industry in their coursework. This includes framing their expertise as solutions to challenges a company may be experiencing, an important skill for the interview process. As most academics have not applied for positions in industry, they may feel ill equipped to help students practice the important skills of framing their expertise within company goals and to discuss the business and financial concepts relevant to careers in scientific industry. Here, we describe an educational activity first developed by a leader in the biotech/pharmaceutical industry that was modified and given educational context by an academic in a class of upper-level undergraduate and graduate students. In this guided inquiry activity, students were instructed to select a start-up company in their field-ideally one to which they intended to apply for a job. Students were empowered by scaffolded hands-on exercises to research the company's scientific focus and finances, and to frame how their expertise could help companies achieve stated goals. Students compiled and delivered their research as an in-class presentation.

Abstract 3732: Characterization of tumor relevant isocitrate dehydrogenase 1 (IDH1)

Isocitrate dehydrogenase 1 (IDH1) is found commonly mutated in grade II and III gliomas and secondary glioblastomas. Normally, IDH1 catalyzes the conversion of isocitrate and NADP+ to α-ketoglutarate and NADPH. However, in mutated form, IDH1 catalyzes α-ketoglutarate and NADPH into NADP+ and D-2-hydroxyglutarate, an oncometabolite. D-2-hydroxyglutarate is toxic to cells and drives a number of pro-tumorigenic pathways. To understand the molecular mechanisms of how mutated IDH1 can lead to brain cancer, we created point mutations in IDH1 to understand the effects on enzyme activity. Site-directed mutagenesis was used to generate mutations at residue 132, the most common site of mutations in cancer patients. Wild type (WT) IDH1 and mutant IDH1 were heterologously expressed in bacteria . Following purification using affinity column chromatography, steady-state kinetic assays were used to compare the catalytic activity of WT versus mutant IDH1. By elucidating the types of mutations at R132 that facilitate D2HG production, we can have a better understanding of how IDH1 mutations affect prognosis and therapeutic response.

This work was funded by a Research Scholar Grant, RSG-19-075-01-TBE, from the American Cancer Society (CDS), National Institutes of Health R00 CA187594 (CDS), National Institutes of Health U54CA132384 (SDSU) & U54CA132379 (UC San Diego), MARC 5T34GM008303 (SDSU), and IMSD 5R25GM058906 (SDSU), and the California Metabolic Research Foundation (SDSU).

Citation Format: Nalani J. Coleman, Ella Thornberg, Danielle Caliger, Lucas Luna, Christal D. Sohl. Characterization of tumor relevant isocitrate dehydrogenase 1 (IDH1) [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3732.

Water Networks and Correlated Motions in Mutant Isocitrate Dehydrogenase 1 (IDH1) Are Critical for Allosteric Inhibitor Binding and Activity

Point mutations in human isocitrate dehydrogenase 1 (IDH1) can drive malignancies, including lower-grade gliomas and secondary glioblastomas, chondrosarcomas, and acute myeloid leukemias. These mutations, which usually affect residue R132, ablate the normal activity of catalyzing the NADP⁺-dependent oxidation of isocitrate to α-ketoglutarate (αKG) while also acquiring a neomorphic activity of reducing αKG to d-2-hydroxyglutarate (D2HG). Mutant IDH1 can be selectively therapeutically targeted due to structural differences that occur in the wild type (WT) versus mutant form of the enzyme, though the full mechanisms of this selectivity are still under investigation. Here we probe the mechanistic features of the neomorphic activity and selective small molecule inhibition through a new lens, employing WaterMap and molecular dynamics simulations. These tools identified a high-energy path of water molecules connecting the inhibitor binding site with the αKG and NADP⁺ binding sites in mutant IDH1. This water path aligns spatially with the α10 helix from WT IDH1 crystal structures. Mutating residues at the termini of this water path specifically disrupted inhibitor binding and/or D2HG production, revealing additional key residues to consider in optimizing druglike molecules against mutant IDH1. Taken together, our findings from molecular simulations and mutant enzyme kinetic assays provide insight into how disrupting water paths through enzyme active sites can impact not only inhibitor potency but also substrate recognition and activity.

Abstract 4381: Investigating the reversible MDH1 catalytic reaction in squamous non-small cell lung cancer

Non-small-cell lung cancer (NSCLC) is a leading cause of death in men and women worldwide. NSCLC relies on glycolysis to support cell proliferation and to synthesize proteins, nucleic acids and lipids. Malate dehydrogenase (MDH) is an enzyme critical for replenishing intermediates in the tricarboxylic acid (TCA) cycle and for helping to drive glycolysis through its reversible coenzyme system NAD+/NADH by interconverting oxaloacetate and malate. The cytosolic form of MDH (MDH1) has been reported to be amplified in squamous cell NSCLC, which has high glycolytic activity. Elucidation of the role of MDH1 in altering cellular metabolism to drive tumorigenic progression has yet to be determined. Thus, our interest lies in identifying the chemical and cellular characteristics required to drive the reversible MDH1 catalytic reaction in squamous cell NSCLC. Previously, we have shown that metabolic dehydrogenases are regulated by microenvironmental changes in pH and oxidation levels, and we are currently extending this work to MDH1. We are also evaluating various cell lines to gain a better understanding of how the MDH1 reaction can be altered to decrease the tumorigenic capacity of squamous cell NSCLC. Experimental knowledge obtained has the potential to increase our understanding of the MDH1 reaction and educate the design of future anti-tumor therapeutics.

Citation Format: Joi L. Weeks, Grace Wells, Sati Alexander, Christian Metallo, Christal D. Sohl. Investigating the reversible MDH1 catalytic reaction in squamous non-small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4381.

Activity of Selected Nucleoside Analogue ProTides against Zika Virus in Human Neural Stem Cells

Zika virus (ZIKV), an emerging flavivirus that causes neurodevelopmental impairment to fetuses and has been linked to Guillain-Barré syndrome continues to threaten global health due to the absence of targeted prophylaxis or treatment. Nucleoside analogues are good examples of efficient anti-viral inhibitors, and prodrug strategies using phosphate masking groups (ProTides) have been employed to improve the bioavailability of ribonucleoside analogues. Here, we synthesized and tested a small library of 13 ProTides against ZIKV in human neural stem cells. Strong activity was observed for 2′-C-methyluridine and 2′-C-ethynyluridine ProTides with an aryloxyl phosphoramidate masking group. Substitution of a 2-(methylthio) ethyl phosphoramidate for the aryloxyl phosphoramidate ProTide group of 2′-C-methyluridine completely abolished antiviral activity of the compound. The aryloxyl phosphoramidate ProTide of 2′-C-methyluridine outperformed the hepatitis C virus (HCV) drug sofosbuvir in suppression of viral titers and protection from cytopathic effect, while the former compound’s triphosphate active metabolite was better incorporated by purified ZIKV NS5 polymerase over time. These findings suggest both a nucleobase and ProTide group bias for the anti-ZIKV activity of nucleoside analogue ProTides in a disease-relevant cell model.

The FGFR1 V561M Gatekeeper Mutation Drives AZD4547 Resistance through STAT3 Activation and EMT

FGFR1 has been implicated in numerous cancer types including squamous cell lung cancer, a subset of non-small cell lung cancer with a dismal 5-year survival rate. Small-molecule inhibitors targeting FGFR1 are currently in clinical trials, with AZD4547 being one of the furthest along; however, the development of drug resistance is a major challenge for targeted therapies. A prevalent mechanism of drug resistance in kinases occurs through mutation of the gatekeeper residue, V561M in FGFR1; however, mechanisms underlying V561M resistance to AZD4547 are not fully understood. Here, the cellular consequences of the V561M gatekeeper mutation were characterized, and it was found that although AZD4547 maintains nanomolar affinity for V561M FGFR1, based on in vitro binding assays, cells expressing V561M demonstrate dramatic resistance to AZD4547 driven by increased STAT3 activation downstream of V561M FGFR1. The data reveal that the V561M mutation biases cells toward a more mesenchymal phenotype, including increased levels of proliferation, migration, invasion, and anchorage-independent growth, which was confirmed using CyTOF, a novel single-cell analysis tool. Using shRNA knockdown, loss of STAT3 restored sensitivity of cancer cells expressing V561M FGFR1 to AZD4547. Thus, the data demonstrate that combination therapies including FGFR and STAT3 may overcome V561M FGFR1-driven drug resistance in the clinic. IMPLICATIONS: The V561M FGFR1 gatekeeper mutation leads to devastating drug resistance through activation of STAT3 and the epithelial-mesenchymal transition; this study demonstrates that FGFR1 inhibitor sensitivity can be restored upon STAT3 knockdown.

Molecular mechanisms of isocitrate dehydrogenase 1 (IDH1) mutations identified in tumors: The role of size and hydrophobicity at residue 132 on catalytic efficiency

Isocitrate dehydrogenase 1 (IDH1) catalyzes the reversible NADP⁺-dependent conversion of isocitrate (ICT) to α-ketoglutarate (αKG) in the cytosol and peroxisomes. Mutations in IDH1 have been implicated in >80% of lower grade gliomas and secondary glioblastomas and primarily affect residue 132, which helps coordinate substrate binding. However, other mutations found in the active site have also been identified in tumors. IDH1 mutations typically result in a loss of catalytic activity, but many also can catalyze a new reaction, the NADPH-dependent reduction of αKG to d-2-hydroxyglutarate (D2HG). D2HG is a proposed oncometabolite that can competitively inhibit αKG-dependent enzymes. Some kinetic parameters have been reported for several IDH1 mutations, and there is evidence that mutant IDH1 enzymes vary widely in their ability to produce D2HG. We report that most IDH1 mutations identified in tumors are severely deficient in catalyzing the normal oxidation reaction, but that D2HG production efficiency varies among mutant enzymes up to ∼640-fold. Common IDH1 mutations have moderate catalytic efficiencies for D2HG production, whereas rarer mutations exhibit either very low or very high efficiencies. We then designed a series of experimental IDH1 mutants to understand the features that support D2HG production. We show that this new catalytic activity observed in tumors is supported by mutations at residue 132 that have a smaller van der Waals volume and are more hydrophobic. We report that one mutation can support both the normal and neomorphic reactions. These studies illuminate catalytic features of mutations found in the majority of patients with lower grade gliomas.

Abstract 2436: Understanding the molecular mechanism of targeted kinase inhibitor resistance mediated by the FGFR1 gatekeeper mutation

Human fibroblast growth factor receptors (FGFRs) can serve as drivers of oncogenesis and are a target of several inhibitors in clinical trials. Unfortunately, resistance severely hampers the long-term success of tyrosine kinase inhibitors (TKIs), with mutations at the gatekeeper residue often serving as an early means of resistance to cancer therapy, resulting in tumor progression. Here we show the first crystal structures of the FGFR1 gatekeeper mutation, V561M FGFR1, and use kinetic and structural methods to explore the mechanisms of the catalytic activation we observe for this mutation. We show that while the FGFR1 gatekeeper mutation confers resistance to two inhibitors currently in clinical trials, the extent of change in efficacy differ, indicating two unique resistance profiles. We show that some affinity for one inhibitor is maintained by the V561M mutation due to a flexible linker that allows multiple binding conformations. This is the first example showing the same inhibitor binding in unique ways to the WT and gatekeeper mutant forms of FGFR, highlighting some regions of drug design flexibility both within the binding pocket and in the linker of the inhibitor. Identifying this flexibility in the context of the structural features of an active conformation of V561M FGFR1, the form most kinase inhibitors target, will provide critical insights for designing future inhibitors effective against the FGFR gatekeeper mutations.

Citation Format: Christal D. Sohl, Molly Ryan, BeiBei Luo, Kathleen Frey, Karen S. Anderson. Understanding the molecular mechanism of targeted kinase inhibitor resistance mediated by the FGFR1 gatekeeper mutation. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2436. doi:10.1158/1538-7445.AM2015-2436

Illuminating the molecular mechanisms of tyrosine kinase inhibitor resistance for the FGFR1 gatekeeper mutation: the Achilles' heel of targeted therapy

Human fibroblast growth factor receptors (FGFRs) 1-4 are a family of receptor tyrosine kinases that can serve as drivers of tumorigenesis. In particular, FGFR1 gene amplification has been implicated in squamous cell lung and breast cancers. Tyrosine kinase inhibitors (TKIs) targeting FGFR1, including AZD4547 and E3810 (Lucitanib), are currently in early phase clinical trials. Unfortunately, drug resistance limits the long-term success of TKIs, with mutations at the "gatekeeper" residue leading to tumor progression. Here we show the first structural and kinetic characterization of the FGFR1 gatekeeper mutation, V561M FGFR1. The V561M mutation confers a 38-fold increase in autophosphorylation achieved at least in part by a network of interacting residues forming a hydrophobic spine to stabilize the active conformation. Moreover, kinetic assays established that the V561M mutation confers significant resistance to E3810, while retaining affinity for AZD4547. Structural analyses of these TKIs with wild type (WT) and gatekeeper mutant forms of FGFR1 offer clues to developing inhibitors that maintain potency against gatekeeper mutations. We show that AZD4547 affinity is preserved by V561M FGFR1 due to a flexible linker that allows multiple inhibitor binding modes. This is the first example of a TKI binding in distinct conformations to WT and gatekeeper mutant forms of FGFR, highlighting adaptable regions in both the inhibitor and binding pocket crucial for drug design. Exploiting inhibitor flexibility to overcome drug resistance has been a successful strategy for combatting diseases such as AIDS and may be an important approach for designing inhibitors effective against kinase gatekeeper mutations.

Fluorescence resonance energy transfer studies of DNA polymerase β: the critical role of fingers domain movements and a novel non-covalent step during nucleotide selection

During DNA repair, DNA polymerase β (Pol β) is a highly dynamic enzyme that is able to select the correct nucleotide opposite a templating base from a pool of four different deoxynucleoside triphosphates (dNTPs). To gain insight into nucleotide selection, we use a fluorescence resonance energy transfer (FRET)-based system to monitor movement of the Pol β fingers domain during catalysis in the presence of either correct or incorrect dNTPs. By labeling the fingers domain with ((((2-iodoacetyl)amino)ethyl)amino)naphthalene-1-sulfonic acid (IAEDANS) and the DNA substrate with Dabcyl, we are able to observe rapid fingers closing in the presence of correct dNTPs as the IAEDANS comes into contact with a Dabcyl-labeled, one-base gapped DNA. Our findings show that not only do the fingers close after binding to the correct dNTP, but that there is a second conformational change associated with a non-covalent step not previously reported for Pol β. Further analyses suggest that this conformational change corresponds to the binding of the catalytic metal into the polymerase active site. FRET studies with incorrect dNTP result in no changes in fluorescence, indicating that the fingers do not close in the presence of incorrect dNTP. Together, our results show that nucleotide selection initially occurs in an open fingers conformation and that the catalytic pathways of correct and incorrect dNTPs differ from each other. Overall, this study provides new insight into the mechanism of substrate choice by a polymerase that plays a critical role in maintaining genome stability.

Probing the molecular mechanism of action of the HIV-1 reverse transcriptase inhibitor 4'-ethynyl-2-fluoro-2'-deoxyadenosine (EFdA) using pre-steady-state kinetics

The novel antiretroviral 4'-ethynyl-2-fluoro-2'-deoxyadenosine (EFdA) is a potent nucleoside HIV-1 reverse transcriptase (RT) inhibitor (NRTI). Unlike other FDA-approved NRTIs, EFdA contains a 3'-hydroxyl. Pre-steady-state kinetics showed RT preferred incorporating EFdA-TP over native dATP. Moreover, RT slowly inserted nucleotides past an EFdA-terminated primer, resulting in delayed chain termination with unaffected fidelity. This is distinct from KP1212, another 3'-hydroxyl-containing RT inhibitor considered to promote viral lethal mutagenesis. New mechanistic features of RT inhibition by EFdA are revealed.

Mutations in human DNA polymerase γ confer unique mechanisms of catalytic deficiency that mirror the disease severity in mitochondrial disorder patients

Human mitochondrial DNA polymerase γ (pol γ) is solely responsible for the replication and repair of the mitochondrial genome. Unsurprisingly, alterations in pol γ activity have been associated with mitochondrial diseases such as Alpers syndrome and progressive external ophthalmoplegia. Thus far, predicting the severity of mitochondrial disease based the magnitude of deficiency in pol γ activity has been difficult. In order to understand the relationship between disease severity in patients and enzymatic defects in vitro, we characterized the molecular mechanisms of four pol γ mutations, A957P, A957S, R1096C and R1096H, which have been found in patients suffering from aggressive Alpers syndrome to mild progressive external ophthalmoplegia. The A957P mutant showed the most striking deficiencies in the incorporation efficiency of a correct deoxyribonucleotide triphosphate (dNTP) relative to wild-type pol γ, with less, but still significant incorporation efficiency defects seen in R1096H and R1096C, and only a small decrease in incorporation efficiency observed for A957S. Importantly, this trend matches the disease severity observed in patients very well (approximated as A957P ≫ R1096C ≥ R1096H ≫ A957S, from most severe disease to least severe). Further, the A957P mutation conferred a two orders of magnitude loss of fidelity relative to wild-type pol γ, indicating that a buildup of mitochondrial genomic mutations may contribute to the death in infancy seen with these patients. We conclude that characterizing the unique molecular mechanisms of pol γ deficiency for physiologically important mutant enzymes is important for understanding mitochondrial disease and for predicting disease severity.

Oxidation of dihydrotestosterone by human cytochromes P450 19A1 and 3A4

Dihydrotestosterone is a more potent androgen than testosterone and plays an important role in endocrine function. We demonstrated that, like testosterone, dihydrotestosterone can be oxidized by human cytochrome P450 (P450) 19A1, the steroid aromatase. The products identified include the 19-hydroxy- and 19-oxo derivatives and the resulting Δ(1,10)-, Δ(5,10)-, and Δ(9,10)-dehydro 19-norsteroid products (loss of 19-methyl group). The overall catalytic efficiency of oxidation was ~10-fold higher than reported for 3α-reduction by 3α-hydroxysteroid dehydrogenase, the major enzyme known to deactivate dihydrotestosterone. These and other studies demonstrate the flexibility of P450 19A1 in removing the 1- and 2-hydrogens from 19-norsteroids, the 2-hydrogen from estrone, and (in this case) the 1-, 5β-, and 9β-hydrogens of dihydrotestosterone. Incubation of dihydrotestosterone with human liver microsomes and NADPH yielded the 18- and 19-hydroxy products plus the Δ(1,10)-dehydro 19-nor product identified in the P450 19A1 reaction. The 18- and 19-hydroxylation reactions were attributed to P450 3A4, and 18- and 19-hydroxydihydrotestosterone were identified in human plasma and urine samples. The change in the pucker of the A ring caused by reduction of the Δ(4,5) bond is remarkable in shifting the course of hydroxylation from the 6β-, 2β-, 1β-, and 15β-methylene carbons (testosterone) to the axial methyl groups (18, 19) in dihydrotestosterone and demonstrates the sensitivity of P450 3A4, even with its large active site, to small changes in substrate structure.

Balancing antiviral potency and host toxicity: identifying a nucleotide inhibitor with an optimal kinetic phenotype for HIV-1 reverse transcriptase

Two novel thymidine analogs, 3'-fluoro-3'-deoxythymidine (FLT) and 2',3'-didehydro-3'-deoxy-4'-ethynylthymidine (Ed4T), have been investigated as nucleoside reverse transcriptase inhibitors (NRTIs) for treatment of HIV infection. Ed4T seems very promising in phase II clinical trials, whereas toxicity halted FLT development during this phase. To understand these different molecular mechanisms of toxicity, pre-steady-state kinetic studies were used to examine the interactions of FLT and Ed4T with wild-type (WT) human mitochondrial DNA polymerase γ (pol γ), which is often associated with NRTI toxicity, as well as the viral target protein, WT HIV-1 reverse transcriptase (RT). We report that Ed4T-triphosphate (TP) is the first analog to be preferred over native nucleotides by RT but to experience negligible incorporation by WT pol γ, with an ideal balance between high antiretroviral efficacy and minimal host toxicity. WT pol γ could discriminate Ed4T-TP from dTTP 12,000-fold better than RT, with only an 8.3-fold difference in discrimination being seen for FLT-TP. A structurally related NRTI, 2',3'-didehydro-2',3'-dideoxythymidine, is the only other analog favored by RT over native nucleotides, but it exhibits only a 13-fold difference (compared with 12,000-fold for Ed4T) in discrimination between the two enzymes. We propose that the 4'-ethynyl group of Ed4T serves as an enzyme selectivity moiety, critical for discernment between RT and WT pol γ. We also show that the pol γ mutation R964C, which predisposes patients to mitochondrial toxicity when receiving 2',3'-didehydro-2',3'-dideoxythymidine to treat HIV, produced some loss of discrimination for FLT-TP and Ed4T-TP. These molecular mechanisms of analog incorporation, which are critical for understanding pol γ-related toxicity, shed light on the unique toxicity profiles observed during clinical trials.

Multi-step oxidations catalyzed by cytochrome P450 enzymes: Processive vs. distributive kinetics and the issue of carbonyl oxidation in chemical mechanisms

Catalysis of sequential oxidation reactions is not unusual in cytochrome P450 (P450) reactions, not only in steroid metabolism but also with many xenobiotics. One issue is how processive/distributive these reactions are, i.e., how much do the "intermediate" products dissociate. Our work with human P450s 2E1, 2A6, and 19A1 on this subject has revealed a mixture of systems, surprisingly with a more distributive mechanism with an endogenous substrate (P450 19A1) than for some xenobiotics (P450s 2E1, 2A6). One aspect of this research involves carbonyl intermediates, and the choice of catalytic mechanism is linked to the hydration state of the aldehyde. The non-enzymatic rates of hydration and dehydration of carbonyls are not rapid and whether P450s catalyze the reversible hydration is unknown. If carbonyl hydration and dehydration are slow, the mechanism may be set by the carbonyl hydration status.

Kinetic analysis of the three-step steroid aromatase reaction of human cytochrome P450 19A1

Cytochrome P450 19A1 (P450 19A1), the aromatase, catalyzes the conversion of androgens to estrogens through a sequential three-step reaction, generating 19-hydroxy and 19-aldehyde intermediates en route to the product estrogen. A procedure for the heterologous expression and purification of P450 19A1 in Escherichia coli was developed (k(cat) of 0.06 s(-1) for the conversion of androstenedione to estrone). Binding of the substrate and intermediates show low micromolar dissociation constants and are at least two-step processes. Rates of reduction of the iron were fast in the presence of substrate, either intermediate, or product. P450 19A1 is a distributive rather than a processive enzyme, with the sequential reaction allowing free dissociation of the intermediates as revealed by pulse-chase experiments. Conversion of androstenedione to estrone (under single turnover conditions) generated a progress curve showing changes in the concentrations of the substrate, intermediates, and product. A minimal kinetic model containing the individual rate constants for the steps in P450 19A1 catalysis was developed to globally fit the time course of the overall reaction, the dissociation constants, the two-step ligand binding, the distributive character, the iron-reduction rates, and the steady-state conversion of the 19-hydroxy androstenedione and 19-aldehyde androstenedione intermediates to estrone.

Measurement of cytochrome P450 and NADPH-cytochrome P450 reductase

Cytochrome P450 (P450) enzymes are important in the metabolism of steroids, vitamins, carcinogens, drugs and other compounds. Two of the commonly used assays in this field are the measurements of total P450 and NADPH-P450 reductase in biological preparations. A detailed protocol is presented for the measurement of P450 by its spectral properties, along with a protocol for measuring NADPH-P450 reductase by its NADPH-cytochrome c reduction activity. Each assay can be completed in 5-10 min. Detailed explanations for the rationale of particular sequences in the protocols are provided, along with potential confounding problems.

Chromatographic assays of drug oxidation by human cytochrome P450 3A4

Cytochrome P450 enzymes (P450s) are heme-thiolate mono-oxygenases involved in the oxidation of many endogenous and exogenous substrates. Herein, we describe two protocols for measuring the activity of a key enzyme of drug metabolism, P450 3A4. In this protocol, the substrate is incubated with human liver microsomes, the reaction is quenched, and the substrates and products are extracted and subjected to liquid chromatography (LC) separation and detection. Oxidation of the calcium-channel blocker nifedipine is measured using UV-Vis spectroscopy in-line with high performance liquid chromatography (HPLC). 6beta-Hydroxytestosterone formation from testosterone is measured by HPLC coupled to mass spectrometry (MS). Both of these procedures are rapid, requiring 2 h or less, and can be used to confirm and measure P450 3A4 activity and can also be used as a guide for developing other assays for measuring P450 catalysis. The separation strategy described here is more rapid than many available methods, except when ultra-performance liquid chromatography (UPLC) is used.

High-throughput fluorescence assay of cytochrome P450 3A4

Cytochrome P450 mono-oxygenases (P450s) are the principal enzymes involved in the oxidative metabolism of drugs and other xenobiotics. In this protocol, we describe a fluorescence-based, high-throughput assay for measuring the activity of P450 3A4, one of the key enzymes involved in drug metabolism. The assay involves the oxidative debenzylation of a substituted coumarin, yielding an increase in fluorescence on reaction. The entire procedure can be accomplished in 1 h or less.

Cooperativity of cytochrome P450 1A2: interactions of 1,4-phenylene diisocyanide and 1-isopropoxy-4-nitrobenzene

Homotropic cooperativity of 1-alkoxy-4-nitrobenzene substrates and also their heterotropic cooperative binding interactions with the iron ligand 1,4-phenylene diisocyanide (Ph(NC)2) had been demonstrated previously with rabbit cytochrome P450 (P450) 1A2 [G.P. Miller, F.P. Guengerich, Biochemistry 40 (2001) 7262-7272]. Multiphasic kinetics were observed for the binding of Ph(NC)2 to both ferric and ferrous P450 1A2, including relatively slow steps. Ph(NC)2 induced an apparently rapid change in the circular dichroism spectrum, consistent with a structural change, but had no effect on tryptophan fluorescence. Ph(NC)2 binds the P450 iron in both the ferric and ferrous forms; ferric P450 1A2 was reduced rapidly in the absence of added ligands, and the rate was attenuated when Ph(NC)2 was bound. No oxidation products of Ph(NC)2 were detected. Docking studies with a rabbit P450 1A2 homology model based on the published structure of a human P450 1A2.alpha-naphthoflavone (alphaNF) complex indicated adequate room for a complex with either two 1-isopropoxy-4-nitrobenzene molecules or a combination of one 1-isopropoxy-4-nitrobenzene and one Ph(NC)2; in the case of alphaNF no space for an extra ligand was available. The patterns of homotropic cooperativity seen with 1-alkoxy-4-nitrobenzenes (biphasic plots of v vs. S) differ from those seen with polycyclic hydrocarbons (positive cooperativity), suggesting that only with the latter does the ligand interaction produce improved catalysis. Consistent with this view, Ph(NC)2 inhibited the oxidation of 1-isopropoxy-4-nitrobenzene and other substrates.

Cooperativity in oxidation reactions catalyzed by cytochrome P450 1A2: highly cooperative pyrene hydroxylation and multiphasic kinetics of ligand binding

Rabbit liver cytochrome P450 (P450) 1A2 was found to catalyze the 5,6-epoxidation of alpha-naphthoflavone (alphaNF), 1-hydroxylation of pyrene, and the subsequent 6-, 8-, and other hydroxylations of 1-hydroxy (OH) pyrene. Plots of steady-state rates of product formation versus substrate concentration were hyperbolic for alphaNF epoxidation but highly cooperative (Hill n coefficients of 2-4) for pyrene and 1-OH pyrene hydroxylation. When any of the three substrates (alphaNF, pyrene, 1-OH pyrene) were mixed with ferric P450 1A2 using stopped-flow methods, the changes in the heme Soret spectra were relatively slow and multiphasic. Changes in the fluorescence of all of the substrates were much faster, consistent with rapid initial binding to P450 1A2 in a manner that does not change the heme spectrum. For binding of pyrene to ferrous P450 1A2, the course of the spectra revealed sequential changes in opposite directions, consistent with P450 1A2 being involved in a series of transitions to explain the kinetic multiphasicity as opposed to multiple, slowly interconverting populations of enzyme undergoing the same event at different rates. Models of rabbit P450 1A2 based on a published crystal structure of a human P450 1A2-alphaNF complex show active site space for only one alphaNF or for two pyrenes. The spectral changes observed for binding and hydroxylation of pyrene and 1-OH pyrene could be fit to a kinetic model in which hydroxylation occurs only when two substrates are bound. Elements of this mechanism may be relevant to other cases of P450 cooperativity.

Recombinant Enzymes Overexpressed in Bacteria Show Broad Catalytic Specificity of Human Cytochrome P450 2W1 and Limited Activity of Human Cytochrome P450 2S1

Human cytochromes P450 2S1 and 2W1 have received only limited attention with regard to characterization of function. Both cytochromes P450 have been reported to be overexpressed in human tumors, and cytochrome P450 2S1 is induced by carcinogenic polycyclic hydrocarbons. We report methods for high-level expression and purification of both cytochromes P450 from Escherichia coli, with the goal of establishing function. The level of expression of human cytochrome P450 2W1 achieved using codon optimization for E. coli was 1800 nmol of cytochrome P450 per liter of culture, the highest level achieved in this laboratory to date. Assays with a number of the typical cytochrome P450 substrates showed no detectable activity, including some for which qualitative reports have appeared in the literature. Cytochrome P450 2W1 catalyzed benzphetamine N-demethylation (k(cat), 3.8/min) and arachidonic acid oxidation, albeit at a very low rate (approximately 0.05/min). In a umu genotoxicity screen, cytochrome P450 2W1 catalyzed the activation of several procarcinogens, particularly polycyclic hydrocarbon diols, but cytochrome P450 2S1 did not. The bioactivation of procarcinogens by cytochrome P450 2W1 may be of significance in the context of reports of preferential expression of the enzyme in tumors, in that activation of procarcinogens could lead to the accumulation of mutations and enhance the carcinogenic process.

Synthesis and solid-state molecular structures of nitrosoalkane complexes of iron porphyrins containing methanol, pyridine, and 1-methylimidazole ligands

Nitrosoalkanes belong to the family of C-nitroso compounds and are known to bind to the iron center in heme proteins. We have prepared and characterized a series of new nitrosoalkane heme model complexes of the form (por)Fe(RNO)(L) (por=porphyrinato dianion; R=isopropyl; L=MeOH, pyridine, 1-methylimidazole) by infrared and 1H NMR spectroscopy and X-ray crystallography. Within the set of octaethylporphyrinato (OEP) compounds, the infrared stretching frequencies of the NO groups decrease in the order (OEP)Fe(iPrNO)(MeOH).MeOH (1433 cm-1) > (OEP)Fe(iPrNO)(py) (1429 cm-1) > (OEP)Fe(iPrNO)(1-MeIm) (1423 cm-1), reflecting the increased backdonation of electron density in the 1-methylimidazole derivative. The molecular structures of the compounds as determined by crystallography reveal N-binding of the nitrosoalkane ligands to the formally ferrous metal centers.