Theoretical Computer‐Aided Mutagenic Study on the Triple Green Fluorescent Protein Mutant S65T/H148D/Y145F

Green fluorescent protein (GFP) mutant S65T/H148D has been proposed to host a photocycle that involves an excited‐state proton transfer between the chromophore (Cro) and the Asp148 residue and takes place in less than 50 fs without a measurable kinetic isotope effect. It has been suggested that the interaction between the unsuspected Tyr145 residue and the chromophore is needed for the ultrafast sub‐50 fs rise in fluorescence. To verify this, we have performed a computer‐aided mutagenic study to introduce the additional mutation Y145F, which eliminates this interaction. By means of QM/MM molecular dynamics simulations and time‐dependent density functional theory studies, we have assessed the importance of the Cro–Tyr145 interaction and the solvation of Asp148 and shown that in the triple mutant S65T/H148D/Y145F a significant loss in the ultrafast rise of the Stokes‐shifted fluorescence should be expected.     

An alternate proton acceptor for excited-state proton transfer in green fluorescent protein: rewiring GFP

The neutral form of the chromophore in wild-type green fluorescent protein (wtGFP) undergoes excited-state proton transfer (ESPT) upon excitation, resulting in characteristic green (508 nm) fluorescence. This ESPT reaction involves a proton relay from the phenol hydroxyl of the chromophore to the ionized side chain of E222, and results in formation of the anionic chromophore in a protein environment optimized for the neutral species (the I* state). Reorientation or replacement of E222, as occurs in the S65T and E222Q GFP mutants, disables the ESPT reaction and results in loss of green emission following excitation of the neutral chromophore. Previously, it has been shown that the introduction of a second mutation (H148D) into S65T GFP allows the recovery of green emission, implying that ESPT is again possible. A similar recovery of green fluorescence is also observed for the E222Q/H148D mutant, suggesting that D148 is the proton acceptor for the ESPT reaction in both double mutants. The mechanism of fluorescence emission following excitation of the neutral chromophore in S65T/H148D and E222Q/H148D has been explored through the use of steady state and ultrafast time-resolved fluorescence and vibrational spectroscopy. The data are contrasted with those of the single mutant S65T GFP. Time-resolved fluorescence studies indicate very rapid (< 1 ps) formation of I* in the double mutants, followed by vibrational cooling on the picosecond time scale. The time-resolved IR difference spectra are markedly different to those of wtGFP or its anionic mutants. In particular, no spectral signatures are apparent in the picosecond IR difference spectra that would correspond to alteration in the ionization state of D148, leading to the proposal that a low-barrier hydrogen bond (LBHB) is present between the phenol hydroxyl of the chromophore and the side chain of D148, with different potential energy surfaces for the ground and excited states. This model is consistent with recent high-resolution structural data in which the distance between the donor and acceptor oxygen atoms is < or = 2.4 A. Importantly, these studies indicate that the hydrogen-bond network in wtGFP can be replaced by a single residue, an observation which, when fully explored, will add to our understanding of the various requirements for proton-transfer reactions within proteins.     

Structural forms of green fluorescent protein by quantum mechanics/molecular mechanics calculations

The molecular modeling of structural forms of the green fluorescent protein (GFP) with the Ser65Thr single-site mutation was performed by the quantum mechanics/molecular mechanics (QM/MM) method. Two model systems were constructed based on the crystallographic structure from the Protein Data Bank (PDB entry code 1EMA.) The model systems differ in the initial protonation state of the side chain of the amino acid residue Glu222 near the chromophore. The atomic coordinates of the protein macromolecule corresponding to the equilibrium geometric configurations were determined by total energy minimization using the QM/MM method within the density functional theory approximation PBE0/cc-pVDZ for the quantum subsystem that consists of the chromophore, a water molecule, and the side chains of Arg96, Glu222, and Ser205, and with the parameters of the AMBER force field for the molecular mechanics subsystem. In the analysis of the results, particular attention was given to the hydrogen bond redistribution in the chromophore-containing region of the protein caused by a change in the protonation state of the chromophore. The results obtained from the model containing the initially protonated side chain of Glu222 suggest a new interpretation of the photophysical processes in the green fluorescent protein.     

One‐Pot Bioconversion of l‐Arabinose to l‐Ribulose in an Enzymatic Cascade

This work reports the one‐pot enzymatic cascade that completely converts l ‐arabinose to l ‐ribulose using four reactions catalyzed by pyranose 2‐oxidase (P2O), xylose reductase, formate dehydrogenase, and catalase. As wild‐type P2O is specific for the oxidation of six‐carbon sugars, a pool of P2O variants was generated based on rational design to change the specificity of the enzyme towards the oxidation of l ‐arabinose at the C2‐position. The variant T169G was identified as the best candidate, and this had an approximately 40‐fold higher rate constant for the flavin reduction (sugar oxidation) step, as compared to the wild‐type enzyme. Computational calculations using quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) showed that this improvement is due to a decrease in the steric effects at the axial C4‐OH of l ‐arabinose, which allows a reduction in the distance between the C2‐H and flavin N5, facilitating hydride transfer and enabling flavin reduction.     

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.     

Probing the range of applicability of structure‐ and energy‐adjusted QM/MM link bonds

The hydrogen‐capping method is one of the most popular and widely used coupling‐schemes for quantum mechanics/molecular mechanics (QM/MM)‐molecular dynamics simulations of macromolecular systems. This is mostly due to the fact that it is fairly convenient to implement and parametrize, thus providing an excellent compromise between accuracy and computational effort. In this work, a viable and straight‐forward approach to optimize the placing of the link atom on a suitable distance ratio between the frontier atoms is discussed. To further increase the accuracy, instead of global parameters for all amino acids, different parameter sets for each type of amino acid are derived. The dependency of the link bond parameters on the chemical environment and the used QM‐method is probed to assess the range of applicability of the parametrization. Suitable sets of parameters for RI‐MP2, B3LYP, (RI)‐B3LYP‐D3, and RI‐BLYP‐D3 at triple‐zeta level for all relevant proteinogenic amino acids are presented. Furthermore, the scope and range of the perturbation, stemming from the introduction of link bonds is evaluated through application of the presented QM/MM scheme in calculations of the active site of 15S‐lipoxygenase. © 2015 Wiley Periodicals, Inc.     

Water solvent effects using continuum and discrete models: The nitromethane molecule,CH3NO2

The first three valence transitions of the two nitromethane conformers (CH₃NO₂) are two dark n → π* transitions and a very intense π → π* transition. In this work, these transitions in gas‐phase and solvated in water of both conformers were investigated theoretically. The polarizable continuum model (PCM), two conductor‐like screening (COSMO) models, and the discrete sequential quantum mechanics/molecular mechanics (S‐QM/MM) method were used to describe the solvation effect on the electronic spectra. Time dependent density functional theory (TDDFT), configuration interaction including all single substitutions and perturbed double excitations (CIS(D)), the symmetry‐adapted‐cluster CI (SAC‐CI), the multistate complete active space second order perturbation theory (CASPT2), and the algebraic‐diagrammatic construction (ADC(2)) electronic structure methods were used. Gas‐phase CASPT2, SAC‐CI, and ADC(2) results are in very good agreement with published experimental and theoretical spectra. Among the continuum models, PCM combined either with CASPT2, SAC‐CI, or B3LYP provided good agreement with available experimental data. COSMO combined with ADC(2) described the overall trends of the transition energy shifts. The effect of increasing the number of explicit water molecules in the S‐QM/MM approach was discussed and the formation of hydrogen bonds was clearly established. By including explicitly 24 water molecules corresponding to the complete first solvation shell in the S‐QM/MM approach, the ADC(2) method gives more accurate results as compared to the TDDFT approach and with similar computational demands. The ADC(2) with S‐QM/MM model is, therefore, the best compromise for accurate solvent calculations in a polar environment. © 2015 Wiley Periodicals, Inc.     

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.     

Insights into the different dioxygen activation pathways of methane and toluene monooxygenase hydroxylases

The methane and toluene monooxygenase hydroxylases (MMOH and TMOH, respectively) have almost identical active sites, yet the physical and chemical properties of their oxygenated intermediates, designated P*, H(peroxo), Q, and Q* in MMOH and ToMOH(peroxo) in a subclass of TMOH, ToMOH, are substantially different. We review and compare the structural differences in the vicinity of the active sites of these enzymes and discuss which changes could give rise to the different behavior of H(peroxo) and Q. In particular, analysis of multiple crystal structures reveals that T213 in MMOH and the analogous T201 in TMOH, located in the immediate vicinity of the active site, have different rotatory configurations. We study the rotational energy profiles of these threonine residues with the use of molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) computational methods and put forward a hypothesis according to which T213 and T201 play an important role in the formation of different types of peroxodiiron(III) species in MMOH and ToMOH. The hypothesis is indirectly supported by the QM/MM calculations of the peroxodiiron(III) models of ToMOH and the theoretically computed Mo?ssbauer spectra. It also helps explain the formation of two distinct peroxodiiron(III) species in the T201S mutant of ToMOH. Additionally, a role for the ToMOD regulatory protein, which is essential for intermediate formation and protein functioning in the ToMO system, is advanced. We find that the low quadrupole splitting parameter in the Mo?ssbauer spectrum observed for a ToMOH(peroxo) intermediate can be explained by protonation of the peroxo moiety, possibly stabilized by the T201 residue. Finally, similarities between the oxygen activation mechanisms of the monooxygenases and cytochrome P450 are discussed.     

Ultrafast electronic and vibrational dynamics of stabilized A state mutants of the green fluorescent protein (GFP): Snipping the proton wire

Two blue absorbing and emitting mutants (S65G/T203V/E222Q and S65T at pH 5.5) of the green fluorescent protein (GFP) have been investigated through ultrafast time resolved infra-red (TRIR) and fluorescence spectroscopy. In these mutants, in which the excited state proton transfer reaction observed in wild-type GFP has been blocked, the photophysics are dominated by the neutral A state. It was found that the A* excited state lifetime is short, indicating that it is relatively less stabilised in the protein matrix than the anionic form. However, the lifetime of the A state can be increased through modifications to the protein structure. The TRIR spectra show that a large shifts in protein vibrational modes on excitation of the A state occurs in both these GFP mutants. This is ascribed to a change in H-bonding interactions between the protein matrix and the excited state.     

Supramolecular Nanoscaled Self‐Assembly of Bifunctional Cyclodextrin and Trifunctional Melamine Precursors

A new type of guest–core supramolecular networks via inclusion complexation of nanoscaled building blocks such as bifunctional cyclodextrin (CD) derivatives and trifunctional melamine derivatives were prepared. By using AFM and an adoption of the cryo‐TEM technique under high acceleration voltage the nanoscale supramolecular network structure, nexus units, and CD molecules could be visualized. In addition to the 2‐D ¹H NMR rotating frame Overhauser effect spectroscopy (ROESY) experiments, theoretical studies on the molecular docking of the CDs and the melamine derivative have been conducted to elucidate the thermodynamic properties by the two‐layered integrated molecular orbital and molecular mechanics (ONIOM) method, which combines both quantum mechanics and molecular mechanics (QM/MM) calculations.

     

Fragment‐based quantum mechanical calculation of protein–protein binding affinities

The electrostatically embedded generalized molecular fractionation with conjugate caps (EE‐GMFCC) method has been successfully utilized for efficient linear‐scaling quantum mechanical (QM) calculation of protein energies. In this work, we applied the EE‐GMFCC method for calculation of binding affinity of Endonuclease colicin–immunity protein complex. The binding free energy changes between the wild‐type and mutants of the complex calculated by EE‐GMFCC are in good agreement with experimental results. The correlation coefficient (R) between the predicted binding energy changes and experimental values is 0.906 at the B3LYP/6‐31G*‐D level, based on the snapshot whose binding affinity is closest to the average result from the molecular mechanics/Poisson–Boltzmann surface area (MM/PBSA) calculation. The inclusion of the QM effects is important for accurate prediction of protein–protein binding affinities. Moreover, the self‐consistent calculation of PB solvation energy is required for accurate calculations of protein–protein binding free energies. This study demonstrates that the EE‐GMFCC method is capable of providing reliable prediction of relative binding affinities for protein–protein complexes. © 2018 Wiley Periodicals, Inc.     

Mitochondria from cultured cells derived from normal and thiamine-responsive megaloblastic anemia individuals efficiently import thiamine diphosphate

Background

Within the family of green fluorescent protein (GFP) homologs, one can mark two main groups, specifically, fluorescent proteins (FPs) and non-fluorescent or chromoproteins (CPs). Structural background of differences between FPs and CPs are poorly understood to date.

Results

Here, we applied site-directed and random mutagenesis in order to to transform CP into FP and vice versa. A purple chromoprotein asCP (asFP595) from Anemonia sulcata and a red fluorescent protein DsRed from Discosoma sp. were selected as representatives of CPs and FPs, respectively. For asCP, some substitutions at positions 148 and 165 (numbering in accordance to GFP) were found to dramatically increase quantum yield of red fluorescence. For DsRed, substitutions at positions 148, 165, 167, and 203 significantly decreased fluorescence intensity, so that the spectral characteristics of these mutants became more close to those of CPs. Finally, a practically non-fluorescent mutant DsRed-NF was generated. This mutant carried four amino acid substitutions, specifically, S148C, I165N, K167M, and S203A. DsRed-NF possessed a high extinction coefficient and an extremely low quantum yield (< 0.001). These spectral characteristics allow one to regard DsRed-NF as a true chromoprotein.

Conclusions

We located a novel point in asCP sequence (position 165) mutations at which can result in red fluorescence appearance. Probably, this finding could be applied onto other CPs to generate red and far-red fluorescent mutants. A possibility to transform an FP into CP was demonstrated. Key role of residues adjacent to chromophore's phenolic ring in fluorescent/non-fluorescent states determination was revealed.     

Calculation of wave‐functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. II. Application of the local basis equation

The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009 , 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self‐consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree–Fock–Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self‐consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave‐function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem. © 2013 Wiley Periodicals, Inc.     

Hidden electronic excited state of enhanced green fluorescent protein

Precise two-photon absorption spectra of the green fluorescent protein (GFP) and the mutants sapphire-GFP (T203I) and enhanced GFP (S65T/F64L), as well as a model compound for the chromophore, 4'-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) were measured by multiplex two-photon absorption spectroscopy. The observed TPA bands of the anionic forms of enhanced GFP and HBDI were significantly shifted to the higher energy compared with the lowest-energy bands in one-photon absorption spectra. This result indicated the existence of a hidden electronic excited state in the vicinity of the lowest excited singlet (S1) state of the anionic form of the GFP chromophore, which is the origin of the blue shift of the two-photon absorption spectra as well as two-photon fluorescence excitation spectra.     

Assessing the validity of DLPNO-CCSD(T) in the calculation of activation and reaction energies of ubiquitous enzymatic reactions

The domain-based local pair natural orbital coupled-cluster with single, double, and perturbative triples excitation (DLPNO-CCSD(T)) method was employed to portray the activation and reaction energies of four ubiquitous enzymatic reactions, and its performance was confronted to CCSD(T)/complete basis set (CBS) to assess its accuracy and robustness in this specific field. The DLPNO-CCSD(T) results were also confronted to those of a set of density functionals (DFs) to understand the benefit of implementing this technique in enzymatic quantum mechanics/molecular mechanics (QM/MM) calculations as a second QM component, which is often treated with DF theory (DFT). On average, the DLPNO-CCSD(T)/aug-cc-pVTZ results were 0.51 kcal·mol⁻¹ apart from the canonic CCSD(T)/CBS, without noticeable biases toward any of the reactions under study. All DFs fell short to the DLPNO-CCSD(T), both in terms of accuracy and robustness, which suggests that this method is advantageous to characterize enzymatic reactions and that its use in QM/MM calculations, either alone or in conjugation with DFT, in a two-region QM layer (DLPNO-CCSD(T):DFT), should enhance the quality and faithfulness of the results.     

Comment on “A stationary‐wave model of enzyme catalysis” by Carlo Canepa

Energy barriers for enzyme‐catalyzed reactions calculated with quantum mechanics/molecular mechanics (QM/MM) and empirical valence bond (EVB) methods can be in excellent agreement with activation energies derived from experiments, supporting the applicability of transition state theory for enzymic reactions. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Whether proton transition to the triphosphate tail of ATP occurs at protein kinase environment: A Car‐Parrinello ab initio molecular dynamics study

Protonation state of the triphosphate tail of ATP (adenosine triphosphate) in protein environment is a fundamental issue, which has significant impact on the mechanism investigation of biochemical processes with ATP involved. Proton transition from surroundings (water molecule coordinating to magnesium, H_W; amino group of Lys, H_L) to the ATP tail in the catalytic core of protein kinase found recently disproved the commonly accepted deprotonation state of ATP tail. In this account, Car‐Parrinello ab initio molecular dynamics (CP‐AIMD) method has been employed to examine whether the proton transition occurs. To provide a comparison basis for the dynamics simulations, static quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM) calculations have also been carried out. Consistent results have been obtained that complete transition of hydrogen from the surroundings to the triphosphate tail of ATP is not allowed. The most dominant conformations correspond to the ones with H_W bonding to O(W) and H‐bonding to O(ATP), [O(W)‐H_W···O(ATP)], H_L bonding to N(Lys) and H‐bonding to O(ATP), [N(Lys)‐H_L···O(ATP)]. Metastable structures with one hydrogen atom bonding with two heavy atoms (hydrogen acceptors) were also located by our dynamic simulations. This bonding mode can satisfy the hungering for hydrogen of the two heavy atoms simultaneously. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008