Short Hydrogen Bonds and Proton Delocalization in Green Fluorescent Protein (GFP)

Innovative indenone-derivative colorimetric fluorescent probe: A approach for copper ion detection in water

Copper (Cu2+) is a metal chemical element closely related to human life and is widely used in many fields. However, with the discharge of copper wastewater, the water quality will be seriously affected, leading to excessive intake of Cu2+ and a variety of diseases. Hence, there is a pressing need for an effective detection method for Cu2+ in aqueous environments. Leveraging the remarkable attributes of GFP chromophores and indenone derivatives, we have created a novel colorimetric fluorescent probe P-Cu2+, tailored for efficient copper ion detection. The addition of Cu2+ causes the solution to visibly change from colorless to a pronounced yellow, enabling naked-eye detection and offering promise for real sample analysis.

The Imperative Use of Bacillus Consortium and Quercetin Contributes to Suppress Fusarium Wilt Disease by Direct Antagonism and Induced Resistance

Fusarium wilt diseases severely influence the growth and productivity of numerous crop plants. The consortium of antagonistic rhizospheric Bacillus strains and quercetin were evaluated imperatively as a possible remedy to effectively manage the Fusarium wilt disease of tomato plants. The selection of Bacillus strains was made based on in-vitro antagonistic bioassays against Fusarium oxysporum f.sp. lycoprsici (FOL). Quercetin was selected after screening a library of phytochemicals during in-silico molecular docking analysis using tomato LysM receptor kinases "SILKY12" based on its dual role in symbiosis and plant defense responses. After the selection of test materials, pot trials were conducted where tomato plants were provided consortium of Bacillus strains as soil drenching and quercetin as a foliar spray in different concentrations. The combined application of consortium (Bacillus velezensis strain BS6, Bacillus thuringiensis strain BS7, Bacillus fortis strain BS9) and quercetin (1.0 mM) reduced the Fusarium wilt disease index up to 69%, also resulting in increased plant growth attributes. Likewise, the imperative application of the Bacillus consortium and quercetin (1.0 mM) significantly increased total phenolic contents and activities of the enzymes of the phenylpropanoid pathway. Non-targeted metabolomics analysis was performed to investigate the perturbation in metabolites. FOL pathogen negatively affected a range of metabolites including carbohydrates, amino acids, phenylpropanoids, and organic acids. Thereinto, combined treatment of Bacillus consortium and quercetin (1.0 mM) ameliorated the production of different metabolites in tomato plants. These findings prove the imperative use of Bacillus consortium and quercetin as an effective and sustainable remedy to manage Fusarium wilt disease of tomato plants and to promote the growth of tomato plants under pathogen stress conditions.

Free-Base Corrole Anion

Free-base corroles have long been known to be acidic, readily undergoing deprotonation by mild bases and in polar solvents. The conjugate base, however, has not been structurally characterized until now. Presented here is a first crystal structure of a free-base corrole anion, derived from tris(p-cyanophenyl)corrole, as the tetrabuylammonium salt. The low-temperature (100 K) structure reveals localized hydrogens on a pair of opposite pyrrole nitrogens. DFT calculations identify such a structure as the global minimum but also point to two cis tautomers only 4-7 kcal/mol above the ground state. In terms of free energy, however, the cis tautomers are above or essentially flush with the trans-to-cis barrier so the cis tautomers are unlikely to exist or be observed as true intermediates. Thus, the hydrogen bond within each dipyrrin unit on either side of the molecular pseudo-C2 axis through C10 (i.e., between pyrrole rings A and B or between C and D) qualifies as or closely approaches a low-barrier hydrogen bond. Proton migration across the pseudo-C2 axis entails much higher activation energies >20 kcal/mol, reflecting the relative rigidity of the molecule along the C1-C19 pyrrole-pyrrole linkage.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Symmetry of Hydrogen Bonds: Application of NMR Method of Isotopic Perturbation and Relevance of Solvatomers

Short, strong, symmetric, low-barrier hydrogen bonds (H-bonds) are thought to be of special significance. We have been searching for symmetric H-bonds by using the NMR technique of isotopic perturbation. Various dicarboxylate monoanions, aldehyde enols, diamines, enamines, acid-base complexes, and two sterically encumbered enols have been investigated. Among all of these, we have found only one example of a symmetric H-bond, in nitromalonamide enol, and all of the others are equilibrating mixtures of tautomers. The nearly universal lack of symmetry is attributed to the presence of these H-bonded species as a mixture of solvatomers, meaning isomers (or stereoisomers or tautomers) that differ in their solvation environment. The disorder of solvation renders the two donor atoms instantaneously inequivalent, whereupon the hydrogen attaches to the less well solvated donor. We therefore conclude that there is no special significance to short, strong, symmetric, low-barrier H-bonds. Moreover, they have no heightened stability or else they would have been more prevalent.

Chemical Features and Machine Learning Assisted Predictions of Protein-Ligand Short Hydrogen Bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 A closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

A novel fluorescent probe with a pyrazolo[4,3-c]quinoline core selectively recognizes c-MYC promoter G-quadruplexes

An aptamer combined with a c-MYC-selective fluorophore could work as the fluorescent core of nucleic acid mimics of fluorescent proteins to locate and image functional biomolecules in cells.

Excitation ratiometric chloride sensing in a standalone yellow fluorescent protein is powered by the interplay between proton transfer and conformational reorganization

Natural and laboratory-guided evolution has created a rich diversity of fluorescent protein (FP)-based sensors for chloride (Cl⁻). To date, such sensors have been limited to the Aequorea victoria green fluorescent protein (avGFP) family, and fusions with other FPs have unlocked ratiometric imaging applications. Recently, we identified the yellow fluorescent protein from jellyfish Phialidium sp. (phiYFP) as a fluorescent turn-on, self-ratiometric Cl⁻ sensor. To elucidate its working mechanism as a rare example of a single FP with this capability, we tracked the excited-state dynamics of phiYFP using femtosecond transient absorption (fs-TA) spectroscopy and target analysis. The photoexcited neutral chromophore undergoes bifurcated pathways with the twisting-motion-induced nonradiative decay and barrierless excited-state proton transfer. The latter pathway yields a weakly fluorescent anionic intermediate , followed by the formation of a red-shifted fluorescent state that enables the ratiometric response on the tens of picoseconds timescale. The redshift results from the optimized π–π stacking between chromophore Y66 and nearby Y203, an ultrafast molecular event. The anion binding leads to an increase of the chromophore pK_a and ESPT population, and the hindrance of conversion. The interplay between these two effects determines the turn-on fluorescence response to halides such as Cl⁻ but turn-off response to other anions such as nitrate as governed by different binding affinities. These deep mechanistic insights lay the foundation for guiding the targeted engineering of phiYFP and its derivatives for ratiometric imaging of cellular chloride with high selectivity.

We discovered an interplay between proton transfer and conformational reorganization that powers a standalone fluorescent-protein-based excitation-ratiometric biosensor for chloride imaging.

Structural Origins of Altered Spectroscopic Properties upon Ligand Binding in Proteins Containing a Fluorescent Noncanonical Amino Acid

Fluorescent noncanonical amino acids (fNCAAs) could serve as starting points for the rational design of protein-based fluorescent sensors of biological activity. However, efforts toward this goal are likely hampered by a lack of atomic-level characterization of fNCAAs within proteins. Here, we describe the spectroscopic and structural characterization of five streptavidin mutants that contain the fNCAA l-(7-hydroxycoumarin-4-yl)ethylglycine (7-HCAA) at sites proximal to the binding site of its substrate, biotin. Many of the mutants exhibited altered fluorescence spectra in response to biotin binding, which included both increases and decreases in fluorescence intensity as well as red- or blue-shifted emission maxima. Structural data were also obtained for three of the five mutants. The crystal structures shed light on interactions between 7-HCAA and functional groups, contributed either by the protein or by the substrate, that may be responsible for the observed changes in the 7-HCAA spectra. These data could be used in future studies aimed at the rational design of fluorescent, protein-based sensors of small molecule binding or dissociation.

Controlling Light-Induced Proton Transfer from the GFP Chromophore

Green Fluorescent Protein (GFP) is known to undergo excited-state proton transfer (ESPT). Formation of a short H-bond favors ultrafast ESPT in GFP-like proteins, such as the GFP S65T/H148D mutant, but the detailed mechanism and its quantum nature remain to be resolved. Here we study in vacuo, light-induced proton transfer from the GFP chromophore in hydrogen-bonded complexes with two anionic proton acceptors, I^- and deprotonated trichloroacetic acid (TCA^- ). We address the role of the strong H-bond and the quantum mechanical proton-density distribution in the excited state, which determines the proton-transfer probability. Our study shows that chemical modifications to the molecular network drastically change the proton-transfer probability and it can become strongly wavelength dependent. The proton-transfer branching ratio is found to be 60 % for the TCA complex and 10 % for the iodide complex, being highly dependent on the photon energy in the latter case. Using high-level ab initio calculations, we show that light-induced proton transfer takes place in S₁ , revealing intrinsic photoacid properties of the isolated GFP chromophore in strongly bound H-bonded complexes. ESPT is found to be very sensitive to the topography of the highly anharmonic potential in S₁ , depending on the quantum-density distribution upon vibrational excitation. We also show that the S₁ potential-energy surface, and hence excited-state proton transfer, can be controlled by altering the chromophore microenvironment.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Switching between Ultrafast Pathways Enables a Green-Red Emission Ratiometric Fluorescent-Protein-Based Ca2+ Biosensor

Ratiometric indicators with long emission wavelengths are highly preferred in modern bioimaging and life sciences. Herein, we elucidated the working mechanism of a standalone red fluorescent protein (FP)-based Ca2+ biosensor, REX-GECO1, using a series of spectroscopic and computational methods. Upon 480 nm photoexcitation, the Ca2+-free biosensor chromophore becomes trapped in an excited dark state. Binding with Ca2+ switches the route to ultrafast excited-state proton transfer through a short hydrogen bond to an adjacent Glu80 residue, which is key for the biosensor's functionality. Inspired by the 2D-fluorescence map, REX-GECO1 for Ca2+ imaging in the ionomycin-treated human HeLa cells was achieved for the first time with a red/green emission ratio change (ΔR/R0) of ~300%, outperforming many FRET- and single FP-based indicators. These spectroscopy-driven discoveries enable targeted design for the next-generation biosensors with larger dynamic range and longer emission wavelengths.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

A Self-Assembled Rhombohedral DNA Crystal Scaffold with Tunable Cavity Sizes and High-Resolution Structural Detail

DNA is an ideal molecule for the construction of 3D crystals with tunable properties owing to its high programmability based on canonical Watson-Crick base pairing, with crystal assembly in all three dimensions facilitated by immobile Holliday junctions and sticky end cohesion. Despite the promise of these systems, only a handful of unique crystal scaffolds have been reported. Herein, we describe a new crystal system with a repeating sequence that mediates the assembly of a 3D scaffold via a series of Holliday junctions linked together with complementary sticky ends. By using an optimized junction sequence, we could determine a high-resolution (2.7 Å) structure containing R3 crystal symmetry, with a slight subsequent improvement (2.6 Å) using a modified sticky-end sequence. The immobile Holliday junction sequence allowed us to produce crystals that provided unprecedented atomic detail. In addition, we expanded the crystal cavities by 50 % by adding an additional helical turn between junctions, and we solved the structure to 4.5 Å resolution by molecular replacement.

Structural and spectroscopic characterization of photoactive yellow protein and photoswitchable fluorescent protein constructs containing heavy atoms

Photo-induced structural rearrangements of chromophore-containing proteins are essential for various light-dependent signaling pathways and optogenetic applications. Ultrafast structural and spectroscopic methods have offered insights into these structural rearrangements across many timescales. However, questions still remain about exact mechanistic details, especially regarding photoisomerization of the chromophore within these proteins femtoseconds to picoseconds after photoexcitation. Instrumentation advancements for time-resolved crystallography and ultrafast electron diffraction provide a promising opportunity to study these reactions, but achieving enough signal-to-noise is a constant challenge. Here we present four new photoactive yellow protein constructs and one new fluorescent protein construct that contain heavy atoms either within or around the chromophore and can be expressed with high yields. Structural characterization of these constructs, most at atomic resolution, show minimal perturbation caused by the heavy atoms compared to wild-type structures. Spectroscopic studies report the effects of the heavy atom identity and location on the chromophore's photophysical properties. None of the substitutions prevent photoisomerization, although certain rates within the photocycle may be affected. Overall, these new proteins containing heavy atoms are ideal samples for state-of-theart time-resolved crystallography and electron diffraction experiments to elucidate crucial mechanistic information of photoisomerization.

Dissecting Optical Response and Molecular Structure of Fluorescent Proteins With Non-canonical Chromophores

Tracking the structural dynamics of fluorescent protein chromophores holds the key to unlocking the fluorescence mechanisms in real time and enabling rational design principles of these powerful and versatile bioimaging probes. By combining recent chemical biology and ultrafast spectroscopy advances, we prepared the superfolder green fluorescent protein (sfGFP) and its non-canonical amino acid (ncAA) derivatives with a single chlorine, bromine, and nitro substituent at the ortho site to the phenolate oxygen of the embedded chromophore, and characterized them using an integrated toolset of femtosecond transient absorption and tunable femtosecond stimulated Raman spectroscopy (FSRS), aided by quantum calculations of the vibrational normal modes. A dominant vibrational cooling time constant of ~4 and 11 ps is revealed in Cl-GFP and Br-GFP, respectively, facilitating a ~30 and 12% increase of the fluorescent quantum yield vs. the parent sfGFP. Similar time constants were also retrieved from the transient absorption spectra, substantiating the correlated electronic and vibrational motions on the intrinsic molecular timescales. Key carbon-halogen stretching motions coupled with phenolate ring motions of the deprotonated chromophores at ca. 908 and 890 cm-1 in Cl-GFP and Br-GFP exhibit enhanced activities in the electronic excited state and blue-shift during a distinct vibrational cooling process on the ps timescale. The retrieved structural dynamics change due to targeted site-specific halogenation of the chromophore thus provides an effective means to design new GFP derivatives and enrich the bioimaging probe toolset for life and medical sciences.

Quantum kinetic energy and isotope fractionation in aqueous ionic solutions

At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na⁺, Cl^- and HPO₄^2- ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.

Electrostatic control of photoisomerization pathways in proteins

Rotation around a specific bond after photoexcitation is central to vision and emerging opportunities in optogenetics, super-resolution microscopy, and photoactive molecular devices. Competing roles for steric and electrostatic effects that govern bond-specific photoisomerization have been widely discussed, the latter originating from chromophore charge transfer upon excitation. We systematically altered the electrostatic properties of the green fluorescent protein chromophore in a photoswitchable variant, Dronpa2, using amber suppression to introduce electron-donating and electron-withdrawing groups to the phenolate ring. Through analysis of the absorption (color), fluorescence quantum yield, and energy barriers to ground- and excited-state isomerization, we evaluate the contributions of sterics and electrostatics quantitatively and demonstrate how electrostatic effects bias the pathway of chromophore photoisomerization, leading to a generalized framework to guide protein design.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

Unveiling coupled electronic and vibrational motions of chromophores in condensed phases

The quest for capturing molecular movies of functional systems has motivated scientists and engineers for decades. A fundamental understanding of electronic and nuclear motions, two principal components of the molecular Schrödinger equation, has the potential to enable the de novo rational design for targeted functionalities of molecular machines. We discuss the development and application of a relatively new structural dynamics technique, femtosecond stimulated Raman spectroscopy with broadly tunable laser pulses from the UV to near-IR region, in tracking the coupled electronic and vibrational motions of organic chromophores in solution and protein environments. Such light-sensitive moieties hold broad interest and significance in gaining fundamental knowledge about the intramolecular and intermolecular Hamiltonian and developing effective strategies to control macroscopic properties. Inspired by recent experimental and theoretical advances, we focus on the in situ characterization and spectroscopy-guided tuning of photoacidity, excited state proton transfer pathways, emission color, and internal conversion via a conical intersection.

Structural Evidence of Photoisomerization Pathways in Fluorescent Proteins

Double-bond photoisomerization in molecules such as the green fluorescent protein (GFP) chromophore can occur either via a volume-demanding one-bond-flip pathway or via a volume-conserving hula-twist pathway. Understanding the factors that determine the pathway of photoisomerization would inform the rational design of photoswitchable GFPs as improved tools for super-resolution microscopy. In this communication, we reveal the photoisomerization pathway of a photoswitchable GFP, rsEGFP2, by solving crystal structures of cis and trans rsEGFP2 containing a monochlorinated chromophore. The position of the chlorine substituent in the trans state breaks the symmetry of the phenolate ring of the chromophore and allows us to distinguish the two pathways. Surprisingly, we find that the pathway depends on the arrangement of protein monomers within the crystal lattice: in a looser packing, the one-bond-flip occurs, whereas, in a tighter packing (7% smaller unit cell size), the hula-twist occurs.

Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins

Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.

Role of Hydrogen Bonding in Green Fluorescent Protein-like Chromophore Emission

The fluorescence emission from green fluorescent protein (GFP) is known to be heavily influenced by hydrogen bonding between the core fluorophore and the surrounding side chains or water molecules. Yet how to utilize this feature for modulating the fluorescence of GFP chromophore or GFP-like fluorophore still remains elusive. Here we present theoretical calculations to predict how hydrogen bonding could influence the excited states of the GFP-like fluorophores. These studies provide both a new perspective for understanding the photophysical properties of GFP as well as a solid basis for the rational design of GFP-based fluorophores.

Unraveling the structural and chemical features of biological short hydrogen bonds

Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.

The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.

Perturbation of Short Hydrogen Bonds in Photoactive Yellow Protein via Noncanonical Amino Acid Incorporation

Photoactive yellow protein (PYP) is a small photoreceptor protein that has two unusually short hydrogen bonds between the deprotonated p-coumaric acid chromophore and two amino acids, a tyrosine and a glutamic acid. This has led to considerable debate as to whether the glutamic acid-chromophore hydrogen bond is a low barrier hydrogen bond, with conflicting results in the literature. We have modified the p Ka of the tyrosine by amber suppression and of the chromophore by chemical substitution. X-ray crystal structures of these modified proteins are nearly identical to the wild-type protein, so the heavy atom distance between proton donor and acceptor is maintained, even though these modifications change the relative proton affinity between donor and acceptor. Despite a considerable change in relative proton affinity, the NMR chemical shifts of the hydrogen-bonded protons are only moderately affected. QM/MM calculations were used to explore the protons' potential energy surface and connect the calculated proton position with empirically measured proton chemical shifts. The results are inconsistent with a low barrier hydrogen bond but in all cases are consistent with a localized proton, suggesting an ionic hydrogen bond rather than a low barrier hydrogen bond.

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.

Enhancing the Anti-Solvatochromic Two-Photon Fluorescence for Cirrhosis Imaging by Forming a Hydrogen-Bond Network

Two-photon imaging is an emerging tool for biomedical research and clinical diagnostics. Electron donor-acceptor (D-A) type molecules are the most widely employed two-photon scaffolds. However, current D-A type fluorophores suffer from solvatochromic quenching in aqueous biological samples. To address this issue, we devised a novel class of D-A type green fluorescent protein (GFP) chromophore analogues that form a hydrogen-bond network in water to improve the two-photon efficiency. Our design results in two-photon chalcone (TPC) dyes with 0.80 quantum yield and large two-photon action cross section (210 GM) in water. This strategy to form hydrogen bonds can be generalized to design two-photon materials with anti-solvatochromic fluorescence. To demonstrate the improved in vivo imaging, we designed a sulfide probe based on TPC dyes and monitored endogenous H₂ S generation and scavenging in the cirrhotic rat liver for the first time.

Crystal Structure of Green Fluorescent Protein Clover and Design of Clover-Based Redox Sensors

We have determined the crystal structure of Clover, one of the brightest fluorescent proteins, and found that its T203H/S65G mutations relative to wild-type GFP lock the critical E222 side chain in a fixed configuration that mimics the major conformer of that in EGFP. The resulting equilibrium shift to the predominantly deprotonated chromophore increases the extinction coefficient (EC), opposes photoactivation, and is responsible for the bathochromic shift. Clover's brightness can further be attributed to a π-π stacking interaction between H203 and the chromophore. Consistent with these observations, the Clover G65S mutant reversed the equilibrium shift, dramatically decreased the EC, and made Clover photoactivatable under conditions that activated photoactivatable GFP. Using the Clover structure, we rationally engineered a non-photoactivatable redox sensor, roClover1, and determined its structure as well as that of its parental template, roClover0.1. These high-resolution structures provide deeper insights into structure-function relationships in GFPs and may aid the development of excitation-improved ratiometric biosensors.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

Unveiling Structural Motions of a Highly Fluorescent Superphotoacid by Locking and Fluorinating the GFP Chromophore in Solution

Superphotoacidity involves ultrafast proton motions implicated in numerous chemical and biological processes. We used conformational locking and strategic addition of electron-withdrawing substituents to synthesize a new GFP chromophore analogue: p-HO-3,5-diF-BDI:BF₂ (diF). It is highly fluorescent and exhibits excited-state proton transfer (ESPT) in various solvents, placing it among the strongest photoacids. Tunable femtosecond stimulated Raman spectroscopy with unique resonance conditions and transient absorption are complementarily employed to elucidate the structural basis for superphotoacidity. We reveal a multistep ESPT reaction from diF to methanol with an initial proton dissociation on the ∼600 fs time scale that forms a charge-separated state, stabilized by solvation, and followed by a diffusion-controlled proton transfer on the ∼350 ps time scale. A ∼1580 cm^-1 phenolic ring motion is uncovered to accompany ESPT before 1 ps. This study provides a vivid movie of the photoinduced proton dissociation of a superphotoacid with bright fluorescence, effectively bridging fundamental mechanistic insights to precise control of macroscopic functions.

Directed Evolution of Key Residues in Fluorescent Protein Inverses the Polarity of Voltage Sensitivity in the Genetically Encoded Indicator ArcLight

Genetically encoded calcium indicators (GECIs) produce unprecedentedly large signals that have enabled routine optical recording of single neuron activity in vivo in rodent brain. Genetically encoded voltage indicators (GEVIs) offer a more direct measure of neuronal electrical status, however the signal-to-noise characteristics and signal polarity of the probes developed to date have precluded routine use in vivo. We applied directed evolution to target modulable areas of the fluorescent protein in GEVI ArcLight to create the first GFP-based GEVI (Marina) that exhibits a ΔF/ΔV with a positive slope relationship. We found that only three rounds of site-directed mutagenesis produced a family of “brightening” GEVIs with voltage sensitivities comparable to that seen in the parent probe ArcLight. This shift in signal polarity is an essential first step to producing voltage indicators with signal-to-noise characteristics comparable to GECIs to support widespread use in vivo.

Nuclear quantum effects in a HIV/cancer inhibitor: The case of ellipticine

Ellipticine is a natural product that is currently being actively investigated for its inhibitory cancer and HIV properties. Here we use path-integral molecular dynamics coupled with excited state calculations to characterize the role of nuclear quantum effects on the structural and electronic properties of ellipticine in water, a common biological solvent. Quantum effects collectively enhance the fluctuations of both light and heavy nuclei of the covalent and hydrogen bonds in ellipticine. In particular, for the ellipticine-water system, where the proton donor and acceptor have different proton affinities, we find that nuclear quantum effects (NQEs) strengthen both the strong and the weak H bonds. This is in contrast to what is observed for the cases where the proton affinity of the donors and acceptors is same. These structural fluctuations cause a significant red-shift in the absorption spectra and an increase in the broadening, bringing it into closer agreement with the experiments. Our work shows that nuclear quantum effects alter both qualitatively and quantitatively the optical properties of this biologically relevant system and highlights the importance of the inclusion of these effects in the microscopic understanding of their optical properties. We propose that isotopic substitution will produce a blue shift and a reduction in the broadening of the absorption peak.

Short Carboxylic Acid-Carboxylate Hydrogen Bonds Can Have Fully Localized Protons

Short hydrogen bonds (H-bonds) have been proposed to play key functional roles in several proteins. The location of the proton in short H-bonds is of central importance, as proton delocalization is a defining feature of low-barrier hydrogen bonds (LBHBs). Experimentally determining proton location in H-bonds is challenging. Here, bond length analysis of atomic (1.15-0.98 Å) resolution X-ray crystal structures of the human protein DJ-1 and its bacterial homologue, YajL, was used to determine the protonation states of H-bonded carboxylic acids. DJ-1 contains a buried, dimer-spanning 2.49 Å H-bond between Glu15 and Asp24 that satisfies standard donor-acceptor distance criteria for a LBHB. Bond length analysis indicates that the proton is localized on Asp24, excluding a LBHB at this location. However, similar analysis of the Escherichia coli homologue YajL shows both residues may be protonated at the H-bonded oxygen atoms, potentially consistent with a LBHB. A Protein Data Bank-wide screen identifies candidate carboxylic acid H-bonds in approximately 14% of proteins, which are typically short [⟨d_O-O⟩ = 2.542(2) Å]. Chemically similar H-bonds between hydroxylated residues (Ser/Thr/Tyr) and carboxylates show a trend of lengthening O-O distance with increasing H-bond donor pK_a. This trend suggests that conventional electronic effects provide an adequate explanation for short, charge-assisted carboxylic acid-carboxylate H-bonds in proteins, without the need to invoke LBHBs in general. This study demonstrates that bond length analysis of atomic resolution X-ray crystal structures provides a useful experimental test of certain candidate LBHBs.

The Low Barrier Hydrogen Bond in the Photoactive Yellow Protein: A Vacuum Artifact Absent in the Crystal and Solution

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large ¹H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK_a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical ¹H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.

A Critical Test of the Electrostatic Contribution to Catalysis with Noncanonical Amino Acids in Ketosteroid Isomerase

The vibrational Stark effect (VSE) has been used to measure the electric field in the active site of ketosteroid isomerase (KSI). These measured fields correlate with ΔG(⧧) in a series of conventional mutants, yielding an estimate for the electrostatic contribution to catalysis (Fried et al. Science 2014, 346, 1510-1513). In this work we test this result with much more conservative variants in which individual Tyr residues in the active site are replaced by 3-chlorotyrosine via amber suppression. The electric fields sensed at the position of the carbonyl bond involved in charge displacement during catalysis were characterized using the VSE, where the field sensitivity has been calibrated by vibrational Stark spectroscopy, solvatochromism, and MD simulations. A linear relationship is observed between the electric field and ΔG(⧧) that interpolates between wild-type and more drastic conventional mutations, reinforcing the evaluation of the electrostatic contribution to catalysis in KSI. A simplified model and calculation are developed to estimate changes in the electric field accompanying changes in the extended hydrogen-bond network in the active site. The results are consistent with a model in which the O-H group of a key active site tyrosine functions by imposing a static electrostatic potential onto the carbonyl bond. The model suggests that the contribution to catalysis from the active site hydrogen bonds is of similar weight to the distal interactions from the rest of the protein. A similar linear correlation was also observed between the proton affinity of KSI's active site and the catalytic rate, suggesting a direct connection between the strength of the H-bond and the electric field it exerts.

Dissecting Proton Delocalization in an Enzyme's Hydrogen Bond Network with Unnatural Amino Acids

Extended hydrogen bond networks are a common structural motif of enzymes. A recent analysis proposed quantum delocalization of protons as a feature present in the hydrogen bond network spanning a triad of tyrosines (Y(16), Y(32), and Y(57)) in the active site of ketosteroid isomerase (KSI), contributing to its unusual acidity and large isotope shift. In this study, we utilized amber suppression to substitute each tyrosine residue with 3-chlorotyrosine to test the delocalization model and the proton affinity balance in the triad. X-ray crystal structures of each variant demonstrated that the structure, notably the O-O distances within the triad, was unaffected by 3-chlorotyrosine substitutions. The changes in the cluster's acidity and the acidity's isotope dependence in these variants were assessed via UV-vis spectroscopy and the proton sharing pattern among individual residues with (13)C nuclear magnetic resonance. Our data show pKa detuning at each triad residue alters the proton delocalization behavior in the H-bond network. The extra stabilization energy necessary for the unusual acidity mainly comes from the strong interactions between Y(57) and Y(16). This is further enabled by Y(32), which maintains the right geometry and matched proton affinity in the triad. This study provides a rich picture of the energetics of the hydrogen bond network in enzymes for further model refinement.