Ultrafast photochemistry of bilirubin-albumin complex reveals the molecular origin of jaundice neonatal blue-green phototherapy

Technical Summary: Ultrafast Photochemistry of Bilirubin-Albumin Complex

Problem Statement
Neonatal jaundice, caused by excessive bilirubin (BR) accumulation, is treated clinically with phototherapy. While blue-green light (460–490 nm) is clinically more effective than shorter-wavelength blue light (450 nm) at clearing bilirubin, the molecular origin of this wavelength dependence has remained unclear. Existing explanations rely on tissue optics (e.g., skin penetration, hemoglobin absorption), but in vitro studies of the bilirubin-human serum albumin (BR-HSA) complex suggest an intrinsic photochemical mechanism. The central challenge is to define how the BR-HSA complex, specifically within the high-affinity binding site, directs the initial photochemical dynamics to favor the formation of lumirubin (LR)—the excretable structural isomer—over reversible configurational isomers. A critical unresolved question is how the complex accesses the productive 4E,15Z-BR intermediate required for LR formation, given that the overall quantum yield for LR is low (0.3%).

Methodology
The authors employed a multi-modal approach combining ultrafast spectroscopy with theoretical modeling to resolve the initial photochemical dynamics of the BR-HSA complex under physiologically relevant conditions (mass ratio ~8 mg/g, mimicking severe hyperbilirubinemia).

• Experimental Techniques:
  • Femtosecond Transient Absorption (fs-TA): Used to track excited-state relaxation and ground-state recovery across excitation wavelengths (400–500 nm). Decay-associated difference spectra (DADS) and global fitting were used to resolve kinetic components.
  • Microsecond Transient Absorption (µs-TA): Used to identify long-lived photoproducts and distinguish between reversible isomers and stable species.
  • Femtosecond Stimulated Raman Spectroscopy (FSRS): Utilized to probe excited-state structural dynamics and vibrational signatures with high temporal resolution. Isotope effects (H₂O vs. D₂O) were analyzed to identify proton transfer events.
  • Sample Preparation: BR-HSA complexes were prepared at a molar ratio of 1.1:1 (HSA:BR) to ensure exclusive binding to the primary high-affinity site.
• Theoretical Modeling:
  • Molecular Dynamics (MD) Simulations: Performed to reconstruct the dynamic hydration environment and binding poses of BR within candidate HSA subdomains (IB and IIA).
  • Time-Dependent Density Functional Theory (TD-DFT): Used on cluster models derived from MD snapshots to calculate excited-state properties, electron-hole distributions, and circular dichroism (CD) spectra to validate binding site models.
  • DFT Harmonic Vibrational Analysis: Used to assign Raman bands to specific molecular motions (e.g., HOOP modes, N-H bending).

Key Contributions and Results

1. Identification of the Asymmetric Binding Pocket (Subdomain IIA):
   The study resolves a long-standing debate regarding the BR binding site in HSA. By comparing MD/TD-DFT results for the crystallographic IB site against the proposed IIA site, the authors demonstrate that only the subdomain IIA model accounts for the experimental photophysics.

   • The IIA pocket provides a hydrophilic micro-environment for the AB moiety (via residues R218, R222, R257, K195, and transient water) necessary for proton transfer, whereas the IB site is hydrophobic.
   • TD-DFT calculations on the IIA model correctly predict the low-energy excited state (S₁) to be localized on the AB moiety and the high-energy state (S₂) on the CD moiety, matching the excitonic splitting observed in fs-TA.
   • Calculated circular dichroism spectra for the IIA model match the experimental P-type enantiomer signature, unlike the IB model.

2. Wavelength-Dependent Exciton-Selective Branching:
   fs-TA reveals that HSA breaks the kinetic equivalence of the two excitonic states, creating distinct photochemical pathways:

   • Blue Light (High-Energy Excitation, ~400 nm): Primarily populates the high-energy exciton localized on the CD moiety. This drives Z→E isomerization on the CD side, yielding the long-lived, non-metabolizable 4Z*,15*E-BR** photoproduct. This pathway shows weak kinetic isotope effects (KIE), indicating minimal proton transfer involvement.
   • Blue-Green Light (Low-Energy Excitation, ~500 nm): Selectively populates the low-energy exciton localized on the AB moiety. This pathway is characterized by pronounced proton transfer and is the precursor to lumirubin.

3. Mechanism of Proton-Coupled Isomerization on the AB Moiety:
   FSRS and isotope effect studies elucidate the ultrafast dynamics of the productive AB-side pathway:

   • Coupled Dynamics: Upon 500 nm excitation, the AB moiety undergoes a coupled cascade of Z→E isomerization and proton transfer.
   • Kinetic Steps:
     • ~120 fs: Ultrafast C4=C5 torsion drives the system through an S₁/S₀ conical intersection, forming a vibrationally hot 4E*,15*Z-like ground-state intermediate**. This step involves the collapse of the initial hydrogen bond network (N-H···COO⁻).
     • ~1.1 ps: Completion of proton dissociation from the B-ring N-H group to the COO⁻ acceptor, yielding a deprotonated 4E,15Z species. This step shows strong KIE (values ~1.75–2.09).
     • ~20 ps: Slower skeletal relaxation of the deprotonated framework.
   • Precursor State: The resulting deprotonated 4E,15Z species is structurally and electronically preorganized for photocyclization. The deprotonation enhances nucleophilicity on the B-ring, while the isomerization deplanarizes the chromophore, concentrating electron density to facilitate ring closure.

4. Origin of Blue-Green Phototherapy Superiority:
   The study establishes that blue-green light is superior not merely because it penetrates tissue better, but because it selectively accesses the AB-side pathway that leads to lumirubin.

   • Blue light predominantly generates the dead-end 4Z,15E isomer.
   • Blue-green light efficiently generates the deprotonated 4E,15Z intermediate, which is the necessary precursor for the irreversible formation of lumirubin.
   • Additionally, blue-green light reduces the re-excitation and reverse photoconversion of the formed lumirubin, which absorbs strongly in the blue region.

Significance and Claims
The paper claims to provide the molecular origin of the wavelength dependence in neonatal phototherapy. It posits that the HSA binding pocket creates an asymmetric reaction environment that differentiates the photochemistry of the two dipyrrinone subunits. The superior efficacy of blue-green light is attributed to its ability to selectively excite the low-energy exciton on the AB moiety, triggering a proton-transfer-coupled isomerization that yields a deprotonated 4E,15Z intermediate preorganized for lumirubin formation.

The authors state that these findings offer a mechanistic basis for optimizing next-generation phototherapy devices. By understanding that the productive pathway requires specific excitation of the AB moiety and involves proton transfer, future device designs can be rationally tuned to maximize the population of this specific excitonic state while minimizing the generation of non-metabolizable isomers and the depletion of the final product. The work does not propose new clinical protocols but provides the physical rationale for existing clinical guidelines recommending blue-green light.