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.

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.

Quantum-based Modeling of Protein-ligand Interaction: The Complex of RutA with Uracil and Molecular Oxygen.

Modern quantum-based methods are employed to model interaction of the flavin-dependent enzyme RutA with the uracil and oxygen molecules. This complex presents the structure of reactants for the chain of chemical reactions of monooxygenation in the enzyme active site, which is important in drug metabolism. In this case, application of quantum-based approaches is an essential issue, unlike conventional modeling of protein-ligand interaction with force fields using molecular mechanics and classical molecular dynamics methods. We focus on two difficult problems to characterize the structure of reactants in the RutA-FMN-O2 -uracil complex, where FMN stands for the flavin mononucleotide species. First, location of a small O2 molecule in the triplet spin state in the protein cavities is required. Second, positions of both ligands, O2 and uracil, must be specified in the active site with a comparable accuracy. We show that the methods of molecular dynamics with the interaction potentials of quantum mechanics/molecular mechanics theory (QM/MM MD) allow us to characterize this complex and, in addition, to surmise possible reaction mechanism of uracil oxygenation by RutA.

How Reproducible Are QM/MM Simulations? Lessons from Computational Studies of the Covalent Inhibition of the SARS-CoV-2 Main Protease by Carmofur.

This work explores the level of transparency in reporting the details of computational protocols that is required for practical reproducibility of quantum mechanics/molecular mechanics (QM/MM) simulations. Using the reaction of an essential SARS-CoV-2 enzyme (the main protease) with a covalent inhibitor (carmofur) as a test case of chemical reactions in biomolecules, we carried out QM/MM calculations to determine the structures and energies of the reactants, the product, and the transition state/intermediate using analogous QM/MM models implemented in two software packages, NWChem and Q-Chem. Our main benchmarking goal was to reproduce the key energetics computed with the two packages. Our results indicate that quantitative agreement (within the numerical thresholds used in calculations) is difficult to achieve. We show that rather minor details of QM/MM simulations must be reported in order to ensure the reproducibility of the results and offer suggestions toward developing practical guidelines for reporting the results of biosimulations.

Protonation States of Molecular Groups in the Chromophore-Binding Site Modulate Properties of the Reversibly Switchable Fluorescent Protein rsEGFP2.

The role of protonation states of the chromophore and its neighboring amino acid side chains of the reversibly switching fluorescent protein rsEGFP2 upon photoswitching is characterized by molecular modeling methods. Numerous conformations of the chromophore-binding site in computationally derived model systems are obtained using the quantum chemistry and QM/MM approaches. Excitation energies are computed using the extended multiconfigurational quasidegenerate perturbation theory (XMCQDPT2). The obtained structures and absorption spectra allow us to provide an interpretation of the observed structural and spectral properties of rsEGFP2 in the active ON and inactive OFF states. The results demonstrate that in addition to the dominating anionic and neutral forms of the chromophore, the cationic and zwitterionic forms may participate in the photoswitching of rsEGFP2. Conformations and protonation forms of the Glu223 and His149 side chains in the chromophore-binding site play an essential role in stabilizing specific protonation forms of the chromophore.

Stalling chromophore synthesis of the fluorescent protein Venus reveals the molecular basis of the final oxidation step.

Fluorescent proteins (FPs) have revolutionised the life sciences, but the mechanism of chromophore maturation is still not fully understood. Here we show that incorporation of a photo-responsive non-canonical amino acid within the chromophore stalls maturation of Venus, a yellow FP, at an intermediate stage; a crystal structure indicates the presence of O2 located above a dehydrated enolate form of the imidazolone ring, close to the strictly conserved Gly67 that occupies a twisted conformation. His148 adopts an "open" conformation so forming a channel that allows O2 access to the immature chromophore. Absorbance spectroscopy supported by QM/MM simulations suggests that the first oxidation step involves formation of a hydroperoxyl intermediate in conjunction with dehydrogenation of the methylene bridge. A fully conjugated mature chromophore is formed through release of H2O2, both in vitro and in vivo. The possibility of interrupting and photochemically restarting chromophore maturation and the mechanistic insights open up new approaches for engineering optically controlled fluorescent proteins.

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.

Interplay between Locally Excited and Charge Transfer States Governs the Photoswitching Mechanism in the Fluorescent Protein Dreiklang.

We present the results of high-level electronic structure and dynamics simulations of the photoactive protein Dreiklang. With the goal of understanding the details of the Dreiklang photocycle, we carefully characterize the excited states of the ON- and OFF-forms of Dreiklang. The key finding of our study is the existence of a low-lying excited state of a charge-transfer character in the neutral ON form and that population of this state, which is nearly isoenergetic with the locally excited bright state, initiates a series of steps that ultimately lead to the formation of the hydrated dark chromophore (OFF state). These results allow us to refine the mechanistic picture of Dreiklang's photocycle and photoactivation.

Structure of the Brain N-Acetylaspartate Biosynthetic Enzyme NAT8L Revealed by Computer Modeling.

We report the results of computational modeling of a three-dimensional all-atom structure of the membrane-associated protein encoded by the NAT8L gene, aspartate N-acetyltransferase, which is essential for brain synthesis of N-acetyl-l-aspartate (NAA). The lack of experimentally derived three-dimensional structures of NAT8L poses one of the obstacles in studies of the mechanism of NAA formation and understanding the precise role of NAA in neurological disorders. We apply a computational protocol employing the contact map prediction, ab initio folding, homology modeling, and refinement to obtain a structure of NAT8L with the aspartate and acetyl coenzyme A cofactors in the protein molecule. To verify the computational protocol, we check its predictive power by reproducing the crystal structure of a related N-acetyltransferase domain, specifically, that from the bacterial N-acetylglutamate synthase. We show that the constructed NAT8L model correlates with structural features of the protein revealed in rare experimental studies. The obtained structure of the enzyme active site with the trapped reactants suggests a mechanism of the acetyl transfer upon NAA formation.

Computational Modeling Reveals the Mechanism of Fluorescent State Recovery in the Reversibly Photoswitchable Protein Dreiklang.

The unique properties of the photoswitchable protein Dreiklang are attributed to a reversible hydration/dehydration reaction at the imidazolinone ring of the chromophore. Recovery of the fluorescent state, which is associated with a chemical reaction of the chromophore's dehydration, is an important part of the photocycle of this protein. Here we characterize the fluorescent (ON) and nonfluorescent (OFF) states of Dreiklang and simulate the thermal recovery reaction OFF → ON using computational approaches. By using molecular modeling methods including the quantum mechanics/molecular mechanics (QM/MM) technique, we characterize the structures and spectra of the ON- and OFF-states. The results are consistent with available experimental data. The computed reaction profile explains the observed recovery reaction and clarifies the mechanism of chemical transformations in the chromophore-containing pocket in Dreiklang.

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.

Computational Challenges in Modeling of Representative Bioimaging Proteins: GFP-Like Proteins, Flavoproteins, and Phytochromes.

Remarkable success in engineering novel efficient biomarkers based on fluorescent and photoactive proteins provokes a question of whether computational modeling of their properties can contribute to this important field. In this Feature Article, we analyze selected papers devoted to computer simulations of three types of photoactive systems: the green fluorescent protein and its derivatives, the flavin-binding proteins, and the phytochrome domains. The main emphasis is on structures, optical spectra, and chemical reactions in the chromophore-containing pockets. Quantum chemistry, quantum mechanics/molecular mechanics, and molecular dynamics methods are effective tools in these simulations. We highlight both the success stories and the persisting challenges, discussing the ways of elevating theoretical approaches to the level of testable predictions.

Computer-designed active human butyrylcholinesterase double mutant with a new catalytic triad.

A computer-designed mutant of human butyrylcholinesterase (BChE), N322E/E325G, with a novel catalytic triad was made. The catalytic triad of the wild-type enzyme (S198·H438·E325) was replaced by S198·H438·N322E in silico. Molecular dynamics for 1.5 µs and Markov state model analysis showed that the new catalytic triad should be operative in the mutant enzyme, suggesting functionality. QM/MM modeling performed for the reaction of wild-type BChE and double mutant with echothiophate showed high reactivity of the mutant towards the organophosphate. A truncated monomeric (L530 stop) double mutant was expressed in Expi293 cells. Non-purified transfected cell culture medium was analyzed. Polyacrylamide gel electrophoresis under native conditions followed by activity staining with BTC as the substrate provided evidence that the monomeric BChE mutant was active. Inhibition of the double mutant by echothiophate followed by polyacrylamide gel electrophoresis and activity staining showed that this enzyme slowly self-reactivated. However, because Expi293 cells secrete an endogenous BChE tetramer and several organophosphate-reacting enzymes, catalytic parameters and self-reactivation constants after phosphorylation of the new mutant were not determined in the crude cell culture medium. The study shows that the computer-designed double mutant (N322E/E325G) with a new catalytic triad (S198·H438·N322E) is a suitable template for design of novel active human BChE mutants that display an organophosphate hydrolase activity.

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.

Optimization of Cholinesterase-Based Catalytic Bioscavengers Against Organophosphorus Agents.

Organophosphorus agents (OPs) are irreversible inhibitors of acetylcholinesterase (AChE). OP poisoning causes major cholinergic syndrome. Current medical counter-measures mitigate the acute effects but have limited action against OP-induced brain damage. Bioscavengers are appealing alternative therapeutic approach because they neutralize OPs in bloodstream before they reach physiological targets. First generation bioscavengers are stoichiometric bioscavengers. However, stoichiometric neutralization requires administration of huge doses of enzyme. Second generation bioscavengers are catalytic bioscavengers capable of detoxifying OPs with a turnover. High bimolecular rate constants (kcat/Km > 106 M-1min-1) are required, so that low enzyme doses can be administered. Cholinesterases (ChE) are attractive candidates because OPs are hemi-substrates. Moderate OP hydrolase (OPase) activity has been observed for certain natural ChEs and for G117H-based human BChE mutants made by site-directed mutagenesis. However, before mutated ChEs can become operational catalytic bioscavengers their dephosphylation rate constant must be increased by several orders of magnitude. New strategies for converting ChEs into fast OPase are based either on combinational approaches or on computer redesign of enzyme. The keystone for rational conversion of ChEs into OPases is to understand the reaction mechanisms with OPs. In the present work we propose that efficient OP hydrolysis can be achieved by re-designing the configuration of enzyme active center residues and by creating specific routes for attack of water molecules and proton transfer. Four directions for nucleophilic attack of water on phosphorus atom were defined. Changes must lead to a novel enzyme, wherein OP hydrolysis wins over competing aging reactions. Kinetic, crystallographic, and computational data have been accumulated that describe mechanisms of reactions involving ChEs. From these studies, it appears that introducing new groups that create a stable H-bonded network susceptible to activate and orient water molecule, stabilize transition states (TS), and intermediates may determine whether dephosphylation is favored over aging. Mutations on key residues (L286, F329, F398) were considered. QM/MM calculations suggest that mutation L286H combined to other mutations favors water attack from apical position. However, the aging reaction is competing. Axial direction of water attack is not favorable to aging. QM/MM calculation shows that F329H+F398H-based multiple mutants display favorable energy barrier for fast reactivation without aging.

Improving the Design of the Triple-Decker Motif in Red Fluorescent Proteins.

We characterize computationally a red fluorescent protein (RFP) with the chromophore (Chro) sandwiched between two aromatic tyrosine rings in a triple-decker motif. According to the original proposal [ J. Phys. Chem. Lett. 2013 , 4 , 1743 ], such a tyrosine-chromophore-tyrosine π-stacked construct can be accommodated in the green fluorescent protein (GFP). A recent study [ ACS Chem. Biol. 2016 , 11 , 508 ] attempted to realize the triple-decker motif and obtained an RFP variant called mRojoA-VYGV with two tyrosine residues surrounding the chromophore. The crystal structure showed that only a tyrosine-chromophore pair was involved in π-stacking, whereas the second tyrosine was oriented perpendicularly, edge-to-face with respect to the chromophore. We propose a more promising variant of this RFP with a perfect triple-decker unit achieved by introducing additional mutations in mRojoA-VYGV. The structures and optical properties of model proteins based on the structures of mCherry and mRojoA are characterized computationally by QM(DFT)/MM. The electronic transitions in the protein-bound chromophores are computed by high-level quantum chemical methods. According to our calculations, the triple-decker chromophore unit in the new RFP variant is stable within the protein and its optical bands are red-shifted with respect to the parent mCherry and mRojoA species.

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.

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.

Analysis of proton wires in the enzyme active site suggests a mechanism of c-di-GMP hydrolysis by the EAL domain phosphodiesterases.

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c-di-GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme-substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c-di-GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c-di-GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c-di-GMP to promote the O3'P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670-1680. © 2016 Wiley Periodicals, Inc.