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.

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.

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.

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.

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.

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.

Mechanisms of ATP to cAMP Conversion Catalyzed by the Mammalian Adenylyl Cyclase: A Role of Magnesium Coordination Shells and Proton Wires.

We report a mechanism of adenosine triphosphate (ATP) to cyclic adenosine monophosphate (cAMP) conversion by the mammalian type V adenylyl cyclase revealed in molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) simulations. We characterize a set of computationally derived enzyme-substrate (ES) structures showing an important role of coordination shells of magnesium ions in the solvent accessible active site. In the lowest energy ES conformation, the coordination shell of MgA2+ does not include the Oδ1 atom of the conserved Asp440 residue. Starting from this conformation, a one-step reaction mechanism is characterized that includes proton transfer from the ribose O3'H3' group in ATP to Asp440 via a shuttling water molecule concerted with PA-O3A bond cleavage and O3'-PA bond formation. The energy profile of this route is consistent with the observed reaction kinetics. The computed energy profiles initiated from higher energy ES complexes are characterized by larger energy expenses to complete the reaction. Consistent with experimental data, we show that the Asp440Ala mutant of the enzyme should exhibit a reduced but retained activity. All considered reaction pathways include proton wires from the O3'H3' group via shuttling water molecules.

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.