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

Metabolic Functions of Biliverdin IXβ Reductase in Redox-Regulated Hematopoietic Cell Fate

Cytoprotective heme oxygenases derivatize heme to generate carbon monoxide, ferrous iron, and isomeric biliverdins, followed by rapid NAD(P)H-dependent biliverdin reduction to the antioxidant bilirubin. Recent studies have implicated biliverdin IXβ reductase (BLVRB) in a redox-regulated mechanism of hematopoietic lineage fate restricted to megakaryocyte and erythroid development, a function distinct and non-overlapping from the BLVRA (biliverdin IXα reductase) homologue. In this review, we focus on recent progress in BLVRB biochemistry and genetics, highlighting human, murine, and cell-based studies that position BLVRB-regulated redox function (or ROS accumulation) as a developmentally tuned trigger that governs megakaryocyte/erythroid lineage fate arising from hematopoietic stem cells. BLVRB crystallographic and thermodynamic studies have elucidated critical determinants of substrate utilization, redox coupling and cytoprotection, and have established that inhibitors and substrates bind within the single-Rossmann fold. These advances provide unique opportunities for the development of BLVRB-selective redox inhibitors as novel cellular targets that retain potential for therapeutic applicability in hematopoietic (and other) disorders.

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.

Repositioning Food and Drug Administration-Approved Drugs for Inhibiting Biliverdin IXβ Reductase B as a Novel Thrombocytopenia Therapeutic Target

Biliverdin IXβ reductase B (BLVRB) has recently been proposed as a novel therapeutic target for thrombocytopenia through its reactive oxygen species (ROS)-associated mechanism. Thus, we aim at repurposing drugs as new inhibitors of BLVRB. Based on IC₅₀ (<5 μM), we have identified 20 compounds out of 1496 compounds from the Food and Drug Administration (FDA)-approved library and have clearly mapped their binding sites to the active site. Furthermore, we show the detailed BLVRB-binding modes and thermodynamic properties (ΔH, ΔS, and K_D) with nuclear magnetic resonance (NMR) and isothermal titration calorimetry together with complex structures of eight water-soluble compounds. We anticipate that the results will serve as a novel platform for further in-depth studies on BLVRB effects for related functions such as ROS accumulation and megakaryocyte differentiation, and ultimately treatments of platelet disorders.

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.

Biochemical characterization of biliverdins IXβ/δ generated by a selective heme oxygenase

The pro-oxidant effect of free heme (Fe2+-protoporphyrin IX) is neutralized by phylogenetically-conserved heme oxygenases (HMOX) that generate carbon monoxide, free ferrous iron, and biliverdin (BV) tetrapyrrole(s), with downstream BV reduction by non-redundant NADPH-dependent BV reductases (BLVRA and BLVRB) that retain isomer-restricted functional activity for bilirubin (BR) generation. Regioselectivity for the heme α-meso carbon resulting in predominant BV IXα generation is a defining characteristic of canonical HMOXs, thereby limiting generation and availability of BVs IXβ, IXδ, and IXγ as BLVRB substrates. We have now exploited the unique capacity of the Pseudomonas aeruginosa (P. aeruginosa) hemO/pigA gene for focused generation of isomeric BVs (IXβ and IXδ). A scalable system followed by isomeric separation yielded highly pure samples with predicted hydrogen-bonded structure(s) as documented by 1H NMR spectroscopy. Detailed kinetic studies established near-identical activity of BV IXβ and BV IXδ as BLVRB-selective substrates, with confirmation of an ordered sequential mechanism of BR/NADP+ dissociation. Halogenated xanthene-based compounds previously identified as BLVRB-targeted flavin reductase inhibitors displayed comparable inhibition parameters using BV IXβ as substrate, documenting common structural features of the cofactor/substrate-binding pocket. These data provide further insights into structure/activity mechanisms of isomeric BVs as BLVRB substrates, with potential applicability to further dissect redox-regulated functions in cytoprotection and hematopoiesis.

Heme-Derived Metabolic Signals Dictate Immune Responses

Heme is one of the most abundant molecules in the body acting as the functional core of hemoglobin/myoglobin involved in the O₂/CO₂ carrying in the blood and tissues, redox enzymes and cytochromes in mitochondria. However, free heme is toxic and therefore its removal is a significant priority for the host. Heme is a well-established danger-associated molecular pattern (DAMP), which binds to toll-like receptor 4 (TLR4) to induce immune responses. Heme-derived metabolites including the bile pigments, biliverdin (BV) and bilirubin (BR), were first identified as toxic drivers of neonatal jaundice in 1800 but have only recently been appreciated as endogenous drivers of multiple signaling pathways involved in protection from oxidative stress and regulators of immune responses. The tissue concentration of heme, BV and BR is tightly controlled. Heme oxygenase-1 (HO-1, encoded by HMOX1) produces BV by heme degradation, while biliverdin reductase-A (BLVR-A) generates BR by the subsequent conversion of BV. BLVR-A is a fascinating protein that possesses a classical protein kinase domain, which is activated in response to BV binding to its enzymatic site and initiates the downstream mitogen-activated protein kinases (MAPK) and phosphatidylinositol 3-kinase (PI3K) pathways. This links BLVR-A activity to cell growth and survival pathways. BLVR-A also contains a bZip DNA binding domain and a nuclear export sequence (NES) and acts as a transcription factor to regulate the expression of immune modulatory genes. Here we will discuss the role of heme-related immune response and the potential for targeting the heme system for therapies directed toward hepatitis and cancer.

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.

Heme degradation enzyme biliverdin IXβ reductase is required for stem cell glutamine metabolism

Bioenergetic requirements of hematopoietic stem cells and pluripotent stem cells (PSCs) vary with lineage fate, and cellular adaptations rely largely on substrate (glucose/glutamine) availability and mitochondrial function to balance tricarboxylic acid (TCA)-derived anabolic and redox-regulated antioxidant functions. Heme synthesis and degradation converge in a linear pathway that utilizes TCA cycle-derived carbon in cataplerotic reactions of tetrapyrrole biosynthesis, terminated by NAD(P)H-dependent biliverdin reductases (IXα, BLVRA and IXβ, BLVRB) that lead to bilirubin generation and cellular antioxidant functions. We now demonstrate that PSCs with targeted deletion of BLVRB display physiologically defective antioxidant activity and cellular viability, associated with a glutamine-restricted defect in TCA entry that was computationally predicted using gene/metabolite topological network analysis and subsequently validated by bioenergetic and isotopomeric studies. Defective BLVRB-regulated glutamine utilization was accompanied by exaggerated glycolytic accumulation of the rate-limiting hexokinase reaction product glucose-6-phosphate. BLVRB-deficient embryoid body formation (a critical size parameter of early lineage fate potential) demonstrated enhanced sensitivity to the pentose phosphate pathway (PPP) inhibitor 6-aminonicotinamide with no differences in the glycolytic pathway inhibitor 2-deoxyglucose. These collective data place heme catabolism in a crucial pathway of glutamine-regulated bioenergetic metabolism and suggest that early stages of lineage fate potential require glutamine anaplerotic functions and an intact PPP, which are, in part, regulated by BLVRB activity. In principle, BLVRB inhibition represents an alternative strategy for modulating cellular glutamine utilization with consequences for cancer and hematopoietic metabolism.

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.