Computational and experimental studies on the catalytic mechanism of biliverdin-IXbeta reductase

Utility of ToxTracker in animal alternative testing strategy for fragrance materials

To determine the utility of the ToxTracker assay in animal alternative testing strategies, the genotoxic potential of four fragrance materials (2-octen-4-one, lauric aldehyde, veratraldehyde, and p-methoxy cinnamaldehyde) were tested in the ToxTracker assay. These materials have been previously evaluated in an in vitro as well as in vivo micronucleus assay, conducted as per OECD guidelines. In addition to these studies, reconstructed human skin micronucleus studies were conducted on all four materials. All four materials were positive in an in vitro micronucleus assay but were negative in both in vivo and 3D skin micronucleus assays. The ToxTracker assay, in combination with in silico methods to predict metabolism was used to identify mechanisms for the misleading positive outcomes observed in the in vitro micronucleus assays. The results show that the ToxTracker assay, in conjunction with in silico predictions, can provide the information needed to aid in the identification of an appropriate animal alternative follow-up assay, for substances with positive results in the standard in vitro test battery. Thus, the ToxTracker assay is a valuable tool to identify the genotoxic potential of fragrance materials and can aid with replacing animal-based follow-up testing with appropriate animal alternative assay(s).

Quantitative interpretation of ToxTracker dose-response data for potency comparisons and mode-of-action determination

ToxTracker is an in vitro mammalian stem cell-based reporter assay that detects activation of specific cellular signaling pathways (DNA damage, oxidative stress, and/or protein damage) upon chemical exposure using flow cytometry. Here we used quantitative methods to empirically analyze historical control data, and dose-response data across a wide range of reference chemicals. First, we analyzed historical control data to define a fold-change threshold for identification of a significant positive response. Next, we used the benchmark dose (BMD) combined-covariate approach for potency ranking of a set of more than 120 compounds; the BMD values were used for comparative identification of the most potent inducers of each reporter. Lastly, we used principal component analysis (PCA) to investigate functional and statistical relationships between the ToxTracker reporters. The PCA results, based on the BMD results for all substances, indicated that the DNA damage (Rtkn, Bscl2) and p53 (Btg2) reporters are functionally complementary and indicative of genotoxic stress. The oxidative stress (Srxn1 and Blvrb) and protein stress (Ddit3) reporters are independent indicators of cellular stress, and essential for toxicological profiling using the ToxTracker assay. Overall, dose-response modeling of multivariate ToxTracker data can be used for potency ranking and mode-of-action determination. In the future, IVIVE (in vitro to in vivo extrapolation) methods can be employed to determine in vivo AED (administered equivalent dose) values that can in turn be used for human health risk assessment.

A tiered approach to investigate the inhalation toxicity of cobalt substances. Tier 2 b: Reactive cobalt substances induce oxidative stress in ToxTracker and activate hypoxia target genes

Cobalt metal and cobalt sulfate are carcinogenic in rodents following inhalation exposure. The pre-carcinogenic effects associated with exposure to these cobalt substances include oxidative stress and genotoxicity. Some, but not all, cobalt substances induce in vitro clastogenicity or an increase in micronuclei. As a result, these substances are classified genotoxic carcinogens, having major impacts on their risk assessment, e.g. assumption of a non-thresholded dose response. Here, we investigated the potential of nine cobalt substances to cause genotoxicity and oxidative stress using the ToxTracker assay, with an extension to measure biomarkers of hypoxia. None of the nine tested substances activated the DNA damage markers in ToxTracker, and five substances activated the oxidative stress response reporters. The same five substances also activated the expression of several hypoxia target genes. Consistent with the lower tier of testing found in the preceding paper of this series, these compounds can be grouped based on their ability to release bioavailable cobalt ion and to trigger subsequent key events.

The in vitro ToxTracker and Aneugen Clastogen Evaluation extension assay as a tool in the assessment of relative genotoxic potential of e-liquids and their aerosols

In vitro (geno)toxicity assessment of electronic vapour products (EVPs), relative to conventional cigarette, currently uses assays, including the micronucleus and Ames tests. Whilst informative on induction of a finite endpoint and relative risk posed by test articles, such assays could benefit from mechanistic supplementation. The ToxTracker and Aneugen Clastogen Evaluation analysis can indicate the activation of reporters associated with (geno)toxicity, including DNA damage, oxidative stress, the p53-related stress response and protein damage. Here, we tested for the different effects of a selection of neat e-liquids, EVP aerosols and Kentucky reference 1R6F cigarette smoke samples in the ToxTracker assay. The assay was initially validated to assess whether a mixture of e-liquid base components, propylene glycol (PG) and vegetable glycerine (VG) had interfering effects within the system. This was achieved by spiking three positive controls into the system with neat PG/VG or phosphate-buffered saline bubbled (bPBS) PG/VG aerosol (nicotine and flavour free). PG/VG did not greatly affect responses induced by the compounds. Next, when compared to cigarette smoke samples, neat e-liquids and bPBS aerosols (tobacco flavour; 1.6% freebase nicotine, 1.6% nicotine salt or 0% nicotine) exhibited reduced and less complex responses. Tested up to a 10% concentration, EVP aerosol bPBS did not induce any ToxTracker reporters. Neat e-liquids, tested up to 1%, induced oxidative stress reporters, thought to be due to their effects on osmolarity in vitro. E-liquid nicotine content did not affect responses induced. Additionally, spiking nicotine alone only induced an oxidative stress response at a supraphysiological level. In conclusion, the ToxTracker assay is a quick, informative screen for genotoxic potential and mechanisms of a variety of (compositionally complex) samples, derived from cigarettes and EVPs. This assay has the potential for future application in the assessment battery for next-generation (smoking alternative) products, including EVPs.

The multiple roles of coenzyme Q in cellular homeostasis and their relevance for the pathogenesis of coenzyme Q deficiency

Coenzyme Q (CoQ) is a redox active lipid that plays a central role in cellular homeostasis. It was discovered more than 60 years ago because of its role as electron transporter in the mitochondrial respiratory chain. Since then it has become evident that CoQ has many other functions, not directly related to bioenergetics. It is a cofactor of several mitochondrial dehydrogenases involved in the metabolism of lipids, amino acids, and nucleotides, and in sulfide detoxification. It is a powerful antioxidant and it is involved in the control of programmed cell death by modulating both apoptosis and ferroptosis. CoQ deficiency is a clinically and genetically heterogeneous group of disorders characterized by the impairment of CoQ biosynthesis. CoQ deficient patients display defects in cellular bioenergetics, but also in the other pathways in which CoQ is involved. In this review we will focus on the functions of CoQ not directly related to the respiratory chain, and on how their impairment is relevant for the pathophysiology of CoQ deficiency. A better understanding of the complex set of events triggered by CoQ deficiency will allow to design novel approaches for the treatment of this condition.

Divergent erythroid megakaryocyte fates in Blvrb-deficient mice establish non-overlapping cytoprotective functions during stress hematopoiesis

Cytoprotective mechanisms of heme oxygenases function by derivatizing heme to generate carbon monoxide, ferrous iron, and isomeric biliverdins, followed by rapid NAD(P)H-dependent biliverdin reduction to the antioxidant bilirubin using two non-overlapping biliverdin reductases that display biliverdin isomer-restricted redox activity. Although cytoprotective functions of heme oxygenases are widely recognized, concomitant effects of downstream biliverdin reductases remain incomplete. A computational model predicated on murine hematopoietic single-cell transcriptomic data identified Blvrb as a biological driver linked to the tumor necrosis factor stress pathway as a predominant source of variation defining hematopoietic cell heterogeneity. In vivo studies using Blvrb-deficient mice established the dispensable role of Blvrb in steady-state hematopoiesis, although model validation using aged Blvrb-deficient mice established an important cytoprotective function in stress hematopoiesis with dichotomous megakaryocyte-biased hematopoietic recovery. Defective stress erythropoiesis was evident in Blvrb^-/- spleens and in bone marrow erythroid development, occurring in conjunction with defective lipid peroxidation as a marker of oxidant mishandling. Cell autonomous effects on megakaryocyte lineage bias were documented using multipotential progenitor assays. These data provide the first physiological function of murine Blvrb in a non-redundant pathway of stress cytoprotection. Divergent effects on erythroid/megakaryocyte lineage speciation impute a novel redox-regulated mechanism for lineage partitioning.

Structure, dynamics and function of the evolutionarily changing biliverdin reductase B family

Biliverdin reductase B (BLVRB) family members are general flavin reductases critical in maintaining cellular redox with recent findings revealing that BLVRB alone can dictate cellular fate. However, as opposed to most enzymes, the BLVRB family remains enigmatic with an evolutionarily changing active site and unknown structural and functional consequences. Here, we applied a multi-faceted approach that combines X-ray crystallography, NMR and kinetics methods to elucidate the structural and functional basis of the evolutionarily changing BLVRB active site. Using a panel of three BLVRB isoforms (human, lemur and hyrax) and multiple human BLVRB mutants, our studies reveal a novel evolutionary mechanism where coenzyme 'clamps' formed by arginine side chains at two co-evolving positions within the active site serve to slow coenzyme release (Positions 14 and 78). We find that coenzyme release is further slowed by the weaker binding substrate, resulting in relatively slow turnover numbers. However, different BLVRB active sites imposed by either evolution or mutagenesis exhibit a surprising inverse relationship between coenzyme release and substrate turnover that is independent of the faster chemical step of hydride transfer also measured here. Collectively, our studies have elucidated the role of the evolutionarily changing BLVRB active site that serves to modulate coenzyme release and has revealed that coenzyme release is coupled to substrate turnover.

Modulating Enzyme Function via Dynamic Allostery within Biliverdin Reductase B

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic "dials" to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.

Exploring the FMN binding site in the mitochondrial outer membrane protein mitoNEET

MitoNEET is a mitochondrial outer membrane protein that hosts a redox active [2Fe-2S] cluster in the C-terminal cytosolic domain. Increasing evidence has shown that mitoNEET has an essential role in regulating energy metabolism in human cells. Previously, we reported that the [2Fe-2S] clusters in mitoNEET can be reduced by the reduced flavin mononucleotide (FMNH₂) and oxidized by oxygen or ubiquinone-2, suggesting that mitoNEET may act as a novel redox enzyme catalyzing electron transfer from FMNH₂ to oxygen or ubiquinone. Here, we explore the FMN binding site in mitoNEET by using FMN analogs and find that lumiflavin, like FMN, at nanomolar concentrations can mediate the redox transition of the mitoNEET [2Fe-2S] clusters in the presence of flavin reductase and NADH (100 μM) under aerobic conditions. The electron paramagnetic resonance (EPR) measurements show that both FMN and lumiflavin can dramatically change the EPR spectrum of the reduced mitoNEET [2Fe-2S] clusters and form a covalently bound complex with mitoNEET under blue light exposure, suggesting that FMN/lumiflavin has specific interactions with the [2Fe-2S] clusters in mitoNEET. In contrast, lumichrome, another FMN analog, fails to mediate the redox transition of the mitoNEET [2Fe-2S] clusters and has no effect on the EPR spectrum of the reduced mitoNEET [2Fe-2S] clusters under blue light exposure. Instead, lumichrome can effectively inhibit the FMNH₂-mediated reduction of the mitoNEET [2Fe-2S] clusters, indicating that lumichrome may act as a potential inhibitor to block the electron transfer activity of mitoNEET.

Aneugen Versus Clastogen Evaluation and Oxidative Stress-Related Mode-of-Action Assessment of Genotoxic Compounds Using the ToxTracker Reporter Assay

Understanding the mode-of-action (MOA) of genotoxic compounds and differentiating between direct DNA interaction and indirect genotoxicity is crucial for their reliable safety assessment. ToxTracker is a stem cell-based reporter assay that detects activation of various cellular responses that are associated with genotoxicity and cancer. ToxTracker consists of 6 different GFP reporter cell lines that can detect the induction of DNA damage, oxidative stress, and protein damage in a single test. The assay can thereby provide insight into the MOA of compounds. Genotoxicity is detected in ToxTracker by activation of 2 independent GFP reporters. Activation of the Bscl2-GFP reporter is associated with induction of DNA adducts and subsequent inhibition of DNA replication and the Rtkn-GFP reporter is activated following the formation of DNA double-strand breaks. Here, we show that the differential activation of these 2 genotoxicity reporters could be used to further differentiate between a DNA reactive and clastogenic or a non-DNA-reactive aneugenic MOA of genotoxic compounds. For further classification of aneugenic and clastogenic compounds, the ToxTracker assay was extended with cell cycle analysis and aneuploidy assessment. The extension was validated using a selection of 16 (genotoxic) compounds with a well-established MOA. Furthermore, indirect genotoxicity related to the production of reactive oxygen species was investigated using the DNA damage and oxidative stress ToxTracker reporters in combination with different reactive oxygen species scavengers. With these new extensions, ToxTracker was able to accurately classify compounds as genotoxic or nongenotoxic and could discriminate between DNA-reactive compounds, aneugens, and indirect genotoxicity caused by oxidative stress.

Biliverdin Reductase B Dynamics Are Coupled to Coenzyme Binding

Biliverdin reductase B (BLVRB) is a newly identified cellular redox regulator that catalyzes the NADPH-dependent reduction of multiple substrates. Through mass spectrometry analysis, we identified high levels of BLVRB in mature red blood cells, highlighting the importance of BLVRB in redox regulation. The BLVRB conformational changes that occur during conezyme/substrate binding and the role of dynamics in BLVRB function, however, remain unknown. Through a combination of NMR, kinetics, and isothermal titration calorimetry studies, we determined that BLVRB binds its coenzyme 500-fold more tightly than its substrate. While the active site of apo BLVRB is highly dynamic on multiple timescales, active site dynamics are largely quenched within holo BLVRB, in which dynamics are redistributed to other regions of the enzyme. We show that a single point mutation of Arg78➔Ala leads to both an increase in active site micro-millisecond motions and an increase in the microscopic rate constants of coenzyme binding. This demonstrates that altering BLVRB active site dynamics can directly cause a change in functional characteristics. Our studies thus address the solution behavior of apo and holo BLVRB and identify a role of enzyme dynamics in coenzyme binding.

In silico and crystallographic studies identify key structural features of biliverdin IXβ reductase inhibitors having nanomolar potency

Heme cytotoxicity is minimized by a two-step catabolic reaction that generates biliverdin (BV) and bilirubin (BR) tetrapyrroles. The second step is regulated by two non-redundant biliverdin reductases (IXα (BLVRA) and IXβ (BLVRB)), which retain isomeric specificity and NAD(P)H-dependent redox coupling linked to BR's antioxidant function. Defective BLVRB enzymatic activity with antioxidant mishandling has been implicated in metabolic consequences of hematopoietic lineage fate and enhanced platelet counts in humans. We now outline an integrated platform of in silico and crystallographic studies for the identification of an initial class of compounds inhibiting BLVRB with potencies in the nanomolar range. We found that the most potent BLVRB inhibitors contain a tricyclic hydrocarbon core structure similar to the isoalloxazine ring of flavin mononucleotide and that both xanthene- and acridine-based compounds inhibit BLVRB's flavin and dichlorophenolindophenol (DCPIP) reductase functions. Crystallographic studies of ternary complexes with BLVRB-NADP⁺-xanthene-based compounds confirmed inhibitor binding adjacent to the cofactor nicotinamide and interactions with the Ser-111 side chain. This residue previously has been identified as critical for maintaining the enzymatic active site and cellular reductase functions in hematopoietic cells. Both acridine- and xanthene-based compounds caused selective and concentration-dependent loss of redox coupling in BLVRB-overexpressing promyelocytic HL-60 cells. These results provide promising chemical scaffolds for the development of enhanced BLVRB inhibitors and identify chemical probes to better dissect the role of biliverdins, alternative substrates, and BLVRB function in physiologically relevant cellular contexts.

The mitochondrial outer membrane protein mitoNEET is a redox enzyme catalyzing electron transfer from FMNH2 to oxygen or ubiquinone

Increasing evidence suggests that mitoNEET, a target of the type II diabetes drug pioglitazone, is a key regulator of energy metabolism in mitochondria. MitoNEET is anchored to the mitochondrial outer membrane via its N-terminal α helix domain and hosts a redox-active [2Fe-2S] cluster in its C-terminal cytosolic region. The mechanism by which mitoNEET regulates energy metabolism in mitochondria, however, is not fully understood. Previous studies have shown that mitoNEET specifically interacts with the reduced flavin mononucleotide (FMNH₂) and that FMNH₂ can quickly reduce the mitoNEET [2Fe-2S] clusters. Here we report that the reduced mitoNEET [2Fe-2S] clusters can be readily oxidized by oxygen. In the presence of FMN, NADH, and flavin reductase, which reduces FMN to FMNH₂ using NADH as the electron donor, mitoNEET mediates oxidation of NADH with a concomitant reduction of oxygen. Ubiquinone-2, an analog of ubiquinone-10, can also oxidize the reduced mitoNEET [2Fe-2S] clusters under anaerobic or aerobic conditions. Compared with oxygen, ubiquinone-2 is more efficient in oxidizing the mitoNEET [2Fe-2S] clusters, suggesting that ubiquinone could be an intrinsic electron acceptor of the reduced mitoNEET [2Fe-2S] clusters in mitochondria. Pioglitazone or its analog NL-1 appears to inhibit the electron transfer activity of mitoNEET by forming a unique complex with mitoNEET and FMNH₂ The results suggest that mitoNEET is a redox enzyme that may promote oxidation of NADH to facilitate enhanced glycolysis in the cytosol and that pioglitazone may regulate energy metabolism in mitochondria by inhibiting the electron transfer activity of mitoNEET.

A substrate-bound structure of cyanobacterial biliverdin reductase identifies stacked substrates as critical for activity

Biliverdin reductase catalyses the last step in haem degradation and produces the major lipophilic antioxidant bilirubin via reduction of biliverdin, using NAD(P)H as a cofactor. Despite the importance of biliverdin reductase in maintaining the redox balance, the molecular details of the reaction it catalyses remain unknown. Here we present the crystal structure of biliverdin reductase in complex with biliverdin and NADP+. Unexpectedly, two biliverdin molecules, which we designated the proximal and distal biliverdins, bind with stacked geometry in the active site. The nicotinamide ring of the NADP+ is located close to the reaction site on the proximal biliverdin, supporting that the hydride directly attacks this position of the proximal biliverdin. The results of mutagenesis studies suggest that a conserved Arg185 is essential for the catalysis. The distal biliverdin probably acts as a conduit to deliver the proton from Arg185 to the proximal biliverdin, thus yielding bilirubin.

Enzymatic Activity and Thermodynamic Stability of Biliverdin IXβ Reductase Are Maintained by an Active Site Serine

Biliverdin reductase IXβ (BLVRB) is a crucial enzyme in heme metabolism. Recent studies in humans have identified a loss-of-function mutation (Ser111Leu) that unmasks a fundamentally important role in hematopoiesis. We have undertaken experimental and thermodynamic modeling studies to provide further insight into the role of the cofactor in substrate accessibility and protein folding properties regulating BLVRB catalytic mechanisms. Site-directed mutagenesis with molecular dynamic (MD) simulations establish the critical role of NAD(P)H-dependent conformational changes on substrate accessibility by forming the "hydrophobic pocket", along with identification of a single key residue (Arg35) modulating NADPH/NADH selectivity. Loop80 and Loop120 block the hydrophobic substrate binding pocket in apo BLVRB (open), whereas movement of these structures after cofactor binding results in the "closed" (catalytically active) conformation. Both enzymatic activity and thermodynamic stability are affected by mutation(s) involving Ser111, which is located in the core of the BLVRB active site. This work 1) elucidates the crucial role of Ser111 in enzymatic catalysis and thermodynamic stability by active site hydrogen bond network; 2) defines a dynamic model for apo BLVRB extending beyond the crystal structure of the binary BLVRB/NADP⁺ complex; 3) provides a structural basis for the "encounter" and "equilibrium" states of the binary complex, which are regulated by NAD(P)H.

Flavin nucleotides act as electron shuttles mediating reduction of the [2Fe-2S] clusters in mitochondrial outer membrane protein mitoNEET

MitoNEET, a primary target of type II diabetes drug pioglitazone, has an essential role in regulating energy metabolism, iron homeostasis, and production of reactive oxygen species in mitochondria. Structurally, mitoNEET is anchored to the mitochondrial outer membrane via its N-terminal transmembrane α-helix. The C-terminal cytosolic domain of mitoNEET hosts a redox active [2Fe-2S] cluster via three cysteine and one histidine residues. Here we report that the reduced flavin nucleotides can rapidly reduce the mitoNEET [2Fe-2S] clusters under anaerobic or aerobic conditions. In the presence of NADH and flavin reductase, 1 molecule of flavin nucleotide is sufficient to reduce about 100 molecules of the mitoNEET [2Fe-2S] clusters in 4min under aerobic conditions. The electron paramagnetic resonance (EPR) measurements show that flavin mononucleotide (FMN), but not flavin adenine dinucleotide (FAD), has a specific interaction with mitoNEET. Molecular docking models further reveal that flavin mononucleotide binds mitoNEET at the region between the N-terminal transmembrane α-helix and the [2Fe-2S] cluster binding domain. The closest distance between the [2Fe-2S] cluster and the bound flavin mononucleotide in mitoNEET is about 10Å, which could facilitate rapid electron transfer from the reduced flavin nucleotide to the [2Fe-2S] cluster in mitoNEET. The results suggest that flavin nucleotides may act as electron shuttles to reduce the mitoNEET [2Fe-2S] clusters and regulate mitochondrial functions in human cells.

Hydrophobic Shielding Drives Catalysis of Hydride Transfer in a Family of F420H2-Dependent Enzymes

A family of flavin/deazaflavin-dependent oxidoreductases (FDORs) from mycobacteria has been recently characterized and found to play a variety of catalytic roles, including the activation of prodrugs such as the candidate anti-tuberculosis drug pretomanid (PA-824). However, our understanding of the catalytic mechanism used by these enzymes is relatively limited. To address this, we have used a combination of quantum mechanics and molecular dynamics calculations to study the catalytic mechanism of the activation of pretomanid by the deazaflavin-dependent nitroreductase (Ddn) from Mycobacterium tuberculosis. The preferred pathway involves an initial hydride transfer step from the deprotonated cofactor (i.e., F₄₂₀H^-), with subsequent protonation, before a series of spontaneous intramolecular reactions to form the final reactive nitrogen species. The most likely proton source is a hydroxonium ion within the solvent accessible active site. Intriguingly, catalysis of the rate-determining hydride transfer step is aided by three tyrosine residues that form a hydrophobic barrier around the active site that, upon reaction, is then disrupted to allow increased water accessibility to facilitate the subsequent proton transfer step. The catalytic mechanism we propose is consistent with previous experimental observations of the Ddn enzyme and will inform the design of improved prodrugs in the future.

Rv2074 is a novel F420 H2 -dependent biliverdin reductase in Mycobacterium tuberculosis

Bilirubin is a potent antioxidant that is produced from the reduction of the heme degradation product biliverdin. In mammalian cells and Cyanobacteria, NADH/NADPH-dependent biliverdin reductases (BVRs) of the Rossmann-fold have been shown to catalyze this reaction. Here, we describe the characterization of Rv2074 from Mycobacterium tuberculosis, which belongs to a structurally and mechanistically distinct family of F420 H2 -dependent BVRs (F-BVRs) that are exclusively found in Actinobacteria. We have solved the crystal structure of Rv2074 bound to its cofactor, F420 , and used this alongside molecular dynamics simulations, site-directed mutagenesis and NMR spectroscopy to elucidate its catalytic mechanism. The production of bilirubin by Rv2074 could exploit the anti-oxidative properties of bilirubin and contribute to the range of immuno-evasive mechanisms that have evolved in M. tuberculosis to allow persistent infection.

Methylene Homologues of Artemisone: An Unexpected Structure-Activity Relationship and a Possible Implication for the Design of C10-Substituted Artemisinins

We sought to establish if methylene homologues of artemisone are biologically more active and more stable than artemisone. The analogy is drawn with the conversion of natural O- and N-glycosides into more stable C-glycosides that may possess enhanced biological activities and stabilities. Dihydroartemisinin was converted into 10β-cyano-10-deoxyartemisinin that was hydrolyzed to the α-primary amide. Reduction of the β-cyanide and the α-amide provided the respective methylamine epimers that upon treatment with divinyl sulfone gave the β- and α-methylene homologues, respectively, of artemisone. Surprisingly, the compounds were less active in vitro than artemisone against P. falciparum and displayed no appreciable activity against A549, HCT116, and MCF7 tumor cell lines. This loss in activity may be rationalized in terms of one model for the mechanism of action of artemisinins, namely the cofactor model, wherein the presence of a leaving group at C10 assists in driving hydride transfer from reduced flavin cofactors to the peroxide during perturbation of intracellular redox homeostasis by artemisinins. It is noted that the carba analogue of artemether is less active in vitro than the O-glycoside parent toward P. falciparum, although extrapolation of such activity differences to other artemisinins at this stage is not possible. However, literature data coupled with the leaving group rationale suggest that artemisinins bearing an amino group attached directly to C10 are optimal compounds.

The Extended ToxTracker Assay Discriminates Between Induction of DNA Damage, Oxidative Stress, and Protein Misfolding

Chemical exposure of cells may damage biomolecules, cellular structures, and organelles thereby jeopardizing cellular homeostasis. A multitude of defense mechanisms have evolved that can recognize specific types of damaged molecules and will initiate distinct cellular programs aiming to remove the damage inflicted and prevent cellular havoc. As a consequence, quantitative assessment of the activity of the cellular stress responses may serve as a sensitive reporter for the induction of specific types of damage. We have previously developed the ToxTracker assay, a mammalian stem cell-based genotoxicity assay employing two green fluorescent protein reporters specific for DNA damage and oxidative stress. We have now expanded the ToxTracker assay with an additional four reporter cell lines to include monitoring of additional stress signaling pathways. This panel of six green fluorescent protein reporters is able to discriminate between different primary reactivity of chemicals being their ability to react with DNA and block DNA replication, induce oxidative stress, activate the unfolded protein response, or cause a general P53-dependent cellular stress response. Extensive validation using the compound library suggested by the European Centre for the Validation of Alternative Methods (ECVAM) and a large panel of reference chemicals shows that the ToxTracker assay has an outstanding sensitivity and specificity. In addition, we developed Toxplot, a dedicated software tool for automated data analysis and graphical representation of the test results. Rapid and reliable identification by the ToxTracker assay of specific biological reactivity can significantly improve in vitro human hazard assessment of chemicals.

Interactions between artemisinins and other antimalarial drugs in relation to the cofactor model--a unifying proposal for drug action

Artemisinins are proposed to act in the malaria parasite cytosol by oxidizing dihydroflavin cofactors of redox-active flavoenzymes, and under aerobic conditions by inducing their autoxidation. Perturbation of redox homeostasis coupled with the generation of reactive oxygen species (ROS) ensues. Ascorbic acid-methylene blue (MB), N-benzyl-1,4-dihydronicotinamide (BNAH)-MB, BNAH-lumiflavine, BNAH-riboflavin (RF), and NADPH-FAD-E. coli flavin reductase (Fre) systems at pH 7.4 generate leucomethylene blue (LMB) and reduced flavins that are rapidly oxidized in situ by artemisinins. These oxidations are inhibited by the 4-aminoquinolines piperaquine (PPQ), chloroquine (CQ), and others. In contrast, the arylmethanols lumefantrine, mefloquine (MFQ), and quinine (QN) have little or no effect. Inhibition correlates with the antagonism exerted by 4-aminoquinolines on the antimalarial activities of MB, RF, and artemisinins. Lack of inhibition correlates with the additivity/synergism between the arylmethanols and artemisinins. We propose association via π complex formation between the 4-aminoquinolines and LMB or the dihydroflavins; this hinders hydride transfer from the reduced conjugates to the artemisinins. The arylmethanols have a decreased tendency to form π complexes, and so exert no effect. The parallel between chemical reactivity and antagonism or additivity/synergism draws attention to the mechanism of action of all drugs described herein. CQ and QN inhibit the formation of hemozoin in the parasite digestive vacuole (DV). The buildup of heme-Fe(III) results in an enhanced efflux from the DV into the cytosol. In addition, the lipophilic heme-Fe(III) complexes of CQ and QN that form in the DV are proposed to diffuse across the DV membrane. At the higher pH of the cytosol, the complexes decompose to liberate heme-Fe(III) . The quinoline or arylmethanol reenters the DV, and so transfers more heme-Fe(III) out of the DV. In this way, the 4-aminoquinolines and arylmethanols exert antimalarial activities by enhancing heme-Fe(III) and thence free Fe(III) concentrations in the cytosol. The iron species enter into redox cycles through reduction of Fe(III) to Fe(II) largely mediated by reduced flavin cofactors and likely also by NAD(P)H-Fre. Generation of ROS through oxidation of Fe(II) by oxygen will also result. The cytotoxicities of artemisinins are thereby reinforced by the iron. Other aspects of drug action are emphasized. In the cytosol or DV, association by π complex formation between pairs of lipophilic drugs must adversely influence the pharmacokinetics of each drug. This explains the antagonism between PPQ and MFQ, for example. The basis for the antimalarial activity of RF mirrors that of MB, wherein it participates in redox cycling that involves flavoenzymes or Fre, resulting in attrition of NAD(P)H. The generation of ROS by artemisinins and ensuing Fenton chemistry accommodate the ability of artemisinins to induce membrane damage and to affect the parasite SERCA PfATP6 Ca(2+) transporter. Thus, the effect exerted by artemisinins is more likely a downstream event involving ROS that will also be modulated by mutations in PfATP6. Such mutations attenuate, but cannot abrogate, antimalarial activities of artemisinins. Overall, parasite resistance to artemisinins arises through enhancement of antioxidant defense mechanisms.

Molecular modeling to provide insight into the substrate binding and catalytic mechanism of human biliverdin-IXα reductase

Human biliverdin-IXα reductase (hBVR-A) catalyzes the conversion of biliverdin-IXα to bilirubin-IXα in the last step of heme degradation and is a key enzyme in regulating a wide range of cellular responses. Though the X-ray structure of hBVR-A is available including cofactor, a crystal structure with a bound substrate would be even more useful as a starting point for protein-structure-based inhibitor design, but none have been reported. The present study employed induced fit docking (IFD) to study the substrate binding modes to hBVR-A of biliverdin-IXα and four analogues. The proposed substrate binding modes were examined further by performing molecular dynamics (MD) simulations followed by molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. The predicted binding free energies for the five biliverdin-IXα analogues match well with the relative potency of their reported experimental binding affinities, supporting that the proposed binding modes are reasonable. Furthermore, the ternary complex structure of hBVR-A binding with biliverdin-IXα and the electron donor cofactor NADPH obtained from MD simulations was exploited to investigate the catalytic mechanism, by calculating the reaction energy profile using the quantum mechanics/molecular mechanics (QM/MM) method. On the basis of our calculations, the energetically preferred pathway consists of an initial protonation of the pyrrolic nitrogen on the biliverdin substrate followed by hydride transfer to yield the reduction product. This conclusion is consistent with a previous mechanistic study on human biliverdin IXβ reductase (hBVR-B).

A partial convergence in action of methylene blue and artemisinins: antagonism with chloroquine, a reversal with verapamil, and an insight into the antimalarial activity of chloroquine

Artemisinins rapidly oxidize leucomethylene blue (LMB) to methylene blue (MB); they also oxidize dihydroflavins such as the reduced conjugates RFH₂ of riboflavin (RF), and FADH₂ of the cofactor flavin adenine dinucleotide (FAD), to the corresponding flavins. Like the artemisinins, MB oxidizes FADH₂, but unlike artemisinins, it also oxidizes NAD(P)H. Like MB, artemisinins are implicated in the perturbation of redox balance in the malaria parasite by interfering with parasite flavoenzyme disulfide reductases. The oxidation of LMB by artemisinin is inhibited by chloroquine (CQ), an inhibition that is abruptly reversed by verapamil (VP). CQ also inhibits artemisinin-mediated oxidation of RFH₂ generated from N-benzyl-1,4-dihydronicotinamide (BNAH)-RF, or FADH₂ generated from NADPH or NADPH-Fre, an effect that is also modulated by verapamil. The inhibition likely proceeds by the association of LMB or dihydroflavin with CQ, possibly involving donor-acceptor or π complexes that hinder oxidation by artemisinin. VP competitively associates with CQ, liberating LMB or dihydroflavin from their respective CQ complexes. The observations explain the antagonism between CQ-MB and CQ-artemisinins in vitro, and are reconcilable with CQ perturbing intraparasitic redox homeostasis. They further suggest that a VP-CQ complex is a means by which VP reverses CQ resistance, wherein such a complex is not accessible to the putative CQ-resistance transporter (PfCRT).

Structure and mechanism of sanguinarine reductase, an enzyme of alkaloid detoxification

Sanguinarine reductase is a plant enzyme that prevents the cytotoxic effects of benzophenanthridine alkaloids, which are the main phytoalexins of Papaveraceae. The enzyme catalyzes the reduction of sanguinarine, the most toxic benzophenanthridine, which re-enters the cytoplasm after its primary accumulation in the cell wall region has reached a threshold concentration. We present the sequence of the gene and protein of sanguinarine reductase isolated from cell cultures of Eschscholzia californica. High sequence similarities indicate that the enzyme evolved from a plant-specific branch of the ubiquitous Rossmann fold NAD(P)H/NAD(P)(+) binding reductases, with NADP-dependent epimerases or hydroxysteroid reductases as the most likely ancestors. Based on the x-ray structure of a close homolog, a three-dimensional model of the spatial conformation and catalytic site of sanguinarine reductase was established and used for in silico screening of known three-dimensional structures. Surprisingly, the enzyme shares high structural similarity with enzymes of human and bacterial origin, which have similar functions as the plant homologs but bear little amino acid sequence similarity. Using site-directed mutagenesis, a series of recombinant enzymes was generated and assayed to reveal the impact of individual amino acids and peptides in the catalytic process. It appears that relatively few innovations were required to generate this selective catalyst for alkaloid detoxication, notably an insertion of 13 amino acids and the generation of a novel catalytic triad of Cys-Asp-His were sufficient.

The effect of pH on the initial rate kinetics of the dimeric biliverdin-IXalpha reductase from the cyanobacterium Synechocystis PCC6803

Biliverdin-IXalpha reductase from Synechocystis PCC6803 (sBVR-A) is a stable dimer and this behaviour is observed under a range of conditions. This is in contrast to all other forms of BVR-A, which have been reported to behave as monomers, and places sBVR-A in the dihydrodiol dehydrogenase/N-terminally truncated glucose-fructose oxidoreductase structural family of dimers. The cyanobacterial enzyme obeys an ordered steady-state kinetic mechanism at pH 5, with NADPH being the first to bind and NADP(+) the last to dissociate. An analysis of the effect of pH on k(cat) with NADPH as cofactor reveals a pK of 5.4 that must be protonated for effective catalysis. Analysis of the effect of pH on k(cat)/K(m)(NADPH) identifies pK values of 5.1 and 6.1 in the free enzyme. Similar pK values are identified for biliverdin binding to the enzyme-NADPH complex. The lower pK values in the free enzyme (pK 5.1) and enzyme-NADPH complex (pK 4.9) are not evident when NADH is the cofactor, suggesting that this ionizable group may interact with the 2'-phosphate of NADPH. His84 is implicated as a crucial residue for sBVR-A activity because the H84A mutant has less than 1% of the activity of the wild-type and exhibits small but significant changes in the protein CD spectrum. Binding of biliverdin to sBVR-A is conveniently monitored by following the induced CD spectrum for biliverdin. Binding of biliverdin to wild-type sBVR-A induces a P-type spectrum. The H84A mutant shows evidence for weak binding of biliverdin and appears to bind a variant of the P-configuration. Intriguingly, the Y102A mutant, which is catalytically active, binds biliverdin in the M-configuration.