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

SIRT5 mutants reveal the role of conserved asparagine and glutamine residues in the NAD+-binding pocket

SIRT5, one of the mammalian sirtuins, specifically recognizes succinyl-lysine residues on proteins and catalyzes the desuccinylation reaction. In this study, we characterized SIRT5 mutants with hydrophobic amino acid substitutions at Q140 and N141, in addition to the catalytic residue H158, known as an active site residue, by the Michaelis-Menten analysis and X-ray crystallography. Kinetic analysis showed that the catalytic efficiency (kcat/Km) of the Q140L and N141V mutants decreased to 0.02 times and 0.0038 times that of the wild-type SIRT5, respectively, with the activity of the N141V mutant becoming comparable to that of the H158M mutant. Our findings indicate that N141 contributes significantly to the desuccinylation reaction.

Theoretical investigation of tautomerism of 2- and 4-pyridones: origin, substituent and solvent effects

Computational investigation at the BHandHLYP/6-311+G(d,p) level of theory of the gas-phase tautomerism of 2- and 4-pyridones confirmed the slight prevalence of lactim in the case of the former, but its dominance in the case of the latter, as shown previously. Examination of aromaticity by using HOMA, EDDB, NBOdel, NICS and AICD led to the conclusion that tautomerization of 4-pyridone results in greater aromaticity gain. It is also driven by the Pauli repulsion relief, which was revealed by the tautomerization energy decomposition analysis. By contrast, in the case of 2-pyridone, lactim is favoured by orbital and electrostatic interactions and disfavoured by the Pauli repulsion. Aromaticity gain in this case is smaller. The position of the tautomeric equilibrium can be modulated by substituent inductive effects (Cl and F), inductive and resonance effects (NH2 and NO2), hydrogen bonding (NO2), and medium polarity, the increase of which increases lactam population.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Photophysics of the Blue Light Using Flavin Domain

Light activated proteins are at the heart of photobiology and optogenetics, so there is wide interest in understanding the mechanisms coupling optical excitation to protein function. In addition, such light activated proteins provide unique insights into the real-time dynamics of protein function. Using pump-probe spectroscopy, the function of a photoactive protein can be initiated by a sub-100 fs pulse of light, allowing subsequent protein dynamics to be probed from femtoseconds to milliseconds and beyond. Among the most interesting photoactive proteins are the blue light using flavin (BLUF) domain proteins, which regulate the response to light of a wide range of bacterial and some euglenoid processes. The photosensing mechanism of BLUF domains has long been a subject of debate. In contrast to other photoactive proteins, the electronic and nuclear structure of the chromophore (flavin) is the same in dark- and light-adapted states. Thus, the driving force for photoactivity is unclear.To address this question requires real-time observation of both chromophore excited state processes and their effect on the structure and dynamics of the surrounding protein matrix. In this Account we describe how time-resolved infrared (IR) experiments, coupled with chemical biology, provide important new insights into the signaling mechanism of BLUF domains. IR measurements are sensitive to changes in both chromophore electronic structure and protein hydrogen bonding interactions. These contributions are resolved by isotope labeling of the chromophore and protein separately. Further, a degree of control over BLUF photochemistry is achieved through mutagenesis, while unnatural amino acid substitution allows us to both fine-tune the photochemistry and time resolve protein dynamics with spatial resolution.Ultrafast studies of BLUF domains reveal non-single-exponential relaxation of the flavin excited state. That relaxation leads within one nanosecond to the original flavin ground state bound in a modified hydrogen-bonding network, as seen in transient and steady-state IR spectroscopy. The change in H-bond configuration arises from formation of an unusual enol (imine) form of a critical glutamine residue. The dynamics observed, complemented by quantum mechanical calculations, suggest a unique sequential electron then double proton transfer reaction as the driving force, followed by rapid reorganization in the binding site and charge recombination. Importantly, studies of several BLUF domains reveal an unexpected diversity in their dynamics, although the underlying structure appears highly conserved. It is suggested that this diversity reflects structural dynamics in the ground state at standard temperature, leading to a distribution of structures and photochemical outcomes. Time resolved IR measurements were extended to the millisecond regime for one BLUF domain, revealing signaling state formation on the microsecond time scale. The mechanism involves reorganization of a β-sheet connected to the chromophore binding pocket via a tryptophan residue. The potential of site-specific labeling amino acids with IR labels as a tool for probing protein structural dynamics was demonstrated.In summary, time-resolved IR studies of BLUF domains (along with related studies at visible wavelengths and quantum and molecular dynamics calculations) have resolved the photoactivation mechanism and real-time dynamics of signaling state formation. These measurements provide new insights into protein structural dynamics and will be important in optimizing the potential of BLUF domains in optobiology.

Direct Observation of Ultrafast Proton Rocking in the BLUF Domain

We present direct observation of ultrafast proton rocking in the central motif of a BLUF domain protein scaffold. The mutant design has taken considerations of modulating the proton-coupled electron transfer (PCET) driving forces by replacing Tyr in the original motif with Trp, as well as of removing the interference of a competing electron transfer pathway. Using femtosecond pump-probe spectroscopy and detailed kinetics analysis, we resolved an electron-transfer-coupled Grotthuss-type forward and reversed proton rocking along the FMN-Gln-Trp proton relay chain. The rates of forward and reversed proton transfer are determined to be very close, namely 51 ps vs 52 ps. The kinetic isotope effect (KIE) constants associated with the forward and reversed proton transfer are 3.9 and 5.3, respectively. The observation of ultrafast proton rocking is not only a crucial step towards revealing the nature of proton relay in BLUF domain, but also provides a new paradigm of proton transfer in proteins for theoretical investigations.

Asparagine Tautomerization in Glycosyltransferase Catalysis. The Molecular Mechanism of Protein O-Fucosyltransferase 1

O-glycosylation is a post-translational protein modification essential to life. One of the enzymes involved in this process is protein O-fucosyltransferase 1 (POFUT1), which fucosylates threonine or serine residues within a specific sequence context of epidermal growth factor-like domains (EGF-LD). Unlike most inverting glycosyltransferases, POFUT1 lacks a basic residue in the active site that could act as a catalytic base to deprotonate the Thr/Ser residue of the EGF-LD acceptor during the chemical reaction. Using quantum mechanics/molecular mechanics (QM/MM) methods on recent crystal structures, as well as mutagenesis experiments, we uncover the enzyme catalytic mechanism, revealing that it involves proton shuttling through an active site asparagine, conserved among species, which undergoes tautomerization. This mechanism is consistent with experimental kinetic analysis of Caenorhabditis elegans POFUT1 Asn43 mutants, which ablate enzyme activity even if mutated to Asp, the canonical catalytic base in inverting glycosyltransferases. These results will aid inhibitor development for Notch-associated O-glycosylation disorders.

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.

Detecting Individual Bond Switching within Amides in a Tunneling Junction

Amides are essential in the chemistry of life. Detecting the chemical bond states within amides could unravel the nature of amide stabilization and planarity, which is critical to the structure and reactivity of such molecules. Yet, so far, no work has been reported to detect or measure the bond changes at the single-molecule level within amides. Here, we show that a transition between single and double bonds between N and C atoms in an amide can be monitored in real time in a nanogap between gold electrodes via the generation of distinctive conductance features. Density functional theory simulations show that the switching between amide isomers proceeds via a proton transfer process facilitated by a water molecule bridge, and the resulting molecular junctions display bimodal conductance states with a difference as much as nine times.

Direct detection of coupled proton and electron transfers in human manganese superoxide dismutase

Human manganese superoxide dismutase is a critical oxidoreductase found in the mitochondrial matrix. Concerted proton and electron transfers are used by the enzyme to rid the mitochondria of O₂^•−. The mechanisms of concerted transfer enzymes are typically unknown due to the difficulties in detecting the protonation states of specific residues and solvent molecules at particular redox states. Here, neutron diffraction of two redox-controlled manganese superoxide dismutase crystals reveal the all-atom structures of Mn³⁺ and Mn²⁺ enzyme forms. The structures deliver direct data on protonation changes between oxidation states of the metal. Observations include glutamine deprotonation, the involvement of tyrosine and histidine with altered pK_as, and four unusual strong-short hydrogen bonds, including a low barrier hydrogen bond. We report a concerted proton and electron transfer mechanism for human manganese superoxide dismutase from the direct visualization of active site protons in Mn³⁺ and Mn²⁺ redox states.

Human manganese superoxide dismutase (MnSOD) is an oxidoreductase that uses concerted proton and electron transfers to reduce the levels of superoxide radicals in mitochondria, but mechanistic insights into this process are limited. Here, the authors report neutron crystal structures of Mn³⁺SOD and Mn²⁺SOD, revealing changes in the protonation states of key residues in the enzyme active site during the redox cycle.

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.

Formation of an unusual glutamine tautomer in a blue light using flavin photocycle characterizes the light-adapted state

Blue light using flavin (BLUF) photoreceptor proteins are critical for many light-activated biological processes and are promising candidates for optogenetics because of their modular nature and long-range signaling capabilities. Although the photocycle of the Slr1694 BLUF domain has been characterized experimentally, the identity of the light-adapted state following photoexcitation of the bound flavin remains elusive. Herein hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations of this photocycle provide a nonequilibrium dynamical picture of a possible mechanism for the formation of the light-adapted state. Photoexcitation of the flavin induces a forward proton-coupled electron transfer (PCET) process that leads to the formation of an imidic acid tautomer of Gln50. The calculations herein show that the subsequent rotation of Gln50 allows a reverse PCET process that retains this tautomeric form. In the resulting purported light-adapted state, the glutamine tautomer forms a hydrogen bond with the flavin carbonyl group. Additional ensemble-averaged QM/MM calculations of the dark-adapted and purported light-adapted states demonstrate that the light-adapted state with the imidic acid glutamine tautomer reproduces the experimentally observed spectroscopic signatures. Specifically, the calculations reproduce the red shifts in the flavin electronic absorption and carbonyl stretch infrared spectra in the light-adapted state. Further hydrogen-bonding analyses suggest the formation of hydrogen-bonding interactions between the flavin and Arg65 in the light-adapted state, providing a plausible explanation for the experimental observation of faster photoinduced PCET in this state. These characteristics of the light-adapted state may also be essential for the long-range signaling capabilities of this photoreceptor protein.

Ag(I)-Catalyzed Addition of Cyclopropenones and Nitrones to Access Imides

An unprecedented Ag-catalyzed addition reaction of cyclopropenones and nitrones to access imides was developed. Sequential C-C bond cleavage, N-O bond cleavage, and Mumm rearrangement were uncovered in this process. This protocol exhibited high efficiency, regioselectivity, good yields, and a broad tolerance of various functional groups.

Identification and in silico structural and functional analysis of a trypsin-like protease from shrimp Macrobrachium carcinus

Macrobrachium carcinus (Linnaeus, 1758) is a species of freshwater shrimp widely distributed from Florida southwards to southern Brazil, including southeast of Mexico. In the present work, we identified a putative trypsin-like protease cDNA fragment of 736 nucleotides from M. carcinus hepatopancreas tissue by the 3'RACE technique and compared the deduced amino acid sequence to other trypsin-related proteases to describe its structure and function relationship. The bioinformatics analyses showed that the deduced amino acid sequence likely corresponds to a trypsin-like protease closely related to brachyurins, which comprise a subset of serine proteases with collagenolytic activity found in crabs and other crustacea. The M. carcinus trypsin-like protease sequence showed a global sequence identity of 94% with an unpublished trypsin from Macrobrachium rosenbergii (GenBank accession no. AMQ98968), and only 57% with Penaeus vannamei trypsin (GenBank accession no. CAA60129). A detailed analysis of the amino acid sequence revealed specific differences with crustacean trypsins, such as the sequence motif at the beginning of the mature protein, activation mechanism of the corresponding zymogen, amino acid residues of the catalytic triad and residues responsible for substrate specificity.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

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.

Hybrid QM/MM vs Pure MM Molecular Dynamics for Evaluating Water Distribution within p21N-ras and the Resulting GTP Electronic Density

p21^ras protein activity, regulated by GTP hydrolysis, constitutes an active field of research for the development of cancer targeted therapies that would concern ∼30% of human tumors to which specific mutations have been associated. Indeed, the catalyzing mechanisms provided by the protein environment during GTP hydrolysis and how they are impaired by specific mutations remain to be fully elucidated. In this article, we present results from molecular mechanics (MM) molecular dynamics (MD) simulations and density functional theory (DFT) calculations carried out for wild-type p21 ^N-ras and six Gln 61 mutants. In the first part, we present the water distribution within the active site of the wild-type protein according to MM MD. Significant differences are observed when comparing the results to the previous distribution assessed through quantum mechanics/molecular mechanics (QM/MM) MD. Such method-dependent results highlight the importance of accounting for the electrostatic coupling between the protein complex and the solvent molecules in identifying hydration sites. In the second part, we present the results from DFT calculations performed to determine the electronic distribution of the GTP ligand, considering the wild-type active site arrangement according to both classical and hybrid approaches. Only in the QM/MM-based configuration is the ligand electronic density similar to that of a GDP-like state observed experimentally. For this reason, in the last set of calculations carried out for p21 ^N-ras Gln 61 mutants, only the active site structural conformations obtained through hybrid MD are considered. Through the analysis of the GTP electronic density, we conclude that the wild-type active site arrangement according to QM/MM MD is closer to a catalytically efficient conformation of the protein than the arrangement according to MM MD. Hence, water distribution according to the hybrid approach must correspond to the optimal placement of solvent in the active site. Within all of the studied Gln 61 substituted proteins, p21^ras major catalyzing effect, which consists of stabilizing a more GDP-like state, is lost.

Reactivity of Gold(I) Monocarbene Complexes with Protein Targets: A Theoretical Study

Neutral N–heterocyclic carbene gold(I) compounds such as IMeAuCl are widely used both in homogeneous catalysis and, more recently, in medicinal chemistry as promising antitumor agents. In order to shed light on their reactivity with protein side chains, we have carried out density functional theory (DFT) calculations on the thermodynamics and kinetics of their reactions with water and various nucleophiles as a model of plausible protein binding sites such as arginine, aspartic acid, asparagine, cysteine, glutamic acid, glutamine, histidine, lysine, methionine, selenocysteine, and the N-terminal group. In agreement with recent experimental data, our results suggest that IMeAuCl easily interacts with all considered biological targets before being hydrated—unless sterically prevented—and allows the establishment of an order of thermodynamic stability and of kinetic reactivity for its binding to protein residues.