Reaction mechanism of photoinduced decarboxylation of the photoactivatable green fluorescent protein: an ONIOM(QM:MM) study

All-Atom Quantum Mechanical Calculation of the Second-Harmonic Generation of Fluorescent Proteins

Fluorescent proteins (FPs) are biotags of choice for second-harmonic imaging microscopy (SHIM). Because of their large size, computing their second-harmonic generation (SHG) response represents a great challenge for quantum chemistry. In this contribution, we propose a new all-atom quantum mechanics methodology to compute SHG of large systems. This is now possible because of two recent implementations: the tight-binding GFN2-xTB method to optimize geometries and a related version of the simplified time-dependent density functional theory (sTD-DFT-xTB) to evaluate quadratic response functions. In addition, a new dual-threshold configuration selection scheme is introduced to reduce the computational costs while retaining overall similar accuracy. This methodology was tested to evaluate the SHG of the proteins iLOV and bacteriorhodopsin (bR). In the case of bR, quantitative agreement with respect to experiment was reached for the out-of-resonance low-energy part of the β_HRS frequency dispersion. This work paves the way toward an accurate prediction of the SHG of large structures-a requirement for the design of new and improved SHIM biotags.

Metal-ion promiscuity of microbial enzyme DapE at its second metal-binding site

Metalloenzymes are ubiquitous in nature catalyzing a number of crucial biochemical processes in animal and plant kingdoms. For better adaptation to the relative abundance of different metal ions in different cellular fluids, many of these enzymes exhibit metal promiscuity. The microbial enzyme DapE, an essential enzyme for bacterial growth and survival and a potentially safe target for antibiotics, continues to show enzyme activity when the two zinc ions in its active site are replaced by other transition metal ions. The effect of metal-ion substitution at the second metal-binding site of DapE on its substrate affinity and catalytic efficiency is investigated by QM/MM treatment of the enzyme-substrate complex, by modelling the enzyme with Mn(II), Co(II), Ni(II), or Cu(II) ion in place of Zn(II) at its second metal-binding site, while retaining Zn(II) ion at the first metal-binding site. On the basis of substrate binding energy and activation energy barrier for the chemical catalysis, it is found that Zn-Mn DapE shows poor binding affinity as well as inefficient chemical catalysis. Although Zn-Cu and Zn-Ni DapEs show activation energy barriers comparable to that of wild-type Zn-Zn DapE, their weaker substrate affinity renders these mixed-metal enzymes less efficient. On the other hand, Zn-Co DapE is found to outperform the naturally occurring Zn-Zn DapE, both in terms of substrate affinity and chemical catalysis. The observed metal promiscuity may have played an important role in the survival of bacteria even in those cellular media where Zn ions are in limited supply. Metal nonspecificity in the catalysis of DapE enzyme allows bacteria to thrive in different cellular media.

Room-temperature photo-induced martensitic transformation in a protein crystal

Martensitic transformations are the first-order crystal-to-crystal phase transitions that occur mostly in materials such as steel, alloys and ceramics, thus having many technological applications. These phase transitions are rarely observed in molecular crystals and have not been detected in protein crystals. Reversibly switchable fluorescent proteins are widely used in biotechnology, including super-resolution molecular imaging, and hold promise as candidate biomaterials for future high-tech applications. Here, we report on a reversibly switchable fluorescent protein, Tetdron, whose crystals undergo a photo-induced martensitic transformation at room temperature. Room-temperature X-ray crystallography demonstrates that at equilibrium Tetdron chromophores are all in the trans configuration, with an ∼1:1 mixture of their protonated and deprotonated forms. Irradiation of a Tetdron crystal with 400 nm light induces a martensitic transformation, which results in Tetdron tetramerization at room temperature revealed by X-ray photocrystallography. Crystal and solution spectroscopic measurements provide evidence that the photo-induced martensitic phase transition is coupled with the chromophore deprotonation, but no trans-cis isomerization is detected in the structure of an irradiated crystal. It is hypothesized that protein dynamics assists in the light-induced proton transfer from the chromophore to the bulk solvent and in the ensuing martensitic phase transition. The unique properties of Tetdron may be useful in developing novel biomaterials for optogenetics, data storage and nanotechnology.

Subatomic resolution X-ray structures of green fluorescent protein

Green fluorescent protein (GFP) is a light-emitting protein that does not require a prosthetic group for its fluorescent activity. As such, GFP has become indispensable as a molecular tool in molecular biology. Nonetheless, there has been no subatomic elucidation of the GFP structure owing to the structural polymorphism around the chromophore. Here, subatomic resolution X-ray structures of GFP without the structural polymorphism are reported. The positions of H atoms, hydrogen-bonding network patterns and accurate geometric parameters were determined for the two protonated forms. Compared with previously determined crystal structures and theoretically optimized structures, the anionic chromophores of the structures represent the authentic resonance state of GFP. In addition, charge-density analysis based on atoms-in-molecules theory and noncovalent interaction analysis highlight weak but substantial interactions between the chromophore and the protein environment. Considered with the derived chemical indicators, the lone pair-π interactions between the chromophore and Thr62 should play a sufficient role in maintaining the electronic state of the chromophore. These results not only reveal the fine structural features that are critical to understanding the properties of GFP, but also highlight the limitations of current quantum-chemical calculations.

Modifications of the chromophore of Spinach aptamer based on QM:MM calculations

Spinach aptamer was developed as an RNA analog of the green fluorescent protein. The aptamer interacts with its ligand and modifies its electronic spectrum so that it fluoresces brightly at the wavelength of 501 nm. Song et al. investigated modifications of the ligand in their experimental study and found a molecule emitting at 523 nm upon creating a complex with the Spinach aptamer. The crystal structure of the aptamer in complex with its original ligand has been published, which enabled us to study the system computationally. In this article, we suggest several new modifications of the ligand that shift the emission maximum of the complex to even longer wavelengths. Our results are based on combined quantum mechanical/molecular mechanical calculations with DFT method used for geometry optimization and TD-DFT for calculations of absorption and emission energies.

Photoinduced Chemistry in Fluorescent Proteins: Curse or Blessing?

Photoinduced reactions play an important role in the photocycle of fluorescent proteins from the green fluorescent protein (GFP) family. Among such processes are photoisomerization, photooxidation/photoreduction, breaking and making of covalent bonds, and excited-state proton transfer (ESPT). Many of these transformations are initiated by electron transfer (ET). The quantum yields of these processes vary significantly, from nearly 1 for ESPT to 10^-4-10^-6 for ET. Importantly, even when quantum yields are relatively small, at the conditions of repeated illumination the overall effect is significant. Depending on the task at hand, fluorescent protein photochemistry is regarded either as an asset facilitating new applications or as a nuisance leading to the loss of optical output. The phenomena arising due to phototransformations include (i) large Stokes shifts, (ii) photoconversions, photoactivation, and photoswitching, (iii) phototoxicity, (iv) blinking, (v) permanent bleaching, and (vi) formation of long-lived intermediates. The focus of this review is on the most recent experimental and theoretical work on photoinduced transformations in fluorescent proteins. We also provide an overview of the photophysics of fluorescent proteins, highlighting the interplay between photochemistry and other channels (fluorescence, radiationless relaxation, and intersystem crossing). The similarities and differences with photochemical processes in other biological systems and in dyes are also discussed.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis–trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S₁ keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S₁ and S₀ potential energy surfaces relevant to enol–keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S₁ keto species from approaching the keto S₁/S₀ conical intersections so that most of trajectories are trapped in the S₁ keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S₁ excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis–trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Arginine 66 Controls Dark-State Formation in Green-to-Red Photoconvertible Fluorescent Proteins

Photoactivated localization microscopy (PALM) is a powerful technique to investigate cellular nanostructures quantitatively and dynamically. However, the use of PALM for molecular counting or single-particle tracking remains limited by the propensity of photoconvertible fluorescent protein markers (PCFPs) to repeatedly enter dark states. By designing the single mutants mEos2-A69T and Dendra2-T69A, we completely swapped the blinking behaviors of mEos2 and Dendra2, two popular PCFPs. We combined X-ray crystallography and single-molecule microscopy to show that blinking in mEos2 and Dendra2 is largely controlled by the orientation of arginine 66, a highly conserved residue in Anthozoan PCFPs. The Arg66 side-chain conformation affects the bleaching and the on-to-off transition quantum yields, as well as the fraction of molecules entering long-lived dark states, resulting in widely different apparent blinking behaviors that largely modulate the efficiency of current blinking correction procedures. The present work provides mechanistic insight into the complex photophysics of Anthozoan PCFPs and will facilitate future engineering of bright and low-blinking variants suitable for PALM.

Emission shaping in fluorescent proteins: role of electrostatics and π-stacking

We obtained the fluorescence spectrum of the GFP with trajectory simulations, and revealed the role of the protein sidechains in emission shifts.

A Light-Induced Reaction with Oxygen Leads to Chromophore Decomposition and Irreversible Photobleaching in GFP-Type Proteins

Photobleaching and photostability of proteins of the green fluorescent protein (GFP) family are crucially important for practical applications of these widely used biomarkers. On the basis of simulations, we propose a mechanism for irreversible bleaching in GFP-type proteins under intense light illumination. The key feature of the mechanism is a photoinduced reaction of the chromophore with molecular oxygen (O2) inside the protein barrel leading to the chromophore's decomposition. Using quantum mechanics/molecular mechanics (QM/MM) modeling we show that a model system comprising the protein-bound Chro(-) and O2 can be excited to an electronic state of the intermolecular charge-transfer (CT) character (Chro(•)···O2(-•)). Once in the CT state, the system undergoes a series of chemical reactions with low activation barriers resulting in the cleavage of the bridging bond between the phenolic and imidazolinone rings and disintegration of the chromophore.

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

Chromophore photoreduction in red fluorescent proteins is responsible for bleaching and phototoxicity

Red fluorescent proteins (RFPs) are indispensable tools for deep-tissue imaging, fluorescence resonance energy transfer applications, and super-resolution microscopy. Using time-resolved optical spectroscopy this study investigated photoinduced dynamics of three RFPs, KillerRed, mRFP, and DsRed. In all three RFPs, a new transient absorption intermediate was observed, which decays on a microsecond-millisecond time scale. This intermediate is characterized by red-shifted absorption at 1.68-1.72 eV (λmax = 720-740 nm). On the basis of electronic structure calculations, experimental evidence, and published literature, the chemical nature of the intermediate is assigned to an unusual open-shell dianionic chromophore (dianion-radical) formed via photoreduction. A doubly charged state that is not stable in the isolated (gas phase) chromophore is stabilized by the electrostatic field of the protein. Mechanistic implications for photobleaching, blinking, and phototoxicity are discussed.

The structural and energetic aspects of substrate binding and the mechanism of action of the DapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE) investigated using a hybrid QM/MM method

Substrate binding and the mechanism of action of the DapE-encodedN-succinyl-l,l-diaminopimelic acid desuccinylase (DapE).