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.

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.

Suppression of human alloreactive T cells by linear tetrapyrroles; relevance for transplantation.

The main limitation to successful transplantation is the antigraft response developed by the recipient immune system, and the adverse side effects of immunosuppressive agents which are associated with significant toxicity and counter indications such as infection and cancer. Furthermore, immunosuppressants do little to prevent ischemia-reperfusion injury during the transplantation procedure itself hence there is a growing need to develop novel immunosuppressive drugs specifically aimed at prolonging graft survival. Linear tetrapyrroles derived from the breakdown of mammalian heme have been shown in numerous studies to play a protective role in allograft transplantation and ischemia-reperfusion injury; however, commercial sources of these products have not been approved for use in humans. Plants and algae produce equivalent linear tetrapyrroles called bilins that serve as chromophores in light-sensing. One such marine-derived tetrapyrrole, phycocyanobilin (PCB), shows significant structural similarity to mammalian biliverdin (BV) and may prove to be a safer alternative for use in the clinic if it can exert direct effects on human immune cells. Using a mixed lymphocyte reaction, we quantified the allogeneic responses of recipient cells to donor cells and found that PCB, like BV, effectively suppressed proliferation and proinflammatory cytokine production. In addition, we found that BV and PCB can directly downregulate the proinflammatory responses of both innate dendritic cells and adaptive T cells. We therefore propose that PCB may be an effective therapeutic drug in the clinical setting of transplantation and may also have wider applications in regulating inappropriate inflammation.

BLVRB redox mutation defines heme degradation in a metabolic pathway of enhanced thrombopoiesis in humans.

Human blood cell counts are tightly maintained within narrow physiologic ranges, largely controlled by cytokine-integrated signaling and transcriptional circuits that regulate multilineage hematopoietic specification. Known genetic loci influencing blood cell production account for <10% of platelet and red blood cell variability, and thrombopoietin/cellular myeloproliferative leukemia virus liganding is dispensable for definitive thrombopoiesis, establishing that fundamentally important modifier loci remain unelucidated. In this study, platelet transcriptome sequencing and extended thrombocytosis cohort analyses identified a single loss-of-function mutation (BLVRB(S111L)) causally associated with clonal and nonclonal disorders of enhanced platelet production. BLVRB(S111L) encompassed within the substrate/cofactor [α/ß dinucleotide NAD(P)H] binding fold is a functionally defective redox coupler using flavin and biliverdin (BV) IXß tetrapyrrole(s) and results in exaggerated reactive oxygen species accumulation as a putative metabolic signal leading to differential hematopoietic lineage commitment and enhanced thrombopoiesis. These data define the first physiologically relevant function of BLVRB and implicate its activity and/or heme-regulated BV tetrapyrrole(s) in a unique redox-regulated bioenergetic pathway governing terminal megakaryocytopoiesis; these observations also define a mechanistically restricted drug target retaining potential for enhancing human platelet counts.

Site-directed mutagenesis of mouse glutathione transferase P1-1 unlocks masked cooperativity, introduces a novel mechanism for 'ping pong' kinetic behaviour, and provides further structural evidence for participation of a water molecule in proton abstraction from glutathione.

Mouse liver glutathione transferase P1-1 has three cysteine residues at positions 14, 47 and 169. We have constructed the single, double and triple cysteine to alanine mutants to define the behaviour of all three thiols. We confirm that C47 is the 'fast' thiol (pK 7.4), and define C169 as the alkaline reactive residue with a pK(a) of 8.6. Only a small proportion of C14 is reactive with 5,5'-dithiobis-(2-nitrobenoic acid) (DTNB) at pH 9 in the C47A/C169A double mutant. The native enzyme and the C169A mutant exhibited Michaelis-Menten kinetics, but all other thiol to alanine mutants exhibited sigmoidal kinetics to varying degrees. The C169A mutant exhibited 'ping pong' kinetics, consistent with a mechanism whereby liberation of a proton from a reduced enzyme-glutathione (GSH) complex to form an enzyme-GS(-) (unprotonated) complex is essentially irreversible. Intriguingly, similar behaviour has recently been reported for a mutant of the yeast prion Ure2p. This cooperative behaviour is 'mirrored' in the crystal structure of the C47A mutant, which binds the p-nitrobenzyl moiety of p-nitrobenzyglutathione in distinct orientations in the two crystallographic subunits. The asymmetry seen in this structure for product binding is associated with absence of a water molecule W0 in the standard wild-type conformation of product binding that is clearly identifiable in the new structure, which may represent a structural model for binding of incoming GSH prior to displacement of W0. Elimination of W0 as a hydroxonium ion may be the mechanism for the initial proton extrusion from the active site.

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.

The use of synthetic linear tetrapyrroles to probe the verdin sites of human biliverdin-IXalpha reductase and human biliverdin-IXbeta reductase.

Many vertebrate species express two enzymes that are capable of catalysing the reduction of various isomers of biliverdin. Biliverdin-IXalpha reductase (BVR-A) is most active with its physiological substrate biliverdin-IXalpha, but can also reduce the three other biliverdin isomers IXbeta, IXdelta and IXgamma. Biliverdin-IXbeta reductase (BVR-B) catalyses the reduction of only the IXbeta, IXdelta and IXgamma isomers of biliverdin. Therefore, the activity of BVR-A can be measured using biliverdin-IXalpha as a specific substrate. We now show that the dimethyl esters of biliverdin-IXbeta and biliverdin-IXdelta are substrates for BVR-B, but not for BVR-A. This provides a useful method for specifically assaying the activity of both BVR-A and BVR-B in crude mixtures, using biliverdin-IXalpha for BVR-A and the dimethyl ester of either biliverdin-IXbeta or biliverdin-IXdelta for BVR-B. Human BVR-A has been suggested as a pharmacological target for neonatal jaundice. Because of the absence of a crystal structure with biliverdin bound to BVR-A, we have investigated indirect ways of examining tetrapyrrole binding. In the present study, we report that a number of sterically locked conformers of 18-ethylbiliverdin-IXalpha are substrates for human BVR-A, and discuss the implications for the biliverdin binding site. The oxidation of bilirubin-IXalpha ditaurate to biliverdin-IXalpha ditaurate is also described. We show that biliverdin-IXalpha ditaurate is a substrate for human BVR-A and discuss the possibility of using a competing substrate, which is reduced to a water soluble and excretable rubin, as a prototypic inhibitor of BVR-A.

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

BVR-B (biliverdin-IXbeta reductase) also known as FR (flavin reductase) is a promiscuous enzyme catalysing the pyridine-nucleotide-dependent reduction of a variety of flavins, biliverdins, PQQ (pyrroloquinoline quinone) and ferric ion. Mechanistically it is a good model for BVR-A (biliverdin-IXalpha reductase), a potential pharmacological target for neonatal jaundice and also a potential target for adjunct therapy to maintain protective levels of biliverdin-IXalpha during organ transplantation. In a commentary on the structure of BVR-B it was noted that one outstanding issue remained: whether the mechanism was a concerted hydride transfer followed by protonation of a pyrrolic anion or protonation of the pyrrole followed by hydride transfer. In the present study we have attempted to address this question using QM/MM (quantum mechanics/molecular mechanics) calculations. QM/MM potential energy surfaces show that the lowest energy pathway proceeds with a positively charged pyrrole intermediate via two transition states. These initial calculations were performed with His(153) as the source of the proton. However site-directed mutagenesis studies with both the H153A and the H153N mutant reveal that His(153) is not required for catalytic activity. We have repeated the calculation with a solvent hydroxonium donor and obtain a similar energy landscape indicating that protonation of the pyrrole is the most likely first step followed by hydride transfer and that the required proton may come from bulk solvent. The implications of the present study for the design of inhibitors of BVR-A are discussed.

Activation of biliverdin-IXalpha reductase by inorganic phosphate and related anions.

The effect of pH on the initial-rate kinetic behaviour of BVR-A (biliverdin-IXalpha reductase) exhibits an alkaline optimum with NADPH as cofactor, but a neutral optimum with NADH as cofactor. This has been described as dual cofactor and dual pH dependent behaviour; however, no mechanism has been described to explain this phenomenon. We present evidence that the apparent peak of activity observed at neutral pH with phosphate buffer and NADH as cofactor is an anion-dependent activation, where inorganic phosphate apparently mimics the role played by the 2'-phosphate of NADPH in stabilizing the interaction between NADH and the enzyme. The enzymes from mouse, rat and human all exhibit this behaviour. This behaviour is not seen with BVR-A from Xenopus tropicalis or the ancient cyanobacterial enzyme from Synechocystis PCC 6803, which, in addition to being refractory to activation by inorganic phosphate, are also differentiated by an acid pH optimum with both nicotinamide nucleotides.