Identifying structural and dynamic changes during the Biliverdin Reductase B catalytic cycle

Biliverdin Reductase B (BLVRB) is an NADPH-dependent reductase that catalyzes the reduction of multiple substrates and is therefore considered a critical cellular redox regulator. In this study, we sought to address whether both structural and dynamics changes occur between different intermediates of the catalytic cycle and whether these were relegated to just the active site or the entirety of the enzyme. Through X-ray crystallography, we determined the apo BLVRB structure for the first time, revealing subtle global changes compared to the holo structure and identifying the loss of a critical hydrogen bond that “clamps” the R78-loop over the coenzyme. Amide and Cα chemical shift perturbations were used to identify environmental and secondary structural changes between intermediates, with more distant global changes observed upon coenzyme binding compared to substrate interactions. NMR relaxation rate measurements provided insights into the dynamic behavior of BLVRB during the catalytic cycle. Specifically, the inherently dynamic R78-loop that becomes ordered upon coenzyme binding persists through the catalytic cycle while similar regions experience dynamic exchange. However, the dynamic exchange processes were found to differ through the catalytic cycle with several groups of residues exhibiting similar dynamic responses. Finally, both local and distal structural and dynamic changes occur within BLVRB that are dependent solely on the oxidative state of the coenzyme. Thus, through a comprehensive analysis here, this study revealed structural and dynamic alterations in BLVRB through its catalytic cycle that are not simply relegated to the active site, but instead, are allosterically coupled throughout the enzyme.


Introduction
The development of solution strategies that probe macromolecular motions has revealed intimate relationships between structure and dynamics that are particularly important for enzymes, as they often rely on conformational changes for catalytic function.Some of the most important findings from nuclear magnetic resonance (NMR) relaxation experiments within the last 20 years has been that active sites of enzymes have evolved to undergo motions on timescales similar to their catalytic rates and these motions are allosterically coupled to distal sites (Lisi and Loria, 2017) nearly every enzyme type, such as kinases (Iverson et al., 2020), phosphatases (Whittier et al., 2013;Cui et al., 2019), proteases (Eisenmesser et al., 2015), nucleases (Narayanan et al., 2017;Narayanan et al., 2018), and oxidoreductases that are the focus of this study (Boehr et al., 2006b).Oxidoreductases represent the largest of the six major enzyme families, which in turn comprise one of the largest classes of enzyme folds referred to as short-chain dehydrogenases/reductases (SDRs) (Kavanagh et al., 2008).While such enzymes are widely reliant on their coenzymes for the functional transfer of a hydride (NADPH/NADP + and NADH/ NAD + ), even the oxidation state of the coenzyme may dictate global folding and dynamics.For example, largescale spectral changes induced by a coenzyme's oxidative state have been identified for Rhodospirillum rubrum transhydrogenase (Quirk et al., 1999) and the oxidative state of the coenzyme dictates its stereochemistry bound to ferredoxin (Kean et al., 2017).Identifying the changes between catalytic intermediates is often difficult, as similar X-ray crystal structures can be observed despite largescale spectral changes observed via NMR that do support global changes in solution (Quirk et al., 1999;Sundaresan et al., 2003).One explanation could be that coenzymes or substrates are often soaked into preformed crystals, yet another explanation is that the ground state observed in some crystals may not reflect changes to dynamically sampled conformations in solution.Consistent with dynamically imparted differences, the anisotropic B-factors are different between crystallized forms of NADPH and NADP + bound to ferredoxin (Kean et al., 2017) and changes to conformational sampling in solution have been proposed for catalytically important active site loops within multiple oxidoreductases that include the "D loop" of Rhodospirillum rubrum transhydrogenase (Quirk et al., 1999) and the "M20 loop" of dihydrofolate reductases (DHFR) (Boehr et al., 2006a;Boehr et al., 2006b).In this report, we have monitored both the backbone changes and dynamic changes of a prototypical SDR family member, Biliverdin Reductase B (BLVRB), in order to determine whether there are global structural/dynamic changes through its catalytic cycle and whether its coenzyme's hydride alone may dictate global coupling.
BLVRB is an NADPH-dependent oxidoreductase that has emerged as a critical redox regulator within the last decade, as it catalyzes reduction of multiple substrates (Figure 1A).BLVRB was originally thought to reduce only the biliverdin-β isomer that is less abundant than the biliverdin-α isomer catalyzed by the BLVRA isoform (Cunningham et al., 2000;McDonagh, 2001).However, it was later discovered that BLVRB is identical to flavin reductase (FR) that is a promiscuous enzyme responsible for reduction of flavin substrates (Shalloe et al., 1996), which includes flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN).BLVRB plays such a critical role in redox regulation that this enzyme alone dictates hematopoietic cell fate (Wu et al., 2016).BLVRBmediated redox regulation may explain why its expression is high in cancer cells.For example, BLVRB is a driver of hepatic cancer cell proliferation (Huan et al., 2016) and is a top prognostic factor for fatal prostate cancer (Ramberg et al., 2021).We have also shown that BLVRB is expressed at much higher levels in red blood cells than previously thought where it likely also plays a role in redox regulation (Paukovich et al., 2018).Biochemically, the ratelimiting step in BLVRB activity has remained unknown since its discovery, yet our biochemical and biophysical studies suggest that the slowest step identified thus far, coenzyme release, is coupled to a rate-limiting step (Duff et al., 2020;Redzic et al., 2021) (Figure 1A).Structurally, BLVRB comprises a typical SDR fold with a central βsheet flanked α-helices on each side (Figure 1B).Like other SDR family members, the coenzyme binding site is within the conserved SDR fold itself, while substrate specificity is dictated by the variable C-terminal lobe.X-ray crystallography has revealed that BLVRB R78 forms a hydrogen bond to T12 CO that serves to straddle the coenzyme (PDB accession 1HDO), while NMR has shown that this loop, referred to as the "R78-loop", is highly dynamic in the absence of the oxidized coenzyme (Paukovich et al., 2018).Unfortunately, it has been difficult to identify changes within BLVRB that may occur through its catalytic cycle, as only holo BLVRB has been crystalized with no apparent changes upon substrate binding (Pereira et al., .Substrate binding (S, green hexagon) forms the E:NADPH: S Michaelis-Menten complex (3).Reduction of the substrate results in the product (P, green circle) and oxidized NADP + (red circle) to form the product complex (4).Release of the product then forms the spent coenzyme complex (5), which is released to restart the catalytic cycle.(B) Model of substrate bound ternary complex of BLVRB based on the X-ray crystal structure of BLVRB bound to NADP + and mesobiliverdin IV (PDB accession 1HE3).Biliverdin (green), NADP + (red) and R78 are highlighted.

2001
).Our own X-ray crystallographic studies of a human holo BLVRB mutant and even its lemur homologue that comprises 16 changes has also crystallized identically (Duff et al., 2020).Thus, here we sought to determine whether we could identify structural and dynamic changes within human BLVRB through its catalytic cycle and whether these were relegated to the active site or globally coupled sites.

Apo BLVRB crystallizes in the same conformation as holo BLVRB
In an effort to identify the conformational changes within BLVRB necessary to facilitate coenzyme binding, we have solved the first X-ray crystal structures of apo human BLVRB that has eluded crystallization until now (Figure 2A).Apo BLVRB crystallized within the same space group of P2 1 2 1 2 1 as holo BLVRB (PDB accession 1HDO) at both cryogenic (apo-cryo) and room temperature conditions (apo-RT) (Table 1).At the cofactor binding site, three arginines show clear remodeling compared to the holo state (Figure 2B).In both apo-RT and apo-cryo structures, R35 and R39 are both found in specific conformations, as indicated by clear electron density.Most surprising is that in the apo-cryo structure, the R78-loop (residues 77-80), which exhibits high R1 relaxation rates indicative of disorder in solution (Paukovich et al., 2018), is well modelled in electron density where the R78 sidechain interacts with its crystallographic mate (not shown).Although this suggests that the R78-loop is "closed" and would therefore block binding and release of the coenzyme, the R78 side chain does not form the hydrogen bond to T12 CO that forms the clamp.To determine if this defined conformation of the R78-loop is an artifact of cooling to cryogenic temperatures, the complementary room-temperature (20 °C) structure was compared.This ambient temperature structure indicates that R35 and R39 are found in the same position as the apo-cryo dataset, suggesting these conformational changes are likely important in cofactor binding.In contrast, the R78 sidechain and loop electron density is relatively absent (Figure 2B), indicating flexibility in agreement with previous NMR data (Paukovich et al., 2018) and illustrated here by B-factor comparisons (Figure 2C).Furthermore, this increased flexibility of the R78-loop would allow for coenzyme binding and release, but only as observed within the apo-RT BLVRB structure.Elsewhere, both apo structures are generally similar, but there are significant regions throughout the enzyme which differ when comparing the apo-structures and the previously determined holo-structure (0.8 Å all-atom RMSD).Specifically, loops 34-46, 115-126, and 153-176 show large displacements and conformational changes with some regions having greater than a 2 Å Cα RMSD (Figure 2D).Interestingly, these same regions are inherently dynamic on multiple timescales in solution and exhibit chemical shift perturbations (CSPs) upon coenzyme engagement (Paukovich et al., 2018), which is also further illustrated by our NMR studies below.Thus, these conformational changes observed when comparing both apo BLVRB structures to the holo-structure indicates partial remodeling of BLVRB global and active site structure upon cofactor binding, and furthermore highlight how ambient temperature experiments may better match the flexibility observed within solution NMR studies.To further identify the regions of BLVRB that change during the catalytic cycle and associated changes to motions after cofactor binding, we turned to NMR solution studies.

Secondary structure changes between intermediates reveal global changes
In order to identify changes within BLVRB through multiple stages of its catalytic cycle, we assigned the backbone of all possible stable complexes, which included both the apo BLVRB and bound NADP + complex previously published (Paukovich et al., 2018), the bound NADPH complex here, and the product ternary complex with both NADP + and bilirubin.While BLVRB catalyzes reduction of multiple flavins, their affinities are in the high micromolar, and saturating conditions are not possible.Therefore, the bilirubin product was used instead, as its affinity is within the low micromolar range and would avoid contributions from exchange for relaxation experiments described below.However, commercial preparations of bilirubin used for these studies do comprise mixtures of isoforms that all bind BLVRB (Cunningham et al., 2000).We note that the Michaelis-Menten complex could not be probed here due to the obvious fact that turnover would occur within the NMR tube (with both bound NADPH and biliverdin substrate).Thus, chemical shift assignments have now been expanded to the four possible catalytic intermediates of 1, 2, 4, and 5 (see Figure 1A).
Backbone assignments facilitated the quantitative determination of both amide and Cα CSPs (Figure 3, left).We note that as we are unable to probe the catalytically active Michaelis-Menten complex of BLVRB/NADPH/biliverdin, differences between the prior and following step were calculated instead to identify changes between these intermediates (Figure 3B).In general, while both amide and Cα CSPs are sensitive markers to chemical environments, the backbone Cα chemical shifts are largely dictated by local secondary structure.For example, Cα chemical shifts are routinely compared to their shifts for the same residue within a random peptide (Hafsa and Wishart, 2014), which are called chemical shift propensities (defined as "ΔCα").Propensities larger than that of a random peptide for the same residue indicate α-helical propensity (i.e., ΔCα>0) while Cα chemical shifts that are lower than that of a random peptide indicate β-strand propensity (ΔCα<0).For the case of BLVRB here, most of these propensities remain in the same direction (Supplementary Figure S1), suggesting that the secondary structure remains similar through the catalytic cycle.However, localized changes are observed for several residues.These include R78 that has no propensity in the apo state to largely helical propensity through the rest of the catalytic cycle.Interestingly, S112 adjacent to the coenzyme's hydride is disordered only when bound to NADPH, suggesting localized flexibility that may be important for hydride transfer.Thus, while the secondary structure remains largely the same through the catalytic cycle, localized differences are observed throughout the enzyme that are further described below.
We mapped both the amide and Cα CSPs between intermediates to quantitatively identify regions exhibiting changes (Figure 3, middle and right).Several observations are immediately apparent for the CSPs monitored here between catalytic cycle intermediates.First, the same regions that exhibit amide CSPs overlap well with the regions that exhibit Cα CSPs.This indicates that changes in chemical environment through the catalytic cycle are likely concomitant with subtle changes to local secondary structure.Second, NADPH binding (Figure 3A) and NADP + release (Figure 3D) induces the largest overall CSPs in both magnitude and number of residues.This is contrast to smaller changes in both magnitude and the number of residues induced by the substrate between intermediates (Figures 3B, C) and can be observed quantitatively by the reduced average changes for the substrate when compared to the coenzyme (Figures 3B,C versus Figures 3A,D).Third, many of the CSPs through the catalytic cycle are not immediately in contact with the coenzyme or bilirubin product.This indicates that long range allosteric changes are induced, which is consistent with our recent study that indicated mutation of BLVRB T164 within the C-terminal lobe induces measurable changes to the active site (Redzic et al., 2021).Finally, one of the most interesting findings here may be that the oxidation state of the coenzyme alone induces global CSPs (Figure 3E).Specifically, although the enzyme does not proceed directly between NADPH to NADP + bound forms, a comparison between these two intermediates illustrates that the oxidative state of the coenzyme alone induces changes (Figure 3E).Interestingly, these CSPs are quantitatively on a similar level to the CSPs induced by bilirubin (Figures 3B, C).Moreover, an allosteric network of Cα CSPs includes residues 153-155 and residues 112-116 that reside in the substrate binding pocket.The fact that a single hydride and thus, the oxidative state of the nicotinamide group, can induce such largescale changes may have functional consequences for substrate binding.For example, subtle differences in coenzyme puckering or the dynamics of the coenzyme itself may differ through the catalytic cycle (Sundaresan et al., 2003;Tsybovsky et al., 2013).For BLVRB, many active site residues that straddle the entire coenzyme exhibit Cα CSPs between the oxidation states (Figure 3E), suggesting that subtle structural changes of to the coenzyme's stereochemistry are indeed by coenzyme oxidation.

FIGURE 3
Amide and Cα CSPs between BLVRB catalytic intermediates.The absolute differences for amide CSPs (navy blue) and Cα CSPs (magenta) are shown (left).The same cutoffs were used to map CSPs in order to compare the number of residues affected.Specifically, while the average CSP plus one standard deviation is listed below quantified from these data, a hard cutoff of the lowest value was used for amide CSPs (center, navy blue balls that are larger than 0.78 ppm) and for Cα CSPs (right, magenta balls that are larger than 0.24 ppm).The average CSP plus one standard deviation is show for S111A is 0.45 ppm (dashed line) (E) Amide CSPs for NADPH/NADP + of WT reported in Figure 3E (blue) compared to H153A (green).The average CSP plus one standard deviation is show for H153A is 0.45 ppm (dashed line) (F) Amide CSPs for NADPH/ NADP + of WT reported in Figure 3E (blue) compared to T164S (magenta).The average CSP plus one standard deviation is show for T164S is 0.51 ppm (Continued ) Even more distal regions are also affected by the oxidative state that are simply below the threshold used here for mapping changes, such as residues within the C-terminal lobe that include residues 163-165.Thus, both local and allosteric changes can be monitored through the catalytic cycle with the coenzyme inducing the largest changes but the oxidative state of the coenzyme still modulating both local and allosteric changes.
Mutations that "short-circuit" global hydride coupling alter catalytic function Considering that the oxidative state of the coenzyme alone, and thus a single hydride, modulates changes both within and distant to the active site (Figure 3E), we sought to further explore networks of allostery modulated by the oxidation state of the coenzyme.To further probe such coupling, we screened several BLVRB mutations both within and distal to the BLVRB active site using CSPs.Specifically, we selected three residues that incurred large CSPs between their bound reduced and oxidized forms of the coenzyme, which included S111 (Figure 4A), H153 (Figure 4B), and T164 (Figure 4C).These sites were also chosen as both S111 and H153 are immediately adjacent to the coenzyme's hydride, while previous studies have indicated that T164 is allosterically coupled to the active site (Redzic et al., 2021).As mutation of T164 to a serine found in other mammalian homologues induces larger functional effects that a random mutation to an alanine (Redzic et al., 2021), we chose to explore T164S.Thus, S111A, H153A, and T164S were used to identify the effects to global hydride-coupled networks.
Amide CSPs between NADPH and NADP + bound forms of BLVRB wild type (WT) calculated from 15 N-HSQC spectra were compared to the respective differences for S111A (Figure 4D), H153A (Figure 4E), and T164S (Figure 4F).We note that the amides and several residues immediately adjacent to the mutated active site of S111A and H153A were not identified either due to overlap or severe line-broadening.Some local changes are observed to differ, such as H153A induced changes to neighboring residues C109.However, the most obvious impact of all three mutations was that hydride-mediated changes were similarly diminished for residues within the substrate binding loop of residues 114-121 (Figure 4G).In other words, whether active site residues are mutated, or a residue shown to be involved in allosteric coupling is mutated, their impacts on hydridecoupled networks are similar.We have previously referred to such diminished effects on either CSPs or the dynamics of networks of coupled residues as "short-circuited" (Holliday et al., 2017).Interestingly though, the catalytic effects of substrate turnover are different between these mutations using the FAD substrate (Figure 4H).FAD was used as no commercial sources of the biliverdin-beta isomer are available.Specifically, we compared turnover of these mutations to the BLVRB WT using the standard UV-kinetics assay under saturating conditions of NADPH while changing FAD (Cunningham et al., 2000;Duff et al., 2020).While we have previously presented similar catalytic data (Redzic et al., 2021), here we performed this assay differently in order to sample higher concentrations of FAD (described in Materials and Methods).Reversal of this UV-kinetics assay to vary NADPH is not attainable considering the high concentrations needed to saturate FAD that saturates the UV detector (Paukovich et al., 2018).However, using this assay while varying FAD, all three mutations induced changes to turnover.Both S111A and T164S mutations lead to reductions in maximum velocity and weaker affinities (higher K M values) while H153A has no effect on maximum velocity but increase affinity (lower K M value).Considering that the global CSPs shortcircuited by these three mutants are largely similar, it is interesting to speculate why their imparted catalytic changes differ.Although such differences could be due to an array of perturbations during the catalytic cycle, catalytic differences could also be influenced by the way the oxidized and reduced forms of the coenzyme may lie within the active site of each mutant.In general though, these catalytic differences relative to the WT enzyme indicate that disturbing residues involved in hydride-coupled networks broadly affect catalytic function.
Standard relaxation rates reveal that the dynamic regions are largely similar through the catalytic cycle Enzymes are not static but instead, are often inherently dynamic, which may be particularly important to allow for the conformational changes that ensue during catalytic function.Standard relaxation methods were therefore applied to BLVRB catalytic intermediates to determine how dynamics may change during the catalytic cycle.These include R1 relaxation rates that are sensitive to the ps-ns timescales and R2/R1ρ relaxation rates that are sensitive to the slower μs-ms timescales, which are both extracted from their time dependencies with examples and their fits shown here (Supplementary Figure S2).
R1 relaxation data indicate that the disordered R78-loop in apo BLVRB remains largely ordered upon coenzyme binding through the catalytic cycle (Figure 5A; Supplementary Figure S3).This can be structurally explained by the likely preservation of the hydrogen bond between the R78 side chain and R12 CO observed within the X-ray crystal structure of BLVRB bound to NADP + (Pereira et al., 2001).It should be noted that R1 relaxation rates are strictly a weighted average and thus, FIGURE 4 (Continued) (dashed line) (G) Two views of the positions of each mutation within the holo BLVRB structure and the region of residues 114-121 (yellow) that exhibits similar imparted changes from all three mutations.Mutations are colored the same as D-F and the C4N position that is reduced is shown (yellow sphere).(H) UV-kinetics turnover of FAD is shown for the WT and all three mutants.Specific Vmax values with KM values in parenthesis of each are as follows: WT is 3.3 ± 0.5 μM*s -1 (242 ± 70 μM), H153A is 2.9 ± 0.2 μM*s -1 (129 ± 40 μM), S111A is 2.5 ± 0.10 μM*s -1 (856 ± 300 μM), and T164S is 2.1 ± 0.01 μM*s -1 (171 ± 45 μM).
these findings do not exclude the likely possibility that this loop must further open to facilitate NADPH binding and NADP + release.However, such R1 relaxation rates do imply that once the R78-loop is closed upon NADPH binding that it stays relatively ordered through the catalytic cycle.
Chemical exchange that contributes to R2 relaxation in the μsms timescale are suppressed by the imparted refocusing field within R1ρ relaxation and therefore identify slower motions within BLVRB that persist through the catalytic cycle (Supplementary Figure S4).Specifically, amides that exhibit differences between R2 and R1ρ greater than their respective uncertainties indicate that many of the same residues exhibit exchange that include residues within the R78loop (Figure 5B) and within residues 164-177 (Supplementary Figure S4).This means that many of the same regions of BLVRB are moving within the μs-ms timescale during multiple steps of the catalytic cycle.However, the fact that such motions persist within multiple intermediates within the catalytic cycle does not necessarily imply that their motions are equivalent.Thus, we subjected BLVRB catalytic intermediates to R2-CPMG dispersion experiments and did so at multiple temperatures, as described below.

R2-CPMG dispersion experiments quantify differential dynamics during the BLVRB catalytic cycle
Based on R2/R1ρ relaxation rate comparisons described above that revealed similar regions exhibit μs-ms timescale motions within different catalytic intermediates, we sought to understand whether the underlying motions differ between intermediates.However, our previously published R2-CPMG studies that focused solely on comparisons between the apo and NADP + bound forms at 20 °C revealed several challenges that required alternative strategies (Paukovich et al., 2018).First, exchange contributions with NADP + bound were very small at 20 °C.Such small contributions preclude data collection at multiple static fields that are required to extract the associated biophysical parameters of motions (Kovrigin et al., 2006).Second, our recent studies indicate that much of the exchange dynamics within BLVRB are highly localized but nonetheless partially coupled throughout the enzyme (Redzic et al., 2021).This means that R2-CPMG fits to a two-site exchange are in general confounded by contributions from FIGURE 5 R1, R2, and R1ρ relaxation rates monitored for BLVRB catalytic intermediates.(A) Top: R1 relaxation rates of apo BLVRB (black), BLVRB bound to NADPH (blue), BLVRB bound to NADP + and BR (green), and BLVRB bound to NADP + (red).R1 relaxation rates for residues 60-100 are shown.Bottom: R1 relaxation rates within apo BLVRB greater than one standard deviation mapped onto a model of the ternary complex, which is quenched for all other intermediates.(B) R2 and R1ρ relaxation rates are shown for BLVRB catalytic intermediates with significant differences larger than their uncertainties mapped onto the model structure of the ternary complex for each intermediate (cyan spheres).R2 and R1ρ relaxation rates for residues 60-100 are shown.All data were collected at 20 °C at 900 MHz with all data shown in Supplementary Figure S3, 4.
numerous processes and are therefore an over-simplification.To partially address these challenges, R2-CPMG dispersions were collected at lower temperatures of 0 o C, 5 o C, and 10 o C at 900 MHz where exchange was larger in some intermediates.To estimate the underlying physical parameters that define μs-ms timescale motions within these BLVRB intermediates, amides that exhibited similar responses within each intermediate were simultaneously grouped and fit at multiple temperatures using the Arrhenius dependence (Millet et al., 2000;Vallurupalli and Kay, 2006).(B) Group II residues are mapped onto the ternary complex (magenta spheres) that include L74, T77, C109, H111, and Q126 with L74 and S111 amide R2-CPMG dispersions shown at 10 °C.(C) Group III residues that include I154, G155, T164, I176, Y205 are mapped onto the ternary complex with I154 and Y205 amide dispersions shown at 10 °C.All data were collected at 900 MHz and all R2-CPMG dispersions are shown with their fits in Supplementary Figure S5.
R2-CPMG dispersions collected here reveal that their profiles change for many sites through the catalytic cycle that indicates the associated dynamics also differ (Supplementary Figure S6).Three groups exhibiting similar changes to their R2-CPMG dispersions were identified, which are referred to herein as "Group I", "Group II", and "Group III" (Figure 6).These three groups were then separately

Apo-RT BLVRB (293K)
Apo-cryo BLVRB (cryogenic) fit using Graphical User-friendly Analysis of Relaxation Dispersion Data (GUARDD) (Kleckner and Foster, 2012).The extracted rates of exchange and populations are tabulated here within the three groups at the three temperatures (Table 2).Specifically, amides within Group I exhibit no exchange contributions within apo BLVRB but exhibit similar R2-CPMG dispersions once coenzyme and product are bound (Figure 6A).Several Group I residues make direct contact to the adenine moiety of the adenosine.The opposite is true for amides within Group II (Figure 6B), which exhibit large exchange in apo BLVRB but have this exchange largely quenched within the other intermediates with little to no R2-CPMG dispersion.However, some of these residues do still exhibit small exchange contributions that increase with lower temperatures with both the reduced and oxidized coenzyme that could be fit to relatively fast exchange rates (Table 2).While this indicates that at least some movements within Group II simply shift to faster rates of exchange (and thus, smaller exchange contributions), there is no observable exchange in the presence of the BR product at multiple temperatures (Supplementary Figure S5).Group II residues are scattered closer to the nicotinamide and substrate binding site.Finally, Group III residues exhibit detectable exchange throughout the catalytic cycle with subtle differences between each intermediate (Figure 5C).R2-CPMG dispersion is even different between the oxidized and reduced form of the coenzyme, suggesting that, like CSPs, the oxidative state of the coenzyme is allosterically coupled to other regions of BLVRB.
Interestingly, the associated rates of motions of Group III exhibit the least relative changes as a function of temperature (Table 2).These residues are located within the C-terminal lobe of BLVRB, which also includes residues within 166-175 that remain unobservable (except for G169) within all intermediates.Thus, although many similar regions are moving within multiple intermediates, the specific motions themselves do differ through the catalytic cycle.
To help understand the physical basis of the dynamically sampled conformations monitored by R2-CPMG dispersions, we compared the calculated chemical shift changes from the fits (Δω cpmg ) to the experimental CSPs measured between all possible BLVRB intermediates (Supplementary Figure S6).Specifically, if sampled populations are identical to either the previous conformation within the catalytic cycle (backward sampling) or the next sampled conformer (forward sampling), there should be a correlation between the calculated and experimental chemical shifts.While such comparisons are often complicated by the fact that the sampled chemical environment may differ in the presence of the bound species, we do observe at least several weak correlations for some intermediates.Specifically, several amides within apo BLVRB appear to dynamically sample similar CSPs to their holo forms (Supplementary Figure S6, top).Interestingly, only I176 and Y205 that are both within the C-terminal lobe of BLVRB do not sample conformations with similar CSPs.Beyond apo BLVRB, about half the measured dynamics within NADP + -bound BLVRB dynamically sample similar CSPs to the apo form (Supplementary Figure S6, bottom).Our inability to monitor the Michaelis-Menten complex prevents us from specifically comparing many of the other dynamically sampled conformations (Supplementary Figure S6, middle).Nonetheless, the remaining comparisons do not indicate that dynamics within other intermediates sample similar conformations to those defined by their CSPs.Thus, apo and holo BLVRBs partially sample each other's conformations, which is supported by the apo BLVRB structure solved here that is similar to that of holo BLVRB.

Structural changes between BLVRB catalytic intermediates
This study identified regions of BLVRB that incur conformational changes during different steps of its catalytic cycle.We began by successfully solving X-ray crystal structures of the apo form of BLVRB for the first time at both cryogenic and TABLE 2 Exchange rates and populations determined for three separate groups.Residues within each group were globally fit using 0 °C, 5 °C, and 10 °C data except for apo BLVRB that also included 20 °C.Data were fit using GUARDD (Kleckner and Foster, 2012).
ambient temperatures, which had previously eluded crystallization attempts.The comparison between the holo BLVRB structure and the apo BLVRB structures revealed that several active site residues undergo conformational changes to accommodate the coenzyme, including R35, R39, and R78.One important difference between the apo BLVRB structures was that the apo-cryo structure exhibited a defined density of the R78-loop that would block coenzyme entry and release, while this region was disordered in the apo-RT structure that would facilitate such an interaction.Such comparisons highlight the potential misinterpretation of structures determined by cryogenic crystallography without complementary ambient temperature experiments.Importantly, even in the apo-RT BLVRB structure, the hydrogen bond between R78 and the T12 CO observed within the NADP + -bound holo BLVRB structure was not formed.This is also consistent with the high flexibility of the R78-loop in the apo state (Paukovich et al., 2018).Additional local differences within the crystals between the apo and holo BLVRB structures were observed distant to the active site, illustrating allosteric rearrangements throughout the enzyme.
NMR solution studies that focused on amide and Cα CSPs revealed environmental and secondary structural changes between BLVRB catalytic intermediates both within the active site and distal regions, consistent with allosterically coupled changes.The largest changes were found to occur between BLVRB intermediates upon coenzyme binding, whether reduced or oxidized, compared to more localized changes due to product interactions.This suggests that global rearrangements are induced to a larger degree by coenzyme interactions relative to the product.Interestingly, the oxidative state of the coenzyme, meaning the differences between NADPH-bound and NADP + -bound, exhibited CSPs on a similar scale to those induced by the product.While such changes within the active site could reflect slight differences in the stereochemistry of the coenzyme that are altered by the oxidation state, as has been observed before (Sundaresan et al., 2003;Tsybovsky et al., 2013), distal residues are nonetheless affected.Thus, the oxidative state of the coenzyme induces allosterically coupled changes.Mutations that alter this coupling to the oxidative state, i.e., that alter hydridemediated coupling, also exhibit catalytic differences.We have therefore shown that allosteric coupling to the oxidation state of a reductase's coenzyme is important for function.

Dynamic changes between BLVRB catalytic intermediates
We utilized NMR relaxation rate measurements to reveal the changes that occur within BLVRB during different steps of the catalytic cycle.For the ps-ns timescale, the inherently dynamic R78loop in the apo state remains largely ordered for the other catalytic intermediates.This suggests that once the coenzyme is in place, the R78-loop is stabilized by the hydrogen bond between the side chain of R78 and the T12 CO observed within the NADP + -bound X-ray crystal structure (Pereira et al., 2001).For slower μs-ms timescales, comparisons of R2 relaxation with R1ρ relaxation were initially used to identify sites of dynamic exchange within BLVRB catalytic intermediates.These comparisons revealed that exchange contributions were quenched within similar regions of BLVRB through the catalytic cycle, indicating that the associated dynamic exchange on the μs-ms timescale is relegated to similar regions.By applying R2-CPMG dispersion experiments to further quantify the dynamic exchange within BLVRB at multiple temperatures, we identified groups of residues that exhibit similar patterns, or responses, through the catalytic cycle that were structurally localized to different regions.Thus, despite similar regions moving on the μs-ms timescale within different intermediates, the dynamic processes themselves differ through the catalytic cycle.
Using R2-CPMG dispersions, we identified three groups of residues that exhibit similar μs-ms motions within BLVRB intermediates, which included both the active site as well as allosterically coupled residues.Two of these groups were within the active site (Group I and II) and appeared to undergo a dynamic "switch".Specifically, Group I exhibited no R2-CPMG dispersion in apo BLVRB but did exhibit dispersion throughout the rest of the catalytic cycle and vice-versa for Group II that was inherently dynamic within apo BLVRB but largely quenched with coenzyme and substrate.In contrast, the final group (Group III) was both proximal to the active site but extended well within the C-terminal lobe.The identification of coupled movements within this C-terminal lobe are consistent with our previous studies that have identified one particular site, T164, which modulates coenzyme affinity allosterically (Redzic et al., 2021).Specifically, the movements within this C-terminal region are drastically altered by mutation of T164 and this directly impacts coenzyme affinity.Just as mutation of T164 has previously been shown to alter active site dynamics, here we showed that even the oxidative state of the coenzyme affects the dynamics within these allosterically coupled sites.Specifically, R2-CPMG dispersions are not identical between bound NADPH and NADP + for these residues within Group III.Interestingly, the rate of exchange within this group extracted from group fitting of R2-CPMG dispersions here also revealed the least change with temperature.This observation that BLVRB retains the underlying movements of this C-terminal lobe to similar rates of motions can be compared to other enzymes that have been observed to limit temperature dependencies within catalytically important regions (Schlegel et al., 2009a).For example, a limited temperature dependence may be an evolutionary feature to maintain catalytic efficiencies.The importance of this C-terminal region comprising Group III residues is also consistent with comparisons to other SDR family members.Specifically, the C-terminal lobe of BLVRB is the distinguishing feature among SDR family members that dictates substrate interactions (Kavanagh et al., 2008).
Finally, it is important to note the dynamic differences of BLVRB to another well-known reductase, DHFR.Unlike concerted motions identified for DHFR during its catalytic cycle that enable sampling of either the previous or following catalytic step (Boehr et al., 2006a;Boehr et al., 2006b), we observed only weak evidence of such concerted movements within BLVRB.Specifically, only in the case of either apo BLVRB or holo BLVRB bound to NADP + did we observed a subset of CSPs that were similar to the calculated chemical shift differences derived from R2-CPMG dispersion.While this is at least consistent with our X-ray crystal structures of apo BLVRB that is indeed similar to bound NADP + , the dynamics within the remaining BLVRB intermediates appear not to sample previous or following conformations.Moreover, unlike DHFR where dynamics have been associated with the ratelimiting step of product release (Boehr et al., 2008), BLVRB turnover is much slower than any of the independent steps that include hydride transfer (Duff et al., 2020).This suggests that the dynamic relationship of BLVRB to function may be a collection of dynamic events that influence catalysis rather than a single movement dictating catalysis.Nonetheless, this study revealed dynamic differences within each intermediate through the catalytic cycle within both the active site and allosterically coupled sites.

Protein expression and purification
Expression and purification of recombinant human BLVRB (UniProt P30043) has previously been described (Paukovich et al., 2018).Briefly, BLVRB cloned after an N-terminal 6xHis-tag and thrombin cleavage site was expressed in pET21 was expressed in BL21/DE3 cells at 37 o C and induced at 0.6 ODs (600 nM) with isopropyl b-D-1-thiogalactopyranoside (IPTG).BLVRB mutations were produced by a nested PCR protocol whereby the initial PCR produced a partial fragment with the mutation that utilized this fragment as a primer for the full-length PCR fragment with insertion into NdeI-cleaved pET22.Hi-Fi PCR amplification (NEB) was used for all PCR reactions and In-Fusion (Takara) was used for vector insertion.Unlabeled proteins were expressed in luria broth (LB) for 3 h and labeled proteins were expressed in M9 minimal media for 5 h 2 H, 15 N, 13 C-labeling for backbone assignments was grown in 99.9% D 2 O M9 minimal media supplemented with 15 N ammonium chloride and 13 C glucose facilitated further assignments of residues initially too line broadened. 2H, 15 N-labeling for R2-CPMG dispersion experiments was grown in 99.9% D 2 O M9 minimal media supplemented with 15 N ammonium chloride.BLVRB proteins were all purified denatured via Ni-affinity (Sigma), refolded through dialysis, and further purified via Superose-75 size-exclusion-chromatography (Cytiva) as previously described (Paukovich et al., 2018;Duff et al., 2020;Redzic et al., 2021).
Human glucose-6-phosphate dehydrogenase (G6PD, UniProt P11413) that was used to recycle NADPH, described further below, was also cloned into pET21 for expression at 25 o C and induced at 0.6 ODs (600 nM) with IPTG for 5 h.Soluble protein was purified via Ni-affinity and fractions concentrated for further purification of the tetramer via Superose-200 size-exclusion-chromatography (Cytiva).

Crystallization and structure solution
Crystals were grown by vapor diffusion where protein/mother liquor drops were equilibrated over a 0.25-0.5 M NaCl solution at 25 °C.3 μL of human BLVRB at 18 mg/mL were mixed with 3 µL of crystallization solution consisting of 0.1M MES pH 6.5% and 1%-4% PEG 20000 (Duff et al., 2020).Initial crystals were harvested for seeds of which 0.5 µL were added to fresh drops after 1 day of equilibration.Crystals were cryoprotected in the crystallization solution with the addition of 10% PEG 400.Crystals were either immersed into liquid nitrogen for cryogenic data collection or wrapped in NVH oil for ambient temperature collection.Data was collected at CHESS ID7B2 beamline using CRLs.Data was index and scaled in DIALS (Winter et al., 2018) and merged using AIMLESS (Evans and Murshudov, 2013).BLVRB structures were solved using MOLREP in the CCP4 suite with PDB accession 1HDO as a search model (Bailey, 1994;Pereira et al., 2001;Vagin and Teplyakov, 2010).Real-space refinement, model building, heteroatom and water addition was completed in COOT (Emsley and Cowtan, 2004).Refinement was completed using Phenix.refine using automatic TLS restraints (Afonine et al., 2012).Validation was completed using Molprobity server (http://molprobity.biochem.duke.edu/index.php).

NMR sample preparation, spectroscopy, and data analysis
All BLVRB samples were prepared in 50 mM bis-tris, pH 6.5, 50 mM NaCl, 1 mM DTT at 500 μM protein with 5% D 2 O. Holo BLVRB with NADP + (oxidized coenzyme) comprised 2 mM coenzyme while samples with the coenzyme product also comprised 700 μM bilirubin.For holo BLVRB with NADPH, the recycling system previously exploited was used to extend the life-time of the sample due to the short half-life of the reduced form of this coenzyme (Boehr et al., 2006b).Specifically, 2 mM NADPH was used with 10 μM G6PD and 10 mM G6P.No interaction was observed for G6P with BLVRB, as assessed by 15 N-HSQC spectra.Backbone assignments using a 2 H, 15 N, 13 C-labeled enzymes were collected on a Bruker 600 spectrometer equipped with a cryo-probe at 20 °C for all intermediates.These included an HNCA, HN(co)CA, HNCACB, and a HN(co)CACB.As the Cα resonances are readily obtained from both the HNCA and HN(co)CA, both the HNCACB and CBCA(co)NH were optimized for CB coupling, specifically by setting the delay for CC coupling to 6 m.All 3D assignment spectra utilized non-uniform sampling (NUS) reconstructed to 72 and 96 increments in the nitrogen and carbon dimensions, respectively (Hyberts et al., 2012).Backbone assignments for BLVRB in the presence of NADPH and in the presence of NADP + and bilirubin have been deposited into the Biological Magnetic Resonance Data Bank (BMRB) with accession numbers 51658 and 51657, respectively, and BLVRB alone and in the presence of NADP + have previously been deposited (Paukovich et al., 2018).
All relaxation experiments were conducted on the Rocky Mountain Varian 900 spectrometer equipped with a cryo-probe at 20 o C with R1, R2, and R1ρ relaxation experiments utilizing a standard BioPack sequence and R2-CPMG dispersions were collected at 0 o C, 5 o C, and 10 o C with TROSY selection as previously described (Schlegel et al., 2009b;Holliday et al., 2017).Standard R1, R2, and R1ρ relaxation rates were determined from fits in CCPNmr (Vranken et al., 2005) and R2-CPMG dispersions were fit to extract the associated exchange rates and populations using GUARDD (Kleckner and Foster, 2012).
Cα CSPs atoms were calculated by taking absolute differences between each intermediate while amide CSPs were calculated by the square root of the sum of squares for both nitrogen and protein shifts between intermediates.For hydride-mediated coupling experiments, BLVRB WT, S111A, H153A, and T164S were compared between reduced and oxidized coenzyme.Specifically, either 2 mM NADP + or NADPH with the above-described recycling system were used and amide CSPs calculated accordingly.

UV-kinetic assays
A Varian Cary 50 UV-Vis Spectrophotomer was used with a 1 mm cell and 300 μL total volumes to observe the diminishment in NADPH at 340 nm, as previously described (Cunningham et al., 2000).NADPH concentrations were kept constant at 300 μM, BLVRB and mutants were kept constant at 10 μM, and the FAD substrate was varied as indicated.Each velocity was determined from an average of three kinetic runs with initial velocities calculated for the first 60 s and fit to the Michaelis-Menton equation using GraphPad Prism version 4.0 (GraphPad Software Inc., San Diego, CA).

FIGURE 1
FIGURE 1 BLVRB catalytic cycle and ternary complex.(A)The BLVRB catalytic cycle.Apo BLVRB (1) binds NADPH (blue hexagon) to form E: NADPH (2).Substrate binding (S, green hexagon) forms the E:NADPH: S Michaelis-Menten complex (3).Reduction of the substrate results in the product (P, green circle) and oxidized NADP + (red circle) to form the product complex (4).Release of the product then forms the spent coenzyme complex (5), which is released to restart the catalytic cycle.(B) Model of substrate bound ternary complex of BLVRB based on the X-ray crystal structure of BLVRB bound to NADP + and mesobiliverdin IV (PDB accession 1HE3).Biliverdin (green), NADP + (red) and R78 are highlighted.
FIGURE 2 X-ray crystal structures of apo BLVRB and comparison to holo BLVRB.(A) Comparison of holo BLVRB bound to NADP + (PDB accession 1HDO, white) to apo-cryo BLVRB (blue) and apo-RT BLVRB (cyan).The bound NADP + in holo BLVRB is shown for orientation (red).Dashed circle indicates the region shown further in (B).(B) Comparison of holo BLVRB aligned to apo-cryo BLVRB (light blue) and apo-RT BLVRB (cyan).The bound NADP + in holo BLVRB is shown for orientation (green).2F o -F c electron density map (blue -1.5σ) is present for indicated arginine side chain and residues 77-80 main-chain (mostly absent).While the Arg78 side chain is visible within the apo-cryo BLVRB structure at 2F o -F c electron density of 0.9σ, this side chain position here is unobservable and deduced from the backbone.All other atoms are colored by type: oxygen-red, nitrogen-blue, phosphorous-purple.(C) Normalized B-factors are shown for each structure (1HDO-red).(D) Cα RMSDs are shown between apo BLVRB-cryo/holo BLVRB (blue) and between apo BLVRB-RT/ holo BLVRB (cyan).
FIGURE 3Amide and Cα CSPs between BLVRB catalytic intermediates.The absolute differences for amide CSPs (navy blue) and Cα CSPs (magenta) are shown (left).The same cutoffs were used to map CSPs in order to compare the number of residues affected.Specifically, while the average CSP plus one standard deviation is listed below quantified from these data, a hard cutoff of the lowest value was used for amide CSPs (center, navy blue balls that are larger than 0.78 ppm) and for Cα CSPs (right, magenta balls that are larger than 0.24 ppm).(A) Apo BLVRB versus BLVRB/NADPH with the average plus one standard deviation of 2.2 ppm for amide CSPs (black dashed line) and 0.49 ppm for Cα CSPs (magenta dashed line).(B) BLVRB/NADPH versus BLVRB/ NADP+:BR with the average plus one standard deviation of 0.90 ppm for amide CSPs (black dashed line) and 0.31 ppm for Cα CSPs (magenta dashed line).(C) BLVRB/NADP+:BR versus BLVRB/NADP+ with the average plus one standard deviation of 0.83 ppm for amide CSPs (black dashed line) and 0.24 ppm for Cα CSPs (magenta dashed line).(D) BLVRB/NADP + versus apo BLVRB with the average plus one standard deviation of 1.91 ppm for amide CSPs (black dashed line) and 0.42 ppm for Cα CSPs (magenta dashed line).(E) BLVRB/NADPH versus BLVRB/NADP + with the average plus one standard deviation of 0.78 ppm for amide CSPs (dashed black line) and 0.26 ppm for Cα CSPs (magenta dashed line).

FIGURE 4
FIGURE 4 Mutagenic changes to hydride-mediated coupling and function.All 15 N-HSQCs were collected at 20 °C at 900 MHz and 3D-HNCAs collected at 20 °C at 600 MHz A) The amide of S111A exhibits a significant CSP between reduced/oxidized forms of coenzyme.3D HNCA planes are shown for S111 due to overlap in the 15 N-HSQC spectra.(B) The amide of H153A exhibits a significant CSP between reduced/oxidized forms of coenzyme.3D HNCA planes are shown for H153 due to overlap in the 15 N-HSQC spectra.(C) The amide of T164S exhibits a significant CSP between reduced/oxidized forms of coenzyme, as shown in 15 N-HSQC spectra.(D) Amide CSPs for NADPH/NADP + of WT reported in Figure 3E (blue) compared to S111A (orange).The average CSP plus one standard deviation is show for S111A is 0.45 ppm (dashed line) (E) Amide CSPs for NADPH/NADP + of WT reported in Figure3E(blue) compared to H153A (green).The average CSP plus one standard deviation is show for H153A is 0.45 ppm (dashed line) (F) Amide CSPs for NADPH/ NADP + of WT reported in Figure3E(blue) compared to T164S (magenta).The average CSP plus one standard deviation is show for T164S is 0.51 ppm(Continued ) FIGURE 6R2-CPMG dispersion profiles for BLVRB catalytic intermediates.R2-CPMG dispersions are shown for three groups of amides that exhibit similar trends for apo BLVRB (black), BLVRB bound to NADPH (blue), BLVRB bound to NADP+ and BR (green), and BLVRB bound to NADP + (red).Intensity-based rates are shown as a function of the imparted CPMG field, ν cpmg .(A) Group I residues are mapped onto the ternary complex (orange spheres) that include V34, D36, R39, V52, and D54 with R2-CPMG dispersions of V34 and D54 amides shown at 10 °C.(B) Group II residues are mapped onto the ternary complex (magenta spheres) that include L74, T77, C109, H111, and Q126 with L74 and S111 amide R2-CPMG dispersions shown at 10 °C.(C) Group III residues that include I154, G155, T164, I176, Y205 are mapped onto the ternary complex with I154 and Y205 amide dispersions shown at 10 °C.All data were collected at 900 MHz and all R2-CPMG dispersions are shown with their fits in Supplementary FigureS5.
. Such inherent plasticity includes examples from

TABLE 1 X
-ray data and model statistics for bothApo-RT and Apo-cryo BLVRB structures.