Histidine-assisted reduction of arylnitrenes upon photo-activation of phenyl azide chromophores in GFP-like fluorescent proteins.

The photochemically active sites of the proteins sfGFP66azF and Venus66azF, members of the green fluorescent protein (GFP) family, contain a non-canonical amino acid residue p-azidophenylalanine (azF) instead of Tyr66. The light-induced decomposition of azF at these sites leads to the formation of reactive arylnitrene (nF) intermediates followed by the formation of phenylamine-containing chromophores. We report the first study of the reaction mechanism of the reduction of the arylnitrene intermediates in sfGFP66nF and Venus66nF using molecular modeling methods. The Gibbs energy profiles for the elementary steps of the chemical reaction in sfGFP66nF are computed using molecular dynamics simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. Structures and energies along the reaction pathway in Venus66nF are evaluated using a QM/MM approach. According to the results of the simulations, arylnitrene reduction is coupled with oxidation of the histidine side chain on the His148 residue located near the chromophore.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms.

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

QM/MM Modeling of the Flavin Functionalization in the RutA Monooxygenase.

Oxygenase activity of the flavin-dependent enzyme RutA is commonly associated with the formation of flavin-oxygen adducts in the enzyme active site. We report the results of quantum mechanics/molecular mechanics (QM/MM) modeling of possible reaction pathways initiated by various triplet state complexes of the molecular oxygen with the reduced flavin mononucleotide (FMN) formed in the protein cavities. According to the calculation results, these triplet-state flavin-oxygen complexes can be located at both re-side and si-side of the isoalloxazine ring of flavin. In both cases, the dioxygen moiety is activated by electron transfer from FMN, stimulating the attack of the arising reactive oxygen species at the C4a, N5, C6, and C8 positions in the isoalloxazine ring after the switch to the singlet state potential energy surface. The reaction pathways lead to the C(4a)-peroxide, N(5)-oxide, or C(6)-hydroperoxide covalent adducts or directly to the oxidized flavin, depending on the initial position of the oxygen molecule in the protein cavities.

Modeling Light-Induced Chromophore Hydration in the Reversibly Photoswitchable Fluorescent Protein Dreiklang.

We report the results of a computational study of the mechanism of the light-induced chemical reaction of chromophore hydration in the fluorescent protein Dreiklang, responsible for its switching from the fluorescent ON-state to the dark OFF-state. We explore the relief of the charge-transfer excited-state potential energy surface in the ON-state to locate minimum energy conical intersection points with the ground-state energy surface. Simulations of the further evolution of model systems allow us to characterize the ground-state reaction intermediate tentatively suggested in the femtosecond studies of the light-induced dynamics in Dreiklang and finally to arrive at the reaction product. The obtained results clarify the details of the photoswitching mechanism in Dreiklang, which is governed by the chemical modification of its chromophore.

Modeling Spectral Tuning in Red Fluorescent Proteins Using the Dipole Moment Variation upon Excitation.

We describe a model for spectral tuning in red fluorescent proteins (RFPs) based on the relation between an electronic structure descriptor, the dipole moment variation upon excitation (DMV), and the excitation energy of a protein. This approach aims to overcome the problem of accurate prediction of excitation energies in RFPs, which span a very narrow window of band maxima. The latter roughly corresponds to the energy range of 0.1 eV, which is comparable with typical errors in calculations of the excitation energy by conventional quantum chemistry methods. In this work, we demonstrate a strong quantitative correlation between DMV values, obtained computationally with modest efforts, and excitation energies ΔEex at the experimental excitation band maxima for a series of RFPs with bands between 570 and 605 nm. Protein models are constructed by motifs of the relevant crystal structures, and atomic coordinates are optimized in quantum mechanics/molecular mechanics (QM/MM) calculations with QM-subsystems composed of large chromophore-containing regions. DMV values are evaluated with the electron density computed at the time-dependent density functional theory (TDDFT) level using several functionals and basis sets. We show that the results obtained with the CAM-B3LYP, BHHLYP, and M06-2X functionals demonstrate favorable correlations between DMV and ΔEex with the mean absolute error less than 0.01 eV. Taking into account the solid theoretical grounds of the relation between the DMV and the excitation energy in fluorescent proteins, the described modeling strategy presents a rational tool for spectral tuning in these efficient markers for in vivo imaging.

Light-Induced Change of Arginine Conformation Modulates the Rate of Adenosine Triphosphate to Cyclic Adenosine Monophosphate Conversion in the Optogenetic System Containing Photoactivated Adenylyl Cyclase.

We report the first computational characterization of an optogenetic system composed of two photosensing BLUF (blue light sensor using flavin adenine dinucleotide) domains and two catalytic adenylyl cyclase (AC) domains. Conversion of adenosine triphosphate (ATP) to the reaction products, cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi), catalyzed by ACs initiated by excitation in photosensing domains has emerged in the focus of modern optogenetic applications because of the request in photoregulated enzymes that modulate cellular concentrations of signaling messengers. The photoactivated AC from the soil bacterium Beggiatoa sp. (bPAC) is an important model showing a considerable increase in the ATP to cAMP conversion rate in the catalytic domain after the illumination of the BLUF domain. The 1 µs classical molecular dynamics simulations reveal that the activation of the BLUF domain leading to tautomerization of Gln49 in the chromophore-binding pocket results in switching of the position of the side chain of Arg278 in the active site of AC. Allosteric signal transmission pathways between Gln49 from BLUF and Arg278 from AC were revealed by the dynamical network analysis. The Gibbs energy profiles of the ATP → cAMP + PPi reaction computed using QM(DFT(ωB97X-D3/6-31G**))/MM(CHARMM) molecular dynamics simulations for both Arg278 conformations in AC clarify the reaction mechanism. In the light-activated system, the corresponding arginine conformation stabilizes the pentacoordinated phosphorus of the α-phosphate group in the transition state, thus lowering the activation energy. Simulations of the bPAC system with the Tyr7Phe replacement in the BLUF demonstrate occurrence of both arginine conformations in an equal ratio, explaining the experimentally observed intermediate catalytic activity of the bPAC-Y7F variant as compared with the dark and light states of the wild-type bPAC.

Modeling photophysical properties of the bacteriophytochrome-based fluorescent protein IFP1.4.

An enhanced interest in the phytochrome-based fluorescent proteins is explained by their ability to absorb and emit light in the far-red and infra-red regions particularly suitable for bioimaging. The fluorescent protein IFP1.4 was engineered from the chromophore-binding domain of a bacteriophytochrome in attempts to increase the fluorescence quantum yield. We report the results of simulations of structures in the ground S0 and excited S1 electronic states of IFP1.4 using the methods of quantum chemistry and quantum mechanics/molecular mechanics. We construct different protonation states of the biliverdin (BV) chromophore in the red-absorbing form of the protein by moving protons from the BV pyrrole rings to a suitable acceptor within the system and show that these structures are close in energy but differ by absorption bands. For the first time, we report structures of the minimum energy conical intersection points S1/S0 on the energy surfaces of BV in the protein environment and describe their connection to the local minima in the excited S1 state. These simulations allow us to characterize the deactivation routes in IFP1.4.

Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase.

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the "oxyanion hole", and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials.

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Tuning Electrostatic Gating of Semiconducting Carbon Nanotubes by Controlling Protein Orientation in Biosensing Devices.

The ability to detect proteins through gating conductance by their unique surface electrostatic signature holds great potential for improving biosensing sensitivity and precision. Two challenges are: (1) defining the electrostatic surface of the incoming ligand protein presented to the conductive surface; (2) bridging the Debye gap to generate a measurable response. Herein, we report the construction of nanoscale protein-based sensing devices designed to present proteins in defined orientations; this allowed us to control the local electrostatic surface presented within the Debye length, and thus modulate the conductance gating effect upon binding incoming protein targets. Using a ß-lactamase binding protein (BLIP2) as the capture protein attached to carbon nanotube field effect transistors in different defined orientations. Device conductance had influence on binding TEM-1, an important ß-lactamase involved in antimicrobial resistance (AMR). Conductance increased or decreased depending on TEM-1 presenting either negative or positive local charge patches, demonstrating that local electrostatic properties, as opposed to protein net charge, act as the key driving force for electrostatic gating. This, in turn can, improve our ability to tune the gating of electrical biosensors toward optimized detection, including for AMR as outlined herein.

Proof of concept for poor inhibitor binding and efficient formation of covalent adducts of KRASG12C and ARS compounds.

The use of selective covalent inhibitors with low binding affinity and high reactivity with the target enzyme is a promising way to solve a long-standing problem of the "undruggable" RAS-like proteins. Specifically, compounds of the ARS family that prevent the activation of the GDP-bound G12C mutant of Kirsten RAS (KRAS) are in the focus of recent experimental research. We report the first computational characterization of the entire reaction mechanism of the covalent binding of ARS-853 to the KRASG12C·GDP complex. The application of molecular dynamics, molecular docking and quantum mechanics/molecular mechanics approaches allowed us to model the inhibitor binding to the protein and the chemical reaction of ARS-853 with Cys12 in the enzyme binding site. We estimated a full set of kinetic constants and carried out numerical kinetic analysis of the process. Thus, we were able to compare directly the physicochemical parameters of the reaction obtained in silico and the macroscopic parameters observed in experimental studies. From our computational results, we explain the observed unusual dependence of the rate constant of covalent complex formation, kobs, on the ARS concentration. The latter depends both on the non-covalent binding step with the equilibrium constant, Ki, and on the rate constant of covalent adduct formation, kinact. The calculated ratio kinact/Ki = 213 M-1 s-1 reproduces the corresponding experimental value of 250 ± 30 M-1 s-1 for the interaction of ARS-853 with KRASG12C. Electron density analysis in the reactive region demonstrates that covalent bond formation occurs efficiently according to the Michael addition mechanism, which assumes the activation of the C[double bond, length as m-dash]C bond of ARS-853 by a water molecule and Lys16 in the binding site of KRASG12C. We also refine the kinact and Ki constants of the ARS-107 compound, which shares common features with ARS-853, and show that the decrease in the kinact/Ki ratio in the case of ARS-107 is explained by changes in both Ki and kinact constants.

Dipole Moment Variation Clears Up Electronic Excitations in the π-Stacked Complexes of Fluorescent Protein Chromophores.

We propose a quantitative structure-property relationship (QSPR) model for prediction of spectral tuning in cyan, green, orange, and red fluorescent proteins, which are engineered by motifs of the green fluorescent protein. Protein variants, in which their chromophores are involved in the π-stacking interaction with amino acid residues tyrosine, phenylalanine, and histidine, are prospective markers useful in bioimaging and super-resolution microscopy. In this work, we constructed training sets of the π-stacked complexes of four fluorescent protein chromophores (of the green, orange, red, and cyan series) with various substituted benzenes and imidazoles and tested the use of dipole moment variation upon excitation (DMV) as a descriptor to evaluate the vertical excitation energies in these systems. To validate this approach, we computed and analyzed electron density distributions of the π-stacked complexes and correlated the QSPR predictions with the reference values of the transition energies obtained using the high-level ab initio quantum chemistry methods. According to our results, the use of the DMV descriptor allows one to predict excitation energies in the π-stacked complexes with errors not exceeding 0.1 eV, which makes this model a practically useful tool in the development of efficient fluorescent markers for in vivo imaging.

Dynamical properties of enzyme-substrate complexes disclose substrate specificity of the SARS-CoV-2 main protease as characterized by the electron density descriptors.

A dynamical approach is proposed to discriminate between reactive (rES) and nonreactive (nES) enzyme-substrate complexes taking the SARS-CoV-2 main protease (Mpro) as an important example. Molecular dynamics simulations with the quantum mechanics/molecular mechanics potentials (QM(DFT)/MM-MD) followed by the electron density analysis are employed to evaluate geometry and electronic properties of the enzyme with different substrates along MD trajectories. We demonstrate that mapping the Laplacian of the electron density and the electron localization function provides easily visible images of the substrate activation that allow one to distinguish rES and nES. The computed fractions of reactive enzyme-substrate complexes along MD trajectories well correlate with the findings of recent experimental studies on the substrate specificity of Mpro. The results of our simulations demonstrate the role of the theory level used in QM subsystems for a proper description of the nucleophilic attack of the catalytic cysteine residue in Mpro. The activation of the carbonyl group of a substrate is correctly characterized with the hybrid DFT functional PBE0, whereas the use of a GGA-type PBE functional, that lacks the admixture of the Hartree-Fock exchange fails to describe activation.

Novel flavin-based fluorescent proteins with red-shifted emission bands: a computational study.

The iLOV protein is a promising member of the class of flavin mononucleotide (FMN) based fluorescent proteins (FbFPs). It is becoming a popular tool for bioanalytical applications and bioimaging as a competitor of the well-known green fluorescent protein and its analogues. The main limitation of FbFPs is that all the members have close values of their absorption and emission band maxima. Therefore the upcoming challenge is to introduce novel variants of FbFPs to extend their color palette. We report the results of computational studies of iLOV variants, introducing point mutations and chromophore analogues. We found that point mutations of the apoprotein and substitution of FMN with either 8-amino-FMN or 8-methylamino-FMN lead to the red shift of emission bands up to 100 nm. Substitution with 1-deaza-FMN and the point mutations of the apoprotein result in a set of novel fluorescent proteins with emission bands in the "transparent" window where light readily penetrates through mammalian tissues. Newly suggested FbFPs can be used for multicolor imaging and also as components of FRET pairs.

Diversity of mechanisms in Ras-GAP catalysis of guanosine triphosphate hydrolysis revealed by molecular modeling.

The mechanism of the deceptively simple reaction of guanosine triphosphate (GTP) hydrolysis catalyzed by the cellular protein Ras in complex with the activating protein GAP is an important issue because of the significance of this reaction in cancer research. We show that molecular modeling of GTP hydrolysis in the Ras-GAP active site reveals a diversity of mechanisms of the intrinsic chemical reaction depending on molecular groups at position 61 in Ras occupied by glutamine in the wild-type enzyme. First, a comparison of reaction energy profiles computed at the quantum mechanics/molecular mechanics (QM/MM) level shows that an assignment of the Gln61 side chain in the wild-type Ras either to QM or to MM parts leads to different scenarios corresponding to the glutamine-assisted or the substrate-assisted mechanisms. Second, replacement of Gln61 by the nitro-analog of glutamine (NGln) or by Glu, applied in experimental studies, results in two more scenarios featuring the so-called two-water and the concerted-type mechanisms. The glutamine-assisted mechanism in the wild-type Ras-GAP, in which the conserved Gln61 plays a decisive role, switching between the amide and imide tautomer forms, is consistent with the known experimental results of structural, kinetic and spectroscopy studies. The results emphasize the role of the Ras residue Gln61 in Ras-GAP catalysis and explain the retained catalytic activity of the Ras-GAP complex towards GTP hydrolysis in the Gln61NGln and Gln61Glu mutants of Ras.

Allosteric Control of N-Acetyl-Aspartate Hydrolysis by the Y231C and F295S Mutants of Human Aspartoacylase.

We present the results of molecular modeling of conformational changes in the Y231C and F295S mutants of human aspartoacylase (hAsp), which allow us to propose a mechanism of allosteric regulation of enzyme activity of these protein variants. The hAsp enzyme hydrolyzes one of the most abundant amino acid derivatives in the brain, N-acetyl-aspartate. It is important to understand the reasons for diminishing activity of the mutated enzymes, which is crucial for Canavan disease patients bearing the mutated gene. We explore a model which suggests operation of hAsp in the dimer form with two dynamically inequivalent subunits. Large-scale molecular dynamics simulations reveal that the replacements Y231C and F295S at the periphery of the protein shift the equilibrium between hAsp conformations with the open and closed gates to the enzyme active site buried inside the protein. Application of the dynamical network analysis and the Markov state model approach allows us to strengthen this conclusion and provide a detailed description of dynamically induced structural changes of the protein. The decreased availability of the active site for substrate molecules in the mutated enzymes explains their diminishing activity observed in clinical experiments.

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.

Competition between two cysteines in covalent binding of biliverdin to phytochrome domains.

In this work, we disclose a mechanism of competing chemical reactions of protein assembly for a bacterial phytochrome using modern methods of molecular modeling. The recently designed variant of a near-infrared fluorescent protein miRFP670 shows novel and unexpected features of covalent binding of the biliverdin chromophore to cysteine residues in the phytochrome domains GAF and PAS. Upon protein assembly, biliverdin reacts either with a cysteine from GAF, or with two cysteines, one from GAF and one from PAS. We characterize computationally a model structure of the noncovalently bound biliverdin molecule inside the protein cleft of miRFP670 and model reactions of the covalent binding. Both cysteines, Cys20 (PAS) and Cys253 (GAF), are located close to the electrophile C32 atom of biliverdin and can act as nucleophiles. The nucleophilic attack of Cys253 from the GAF domain results in a single C-S bond formation with an activation energy of 16 kcal mol-1. Another pathway, leading to the biliverdin adduct with two C-S bonds, is characterized by lower energy barriers, less than 11 kcal mol-1. Competition between these reaction pathways explains the experimentally obtained mixture of both adducts. On the basis of our first simulations of covalent BV binding to the phytochrome domains, we propose an approach of a direct experimental validation of the reaction mechanisms using IR vibrational spectroscopy.

Amide-imide tautomerization in the glutamine side chain in enzymatic and photochemical reactions in proteins.

Amide-imide tautomerization presents a pervasive class of chemical transformations in organic chemistry of natural compounds. In this Perspective, we describe two distinctively different protein systems, in which the amide-imide tautomerization in the glutamine side chain takes place in enzymatic or photochemical reactions. First, hydrolysis of guanosine triphosphate (GTP) catalyzed by the Ras-GAP protein complex suggests the occurrence of the imide tautomer of glutamine in reaction intermediates. Second, photoexcitation of flavin-binding protein domains (BLUFs) initiates a chain of reactions in the chromophore-binding pocket, including amide-imide tautomerization of glutamine. Mechanisms of these reactions at the atomic level have been revealed in quantum mechanics/molecular mechanics (QM/MM) simulations. To reinforce conclusions on the critical role of amide-imide tautomerization of glutamine in these reactions we describe results of new quantum chemistry and QM/MM calculations for relevant molecular model systems. We reexamine results of the recent IR spectroscopy studies of BLUF domains, which provide experimental evidences of Gln tautomerization in proteins. We also propose to validate the glutamine-assisted mechanism of enzymatic GTP hydrolysis by using IR spectroscopy in a proper range of wavenumbers.

Molecular Modeling Clarifies the Mechanism of Chromophore Maturation in the Green Fluorescent Protein.

We report the first complete theoretical description of the chain of elementary reactions resulting in chromophore maturation in the green fluorescent protein (GFP). All reaction steps including cyclization, dehydration, and oxidation are characterized at the uniform quantum mechanics/molecular mechanics (QM/MM) computational level using density functional theory in quantum subsystems. Starting from a structure of the wild-type protein with the noncyclized Ser65-Tyr66-Gly67 tripeptide, we modeled cyclization and dehydration reactions. We then added molecular oxygen to the system and modeled the oxidation reaction resulting in the mature protein-bound chromophore. Computationally derived structures of the reaction product and several reaction intermediates agree well with the relevant crystal structures, validating the computational protocol. The highest computed energy barriers at the cyclization-dehydration (17 kcal/mol) and oxidation (21 kcal/mol) steps agree well with the values derived from the kinetics measurements (20.7 and 22.7 kcal/mol, respectively). The simulations provide strong support to the mechanism involving the cyclization-dehydration-oxidation sequence of the chromophore's maturation reactions. The results also establish a solid basis for predictions of maturation mechanisms in other fluorescent proteins.