A concerted mechanism of proton transfer in green fluorescent protein. A theoretical study

Quantum dynamics of excited state proton transfer in green fluorescent protein

Photoexcitation of green fluorescent protein (GFP) triggers long-range proton transfer along a "wire" of neighboring protein residues, which, in turn, activates its characteristic green fluorescence. The GFP proton wire is one of the simplest, most well-characterized models of biological proton transfer but remains challenging to simulate due to the sensitivity of its energetics to the surrounding protein conformation and the possibility of non-classical behavior associated with the movement of lightweight protons. Using a direct dynamics variational multiconfigurational Gaussian wavepacket method to provide a fully quantum description of both electrons and nuclei, we explore the mechanism of excited state proton transfer in a high-dimensional model of the GFP chromophore cluster over the first two picoseconds following excitation. During our simulation, we observe the sequential starts of two of the three proton transfers along the wire, confirming the predictions of previous studies that the overall process starts from the end of the wire furthest from the fluorescent chromophore and proceeds in a concerted but asynchronous manner. Furthermore, by comparing the full quantum dynamics to a set of classical trajectories, we provide unambiguous evidence that tunneling plays a critical role in facilitating the leading proton transfer.

A Not Obvious Correlation Between the Structure of Green Fluorescent Protein Chromophore Pocket and Hydrogen Bond Dynamics: A Choreography From ab initio Molecular Dynamics

The Green Fluorescent Protein (GFP) is a widely studied chemical system both for its large amount of applications and the complexity of the excited state proton transfer responsible of the change in the protonation state of the chromophore. A detailed investigation on the structure of the chromophore environment and the influence of chromophore form (either neutral or anionic) on it is of crucial importance to understand how these factors could potentially influence the protein function. In this study, we perform a detailed computational investigation based on the analysis of ab-initio molecular dynamics simulations, to disentangle the main structural quantities determining the fine balance in the chromophore environment. We found that specific hydrogen bonds interactions directly involving the chromophore (or not), are correlated to quantities, such as the volume of the cavity in which the chromophore is embedded and that it is importantly affected by the chromophore protonation state. The cross-correlation analysis performed on some of these hydrogen bonds and the cavity volume, demonstrates a direct correlation among them and we also identified the ones specifically involved in this correlation. We also found that specific interactions among residues far in the space are correlated, demonstrating the complexity of the chromophore environment and that many structural quantities have to be taken into account to properly describe and understand the main factors tuning the active site of the protein. From an overall evaluation of the results obtained in this work, it is shown that the residues which a priori are perceived to be spectators play instead an important role in both influencing the chromophore environment (cavity volume) and its dynamics (cross-correlations among spatially distant residues).

The mechanism of a green fluorescent protein proton shuttle unveiled in the time-resolved frequency domain by excited state ab initio dynamics

We simulated an excited state proton transfer in green fluorescent protein by excited state ab initio dynamics, and examined the reaction mechanism in both the time and the frequency domain through a multi resolution wavelet analysis. This original approach allowed us, for the first time, to directly compare the trends of photoactivated vibrations to femtosecond stimulated Raman spectroscopy results, and to give an unequivocal interpretation of the role played by low frequency modes in promoting the reaction. We could attribute the main driving force of the reaction to an important photoinduced softening of the ring-ring orientational motion of the chromophore, thus permitting the tightening of the hydrogen bond network and the opening of the reaction pathway. We also found that both the chromophore (in terms of its inter-ring dihedral angle and phenolic C-O and imidazolinone C-N bond distances) and its pocket (in terms of the inter-molecular oxygen's dihedral angle of the chromophore pocket) relaxations are modulated by low frequency (about 120 cm-1) modes involving the oxygen atoms of the network. This is in agreement with the femtosecond Raman spectroscopy findings in the time-frequency domain. Moreover, the rate in proximity to the Franck Condon region involves a picosecond time scale, with a significant influence from fluctuations of nearby hydrogen bonded residues such as His148. This approach opens a new scenario with ab initio simulations as routinely used tools to understand photoreactivity and the results of advanced time resolved spectroscopy techniques.

Role of histidine 148 in stability and dynamics of a highly fluorescent GFP variant

The armory of GFP mutants available to biochemists and molecular biologists is huge. Design and selection of mutants are usually driven by tailored spectroscopic properties, but some key aspects of stability, folding and dynamics of selected GFP variants still need to be elucidated. We have prepared, expressed and characterized three H148 mutants of the highly fluorescent variant GFPmut2. H148 is known to be involved in the H-bonding network surrounding the chromophore, and all the three mutants, H148G, H148R and H148K, show increased pKa values of the chromophore. Only H148G GFPmut2 (Mut2G) gave good expression and purification yields, indicating that position 148 is critical for efficient folding in vivo. The chemical denaturation of Mut2G was monitored by fluorescence emission, absorbance and far-UV circular dichroism spectroscopy. The mutation has little effect on the spectroscopic properties of the protein and on its stability in solution. However, the unfolding kinetics of the protein encapsulated in wet nanoporous silica gels, a system that allows to stabilize conformations that are poorly or only transiently populated in solution, indicate that the unfolding pathway of Mut2G is markedly different from the parent molecule. In particular, encapsulation allowed to identify an unfolding intermediate that retains a native-like secondary structure despite a destructured chromophore environment. Thus, H148 is a critical residue not only for the chromophoric and photodynamic properties, but also for the correct folding of GFP, and its substitution has great impact on expression yields and stability of the mature protein.

Isomerization mechanism of the HcRed fluorescent protein chromophore

To understand how the protein achieves fluorescence, the isomerization mechanism of the HcRed chromophore is studied both under vacuum and in the solvated red fluorescent protein. Quantum mechanical (QM) and quantum mechanical/molecular mechanical (QM/MM) methods are applied both for the ground and the first excited state. The photoinduced processes in the chromophore mainly involve torsions around the imidazolinone-bridge bond (τ) and the phenoxy-bridge bond (φ). Under vacuum, the isomerization of the cis-trans chromophore essentially proceeds by τ twisting, while the radiationless decay requires φ torsion. By contrast, the isomerization of the cis-trans chromophore in HcRed occurs via simultaneous τ and φ twisting. The protein environment significantly reduces the barrier of this hula twist motion compared with vacuum. The excited-state isomerization barrier via the φ rotation of the cis-coplanar conformer in HcRed is computed to be significantly higher than that of the trans-non-coplanar conformer. This is consistent with the experimental observation that the cis-coplanar-conformation of the chromophore is related to the fluorescent properties of HcRed, while the trans-non-planar conformation is weakly fluorescent or non-fluorescent. Our study shows how the protein modifies the isomerization mechanism, notably by interactions involving the nearby residue Ile197, which keeps the chromophore coplanar and blocks the twisting motion that leads to photoinduced radiationless decay.

QM/MM studies of structural and energetic properties of the far-red fluorescent protein HcRed

The far-red fluorescent protein HcRed was investigated using molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Three models of HcRed (anionic chromophore) were considered, differing in the protonation states of nearby Glu residues (A: Glu214 and Glu146 both protonated; B: Glu214 protonated and Glu146 deprotonated; C: Glu214 and Glu146 both deprotonated). SCC-DFTB/MM MD simulations of model B yield good agreement with the available crystallographic data at ambient pH. Bond lengths in the QM region are well reproduced, with a root mean square (rms) deviation between experimental and average MD data of 0.079 A; the chromophore is almost co-planar, which is consistent with experimental observation; and the five hydrogen bonds involving the chromophore are conserved. QM/MM geometry optimizations were performed on representative snapshot structures from the MD simulations for each model. They confirm the structural features observed in the MD simulations. According to the DFT(B3LYP)/MM results, the cis-conformation of the chromophore is more stable than the trans-form by 9.1-12.9 kcal mol(-1) in model B, and by 12.4-19.9 kcal mol(-1) in model C, consistent with the experimental preference for the cis-isomer. However, in model A when both Glu214 and Glu146 are protonated, the stability is inverted with the trans-form being favored. The different protonation states of the titratable active-site residues Glu214 and Glu146 thus critically influence the manner in which the relative stability and degree of planarity of the cis- and trans-conformers vary with pH. Coupled with the known correlation of chromophore conformation with fluorescence efficiency, this work provides a detailed structural basis for the observed phenomenon that red fluorescent proteins such as HcRed, mKate and Rtms5 show bright fluorescence at high pH.

Quantum-chemical calculations of a long proton wire. Application of a harmonic model to analysis of the structure of an ionic defect in a water chain with an excess proton

Quantum-chemical calculations of molecular complexes (NH(3))(3)Zn(2+)...(H(2)O)(n)...NH(3) (C(n), n = 11, 16, 21, and 30) simulating a proton wire donor-water chain-acceptor were carried out. Earlier found periodicity in the length of the O-H bonds in water chain is explained within the framework of a one-component harmonic model. In complexes C(n), the geometry and electronic structure of ionic defect in water chain with an excess proton were studied. Calculations carried out at ab initio (B3LYP/6-31+G**) and semiempirical (PM3) levels of theory predict different patterns of distribution of the O-H bonds lengths and positive charge on the H-bond hydrogen atoms in the region of ionic defect. The obtained data show how a length of water chain and position of a protonated water link in the chain influence the ionic defect structure. To describe the observed structures of ionic defect, the harmonic model was used and the role of parameters of the H-bonded chain was investigated. The performed analysis explains different mechanisms (concerted and stepwise) of proton transfer along the H-bonded chain derived from ab initio and semiempirical calculation schemes.

Structural and Relaxation Effects in Proton Wire Energetics: Model Studies of the Green Fluorescent Protein Photocycle

We present the results of a systematic series of constrained minimum energy pathway calculations on ground state potential energy surfaces, for a cluster model of the proton chain transfer that mediates the photocycle of the green fluorescent protein, as well as for a model including the solvated protein environment. The calculations vary in terms of the types of modes that are assumed to be capable of relaxing in concert with the movement of the protons and the results demonstrate that the nature and extent of dynamical relaxation has a substantive impact on the activation energy for the proton transfer. We discuss the implications of this in terms of currently available dynamical models and chemical rate theories that might be brought to bear on the kinetics of this important example of proton chain transfer in a biological system.

Balance between ultrafast parallel reactions in the green fluorescent protein has a structural origin

The fluorescence photocycle of the green fluorescent protein is functionally dependent on the specific structural protein environment. A direct relationship between equilibrium protein side-chain conformation of glutamate 222 and reactivity is established, particularly the rate of ultrafast proton transfer reactions in the fluorescence photocycle. We show that parallel transformations in the photocycle have a structural origin, and we report on the vibrational properties of responsive amino acids on an ultrafast timescale. Blue excitation of GFP drives two parallel, excited-state deuteron transfer reactions with 10 ps and 75 ps time constants to the buried carboxylic acid side chain of glutamate 222 via a hydrogen-bonding network. Assignment of 1456 cm(-1) and 1441 cm(-1) modes to nu(sym) and assignment of 1564 cm(-1) and 1570 cm(-1) features to nu(asym) of E222 in the 10 ps and 75 ps components, respectively, was possible from the analysis of the transient absorption data of an E222D mutant and was consistent with photoselection measurements. In contrast to the wild-type, measurements of E222D can be described with only one difference spectrum, with the nu(sym) mode at 1435 cm(-1) and the nu(asym) mode at 1567 cm(-1), also correlating a large Deltanu(asym-sym) with slow excited-state proton transfer kinetics. Density Functional Theory calculations and published model compound and theoretical studies relate differences in Deltanu(asym-sym) to the strength and number of hydrogen-bonding interactions that are detected via equilibrium geometry and COO- stretching frequency differences of the carboxylate. The correlation of photocycle kinetics with side-chain conformation of the acceptor suggests that proton transfer from S205 to E222 controls the rate of the overall excited-state proton transfer process, which is consistent with recent theoretical predictions. Photoselection measurements show agreement for localized C=O vibrations of chromophore, Q69, and E222 with Density Functional Theory and ab initio calculations placed in the x-ray geometry and provide their vibrational response in the intermediates in the photocycle.

Variation of spectral, structural, and vibrational properties within the intrinsically fluorescent proteins family: a density functional study

We report an extensive (time-dependent) density functional study of the whole series of the chromophores within the Intrinsically Fluorescent Protein family, in the relevant conformations and protonation states. Over 30 structures are considered, including three newly discovered chromophores (zFP, Kaede, and the Orange variant of DsRed). Ab initio calculations on selected structures are also performed in order to assess the performances of TDDFT along the family. The use of a uniform scheme for all structures allows to establish relationships between the absorption energy and electrostatic, structural or vibrational properties of the chromophores.

Mechanistic aspects of proton chain transfer in the green fluorescent protein : Part II. A comparison of minimal quantum chemical models

In this paper we report the results of extensive quantum chemical reaction pathway calculations for the electronic ground state of several different cluster models that mimic the proton chain transfer path within the green fluorescent protein (GFP). Our principal objective is to establish the robustness with respect to variations in the model of our recent mechanistic inferences for the ground state proton chain transfer [S. Wang and S. C. Smith, J. Phys. Chem. B, 2006, 110, 5084]. Additionally, comparison of our ground state results with the excited state proton transfer (ESPT) study by Vendrell et al. [O. Vendrell, R. Gelabert, M. Moreno and J. M. Lluch, J. Am. Chem. Soc., 2006, 128, 3564] leads to the conclusion that the mechanism of proton chain transfer may be expected to be analogous in ground and excited states, principally because in both cases the loss of the chromophore's phenolic proton contributes strongly to the reaction coordinate only late in the reaction path.

Computational study of the absorption spectra of green fluorescent protein mutants

In this work, we present a theoretical study of the relationship between molecular structure and the red-shift in absorption spectra of S65G and S65T green fluorescent protein (GFP) mutants. To identify the effects of the protein environment, we combined results from molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics calculations to obtain structural properties, and applied time-dependent density functional theory to calculate the excitation energies. By using results from the MD simulations, we were able to provide a systematic analysis of the structural details that may effect the red-shift in the absorption spectra when taking into account temperature effects. Furthermore, a detailed study of hydrogen bonding during the MD simulations demonstrated differences between S65G and S65T, for example, regarding hydrogen bonding with Glu222. An analysis of the absorption spectra for different forms of the chromophore emphasized the dominance of the anionic forms in solution for the S65G and S65T GFP mutants.

Ground-State Structures and Vertical Excitations for the Kindling Fluorescent Protein asFP595

Geometry configurations of a large fraction of the kindling fluorescent protein asFP595 around the chromophore region were optimized by using the effective fragment potential quantum mechanical-molecular mechanical (QM/MM) method. The initial coordinates of heavy atoms were taken from the structure from the Protein Data Bank archive corresponding to the dark-adapted state of the Ala143 --> Gly mutant of asFP595. Optimization of geometry parameters was performed for all internal coordinates in the QM part composed of the chromophore unit and the side chains of His197, Glu215, and Arg92 as well as for positions of effective fragments constituting the MMpart. The structures corresponding to the anion trans, anion cis, and zwitterion trans moieties were considered among various alternatives for the chromophore unit inside the protein matrix. The QM/MM simulations show that the protein environment provides stabilization for the trans-zwitterion isomer compared to the gas-phase conditions. By using the multiconfigurational CASSCF and the time-dependent density functional theory calculations, we estimated positions of spectral bands corresponding to vertical S(0)-S(1) transitions. The results of simulations support the assumption that the dark state of asFP595 corresponds to the anionic or zwitterionic trans-conformation, while the kindled state corresponds to the anionic cis-conformation.

Cyan fluorescent protein: molecular dynamics, simulations, and electronic absorption spectrum

The dynamics and electronic absorption spectrum of enhanced cyan fluorescent protein (ECFP), a mutant of green fluorescent protein (GFP), have been studied by means of a 1 ns molecular dynamics (MD) simulation. The two X-ray conformations A' and B' of ECFP were considered. The chromophore was assumed to be neutral, and all titratable residues were taken in their standard protonation state at neutral pH. The protein was embedded in a box of water molecules (and counterions). The first result is that the two conformations A' and B' are found to be stable all along the simulation. Then, an analysis of the hydrogen-bond networks shows strong differences between the two conformations in the surroundings of the nitrogen atom of the indolic part of the chromophore. This is partly due to the imperfection in the beta barrel near the His148 residue, which allows the access of one solvent molecule inside the protein in conformation A'. Finally, quantum mechanical calculations of the electronic transition energies of the chromophore in the charge cloud of the protein and solvent water molecules were performed using the TDDFT method on 160 snapshots extracted every 5 ps of the MD trajectories. It is found that conformations A' and B' exhibit very similar spectra despite different H-bond networks involving the chromophore. This similarity is related to the weak charge transfer involved in the electronic transition and the weak electrostatic field created by ECFP near the chromophore, within the hypotheses made in the present simulation.

Potential Energy Landscape of the Photoinduced Multiple Proton-Transfer Process in the Green Fluorescent Protein: Classical Molecular Dynamics and Multiconfigurational Electronic Structure Calculations

The green fluorescent protein proton wire operating upon photoexcitation of the internally caged chromophore is investigated by means of classical molecular dynamics and multiconfigurational electronic structure calculations. The structure of the proton wire is studied for the solvated protein, showing that the wire is likely to be found in a configuration ready to operate as soon as the chromophore is photoexcited, and leading to a total of three proton translocations in the vicinity of the chromophore. Multiconfigurational CASSCF and CASPT2 calculations provide a detailed overview of the energy landscape of the proton wire for the ground electronic state S0, the photoactive 1pi pi* state, and the charge-transfer 1pi sigma* state. The results allow discussion of the operation of the wire in terms of the sequence of proton-transfer events and the participation of each electronic state.

Mechanistic Aspects of Proton Chain Transfer: A Computational Study for the Green Fluorescent Protein Chromophore

We explore several models for the ground-state proton chain transfer pathway between the green fluorescent protein chromophore and its surrounding protein matrix, with a view to elucidating mechanistic aspects of this process. We have computed quantum chemically the minimum energy pathways (MEPs) in the ground electronic state for one-, two-, and three-proton models of the chain transfer. There are no stable intermediates for our models, indicating that the proton chain transfer is likely to be a single, concerted kinetic step. However, despite the concerted nature of the overall energy profile, a more detailed analysis of the MEPs reveals clear evidence of sequential movement of protons in the chain. The ground-state proton chain transfer does not appear to be driven by the movement of the phenolic proton off the chromophore onto the neutral water bridge. Rather, this proton is the last of the three protons in the chain to move. We find that the first proton movement is from the bridging Ser205 moiety to the accepting Glu222 group. This is followed by the second proton moving from the bridging water to the Ser205--for our model this is where the barrier occurs. The phenolic proton on the chromophore is hence the last in the chain to move, transferring to a bridging "water" that already has substantial negative charge.