The mechanism of cyclization in chromophore maturation of green fluorescent protein: a theoretical study

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

Power Density Titration of Reversible Photoisomerization of a Fluorescent Protein Chromophore in the Presence of Thermally Driven Barrier Crossing Shown by Quantitative Millisecond Serial Synchrotron X-ray Crystallography

We present millisecond quantitative serial X-ray crystallography at 1.7 Å resolution demonstrating precise optical control of reversible population transfer from Trans-Cis and Cis-Trans photoisomerization of a reversibly switchable fluorescent protein, rsKiiro. Quantitative results from the analysis of electron density differences, extrapolated structure factors, and occupancy refinements are shown to correspond to optical measurements of photoinduced population transfer and have sensitivity to a few percent in concentration differences. Millisecond time-resolved concentration differences are precisely and reversibly controlled through intense continuous wave laser illuminations at 405 and 473 nm for the Trans-to-Cis and Cis-to-Trans reactions, respectively, while the X-ray crystallographic measurement and laser illumination of the metastable Trans chromophore conformation causes partial thermally driven reconversion across a 91.5 kJ/mol thermal barrier from which a temperature jump between 112 and 128 K is extracted.

Combined Structural and Computational Study of the mRubyFT Fluorescent Timer Locked in Its Blue Form

The mRubyFT is a monomeric genetically encoded fluorescent timer based on the mRuby2 fluorescent protein, which is characterized by the complete maturation of the blue form with the subsequent conversion to the red one. It has higher brightness in mammalian cells and higher photostability compared with other fluorescent timers. A high-resolution structure is a known characteristic of the mRubyFT with the red form chromophore, but structural details of its blue form remain obscure. In order to obtain insight into this, we obtained an S148I variant of the mRubyFT (mRubyFT^S148I) with the blocked over time blue form of the chromophore. X-ray data at a 1.8 Å resolution allowed us to propose a chromophore conformation and its interactions with the neighboring residues. The imidazolidinone moiety of the chromophore is completely matured, being a conjugated π-system. The methine bridge is not oxidized in the blue form bringing flexibility to the phenolic moiety that manifests itself in poor electron density. Integration of these data with the results of molecular dynamic simulation disclosed that the OH group of the phenolic moiety forms a hydrogen bond with the side chain of the T163 residue. A detailed comparison of mRubyFT^S148I with other available structures of the blue form of fluorescent proteins, Blue102 and mTagBFP, revealed a number of characteristic differences. Molecular dynamic simulations with the combined quantum mechanic/molecular mechanic potentials demonstrated that the blue form exists in two protonation states, anion and zwitterion, both sharing enolate tautomeric forms of the C=C-O^- fragment. These two forms have similar excitation energies, as evaluated by calculations. Finally, excited state molecular dynamic simulations showed that excitation of the chromophore in both protonation states leads to the same anionic fluorescent state. The data obtained shed light on the structural features and spectral properties of the blue form of the mRubyFT timer.

Yellow and Orange Fluorescent Proteins with Tryptophan-based Chromophores

Rapid development of new microscopy techniques exposed the need for genetically encoded fluorescent tags with special properties. Recent works demonstrated the potential of fluorescent proteins with tryptophan-based chromophores. We applied rational design and random mutagenesis to the monomeric red fluorescent protein FusionRed and found two groups of mutants carrying a tryptophan-based chromophore: with yellow (535 nm) or orange (565 nm) emission. On the basis of the properties of proteins, a model synthetic chromophore, and a computational modeling, we concluded that the presence of a ketone-containing chromophore in different isomeric forms can explain the observed yellow and orange phenotypes.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

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.

Theoretical study on the chemical mechanism of enoyl-CoA hydratase and the form of inhibitor binding

Enoyl-CoA hydratase (ECH) catalyzes the second step in the vital β-oxidation pathway of fatty acid metabolism. This enzyme catalyzes the syn-addition of a water molecule across the double bond of 4-(N,N-dimethylamino) cinnamoyl-CoA (DAC-CoA). In this work, the reaction mechanisms of ECH were investigated using the density functional theory (DFT) methods. The different protonation states in which the important residues Glu164 and Glu144 are either neutral or ionized were considered. Four models of the active site were designed based on the X-ray crystal structure of the enzyme. The calculations gave strong support to the proposed mechanism and confirmed that both Glu164 and Glu144 are in a deprotonated state in the reaction mechanism of ECH. In addition, we constructed a model of the active site with the inhibitor acetoacetyl-CoA based on the crystal structure. Caomparison of the calculated energy barriers showed that binding of the keto-enol form of the inhibitor is more reasonable than that of the di-keto form in the inhibition process. Moreover, acetoacetyl-CoA was found to exhibit a keto-enol tautomerism when it acts as an inhibitor in the reaction. The present theoretical results indicated that both residues Glu164 and Glu144 are unprotonated in ECH with the substrate bound, while only Glu164 is unprotonated when the inhibitor binds ECH.

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.

Theoretical studies of chromophore maturation in the wild-type green fluorescent protein: ONIOM(DFT:MM) investigation of the mechanism of cyclization

The availability of a gene encoding green fluorescence immediately stimulates interest in the puzzle of autocatalytic formation of the green fluorescent protein (GFP) chromophore. Numerous experimental and theoretical studies have indicated that cyclization is the first and most important step in the maturation process of the GFP. In our previous paper based on cluster models [J. Phys. Chem. B2010, 114, 9698-9705], two possible mechanisms have been investigated with the conclusion that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 α-carbon. However, the impact of the protein environment on the reaction mechanism remains to be explored. In this paper, we investigated the two possible mechanisms with inclusion of protein environmental effects by using molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Our calculations reveal no hydrogen bonding network that would facilitate deprotonation of the amide nitrogen of Gly67, although it is the lower energy pathway in the cluster model system. Contrastingly, there is a hydrogen bonding network between Tyr66 α-carbon and Glu222, which is in good agreement with X-ray data. The ONIOM studies show that proton transfer from Tyr66 α-carbon to Glu222 is a long-distance charge transfer process. The charge distribution of the MM region has a significant perturbation to the wave function for the QM region, with the QM energy for the proton transfer product being increased under the influence of the electrostatic protein environment. The barrier for the rate-limiting step in cyclization is quite high, about 40.0 kcal/mol in the case of ONIOM-EE.