Expedient Synthesis and Characterization of π-Extended Luciferins.

Bioluminescence imaging enables the sensitive tracking of cell populations and the visualization of biological processes in living systems. Bioluminescent luciferase/luciferin pairs with far-red and near-infrared emission benefit from the reduced competitive absorption by blood and tissue while also facilitating multiplexing strategies. Luciferins with extended π-systems, such as AkaLumine and recently reported CouLuc-1 and -3, can be used for bioluminescence imaging in this long wavelength regime. Existing synthetic routes to AkaLumine and similar π-extended compounds require a multistep sequence to install the thiazoline heterocycle. Here we detail the development of a two-step strategy for accessing these molecules via a Horner-Wadsworth-Emmons reaction and cysteine condensation sequence from readily available aldehyde starting materials. We detail an improved synthesis of AkaLumine, as well as the corresponding two-carbon homologues, Tri- and Tetra-AkaLumine. We then extended this approach to prepare coumarin- and naphthalene-derived luciferins. These putative luciferins were tested against a panel of luciferases to identify capable emitters. Of these, an easily prepared naphthalene derivative exhibits photon emission on par with that of the broadly used Akaluc/AkaLumine pair with similar emission maxima. Overall, this chemistry provides efficient access to several bioluminescent probes for a variety of imaging applications.

Covalent Inhibition by a Natural Product-Inspired Latent Electrophile.

Strategies to target specific protein cysteines are critical to covalent probe and drug discovery. 3-Bromo-4,5-dihydroisoxazole (BDHI) is a natural product-inspired, synthetically accessible electrophilic moiety that has previously been shown to react with nucleophilic cysteines in the active site of purified enzymes. Here, we define the global cysteine reactivity and selectivity of a set of BDHI-functionalized chemical fragments using competitive chemoproteomic profiling methods. Our study demonstrates that BDHIs capably engage reactive cysteine residues in the human proteome and the selectivity landscape of cysteines liganded by BDHI is distinct from that of haloacetamide electrophiles. Given its tempered reactivity, BDHIs showed restricted, selective engagement with proteins driven by interactions between a tunable binding element and the complementary protein sites. We validate that BDHI forms covalent conjugates with glutathione S-transferase Pi (GSTP1) and peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (PIN1), emerging anticancer targets. BDHI electrophile was further exploited in Bruton's tyrosine kinase (BTK) inhibitor design using a single-step late-stage installation of the warhead onto acrylamide-containing compounds. Together, this study expands the spectrum of optimizable chemical tools for covalent ligand discovery and highlights the utility of 3-bromo-4,5-dihydroisoxazole as a cysteine-reactive electrophile.

Red-Shifted Coumarin Luciferins for Improved Bioluminescence Imaging.

Multicomponent bioluminescence imaging in vivo requires an expanded collection of tissue-penetrant probes. Toward this end, we generated a new class of near-infrared (NIR) emitting coumarin luciferin analogues (CouLuc-3s). The scaffolds were easily accessed from commercially available dyes. Complementary mutant luciferases for the CouLuc-3 analogues were also identified. The brightest probes enabled sensitive imaging in vivo. The CouLuc-3 scaffolds are also orthogonal to popular bioluminescent reporters and can be used for multicomponent imaging applications. Collectively, this work showcases a new set of bioluminescent tools that can be readily implemented for multiplexed imaging in a variety of biological settings.

Modified norcyanines enable ratiometric pH imaging beyond 1000 nm.

Activatable fluorophores with emission beyond 1000 nm have the potential to enable high contrast imaging in complex in vivo settings. However, there are few scaffolds that can be applied to this challenge. Here we detail the synthesis and evaluation of benzo[c,d]indole-substituted norcyanines that enable pH responsive fluorescence imaging in the long wavelength (>1150 nm) range. A key component of our molecular design is the installation of a hydrophilic substituted quaternary amine in the central dihydropyridine ring system. A compound with a C4'-phenyl substituent, but not the C4'-protio homologue, exhibits absorbance maxima of 740 nm and 1130 nm in basic and acidic media, respectively, with evidence of J-aggregate-like properties. These two distinct absorbances enabled ratiometric imaging of probe internalization in a tumor model. Overall, these studies provide a new class of activatable long-wavelength responsive fluorophores with promising photophysical properties.

Targeted multicolor in vivo imaging over 1,000 nm enabled by nonamethine cyanines.

Recent progress has shown that using wavelengths between 1,000 and 2,000 nm, referred to as the shortwave-infrared or near-infrared (NIR)-II range, can enable high-resolution in vivo imaging at depths not possible with conventional optical wavelengths. However, few bioconjugatable probes of the type that have proven invaluable for multiplexed imaging in the visible and NIR range are available for imaging these wavelengths. Using rational design, we have generated persulfonated indocyanine dyes with absorbance maxima at 872 and 1,072 nm through catechol-ring and aryl-ring fusion, respectively, onto the nonamethine scaffold. Multiplexed two-color and three-color in vivo imaging using monoclonal antibody and dextran conjugates in several tumor models illustrate the benefits of concurrent labeling of the tumor and healthy surrounding tissue and lymphatics. These efforts are enabled by complementary advances in a custom-built NIR/shortwave-infrared imaging setup and software package for multicolor real-time imaging.

Protonation-Dependent Sequencing of 5-Formylcytidine in RNA.

Chemical modification of cytidine in noncoding RNAs plays a key role in regulating translation and disease. However, the distribution and dynamics of many of these modifications remain unknown due to a lack of sensitive site-specific sequencing technologies. Here, we report a protonation-dependent sequencing reaction for the detection of 5-formylcytidine (5fC) and 5-carboxycytidine (5caC) in RNA. First, we evaluate how protonation combined with electron-withdrawing substituents alters the molecular orbital energies and reduction of modified cytidine nucleosides, highlighting 5fC and 5caC as reactive species. Next, we apply this reaction to detect these modifications in synthetic oligonucleotides as well as endogenous human transfer RNA (tRNA). Finally, we demonstrate the utility of our method to characterize a patient-derived model of 5fC deficiency, where it enables facile monitoring of both pathogenic loss and exogenous rescue of NSUN3-dependent 5fC within the wobble base of human mitochondrial tRNAMet. These studies showcase the ability of protonation to enhance the reactivity and sensitive detection of 5fC in RNA and more broadly provide a molecular foundation for using optimized sequencing reactions to better understand the role of oxidized RNA cytidine residues in diseases.

Structural insights of the conserved "priming loop" of hepatitis B virus pre-genomic RNA.

Initiation of protein-primed (-) strand DNA synthesis in hepatitis B virus (HBV) requires interaction of the viral polymerase with a cis-acting regulatory signal, designated epsilon (ε), located at the 5'-end of its pre-genomic RNA (pgRNA). Binding of polymerase to ε is also necessary for pgRNA encapsidation. While the mechanistic basis of this interaction remains elusive, mutagenesis studies suggest its internal 6-nt "priming loop" provides an important structural contribution. ε might therefore be considered a promising target for small molecule interventions to complement current nucleoside-analog based anti-HBV therapies. An ideal prerequisite to any RNA-directed small molecule strategy would be a detailed structural description of this important element. Herein, we present a solution NMR structure for HBV ε which, in combination with molecular dynamics and docking simulations, reports on a flexible ligand "pocket", reminiscent of those observed in proteins. We also demonstrate the binding of the selective estrogen receptor modulators (SERMs) Raloxifene, Bazedoxifene, and a de novo derivative to the priming loop.Communicated by Ramaswamy H. Sarma.

Ketone Incorporation Extends the Emission Properties of the Xanthene Scaffold Beyond 1000 nm.

Imaging in the shortwave-infrared region (SWIR, λ = 1000-2500 nm) has the potential to enable deep tissue imaging with high resolution. Critical to the development of these methods is the identification of low molecular weight, biologically compatible fluorescent probes that emit beyond 1000 nm. Exchanging the bridging oxygen atom on the xanthene scaffold (C10' position) with electron withdrawing groups has been shown to lead to significant redshifts in absorbance and emission. Guided by quantum chemistry computational modeling studies, we investigated the installation of a ketone bridge at the C10' position. This simple modification extends the absorbance maxima to 860 nm and the emission beyond 1000 nm, albeit with reduced photon output. Overall, these studies demonstrate that broadly applied xanthene dyes can be extended into the SWIR range.

Impact of Cyanine Conformational Restraint in the Near-Infrared Range.

Appending conformationally restraining ring systems to the cyanine chromophore creates exceptionally bright fluorophores in the visible range. Here, we report the application of this strategy in the near-infrared range through the preparation of the first restrained heptamethine indocyanine. Time-resolved absorption spectroscopy and fluorescence correlation spectroscopy verify that, unlike the corresponding parent unrestrained variant, the restrained molecule is not subject to photoisomerization. Notably, however, the room-temperature emission efficiency and the fluorescence lifetime of the restrained cyanine are not extended relative to the parent cyanine, even in viscous solvents. Thus, in contrast to prior reports, the photoisomerization of heptamethine cyanines does not contribute significantly to the excited-state chemistry of these molecules. We also find that the fluorescence lifetime of the restrained heptamethine cyanine is temperature-insensitive and significantly extended at moderately elevated temperatures relative to the parent cyanine. Finally, computational studies have been used to evaluate the impact of the conformational restraint on atomic and orbital structure across the cyanine series. These studies clarify the role of photoisomerization in the heptamethine cyanine scaffold and demonstrate the dramatic effect of restraint on the temperature sensitivity of these dyes.

Core remodeling leads to long wavelength fluoro-coumarins.

Low molecular weight, uncharged far-red and NIR dyes would be enabling for a range of imaging applications. Rational redesign of the coumarin scaffold leads to Fluoro-Coumarins (FCs), the lowest molecular weight dyes with emission maxima beyond 700, 800, and 900 nm. FCs display large Stokes shifts and high environmental sensitivity, with a 40-fold increase in emission intensity in hydrophobic solvents. Untargeted variants exhibit selective lipid droplet and nuclear staining in live cells. Furthermore, sulfo-lipid derivatization enables active targeting to the plasma membrane. Overall, these studies report a promising platform for the development of biocompatible, context-responsive imaging agents.

Repurposing the Pummerer Rearrangement: Determination of Methionine Sulfoxides in Peptides.

The reversible oxidation of methionine residues in proteins has emerged as a biologically important post-translational modification. However, detection and quantitation of methionine sulfoxide in proteins is difficult. Our aim is to develop a method for specifically derivatizing methionine sulfoxide residues. We report a Pummerer rearrangement of methionine sulfoxide treated sequentially with trimethylsilyl chloride and then 2-mercaptoimidazole or pyridine-2-thiol to produce a dithioacetal product. This derivative is stable to standard mass spectrometry conditions, and its formation identified oxidized methionine residues. The scope and requirements of dithioacetal formation are reported for methionine sulfoxide and model substrates. The reaction intermediates have been investigated by computational techniques and by 13 C NMR spectroscopy. These provide evidence for an α-chlorinated intermediate. The derivatization allows for detection and quantitation of methionine sulfoxide in proteins by mass spectrometry and potentially by immunochemical methods.

Defining the conditional basis of silicon phthalocyanine near-IR ligand exchange.

Bond cleavage reactions initiated by long-wavelength light are needed to extend the scope of the caged-uncaged paradigm into complex physiological settings. Axially unsymmetrical silicon phthalocyanines (SiPcs) undergo efficient release of phenol ligands in a reaction contingent on three factors - near-IR light (690 nm), hypoxia, and a thiol reductant. These studies detail efforts to define the mechanistic basis for this unique conditionally-dependent bond cleavage reaction. Spectroscopic studies provide evidence for the formation of a key phthalocyanine radical anion intermediate formed from the triplet state in a reductant-dependent manner. Computational chemistry studies indicate that phenol ligand solvolysis proceeds through a heptacoordinate silicon transition state and that this solvolytic process is favored following SiPc radical anion formation. These results provide insight regarding the central role that radical anion intermediates formed through photoinduced electron transfer with biological reductants can play in long-wavelength uncaging reactions.

Hybrid Correlation Energy (HyCE): An Approach Based on Separate Evaluations of Internal and External Components.

A novel hybrid correlation energy (HyCE) approach is proposed that determines the total correlation energy via distinct computation of its internal and external components. This approach evolved from two related studies. First, rigorous assessment of the accuracies and size extensivities of a number of electron correlation methods, that include perturbation theory (PT2), coupled-cluster (CC), configuration interaction (CI), and coupled electron pair approximation (CEPA), shows that the CEPA(0) variant of the latter and triples-corrected CC methods consistently perform very similarly. These findings were obtained by comparison to near full CI results for four small molecules and by charting recovered correlation energies for six steadily growing chain systems. Second, by generating valence virtual orbitals (VVOs) and utilizing the CEPA(0) method, we were able to partition total correlation energies into internal (or nondynamic) and external (or dynamic) parts for the aforementioned six chain systems and a benchmark test bed of 36 molecules. When using triple-ζ basis sets it was found that per orbital internal correlation energies were appreciably larger than per orbital external energies and that the former showed far more chemical variation than the latter. Additionally, accumulations of external correlation energies were seen to proceed smoothly, and somewhat linearly, as the virtual space is gradually increased. Combination of these two studies led to development of the HyCE approach, whereby the internal and external correlation energies are determined separately by CEPA(0)/VVO and PT2/external calculations, respectively. When applied to the six chain systems and the 36-molecule benchmark test set it was found that HyCE energies followed closely those of triples-corrected CC and CEPA(0) while easily outperforming MP2 and CCSD. The success of the HyCE approach is more notable when considering that its cost is only slightly more than MP2 and significantly cheaper than the CC approaches.

In Vivo Activation of Duocarmycin-Antibody Conjugates by Near-Infrared Light.

Near-IR photocaging groups based on the heptamethine cyanine scaffold present the opportunity to visualize and then treat diseased tissue with potent bioactive molecules. Here we describe fundamental chemical studies that enable biological validation of this approach. Guided by rational design, including computational analysis, we characterize the impact of structural alterations on the cyanine uncaging reaction. A modest change to the ethylenediamine linker (N,N'-dimethyl to N,N'-diethyl) leads to a bathochromic shift in the absorbance maxima, while decreasing background hydrolysis. Building on these structure-function relationship studies, we prepare antibody conjugates that uncage a derivative of duocarmycin, a potent cytotoxic natural product. The optimal conjugate, CyEt-Pan-Duo, undergoes small molecule release with 780 nm light, exhibits activity in the picomolar range, and demonstrates excellent light-to-dark selectivity. Mouse xenograft studies illustrate that the construct can be imaged in vivo prior to uncaging with an external laser source. Significant reduction in tumor burden is observed following a single dose of conjugate and near-IR light. These studies define key chemical principles that enable the identification of cyanine-based photocages with enhanced properties for in vivo drug delivery.

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study.

The metastable ring structure of the ozone 1(1)A1 ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two (1)A1 states. In the present work, valence correlated energies of the 1(1)A1 state and the 2(1)A1 state were calculated at the 1(1)A1 open minimum, the 1(1)A1 ring minimum, the transition state between these two minima, the minimum of the 2(1)A1 state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1(1)A1 state, the present calculations yield the estimates of (ring minimum-open minimum) ∼45-50 mh and (transition state-open minimum) ∼85-90 mh. For the (2(1)A1-(1)A1) excitation energy, the estimate of ∼130-170 mh is found at the open minimum and 270-310 mh at the ring minimum. At the transition state, the difference (2(1)A1-(1)A1) is found to be between 1 and 10 mh. The geometry of the transition state on the 1(1)A1 surface and that of the minimum on the 2(1)A1 surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

Infrared Spectroscopy and Structure of (NO)(n) Clusters.

Nitrogen oxide clusters (NO)n have been studied in He droplets via infrared depletion spectroscopy and by quantum chemical calculations. The ν1 and ν5 bands of cis-ON-NO dimer have been observed at 1868.2 and 1786.5 cm(-1), respectively. Furthermore, spectral bands of the trimer and tetramer have been located in the vicinity of the corresponding dimer bands in accord with computed frequencies that place NO-stretch bands of dimer, trimer, and tetramer within a few wavenumbers of each other. In addition, a new line at 1878.1 cm(-1) close to the band origin of single molecules was assigned to van der Waals bound dimers of (NO)2, which are stabilized due to the rapid cooling in He droplets. Spectra of larger clusters (n > 5), have broad unresolved features in the vicinity of the dimer bands. Experiments and calculations indicate that trimers consist of a dimer and a loosely bound third molecule, whereas the tetramer consists of two weakly bound dimers.

Reactive Species Involved in the Regioselective Photooxidation of Heptamethine Cyanines.

Heptamethine cyanines are important near-IR fluorophores used in many fluorescence applications. Despite this utility, these molecules are susceptible to light-promoted reactions (photobleaching) involving photochemically generated reactive oxygen species (ROS). Here, we have sought to define key chemical aspects of this nearly inescapable process. Near-IR photolysis of a model heptamethine cyanine leads to the regioselective oxidative cleavage of the cyanine polyene. We report the first quantitative analysis of the major reaction pathway following either photolysis or exposure to candidate ROS. These studies clearly indicate that only singlet oxygen (1O2), and not other feasible ROS, recapitulates the direct photolysis pathway. Computational studies were employed to investigate the regioselectivity of the oxidative cleavage process, and the theoretical ratio is comparable to observed experimental values. These results provide a more complete picture of heptamethine cyanine photooxidation, and provide insight for design of improved compounds for future applications.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion.

An analysis based on the variation principle shows that in the molecules H2 (+), H2, B2, C2, N2, O2, F2, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation.