Isomerization mechanism of the HcRed fluorescent protein chromophore

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).

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.

Generation of bright monomeric red fluorescent proteins via computational design of enhanced chromophore packing

Red fluorescent proteins (RFPs) have found widespread application in chemical and biological research due to their longer emission wavelengths. Here, we use computational protein design to increase the quantum yield and thereby brightness of a dim monomeric RFP (mRojoA, quantum yield = 0.02) by optimizing chromophore packing with aliphatic residues, which we hypothesized would reduce torsional motions causing non-radiative decay. Experimental characterization of the top 10 designed sequences yielded mSandy1 (λ em = 609 nm, quantum yield = 0.26), a variant with equivalent brightness to mCherry, a widely used RFP. We next used directed evolution to further increase brightness, resulting in mSandy2 (λ em = 606 nm, quantum yield = 0.35), the brightest Discosoma sp. derived monomeric RFP with an emission maximum above 600 nm reported to date. Crystallographic analysis of mSandy2 showed that the chromophore p-hydroxybenzylidene moiety is sandwiched between the side chains of Leu63 and Ile197, a structural motif that has not previously been observed in RFPs, and confirms that aliphatic packing leads to chromophore rigidification. Our results demonstrate that computational protein design can be used to generate bright monomeric RFPs, which can serve as templates for the evolution of novel far-red fluorescent proteins.

Local Electric Field Controls Fluorescence Quantum Yield of Red and Far-Red Fluorescent Proteins

Genetically encoded probes with red-shifted absorption and fluorescence are highly desirable for imaging applications because they can report from deeper tissue layers with lower background and because they provide additional colors for multicolor imaging. Unfortunately, red and especially far-red fluorescent proteins have very low quantum yields, which undermines their other advantages. Elucidating the mechanism of nonradiative relaxation in red fluorescent proteins (RFPs) could help developing ones with higher quantum yields. Here we consider two possible mechanisms of fast nonradiative relaxation of electronic excitation in RFPs. The first, known as the energy gap law, predicts a steep exponential drop of fluorescence quantum yield with a systematic red shift of fluorescence frequency. In this case the relaxation of excitation occurs in the chromophore without any significant changes of its geometry. The second mechanism is related to a twisted intramolecular charge transfer in the excited state, followed by an ultrafast internal conversion. The chromophore twisting can strongly depend on the local electric field because the field can affect the activation energy. We present a spectroscopic method of evaluating local electric fields experienced by the chromophore in the protein environment. The method is based on linear and two-photon absorption spectroscopy, as well as on quantum-mechanically calculated parameters of the isolated chromophore. Using this method, which is substantiated by our molecular dynamics simulations, we obtain the components of electric field in the chromophore plane for seven different RFPs with the same chromophore structure. We find that in five of these RFPs, the nonradiative relaxation rate increases with the strength of the field along the chromophore axis directed from the center of imidazolinone ring to the center of phenolate ring. Furthermore, this rate depends on the corresponding electrostatic energy change (calculated from the known fields and charge displacements), in quantitative agreement with the Marcus theory of charge transfer. This result supports the dominant role of the twisted intramolecular charge transfer mechanism over the energy gap law for most of the studied RFPs. It provides important guidelines of how to shift the absorption wavelength of an RFP to the red, while keeping its brightness reasonably high.

Ultrasound-induced emission color and transmittance changes of organogel based on "trans-to-cis" isomerization

Herein, instant and precise control of fluorescent emission color and transmittance could be carried out by ultrasound-promoted gel-to-gel transition of naphthalimide derivatives containing CN unit. It is proved that ultrasound triggered an irreversible and efficient configuration transformation of N1 from "trans to cis" form in gel state, which is stabilized by intermolecular hydrogen bonding interaction and not observed in the solution state.

Far Red Fluorescent Proteins: Where Is the Limit of the Acylimine Chromophore?

The search for new near-infrared probes for fluorescence imaging applications is a rapidly growing field of research. Monomeric fluorescent proteins that autocatalyze their chromophore are the most versatile markers for in vivo applications, but the development of bright far-red fluorescent proteins (RFPs) has proven difficult. In this contribution, we search for the theoretical limit of the red shift and how it can be reached without sacrificing the fluorescence quantum yield. Through extensive excited-state pathway calculations, molecular dynamics sampling, and statistical modeling using QM/MM schemes, we provide a new understanding of the chromophore's photophysics including the role of its acylimine extension, which is the main difference from other families of fluorescent proteins. The excited-state dynamics of the mPlum RFP and its mutants provide an ideal basis due to mPlum's flexible binding pocket and extended dynamic Stokes shift. We found a large number of structural species with red-shifted emission that differ in rotamer states and H-bonds between key amino acid residues in the binding pocket. By analyzing their spectral and structural features, we derive guidelines for future rational genetic design strategies.

Dynamic water bridging and proton transfer at a surface carboxylate cluster of photosystem II

A hydrogen-bond cluster at a negatively-charged protein interface with a bound protein and long-lived waters might be a proton storage site.

The mechanism of oxidation in chromophore maturation of wild-type green fluorescent protein: a theoretical study

DFT calculations suggested that the thermodynamically unfavourable cyclized product was trapped by oxidation.

Catalytic strategy used by the myosin motor to hydrolyze ATP

Myosin is a molecular motor responsible for biological motions such as muscle contraction and intracellular cargo transport, for which it hydrolyzes adenosine 5'-triphosphate (ATP). Early steps of the mechanism by which myosin catalyzes ATP hydrolysis have been investigated, but still missing are the structure of the final ADP·inorganic phosphate (Pi) product and the complete pathway leading to it. Here, a comprehensive description of the catalytic strategy of myosin is formulated, based on combined quantum-classical molecular mechanics calculations. A full exploration of catalytic pathways was performed and a final product structure was found that is consistent with all experiments. Molecular movies of the relevant pathways show the different reorganizations of the H-bond network that lead to the final product, whose γ-phosphate is not in the previously reported HPγO4(2-) state, but in the H2PγO4(-) state. The simulations reveal that the catalytic strategy of myosin employs a three-pronged tactic: (i) Stabilization of the γ-phosphate of ATP in a dissociated metaphosphate (PγO3(-)) state. (ii) Polarization of the attacking water molecule, to abstract a proton from that water. (iii) Formation of multiple proton wires in the active site, for efficient transfer of the abstracted proton to various product precursors. The specific role played in this strategy by each of the three loops enclosing ATP is identified unambiguously. It explains how the precise timing of the ATPase activation during the force generating cycle is achieved in myosin. The catalytic strategy described here for myosin is likely to be very similar in most nucleotide hydrolyzing enzymes.

Thermal isomerization of the chromoprotein asFP595 and its kindling mutant A143G: QM/MM molecular dynamics simulations

Chromoprotein asFP595 and its A143G variant called kindling fluorescent protein (KFP) are among the chronologically first species for which trans-cis chromophore isomerization has been proposed as a driving force of photoswitching. In spite of long-lasting efforts to characterize the route between protein conformations referring to the trans and cis forms of the chromophore, the molecular mechanism of this transformation is still under debate. We report the results of computational studies of the trans-cis isomerization of the anionic and neutral chromophore inside the protein matrices in the ground electronic state for both variants, asFP595 and KFP. Corresponding free energy profiles (potentials of mean force) were evaluated by using molecular dynamics simulations with the quantum mechanical-molecular mechanical (QM/MM) forces. The computed free energy barrier for the cis-trans ground state (thermal) isomerization reaction is about 2 kcal/mol higher in KFP than that in asFP595. These results provide interpretation of experimental studies on thermal relaxation from the light-induced activation of fluorescence of these proteins and correctly show that the A143G mutation in asFP595 noticeably increases the lifetime of the fluorescence species.