Computational Characterization of “Dark” Intermediates in the Ultrafast Deactivation of Photoexcited Bilirubin

Development of Parallel On-the-Fly Crystal Algorithm for Reaction Discovery in Large and Complex Molecular Systems

The parallel on-the-fly Crystal algorithm is a new, efficient global search algorithm for exploring single-state potential energy surfaces and conical intersection seam spaces of a wide range of molecules. Despite major developments, its application to complex molecular systems, especially in the condensed phase, remains challenging due to the high dimensionality of the configurational space. In this work, we address this challenge and extend its applicability to the reaction discovery of large and complex molecular photoswitches in various molecular environments, including in the condensed phase with explicit solvent molecules. This is achieved by performing an explicit exploration of a comparatively large Crystal configurational subspace, while gradually relaxing the remaining degrees of freedom. The new Crystal algorithm is applied to the reaction discovery of bilirubin and donor-acceptor Stenhouse adducts, a next-generation class of molecular photoswitches, in vacuum and in the aqueous solution. To this end, we designed an automated and systematic workflow for Crystal to discover and characterize new minima and low-energy reaction pathways in these challenging and complex systems. Our findings demonstrate the algorithm's effectiveness in quickly exploring the configuration space and uncovering kinetically accessible products, offering new insights into the intricate chemical reactivities of these molecules and the roles of molecular environments on the reaction pathways. The results underscore the promising potential of parallelized global exploration methods for reaction discovery in biomolecular systems.

Bilirubin: Photophysical and photochemical properties, phototherapy, analytical methods of measurement. A short review

Bilirubin, a yellow bile pigment, plays an important role in the body, being a potent antioxidant and having anti-inflammatory, immunomodulatory, cytoprotective, and neuroprotective functions. This makes bilirubin promising as a therapeutic and diagnostic agent in biomedicine. However, excess bilirubin is toxic and should be removed from the body. Bilirubin exhibits photochemical activity, which has been the subject of numerous studies up to now. Such studies are relevant because the bilirubin photochemistry provides the basis for bilirubin removing in phototherapy of neonatal jaundice (neonatal hyperbilirubinemia) and for some therapeutic applications. Furthermore, it can model several elementary processes of molecular photonics. In particular, the bilirubin molecule is capable of ultrafast Z-E photoisomerization and contains two almost identical dipyrromethenone chromophores capable of exciton coupling. The present review considers the data on the photophysical and photochemical properties of bilirubin and ultrafast routes of its phototransformations, as well as its photochemical reactions in phototherapy of neonatal hyperbilirubinemia and the ways to decrease the possible adverse effects of the phototherapy. The main analytical methods of bilirubin measurement in biological systems are also viewed.

Investigation of Ultrafast Configurational Photoisomerization of Bilirubin Using Femtosecond Stimulated Raman Spectroscopy

Phototherapy is an efficient and safe way to reduce high levels of free 4Z,15Z-bilirubin (ZZ-BR) in the serum of newborns. The success of BR phototherapy lies in photoinduced configurational and structural isomerization processes that form excretable isomers. However, the physical picture of photoinduced photoisomerization of ZZ-BR is still unclear. Here, we strategically implement tunable femtosecond stimulated Raman spectroscopy and several time-resolved electronic spectroscopies, assisted by quantum chemical calculations, to dissect the detailed primary configurational isomerization dynamics of free ZZ-BR in organic solvents. The results of this study demonstrate that upon photoexcitation, ultrafast configurational isomerization proceeds by a volume-conserving "hula twist", followed by intramolecular hydrogen-bond distortion and large-scale rotation of the two dipyrrinone halves of the ZZ-BR isomer in a few picoseconds. After that, most of the population recovers back to ZZ-BR, and a very small amount is converted into stable BR isomers via structural isomerization.

Conformational Control of the Photodynamics of a Bilirubin Dipyrrinone Subunit: Femtosecond Spectroscopy Combined with Nonadiabatic Simulations

The photochemistry of bilirubin has been extensively studied due to its importance in the phototherapy of hyperbilirubinemia. In the present work, we investigated the ultrafast photodynamics of a bilirubin dipyrrinone subunit, vinylneoxanthobilirubic acid methyl ester. The photoisomerization and photocyclization reactions of its (E) and (Z) isomers were studied using femtosecond transient absorption spectroscopy and by multireference electronic structure theory, where the nonadiabatic dynamics was modeled with a Landau-Zener surface hopping technique. The following picture has emerged from the combined theoretical and experimental approach. Upon excitation, dipyrrinone undergoes a very fast vibrational relaxation, followed by an internal conversion on a picosecond time scale. The internal conversion leads either to photoisomerization or regeneration of the starting material. Further relaxation dynamics on the order of tens of picoseconds was observed in the ground state. The nonadiabatic simulations revealed a strong conformational control of the photodynamics. The ultrafast formation of a cyclic photochemical product from a less-populated conformer of the studied subunit was predicted by our calculations. We discuss the relevance of the present finding for the photochemistry of native bilirubin. The work has also pointed to the limits of semiclassical nonadiabatic simulations for simulating longer photochemical processes, probably due to the zero-point leakage issue.

Ultrafast spectroscopy of biliverdin dimethyl ester in solution: pathways of excited-state depopulation

Biliverdin is a bile pigment that has a very low fluorescence quantum yield in solution, but serves as a chromophore in far-red fluorescent proteins being developed for bio-imaging. In this work, excited-state dynamics of biliverdin dimethyl ether (BVE) in solvents were investigated using femtosecond (fs) and picosecond (ps) time-resolved absorption and fluorescence spectroscopy. This study is the first fs timescale investigation of BVE in solvents, and therefore revealed numerous dynamics that were not resolved in previous, 200 ps time resolution measurements. Viscosity- and isotope-dependent experiments were performed to identify the contributions of isomerization and proton transfer to the excited-state dynamics. In aprotic solvents, a ∼2 ps non-radiative decay accounts for 95% of the excited-state population loss. In addition, a minor ∼30 ps emissive decay pathway is likely associated with an incomplete isomerization process around the C15[double bond, length as m-dash]C16 double bond that results in a flip of the D-ring. In protic solvents, the dynamics are more complex due to hydrogen bond interactions between solute and solvent. In this case, the ∼2 ps decay pathway is a minor channel (15%), whereas ∼70% of the excited-state population decays through an 800 fs emissive pathway. The ∼30 ps timescale associated with isomerization is also observed in protic solvents. The most significant difference in protic solvents is the presence of a >300 ps timescale in which BVE can decay through an emissive state, in parallel with excited-state proton transfer to the solvent. Interestingly, a small fraction of a luminous species, which we designate lumin-BVE (LBVE), is present in protic solvents.

Ultrafast internal conversion dynamics of bilirubin bound to UnaG and its N57A mutant

Fluorescent proteins (FPs) have become fundamental tools for live cell imaging. Most FPs currently used are members of the green fluorescent protein super-family, but new fluorophores such as bilin-FPs are being developed and optimized. In particular, the UnaG FP incorporates bilirubin (BR) as a chromophore, enhancing its fluorescence quantum yield by three orders of magnitude relative to that in solution. To investigate the mechanism of this dramatic enhancement and provide a basis for further engineering of UnaG and other tetrapyrrole-based fluorophores, we performed picosecond fluorescence and femtosecond transient absorption measurements of BR bound to UnaG and its N57A site-directed mutant. The dynamics of wt-UnaG, which has a fluorescence QY of 0.51, are largely homogeneous, showing an excited state relaxation of ∼200 ps, and a 2.2 ns excited-state lifetime decay with a kinetic isotope effect (KIE) of 1.1 for D2O vs. H2O buffer. In contrast, for UnaG N57A (fluorescence QY 0.01) the results show a large spectral inhomogeneity with excited state decay timescales of 47 and 200 ps and a KIE of 1.4. The non-radiative deactivation of the excited state is limited by proton transfer. The loss of direct hydrogen bonds to the endo-vinyl dipyrrinone moiety of BR leads to high flexibility and structural heterogeneity of UnaG N57A, as seen in the X-ray crystal structure.