Electrostatic control of photoisomerization pathways in proteins

Ultrafast photophysics of the cyan fluorescent protein chromophore in solution

Incorporation of fluorescent proteins (FPs) into biological systems has revolutionised bioimaging and the understanding of cellular processes. Ongoing developments of FPs are driving efforts to characterise the fundamental photoactive unit (chromophore) embedded within the protein. Cyan FP has a blue emitting chromophore and is widely used in Förster resonance energy transfer studies. Here, we probe the ultrafast photophysics of the cyan FP chromophore in solution using time-resolved fluorescence up-conversion and transient absorption spectroscopies. The ultrafast dynamics are characterised by two lifetimes, sub-picosecond τ1 (or τF) associated with loss of the fluorescent Franck-Condon state, and lifetime τ2 on the order of several picoseconds that is linked with cooling of a hot ground state. MRSF-TDDFT calculations show that the relaxed S1 state equilibrium geometry is classified as a partial twisted intramolecular charge-transfer state, and lies close in energy to a conical intersection seam associated with torsion about the central double bond leading to facile internal conversion. The excited state dynamics exhibit only a weak viscosity dependence, consistent with a barrierless and near-volume-conserving non-radiative decay mechanism. Fluorescence lifetimes for the deprotonated anion are twice those for the neutral.

Deciphering ultrafast structural dynamics of the red Kaede chromophore in solution and a photoconvertible fluorescent protein

Photoconvertible fluorescent proteins (pcFPs) have enabled exquisite images of cellular structures due to their genetic encodability and red-shifted emission with high brightness, hence receiving increased traction in the field. However, the red form of Kaede-like pcFPs after photoconversion remains underexplored. We implemented ultrafast electronic and vibrational spectroscopies on the red Kaede chromophore in solution vs the protein pocket of the least-evolved ancestor (LEA, a Kaede-like green-to-red pcFP) to gain crucial insights into the photophysical processes of the chromophore. The measured fluorescence quantum yield (FQY) values were correlated with ultrafast dynamics to reveal that hydrogen-bonding interactions with the solvent can quench the excited-state Kaede in solution. A viscosity-dependent sub-ps decay indicates nonradiative relaxation involving swift chromophore conformational motions. Femtosecond transient absorption and stimulated Raman spectroscopy (FSRS) reveal an additional ∼1 ps decay of the photoconverted red form of LEA that is absent in green LEA before photoconversion. Transient structural dynamics from FSRS elucidate this decay to involve the phenolate and imidazolinone ring twists that are implicated during cis → trans isomerization and on → off photoswitching in phototransformable fluorescent proteins (FPs). Compared to green-emitting species, the FQY of red LEA (∼0.58) and many other red FPs are often reduced, limiting their applications in modern bioimaging techniques. By shining more light on the often overlooked photoconverted form of pcFPs with ultrafast spectroscopies, we envision such essential mechanistic insights to enable a bottom-up approach for rationally improving the brightness of red-emitting LEA and many other controllable bioprobes, including FPs.

Light-Induced Conformational Heterogeneity Induces Positive Photoswitching in Photoconvertible Fluorescent Proteins of the EosFP Family

Green-to-red photoconvertible fluorescent proteins (PCFPs) of the EosFP family are commonly used in ensemble pulse-chase and single-molecule localization or tracking approaches. However, these fluorescent proteins exhibit highly complex photophysical behaviors. In the green-form, recent NMR experiments revealed that mEos4b and other PCFP variants exist in two different conformational states at thermal equilibrium, which limits their effective photoconversion efficiency. Here, we investigate the conformational heterogeneity of mEos4b in the photoconverted red-form, employing a combination of solution NMR, UV-vis spectroscopy and fluorescence imaging. Only a single red population of mEos4b is observed at thermal equilibrium. However, a second population emerges under illumination with 405 or 488 nm light, which slowly decays in the dark or can be swiftly reverted under 561 nm light. This second population manifests itself through a pH-dependent positive photoswitching mechanism that adds to the already characterized negative photoswitching assigned to cis-trans isomerization of the chromophore. Our data indicate that positive photoswitching, instead, results from the light-induced formation of a second fluorescent state with a cis configuration of the chromophore that exhibits a substantially increased pKa. Such a mechanism, suggested to result from rewiring of the H-bonding network around the first amino acid of the chromophore, adds to the panoply of switching scenarios observed in fluorescent proteins. It bears consequences for the spectroscopic characterization of PCFPs, reduces their apparent brightness and generates short-lived off-times perturbing single-molecule localization microscopy applications.

Turning on the Fluorescence from Isolated GFP Chromophore Anions at Cryogenic Temperatures

The chromophore anion derived from the green fluorescent protein is one of the best-studied chromophores in the gas phase, but attempts to measure fluorescence have failed at room temperature. Here we unequivocally show that the chromophore exhibits fluorescence in the gas phase when cooled to low temperatures (<150  K), thereby validating previous hypotheses. The experimental confirmation is enabled by a unique mass-spectroscopy setup, allowing for fluorescence observation near or at the 0-0 transition without inducing heat in the ions upon photon absorption. The low-temperature conditions effectively simulate the restricted motion experienced within the protein, inhibiting internal conversion via a conical intersection along a twist motion coordinate. Fluorescence-excitation experiments at 100 K reveal an absorption-band maximum at 481.6±0.2  nm, while the dispersed fluorescence spectrum shows maximum emission at 483.6±0.5  nm. Remarkably, both values closely resemble those for proteins cooled to 77 K. We estimate that after excitation at the band maximum, radiation is the only pathway back to the ground state. Franck-Condon simulations at the ωB97XD/aug-cc-pVDZ level of theory nicely reproduce the experimental spectra and identify the fluorescent form to be planar, and that an in-plane scissoring mode (80  cm^{-1}) is active for both absorption and emission.

Substituent effects on the photophysics of the kaede chromophore

Kaede is the prototype of the optical highlighter proteins, which are an important subclass of the fluorescent proteins that can be permanently switched from green to red emitting forms by UV irradiation. This transformation has important applications in bioimaging. Optimising brightness, i.e. enhancing fluorescence characteristics, in these proteins is an important objective. At room temperature, the excited state dynamics of the red form of the kaede chromophore are dominated by a broad distribution of conformers with distinct excited state kinetics. Here, we investigate substituent effects on the photophysics of this form of the kaede chromophore. While an electron withdrawing substituent (nitro) red shifts the electronic spectra, the modified chromophores showed no significant solvatochromism. The lack of solvatochromism suggests small changes in permanent dipole moment between ground and excited electronic states, which is consistent with quantum chemical calculations. Ultrafast fluorescence and transient absorption spectroscopy reveal correlations between radiative and nonradiative decay rates of different conformers in the chromophores. The most significant effect of the substituents is to modify the distribution of conformers. The results are discussed in the context of enhancing brightness of optical highlighter proteins.

Fluorescence Enhancement of Nonemissive Monodeprotonated Luteolin in a Poly(vinyl alcohol) Film

Solid polymer matrixes can modulate the electronic states of embedded chromophores and have been widely used in flexible optoelectronic and optical materials. Luteolin is one of the most common natural flavonoids, and its neutral and monodeprotonated forms are nonemissive in aqueous solution induced by ultrafast excited-state proton transfer (ESPT) followed by nonradiative relaxation. In this study, we have incorporated luteolin into poly(vinyl alcohol) (PVA) films and studied their fluorescence behaviors. Neutral and one monodeprotonated luteolin coexist in the PVA film. Weak steady-state fluorescence of neutral luteolin peaking at about 440 nm is observed for the first time. In addition, the monodeprotonated luteolin in PVA film exhibits obvious fluorescence peaking at 500 nm, with a fluorescence quantum yield of as high as 0.4 and a fluorescence lifetime of as long as 2.4 ns. Time-dependent density functional theory calculations have determined that the ESPT of neutral luteolin is barrierless but that of monodeprotonated luteolin needs to surmount a barrier, explaining their distinct emission properties. These results indicate the modulation ability of the PVA film in both ground-state deprotonation and ESPT, broadening the application areas of the solid polymer matrix.

Modulating the pH dependent photophysical properties of green fluorescent protein

The photophysical properties of the β-barrel superfolder green fluorescent protein (sfGFP) arise from the chromophore that forms post-translationally in the interior of the protein. Specifically, the protonation state of the side chain of tyrosine 66 in the chromophore, in addition to the network of hydrogen bonds between the chromophore and surrounding residues, is directly related to the electronic absorbance and emission properties of the protein. The pH dependence of the photophysical properties of this protein were modulated by the genetic, site-specific incorporation of 3-nitro-l-tyrosine (mNO2Y) at site 66 in sfGFP. The altered photophysical properties of this noncanonical amino acid (ncAA) sfGFP construct were assessed by absorbance and fluorescence spectroscopies. Notably, a comparison of the pK a of the 3-nitrophenol side chain of mNO2Y incorporated in the protein relative to the phenol side chain of the tyrosine at site 66 in the native chromophore as well as the pK a of the 3-nitrophenol side chain of the free ncAA were measured and are compared. A structural analysis of the ncAA containing sfGFP construct is presented to yield molecular insight into the origin of the altered absorbance and fluorescence properties of the protein.

Unified Role of the 145th Residue on the Fluorescence Lifetime of Fluorescent Proteins from the Jellyfish Aequorea victoria

Finding a unified fluorescence mechanism is essential to develop and utilize fluorescent proteins appropriately. Here, we report the unified role of the 145th residue on the fluorescence efficiency of fluorescent proteins developed from the jellyfish Aequorea victoria by demonstrating the difference and similarity between two representative fluorescent proteins, enhanced green fluorescent protein (eGFP), and enhanced yellow fluorescent protein (eYFP). We determined the fluorescence lifetimes of the 19 different Y145 mutants of eGFP and eYFP by picosecond time-resolved fluorescence spectroscopy. We found that the effect of the 145th mutation on the fluorescence lifetime is significant for eYFP but moderate for eGFP. We compared known crystal structures to clarify the observed difference between eGFP and eYFP. As a result, we conclude that the efficiency of the steric restriction of the chromophore motion by the 145th side chain is essentially the same for both eGFP and eYFP. Meanwhile, the restriction of the chromophore motion by hydrogen bonds is more pronounced for eGFP than for YFP. Balance of the steric effect and hydrogen bonding controls the lifetime of the Y145 mutants for eGFP and eYFP. Furthermore, the steric restriction is induced by the electrostatic effect; the different 145th residue induces a different electrostatic environment around the chromophore. The finding in this study reasonably explains the reported lifetimes of other fluorescent proteins and allows the prediction of the lifetime of unknown fluorescent proteins from jellyfish.

XFEL Beamline Optical Instrumentation for Ultrafast Science

Free electron lasers operating in the soft and hard X-ray regime provide capabilities for ultrafast science in many areas, including X-ray spectroscopy, diffractive imaging, solution and material scattering, and X-ray crystallography. Ultrafast time-resolved applications in the picosecond, femtosecond, and attosecond regimes are often possible using single-shot experimental configurations. Aside from X-ray pump and X-ray probe measurements, all other types of ultrafast experiments require the synchronized operation of pulsed laser excitation for resonant or nonresonant pumping. This Perspective focuses on the opportunities for the optical control of structural dynamics by applying techniques from nonlinear spectroscopy to ultrafast X-ray experiments. This typically requires the synthesis of two or more optical pulses with full control of pulse and interpulse parameters. To this end, full characterization of the femtosecond optical pulses is also highly desirable. It has recently been shown that two-color and two-pulse femtosecond excitation of fluorescent protein crystals allowed a Tannor-Rice coherent control experiment, performed under characterized conditions. Pulse shaping and the ability to synthesize multicolor and multipulse conditions are highly desirable and would enable XFEL facilities to offer capabilities for structural dynamics. This Perspective will give a summary of examples of the types of experiments that could be achieved, and it will additionally summarize the laser, pulse shaping, and characterization that would be recommended as standard equipment for time-resolved XFEL beamlines, with an emphasis on ultrafast time-resolved serial femtosecond crystallography.

Non-Degenerate Two-Photon Absorption of Fluorescent Protein Chromophores

Two-photon absorption (2PA), where a pair of photons are absorbed simultaneously, is recognized as a potent bioimaging technique, which depends on the quantified 2PA probability, defined as cross-section (σ2PA). The absorbed photons either have equivalent (ω1 = ω2) or different frequencies (ω1 ≠ ω2), where the former is degenerate 2PA (D-2PA) and the latter is nondegenerate 2PA (ND-2PA). ND-2PA is of particular interest since it is a promising imaging technology with flexibility of photon frequencies and enhanced cross sections, however, it remains a relatively unexplored area compared to D-2PA. This work utilizes time-dependent density functional theory (TD-DFT) and second-order approximate coupled-cluster with the resolution-of-identity approximation (RI-CC2), for the excitation from S0 to S1, to investigate σD-2PA and σND-2PA of FP chromophore models. Interestingly, comparing CAM-B3LYP with the RI-CC2 computations shows qualitative and, in fact, near quantitative agreement in the computed improvements of σND-2PA for comparable (relative) frequency detunings, despite the known underestimations of 2PA cross sections, for TD-DFT results relative to RI-CC2 values. As expected from the 2-state model, the computed values of σND-2PA are quantitatively larger than σD-2PA, where chromophores with the largest values of σD-2PA show greater potential for σND-2PA improvement. Anionic chromophores demonstrated improvements up to 14%, while substantial enhancements were observed in neutral chromophores with some achieving a 30% increase. This work investigates the ND-2PA photophysical characteristics of FP chromophores and identifies qualitative patterns in the computed properties of ND-2PA relative to D-2PA.

The β-subunit of tryptophan synthase is a latent tyrosine synthase

Aromatic amino acids and their derivatives are diverse primary and secondary metabolites with critical roles in protein synthesis, cell structure and integrity, defense and signaling. All de novo aromatic amino acid production relies on a set of ancient and highly conserved chemistries. Here we introduce a new enzymatic transformation for L-tyrosine synthesis by demonstrating that the β-subunit of tryptophan synthase-which natively couples indole and L-serine to form L-tryptophan-can act as a latent 'tyrosine synthase'. A single substitution of a near-universally conserved catalytic residue unlocks activity toward simple phenol analogs and yields exclusive para carbon-carbon bond formation to furnish L-tyrosines. Structural and mechanistic studies show how a new active-site water molecule orients phenols for a nonnative mechanism of alkylation, with additional directed evolution resulting in a net >30,000-fold rate enhancement. This new biocatalyst can be used to efficiently prepare valuable L-tyrosine analogs at gram scales and provides the missing chemistry for a conceptually different pathway to L-tyrosine.

Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins

Red fluorescent protein (RFP) variants are highly sought after for in vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pretwisted chromophore conformation which provides easier access to conical intersections. The main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

Photoswitching the fluorescence of nanoparticles for advanced optical applications

The dynamic optical response properties and the distinct features of nanomaterials make photoswitchable fluorescent nanoparticles (PF NPs) attractive candidates for advanced optical applications. Over the past few decades, the design of PF NPs by coupling photochromic and fluorescent motifs at the nanoscale has been actively pursued, and substantial efforts have been made to exploit their potential applications. In this perspective, we critically summarize various design principles for fabricating these PF NPs. Then, we discuss their distinct optical properties from different aspects by highlighting the capability of NPs in fabricating new, robust photoswitch systems. Afterwards, we introduce the pivotal role of PF NPs in advanced optical applications, including sensing, anti-counterfeiting and imaging. Finally, current challenges and future development of PF NPs are briefly discussed.

Elementary Reactions in the Functional Triads of the Blue-Light Photoreceptor BLUF Domain

The blue light using the flavin (BLUF) domain is one of the smallest photoreceptors in nature, which consists of a unique bidirectional electron-coupled proton relay process in its photoactivation reaction cycle. This perspective summarizes our recent efforts in dissecting the photocycle into three elementary processes, including proton-coupled electron transfer (PCET), proton rocking, and proton relay. Using ultrafast spectroscopy, we have determined the temporal sequence, rates, kinetic isotope effects (KIEs), and concertedness of these elementary steps. Our findings provide important implications for illuminating the photoactivation mechanism of the BLUF domain and suggest an engineering platform to characterize intricate reactions involving proton motions that are ubiquitous in nonphotosensitive protein machines.

Chemical control of excited-state reactivity of the anionic green fluorescent protein chromophore

Controlling excited-state reactivity is a long-standing challenge in photochemistry, as a desired pathway may be inaccessible or compete with other unwanted channels. An important example is internal conversion of the anionic green fluorescent protein (GFP) chromophore where non-selective progress along two competing torsional modes (P: phenolate and I: imidazolinone) impairs and enables Z-to-E photoisomerization, respectively. Developing strategies to promote photoisomerization could drive new areas of applications of GFP-like proteins. Motivated by the charge-transfer dichotomy of the torsional modes, we explore chemical substitution on the P-ring of the chromophore as a way to control excited-state pathways and improve photoisomerization. As demonstrated by methoxylation, selective P-twisting appears difficult to achieve because the electron-donating potential effects of the substituents are counteracted by inertial effects that directly retard the motion. Conversely, these effects act in concert to promote I-twisting when introducing electron-withdrawing groups. Specifically, 2,3,5-trifluorination leads to both pathway selectivity and a more direct approach to the I-twisted intersection which, in turn, doubles the photoisomerization quantum yield. Our results suggest P-ring engineering as an effective approach to boost photoisomerization of the anionic GFP chromophore.

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

Capturing excited-state structural snapshots of evolutionary green-to-red photochromic fluorescent proteins

Photochromic fluorescent proteins (FPs) have proved to be indispensable luminous probes for sophisticated and advanced bioimaging techniques. Among them, an interplay between photoswitching and photoconversion has only been observed in a limited subset of Kaede-like FPs that show potential for discovering the key mechanistic steps during green-to-red photoconversion. Various spectroscopic techniques including femtosecond stimulated Raman spectroscopy (FSRS), X-ray crystallography, and femtosecond transient absorption were employed on a set of five related FPs with varying photoconversion and photoswitching efficiencies. A 3-methyl-histidine chromophore derivative, incorporated through amber suppression using orthogonal aminoacyl tRNA synthetase/tRNA pairs, displays more dynamic photoswitching but greatly reduced photoconversion versus the least-evolved ancestor (LEA). Excitation-dependent measurements of the green anionic chromophore reveal that the varying photoswitching efficiencies arise from both the initial transient dynamics of the bright cis state and the final trans-like photoswitched off state, with an exocyclic bridge H-rocking motion playing an active role during the excited-state energy dissipation. This investigation establishes a close-knit feedback loop between spectroscopic characterization and protein engineering, which may be especially beneficial to develop more versatile FPs with targeted mutations and enhanced functionalities, such as photoconvertible FPs that also feature photoswitching properties.

Optical control of ultrafast structural dynamics in a fluorescent protein

The photoisomerization reaction of a fluorescent protein chromophore occurs on the ultrafast timescale. The structural dynamics that result from femtosecond optical excitation have contributions from vibrational and electronic processes and from reaction dynamics that involve the crossing through a conical intersection. The creation and progression of the ultrafast structural dynamics strongly depends on optical and molecular parameters. When using X-ray crystallography as a probe of ultrafast dynamics, the origin of the observed nuclear motions is not known. Now, high-resolution pump-probe X-ray crystallography reveals complex sub-ångström, ultrafast motions and hydrogen-bonding rearrangements in the active site of a fluorescent protein. However, we demonstrate that the measured motions are not part of the photoisomerization reaction but instead arise from impulsively driven coherent vibrational processes in the electronic ground state. A coherent-control experiment using a two-colour and two-pulse optical excitation strongly amplifies the X-ray crystallographic difference density, while it fully depletes the photoisomerization process. A coherent control mechanism was tested and confirmed the wave packets assignment.

Targeting Ultrafast Spectroscopic Insights into Red Fluorescent Proteins

Red fluorescent proteins (RFPs) represent an increasingly popular class of genetically encodable bioprobes and biomarkers that can advance next-generation breakthroughs across the imaging and life sciences. Since the rational design of RFPs with improved functions or enhanced versatility requires a mechanistic understanding of their working mechanisms, while fluorescence is intrinsically an ultrafast event, a suitable toolset involving steady-state and time-resolved spectroscopic techniques has become powerful in delineating key structural features and dynamic steps which govern irreversible photoconverting or reversible photoswitching RFPs, and large Stokes shift (LSS)RFPs. The pertinent cis-trans isomerization and protonation state change of RFP chromophores in their local environments, involving key residues in protein matrices, lead to rich and complicated spectral features across multiple timescales. In particular, ultrafast excited-state proton transfer in various LSSRFPs showcases the resolving power of wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS) in mapping a photocycle with crucial knowledge about the red-emitting species. Moreover, recent progress in noncanonical RFPs with a site-specifically modified chromophore provides an appealing route for efficient engineering of redder and brighter RFPs, highly desirable for bioimaging. Such an effective feedback loop involving physical chemists, protein engineers, and biomedical microscopists will enable future successes to expand fundamental knowledge and improve human health.

Alkylated green fluorescent protein chromophores: dynamics in the gas phase and in aqueous solution

Fluorescent labelling of macromolecular samples, including using the green fluorescent protein (GFP), has revolutionised the field of bioimaging. The ongoing development of fluorescent proteins require a detailed understanding of the photophysics of the biochromophore, and how chemical derivatisation influences the excited state dynamics. Here, we investigate the photophysical properties associated with the S1 state of three alkylated derivatives of the chromophore in GFP, in the gas phase using time-resolved photoelectron imaging, and in water using femtosecond fluorescence upconversion. The gas-phase lifetimes (1.6-10 ps), which are associated with the intrinsic (environment independent) dynamics, are substantially longer than the lifetimes in water (0.06-3 ps), attributed to stabilisation of both twisted intermediate structures and conical intersection seams in the condensed phase. In the gas phase, alkylation on the 3 and 5 positions of the phenyl ring slows the dynamics due to inertial effects, while a 'pre-twist' of the methine bridge through alkylation on the 2 and 6 positions significantly shortens the excited state lifetimes. Formation of a minor, long-lived (≫ 40 ps) excited state population in the gas phase is attributed to intersystem crossing to a triplet state, accessed because of a T1/S1 degeneracy in the so-called P-trap potential energy minimum associated with torsion of the single-bond in the bridging unit connecting to the phenoxide ring. A small amount of intersystem crossing is supported through TD-DFT molecular dynamics trajectories and MS-CASPT2 calculations. No such intersystem crossing occurs in water at T = 300 K or in ethanol at T ≈ 77 K, due to a significantly altered potential energy surface and P-trap geometry.

Theoretical Understanding on the Facilitated Photoisomerization of a Carbonyl Supported Borane System

Boron compound BOMes2 containing an internal B-O bond undergoes highly efficient photoisomerization, followed by sequential structural transformations, resulting in a rare eight-membered B, O-heterocycle (S. Wang, et al. Org. Lett. 2019, 21, 5285-5289). In this work, the detailed reaction mechanisms of such a unique carbonyl-supported tetracoordinate boron system in the first excited singlet (S1 ) state and the ground (S0 ) state were investigated by using the complete active space self-consistent field and its second-order perturbation (MS-CASPT2//CASSCF) method combined with time-dependent density functional theory (TD-DFT). Moreover, an imine-substituted tetracoordinated organic boron system (BNMes2 ) was selected for comparative study to explore the intrinsic reasons for the difference in reactivity between the two types of compounds. Steric factor was found to influence the photoisomerization activity of BNMes2 and BOMes2 . These results rationalize the experimental observations and can provide helpful insights into understanding the excited-state dynamics of heteroatom-doped tetracoordinate organoboron compounds, which facilitates the rational design of boron-based materials with superior photoresponsive performances.

Fluorescent protein lifetimes report densities and phases of nuclear condensates during embryonic stem-cell differentiation

Fluorescent proteins (FP) are frequently used for studying proteins inside cells. In advanced fluorescence microscopy, FPs can report on additional intracellular variables. One variable is the local density near FPs, which can be useful in studying densities within cellular bio-condensates. Here, we show that a reduction in fluorescence lifetimes of common monomeric FPs reports increased levels of local densities. We demonstrate the use of this fluorescence-based variable to report the distribution of local densities within heterochromatin protein 1α (HP1α) in mouse embryonic stem cells (ESCs), before and after early differentiation. We find that local densities within HP1α condensates in pluripotent ESCs are heterogeneous and cannot be explained by a single liquid phase. Early differentiation, however, induces a change towards a more homogeneous distribution of local densities, which can be explained as a liquid-like phase. In conclusion, we provide a fluorescence-based method to report increased local densities and apply it to distinguish between homogeneous and heterogeneous local densities within bio-condensates.

Structural Characterization of Fluorescent Proteins Using Tunable Femtosecond Stimulated Raman Spectroscopy

The versatile functions of fluorescent proteins (FPs) as fluorescence biomarkers depend on their intrinsic chromophores interacting with the protein environment. Besides X-ray crystallography, vibrational spectroscopy represents a highly valuable tool for characterizing the chromophore structure and revealing the roles of chromophore-environment interactions. In this work, we aim to benchmark the ground-state vibrational signatures of a series of FPs with emission colors spanning from green, yellow, orange, to red, as well as the solvated model chromophores for some of these FPs, using wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS) in conjunction with quantum calculations. We systematically analyzed and discussed four factors underlying the vibrational properties of FP chromophores: sidechain structure, conjugation structure, chromophore conformation, and the protein environment. A prominent bond-stretching mode characteristic of the quinoidal resonance structure is found to be conserved in most FPs and model chromophores investigated, which can be used as a vibrational marker to interpret chromophore-environment interactions and structural effects on the electronic properties of the chromophore. The fundamental insights gained for these light-sensing units (e.g., protein active sites) substantiate the unique and powerful capability of wavelength-tunable FSRS in delineating FP chromophore properties with high sensitivity and resolution in solution and protein matrices. The comprehensive characterization for various FPs across a colorful palette could also serve as a solid foundation for future spectroscopic studies and the rational engineering of FPs with diverse and improved functions.

Stimuli-Induced Subconformation Transformation of the PSI-LHCI Protein at Single-Molecule Resolution

Photosynthesis is a very important process for the current biosphere which can maintain such a subtle and stable circulatory ecosystem on earth through the transformation of energy and substance. Even though been widely studied in various aspects, the physiological activities, such as intrinsic structural vibration and self-regulation process to stress of photosynthetic proteins, are still not in-depth resolved in real-time. Herein, utilizing silicon nanowire biosensors with ultrasensitive temporal and spatial resolution, real-time responses of a single photosystem I-light harvesting complex I (PSI-LHCI) supercomplex of Pisum sativum to various conditions, including gradient variations in temperature, illumination, and electric field, are recorded. Under different temperatures, there is a bi-state switch process associated with the intrinsic thermal vibration behavior. When the variations of illumination and the bias voltage are applied, two additional shoulder states, probably derived from the self-conformational adjustment, are observed. Based on real-time monitoring of the dynamic processes of the PSI-LHCI supercomplex under various conditions, it is successively testified to promising nanotechnology for protein profiling and biological functional integration in photosynthesis studies.

Structural origin and rational development of bright red noncanonical variants of green fluorescent protein

The incorporation of noncanonical amino acids (ncAAs) into fluorescent proteins is promising for red-shifting their fluorescence and benefiting tissue imaging with deep penetration and low phototoxicity. However, ncAA-based red fluorescent proteins (RFPs) have been rare. The 3-aminotyrosine modified superfolder green fluorescent protein (aY-sfGFP) represents a recent advance, yet the molecular mechanism for its red-shifted fluorescence remains elusive while its dim fluorescence hinders applications. Herein, we implement femtosecond stimulated Raman spectroscopy to obtain structural fingerprints in the electronic ground state and reveal that aY-sfGFP possesses a GFP-like instead of RFP-like chromophore. Red color of aY-sfGFP intrinsically arises from a unique "double-donor" chromophore structure that raises ground-state energy and enhances charge transfer, notably differing from the conventional conjugation mechanism. We further developed two aY-sfGFP mutants (E222H and T203H) with significantly improved (∼12-fold higher) brightness by rationally restraining the chromophore's nonradiative decay through electronic and steric effects, aided by solvatochromic and fluorogenic studies of the model chromophore in solution. This study thus provides functional mechanisms and generalizable insights into ncAA-RFPs with an efficient route for engineering redder and brighter fluorescent proteins.

Two-photon excited red-green "discoloration" bioprobes for monitoring lipid droplets and lipid droplet-lysosomal autophagy

Lipid droplets (LDs) and their autophagy by lysosomes are closely related to a variety of physiological and pathological conditions. Therefore, identifying and tracking LDs and the dynamic process of autophagy can provide useful information for the diagnostics and treatment of related diseases. However, few organic small molecule-based fluorescent probes can specifically recognize LDs and dynamically track their autophagy process. Herein, we synthesized a "discoloration" fluorescent bioprobe DPABP-BI with distinguishable features including red fluorescence emission (630 nm), large Stokes shift (145 nm), two-photon excitation and outstanding photostability and biocompatibility. In particular, LDs could be specifically identified via the red fluorescence emission of DPABP-BI (colocalization constant of 0.98), while autophagolysosomes could be visualized via the green fluorescence emission of its acid-hydrolyzed product (colocalization constant of 0.90) to track the autophagy dynamic process. In addition, DPABP-BI enabled the specific recognition of fatty substances in zebrafish larvae. In this study, a two-photon excited red light small molecule probe was constructed to identify LDs and track their autophagy dynamic process by changing the fluorescence emission wavelength.

Multiple Photoisomerization Pathways of the Green Fluorescent Protein Chromophore in a Reversibly Photoswitchable Fluorescent Protein: Insights from Quantum Mechanics/Molecular Mechanics Simulations

Herein, we have employed a combined CASPT2//CASSCF approach within the quantum mechanics/molecular mechanics (QM/MM) framework to explore the early time photoisomerization of rsEGFP2 starting from its two OFF trans states, i.e., Trans1 and Trans2. The results show similar vertical excitation energies to the S₁ state in their Franck-Condon regions. Considering the clockwise and counterclockwise rotations of the C11-C9 bond, four pairs of the S₁ excited-state minima and low-lying S₁/S₀ conical intersections were optimized, based on which we determined four S₁ photoisomerization paths that are essentially barrierless to the relevant S₁/S₀ conical intersections leading to efficient excited-state deactivation to the S₀ state. Most importantly, our work first identified multiple photoisomerization and excited-state decay paths, which must be seriously considered in the future. This work not only sheds significant light on the primary trans-cis photoisomerization of rsEGFP2 but also aids in the understanding of the microscopic mechanism of GFP-like RSFPs and the design of novel GFP-like fluorescent proteins.

Delineating Ultrafast Structural Dynamics of a Green-Red Fluorescent Protein for Calcium Sensing

Fluorescent proteins (FPs) are indispensable tools for noninvasive bioimaging and sensing. Measuring the free cellular calcium (Ca²⁺) concentrations in vivo with genetically encodable FPs can be a relatively direct measure of neuronal activity due to the complex signaling role of these ions. REX-GECO1 is a recently developed red-green emission and excitation ratiometric FP-based biosensor that achieves a high dynamic range due to differences in the chromophore response to light excitation with and without calcium ions. Using steady-state electronic measurements (UV/Visible absorption and emission), along with time-resolved spectroscopic techniques including femtosecond transient absorption (fs-TA) and femtosecond stimulated Raman spectroscopy (FSRS), the potential energy surfaces of these unique biosensors are unveiled with vivid details. The ground-state structural characterization of the Ca²⁺-free biosensor via FSRS reveals a more spacious protein pocket that allows the chromophore to efficiently twist and reach a dark state. In contrast, the more compressed cavity within the Ca²⁺-bound biosensor results in a more heterogeneous distribution of chromophore populations that results in multi-step excited state proton transfer (ESPT) pathways on the sub-140 fs, 600 fs, and 3 ps timescales. These results enable rational design strategies to enlarge the spectral separation between the protonated/deprotonated forms and the Stokes shift leading to a larger dynamic range and potentially higher fluorescence quantum yield, which should be broadly applicable to the calcium imaging and biosensor communities.

4-Cyanotryptophan as a Sensitive Fluorescence Probe of Local Electric Field of Proteins

Electrostatic interactions are key determinants of protein structure, dynamics, and function. Since protein electrostatics are nonuniform, assessment of the internal electric fields (EFs) of proteins requires spatial resolution at the amino acid residue level. In this regard, vibrational Stark spectroscopy, in conjunction with various unnatural amino acid-based vibrational probes, has become a common method for site-specific interrogation of protein EFs. However, application of this method is often limited to proteins with relatively high solubility, due to the intrinsically low oscillator strength of vibrational transitions. Therefore, it would be useful to develop an alternative method that can overcome this limitation. To this end, we show that, using solvatochromic study and molecular dynamics simulations, the frequency of maximum emission intensity of the fluorophore of 4-cyanotryptophan (4CN-Trp), 3-methyl-1H-indole-4-carbonitrile, exhibits a linear dependence on the local EF. Since the absorption and emission spectra of 4CN-Trp are easily distinguishable from those of naturally occurring aromatic amino acids, we believe that this linear relationship provides an easier and more sensitive means to determine the local EF of proteins.

Conformational Dynamics of mCherry Variants: A Link between Side-Chain Motions and Fluorescence Brightness

The 3-fold higher brightness of the recently developed mCherry-XL red fluorescent protein (FP) compared to its progenitor, mCherry, is due to a significant decrease in the nonradiative decay rate underlying its increased fluorescence quantum yield. To examine the structural and dynamic role of the four mutations that distinguish the two FPs and closely related variants, we employed microsecond time scale, all-atom molecular dynamics simulations. The simulations revealed that the I197R mutation leads to the formation of multiple hydrogen-bonded contacts and increased rigidity of the β-barrel. In particular, mCherryXL showed reduced nanosecond time scale breathing of the gap between the β7 and β10-strands, which was previously shown to be the most flexible region of mCherry. Together with experimental results, the simulations also reveal steric interactions of residue 161 and a network of hydrogen-bonding interactions of the chromophore with residues at positions 59, 143, and 163 that are critical in perturbing the chromophore electronic structure. Finally, we shed light on the conformational dynamics of the conserved residues R95 and S146, which are hydrogen-bonded to the chromophore, and provide physical insights into the observed photophysics. To the best of our knowledge, this is the first study that evaluates the conformational space for a set of closely related FPs generated by directed evolution.

Photophysics of the red-form Kaede chromophore

The chromophore responsible for colour switching in the optical highlighting protein Kaede has unexpectedly complicated excited state dynamics, which are measured and analysed here. This will inform the development of new imaging proteins.

On the fluorescence enhancement of arch neuronal optogenetic reporters

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

AMOEBA Force Field Trajectories Improve Predictions of Accurate pKa Values of the GFP Fluorophore: The Importance of Polarizability and Water Interactions

Precisely quantifying the magnitude, direction, and biological functions of electric fields in proteins has long been an outstanding challenge in the field. The most widely implemented experimental method to measure such electric fields at a particular residue in a protein has been through changes in pKa of titratable residues. While many computational strategies exist to predict these values, it has been difficult to do this accurately or connect predicted results to key structural or mechanistic features of the molecule. Here, we used experimentally determined pKa values of the fluorophore in superfolder green fluorescent protein (GFP) with amino acid mutations made at position Thr 203 to evaluate the pKa prediction ability of molecular dynamics (MD) simulations using a polarizable force field, AMOEBA. Structure ensembles from AMOEBA were used to calculate pKa values of the GFP fluorophore. The calculated pKa values were then compared to trajectories using a conventional fixed charge force field (Amber03 ff). We found that the position of water molecules included in the pKa calculation had opposite effects on the pKa values between the trajectories from AMOEBA and Amber03 force fields. In AMOEBA trajectories, the inclusion of water molecules within 35 Å of the fluorophore decreased the difference between the predicted and experimental values, resulting in calculated pKa values that were within an average of 0.8 pKa unit from the experimental results. On the other hand, in Amber03 trajectories, including water molecules that were more than 5 Å from the fluorophore increased the differences between the calculated and experimental pKa values. The inaccuracy of pKa predictions determined from Amber03 trajectories was caused by a significant stabilization of the deprotonated chromophore's free energy compared to the result in AMOEBA. We rationalize the cutoffs for explicit water molecules when calculating pKa to better predict the electrostatic environment surrounding the fluorophore buried in GFP. We discuss how the results from this work will assist the prospective prediction of pKa values or other electrostatic effects in a wide variety of folded proteins.

Steric and Electronic Origins of Fluorescence in GFP and GFP-like Proteins

Fluorescent proteins have become routine tools for biological imaging. However, their nanosecond lifetimes on the excited state present computational hurdles to a full understanding of these photoactive proteins. In this work, we simulate approximately 0.5 nanoseconds of ab initio molecular dynamics to elucidate steric and electronic features responsible for fluorescent protein behavior. Using green fluorescent protein (GFP) and Dronpa2─widely used fluorescent proteins with contrasting functionality─as case studies, we leverage previous findings in the gas phase and solution to explore the deactivation mechanisms available to these proteins. Starting with ground-state analyses, we identify steric (the distribution of empty pockets near the chromophore) and electronic (electric fields exerted on chromophore moieties) factors that offer potential avenues for rational design. Picosecond timescale simulations on the excited state reveal that the chromophore can access twisted structures in Dronpa2, while the chromophore is largely confined to planarity in GFP. We couple ab initio multiple spawning (AIMS) and enhanced sampling simulations to discover and characterize conical intersection seams that facilitate internal conversion, which is a rare event in both systems. Our AIMS simulations correctly capture the relative fluorescence profiles of GFP and Dronpa2 within the first few picoseconds, and we attribute the diminished fluorescence intensity of Dronpa2, relative to GFP, to flexible chromophore intermediates on the excited state. Furthermore, we predict that twisted chromophore intermediates produce red-shifted intensities in the Dronpa2 fluorescence spectrum. If confirmed experimentally, this spectroscopic signature would provide valuable insights when screening and developing novel fluorescent proteins.

Directed Evolution of a Bright Variant of mCherry: Suppression of Nonradiative Decay by Fluorescence Lifetime Selections

The approximately linear scaling of fluorescence quantum yield (ϕ) with fluorescence lifetime (τ) in fluorescent proteins (FPs) has inspired engineering of brighter fluorophores based on screening for increased lifetimes. Several recently developed FPs such as mTurquoise2, mScarlet, and FusionRed-MQV which have become useful for live cell imaging are products of lifetime selection strategies. However, the underlying photophysical basis of the improved brightness has not been scrutinized. In this study, we focused on understanding the outcome of lifetime-based directed evolution of mCherry, which is a popular red-FP (RFP). We identified four positions (W143, I161, Q163, and I197) near the FP chromophore that can be mutated to create mCherry-XL (eXtended Lifetime: ϕ = 0.70; τ = 3.9 ns). The 3-fold higher quantum yield of mCherry-XL is on par with that of the brightest RFP to date, mScarlet. We examined selected variants within the evolution trajectory and found a near-linear scaling of lifetime with quantum yield and consistent blue-shifts of the absorption and emission spectra. We find that the improvement in brightness is primarily due to a decrease in the nonradiative decay of the excited state. In addition, our analysis revealed the decrease in nonradiative rate is not limited to the blue-shift of the energy gap and changes in the excited state reorganization energy. Our findings suggest that nonradiative mechanisms beyond the scope of energy-gap models such the Englman-Jortner model are suppressed in this lifetime evolution trajectory.

Molecular origins of the multi-donor strategy in inducing bathochromic shifts and enlarging Stokes shifts of fluorescent proteins

Long-wavelength fluorescent proteins (LWFPs) and LWFP-based sensors are indispensable tools for bioimaging and biosensing applications. However, it remains challenging to develop LWFPs with outstanding brightness and/or sensitivities, largely due to the lack of simple and effective molecular design strategies. Herein, we rationalized the molecular origins of a multi-donor strategy that affords significant bathochromic shifts and large Stokes shifts with minimal structural changes in the resulting protein fluorophores. We analyzed three key factors that affect the spectral properties of these fluorophores, including the (1) substituent position, (2) electron-donating strength, and (3) number of electron-donating groups. We further demonstrated that this simple design strategy is generalizable to various fluorophore families. We expect that this work can provide rational guidelines for developing fluorescent proteins (and small-molecule fluorophores) with long emission wavelengths and large Stokes shifts.

Photoswitchable Fluorescent Proteins: Mechanisms on Ultrafast Timescales

The advancement of super-resolution imaging (SRI) relies on fluorescent proteins with novel photochromic properties. Using light, the reversibly switchable fluorescent proteins (RSFPs) can be converted between bright and dark states for many photocycles and their emergence has inspired the invention of advanced SRI techniques. The general photoswitching mechanism involves the chromophore cis-trans isomerization and proton transfer for negative and positive RSFPs and hydration-dehydration for decoupled RSFPs. However, a detailed understanding of these processes on ultrafast timescales (femtosecond to millisecond) is lacking, which fundamentally hinders the further development of RSFPs. In this review, we summarize the current progress of utilizing various ultrafast electronic and vibrational spectroscopies, and time-resolved crystallography in investigating the on/off photoswitching pathways of RSFPs. We show that significant insights have been gained for some well-studied proteins, but the real-time "action" details regarding the bidirectional cis-trans isomerization, proton transfer, and intermediate states remain unclear for most systems, and many other relevant proteins have not been studied yet. We expect this review to lay the foundation and inspire more ultrafast studies on existing and future engineered RSFPs. The gained mechanistic insights will accelerate the rational development of RSFPs with enhanced two-way switching rate and efficiency, better photostability, higher brightness, and redder emission colors.

Oriented internal electrostatic fields: an emerging design element in coordination chemistry and catalysis

The power of oriented electrostatic fields (ESFs) to influence chemical bonding and reactivity is a phenomenon of rapidly growing interest. The presence of strong ESFs has recently been implicated as one of the most significant contributors to the activity of select enzymes, wherein alignment of a substrate's changing dipole moment with a strong, local electrostatic field has been shown to be responsible for the majority of the enzymatic rate enhancement. Outside of enzymology, researchers have studied the impacts of "internal" electrostatic fields via the addition of ionic salts to reactions and the incorporation of charged functional groups into organic molecules (both experimentally and computationally), and "externally" via the implementation of bulk fields between electrode plates. Incorporation of charged moieties into homogeneous inorganic complexes to generate internal ESFs represents an area of high potential for novel catalyst design. This field has only begun to materialize within the past 10 years but could be an area of significant impact moving forward, since it provides a means for tuning the properties of molecular complexes via a method that is orthogonal to traditional strategies, thereby providing possibilities for improved catalytic conditions and novel reactivity. In this perspective, we highlight recent developments in this area and offer insights, obtained from our own research, on the challenges and future directions of this emerging field of research.

Structure-based analysis and evolution of a monomerized red-colored chromoprotein from the Olindias formosa jellyfish

GFP-like chromoproteins (CPs) with non-fluorescence ability have been used as bioimaging probes. Existing CPs have voids in the optical absorption window which limits their extensibility. The development of new CP color is therefore ongoing. Here, we cloned CPs from the jellyfish, Olindias formosa, and developed a completely non-fluorescent monomeric red CP, R-Velour, with an absorption peak at 528 nm. To analyze the photophysical properties from a structural aspect, we determined the crystal structure of R-Velour at a 2.1 Å resolution. R-Velour has a trans-chromophore similar to the green fluorescence protein, Gamillus, derived from the same jellyfish. However, in contrast to the two coplanar chromophoric rings in Gamillus, R-Velour has a large torsion inducing non-fluorescence property. Through site-directed mutagenesis, we surveyed residues surrounding the chromophore and found a key residue, Ser155, which contributes to the generation of four-color variants with the bathochromic and hypsochromic shift of the absorption peak, ranging from 506 to 554 nm. The recently proposed spectrum shift theory, based on the Marcus-Hush model, supports the spectrum shift of these mutants. These findings may support further development of R-Velour variants with useful absorption characteristics for bioimaging, including fluorescence lifetime imaging and photoacoustic imaging.

Prediction of Fluorophore Brightness in Designed Mini Fluorescence Activating Proteins

The de novo computational design of proteins with predefined three-dimensional structure is becoming much more routine due to advancements both in force fields and algorithms. However, creating designs with functions beyond folding is more challenging. In that regard, the recent design of small beta barrel proteins that activate the fluorescence of an exogenous small molecule chromophore (DFHBI) is noteworthy. These proteins, termed mini fluorescence activating proteins (mFAPs), have been shown to increase the brightness of the chromophore more than 100-fold upon binding to the designed ligand pocket. The design process created a large library of variants with different brightness levels but gave no rational explanation for why one variant was brighter than another. Here, we use quantum mechanics and molecular dynamics simulations to investigate how molecular flexibility in the ground and excited states influences brightness. We show that the ability of the protein to resist dihedral angle rotation of the chromophore is critical for predicting brightness. Our simulations suggest that the mFAP/DFHBI complex has a rough energy landscape, requiring extensive ground-state sampling to achieve converged predictions of excited-state kinetics. While computationally demanding, this roughness suggests that mFAP protein function can be enhanced by reshaping the energy landscape toward conformations that better resist DFHBI bond rotation.

Unifying structural descriptors for biological and bioinspired nanoscale complexes

Biomimetic nanoparticles are known to serve as nanoscale adjuvants, enzyme mimics and amyloid fibrillation inhibitors. Their further development requires better understanding of their interactions with proteins. The abundant knowledge about protein-protein interactions can serve as a guide for designing protein-nanoparticle assemblies, but the chemical and biological inputs used in computational packages for protein-protein interactions are not applicable to inorganic nanoparticles. Analysing chemical, geometrical and graph-theoretical descriptors for protein complexes, we found that geometrical and graph-theoretical descriptors are uniformly applicable to biological and inorganic nanostructures and can predict interaction sites in protein pairs with accuracy >80% and classification probability ~90%. We extended the machine-learning algorithms trained on protein-protein interactions to inorganic nanoparticles and found a nearly exact match between experimental and predicted interaction sites with proteins. These findings can be extended to other organic and inorganic nanoparticles to predict their assemblies with biomolecules and other chemical structures forming lock-and-key complexes.

Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol

In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme

The past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an "excited-state enzyme". Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency

The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

Complexation of Green and Red Kaede Fluorescent Protein Chromophores by a Zwitterion to Probe Electrostatic and Induction Field Effects

The photophysics of green fluorescent protein (GFP) and red Kaede fluorescent protein (rKFP) are defined by the intrinsic properties of the light-absorbing chromophore and its interaction with the protein binding pocket. This work deploys photodissociation action spectroscopy to probe the absorption profiles for a series of synthetic GFP and rKFP chromophores as the bare anions and as complexes with the betaine zwitterion, which is assumed as a model for dipole microsolvation. Electronic structure calculations and energy decomposition analysis using Symmetry-Adapted Perturbation Theory are used to characterize gas-phase structures and complex cohesion forces. The calculations reveal a preponderance for coordination of betaine to the phenoxide deprotonation site predominantly through electrostatic forces. Calculations using the STEOM-DLPNO-CCSD method are able to reproduce absolute and relative vertical excitation energies for the bare anions and anion-betaine complexes. On the other hand, treatment of the betaine molecule with a point-charge model, in which the charges are computed from some common electron density population analysis schemes, show that just electrostatic and point-charge induction interactions are unable to account for the betaine-induced spectral shift. The present methodology could be applied to investigate cluster forces and optical properties in other gas-phase ion-zwitterion complexes.

Fragment-Based Quantum Mechanical Calculation of Excited-State Properties of Fluorescent RNAs

Fluorescent RNA aptamers have been successfully applied to track and tag RNA in a biological system. However, it is still challenging to predict the excited-state properties of the RNA aptamer–fluorophore complex with the traditional electronic structure methods due to expensive computational costs. In this study, an accurate and efficient fragmentation quantum mechanical (QM) approach of the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) scheme was applied for calculations of excited-state properties of the RNA aptamer–fluorophore complex. In this method, the excited-state properties were first calculated with one-body fragment quantum mechanics/molecular mechanics (QM/MM) calculation (the excited-state properties of the fluorophore) and then corrected with a series of two-body fragment QM calculations for accounting for the QM effects from the RNA on the excited-state properties of the fluorophore. The performance of the EE-GMFCC on prediction of the absolute excitation energies, the corresponding transition electric dipole moment (TEDM), and atomic forces at both the TD-HF and TD-DFT levels was tested using the Mango-II RNA aptamer system as a model system. The results demonstrate that the calculated excited-state properties by EE-GMFCC are in excellent agreement with the traditional full-system time-dependent ab initio calculations. Moreover, the EE-GMFCC method is capable of providing an accurate prediction of the relative conformational excited-state energies for different configurations of the Mango-II RNA aptamer system extracted from the molecular dynamics (MD) simulations. The fragmentation method further provides a straightforward approach to decompose the excitation energy contribution per ribonucleotide around the fluorophore and then reveals the influence of the local chemical environment on the fluorophore. The applications of EE-GMFCC in calculations of excitation energies for other RNA aptamer–fluorophore complexes demonstrate that the EE-GMFCC method is a general approach for accurate and efficient calculations of excited-state properties of fluorescent RNAs.