Search

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol -1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway., (© 2022 Wiley-VCH GmbH.)

To protect their intracellular proteins, extremophile microorganisms synthesize molecules called compatible solutes. These molecules are the result of the attachment of a small negatively charged molecule to a sugar molecule. It has been found that these molecules, not only protect the microorganism against osmotic stress but also against other extreme conditions. They can also confer protection against extreme conditions to isolated enzymes from different organisms making them an exciting prospect for potential biotechnological applications. One of the most widespread compatible solute in hyperthermophile organisms is the molecule 2-O-α-D-mannosyl-D-glycerate (MG). In addition to confer protection to proteins against extreme conditions, MG was found to prevent Alzheimer's β-amyloid aggregation and reduce α-synuclein fibril formation in Parkinson's disease. In this work we studied, using computational methods, the catalytic mechanism of the synthesis of MG by the enzyme mannosylglycerate synthase (MGS) from the thermophilic bacteria Rhodothermus marinus., (© 2021 Wiley-VCH GmbH.)

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg 2+ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The enhancement of the catalytic power of enzymes is a subject of enormous interest both for science and for industry. The latter, in particular, due to the vast applications enzymes can have in industrial processes, for instance in the desulfurization of crude oil, which is mandatory by law in many developed countries and is currently performed using costly chemical processes. In this work we sought to enhance the turnover rate of DszD from Rhodococcus erythropolis, a NADH-FMN oxidoreductase responsible for supplying FMNH 2 to DszA and DszC in the biodesulfurization process of crude oil, the 4S pathway. For this purpose, we replaced the wild type spectator residue of the rate limiting step of the reduction of FMN to FMNH 2 , a process catalysed by DszD and known to play an important role in the reaction energy profile. As replacements, we used all the naturally occurring amino acids, one at a time, using computational methodologies, and repeated the above-mentioned reaction with each mutant. To calculate the different free energy profiles, one for each mutated model, we applied quantum mechanics/molecular mechanics (QM/MM) methods within an ONIOM scheme. The free energy barriers obtained varied between 15.1 and 29.9 kcal mol -1 . Multiple factors contributed to the different ΔG values. The most relevant were electrostatic interactions and the induction of a favourable alignment between substrate and cofactor. These results confirm the great potential that chirurgic mutations have for increasing the catalytic power of DszD in relation to the wild type (wt) enzyme., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Cyclic peptides are a class of compounds with high therapeutic potential, possessing bioactivities including antitumor and antiviral (including anti-HIV). Despite their desirability, efficient design and production of these compounds has not been achieved to date. The catalytic mechanism of patellamide macrocyclization by the PatG macrocyclase domain has been computationally investigated by using quantum mechanics/molecular mechanics methodology, specifically ONIOM(M06/6-311++G(2d,2p):ff94//B3LYP/6-31G(d):ff94). The mechanism proposed herein begins with a proton transfer from Ser783 to His 618 and from the latter to Asp548. Nucleophilic attack of Ser783 on the substrate leads to the formation of an acyl-enzyme covalent complex. The leaving group Ala-Tyr-Asp-Gly (AYDG) of the substrate is protonated by the substrate's N terminus, leading to the breakage of the P1-P1' bond. Finally, the substrate's N terminus attacks the P1 residue, decomposing the acyl-enzyme complex forming the macrocycle. The formation and decomposition of the acyl-enzyme complex have the highest activation free energies (21.1 kcal mol(-1) and 19.8 kcal mol(-1) respectively), typical of serine proteases. Understanding the mechanism behind the macrocyclization of patellamides will be important to the application of the enzymes in the pharmaceutical and biotechnological industries., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

This paper is devoted to the understanding of the reaction mechanism of mycobacterium tuberculosis glutamine synthetase (mtGS) with atomic detail, using computational quantum mechanics/molecular mechanics (QM/MM) methods at the ONIOM M06-D3/6-311++G(2d,2p):ff99SB//B3LYP/6-31G(d):ff99SB level of theory. The complete reaction undergoes a three-step mechanism: the spontaneous transfer of phosphate from ATP to glutamate upon ammonium binding (ammonium quickly loses a proton to Asp54), the attack of ammonia on phosphorylated glutamate (yielding protonated glutamine), and the deprotonation of glutamine by the leaving phosphate. This exothermic reaction has an activation free energy of 21.5 kcal mol(-1) , which is consistent with that described for Escherichia coli glutamine synthetase (15-17 kcal mol(-1) ). The participating active site residues have been identified and their role and energy contributions clarified. This study provides an insightful atomic description of the biosynthetic reaction that takes place in this enzyme, opening doors for more accurate studies for developing new anti-tuberculosis therapies., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In order to elucidate the catalytic mechanism of the Mn-Mn containing serine/threonine protein phosphatase 5 (PP5), we present a density functional theory study with a cluster model approach. According to our results, the reaction occurs through an in-line concerted transition state with an energy of 15.8 kcal mol(-1) , and no intermediates are formed. The important role played by His304 and Asp274 as stabilizers of the leaving group has been shown, whereas the role played by the metal ions seems to be mostly electrostatic. The indispensable requirement of having a neutral active center has been demonstrated by testing different protonation states of the cluster model. We have shown also the importance of describing properly the electronic configuration of the Mn-Mn binuclear centers., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Atomic-level portrait: The mechanism of the reaction catalyzed by the puzzling enzyme farnesyltransferase is elucidated by using computational methods, allowing the obtainment of the first real detailed atomistic quantum-chemical transition-state structure (see figure) for the reaction catalyzed by this enzyme. The results obtained provide an atomic-level framework for the design of more potent and specific inhibitors for this important enzyme.

Glutathione transferases are enzymes of the cellular detoxification system that metabolize a vast spectrum of xenobiotic and endobiotic toxic compounds. They are homodimers or heterodimers and each monomer has an active center composed of a G-site in which glutathione (GSH) binds and an H-site for the electrophilic substrate. When GSH binds to the G-site, the pKa value of its thiol group drops by 2.5 units; this promotes its deprotonation and, therefore, produces a strong nucleophilic thiolate that is able to react with the electrophilic substrate. The mechanism behind the deprotonation of the thiol group is still unknown. Some studies point to the fact that the GSH glutamyl alpha-carboxylate group is essential for GSH activation, whereas others indicate the importance of the active-center water molecules. On the basis of QM/MM calculations, we propose a mechanism of GSH activation in which a water molecule, acting as a bridge, is able to assist in the transfer of the proton from the GSH thiol group to the GSH glutamyl alpha-carboxylate group, after an initial GSH conformational rearrangement. We calculated the potential of mean force of this GSH structural rearrangement that would be necessary for the approach of both groups and we then performed a QM/MM ONIOM scan of water-assisted proton transfer. The overall free-energy barrier for the process is consistent with experimental studies of the enzyme kinetics.

This paper focuses on the inhibition of ribonucleotide reductase (RNR) by gemcitabine, 2',2'-difluoro-2'-deoxycytidine (dFdC), a deoxycytidine analogue that is a very active drug against solid tumors and is currently commercialized as gemzar. RNR inactivation is reductant-dependent and occurs in a very different way from that of other known substrate analogues. In the presence of reductants monomer R1 of RNR is inhibited, whereas in the absence of reductants the radical is lost and monomer R2 is inhibited. As inside the cell reductants are available, it is likely that R1 inactivation is the most favorable mechanism responsible for drug cytotoxicity. This inhibition pathway has been unknown to date, but we have conducted a theoretical study that has led us to the first proposal of a mechanism for RNR inhibition by dFdC in the presence of reductants.

The mechanism for thiol/disulfide exchange has been studied with high-level theoretical calculations. Free energies, transition structures, charge densities, and solvent effects along the reaction pathway have been determined for the first time. Mechanistic results agree with experimental data, and support the idea that the thiolate is the reacting species and that the reaction indeed proceeds through an uncomplicated S(N)2 transition state. The transition structures have the charge density evenly concentrated in the attacking and leaving sulfur atoms. The charge densities allow us to rationalize the solvent effects. As transition structures have the charge density more widely distributed than reactants, hydrophobic environments catalyze the reaction. The effect can be so dramatic that disulfide exchange inside the active site of ribonucleotide reductase is estimated to be catalyzed 10(3) times faster than the reaction in water. It was also found that attack by thiol is much faster than previously assumed, if mediated through water chains. Although the present results, as well as experimental data, still suggest that thiolate is the main reaction species, water-mediated thiol attack is almost kinetically competitive, and can eventually become competitive under specific experimental conditions.

Several 2'-substituted-2'-deoxyribonucleotides are potent time-dependent inactivators of the enzyme ribonucleotide reductase (RNR), which function by destructing its essential tyrosil radical and/or by performing covalent addition to the enzyme. The former leads to inhibition of the R2 dimer of RNR and the latter to inhibition of the R1 dimer. Efforts to elucidate the mechanism of inhibition have been undertaken in the last decades, and a general mechanistic scheme has emerged. Accordingly, two alternative pathways lead either to the inhibition of R1 or R2, for which the 2'-chloro-2'-deoxynucleotides serve as the model for the inhibition of R1 and the 2'-azido-2'-deoxynucleotides the model for the inhibition of R2. However, the underlying reason for the different behavior of the inhibitors has remained unknown until now. Moreover, a fundamental mechanistic alternative has been proposed, based on results from biomimetic reactions, in which the 2'-substituents would be eliminated as radicals, and not as anions, as previously assumed. This would lead to further reactions not predicted by the existing mechanistic scheme. To gain a better understanding we have performed high-level theoretical calculations on the active site of RNR. Results from this work support the general Stubbe's paradigm, although some changes to that mechanism are necessary. In addition, a rational explanation of the factors that determine which of the dimers (R1 or R2) will be inactivated is provided for the first time. It has been demonstrated also that the 2'-substituents are indeed eliminated as anions, and not as radicals. Biomimetic experiments have led to different results because they lack a basic group capable of deprotonating the 3'-HO group of the substrate. It has been found here that the chemical character of the leaving group (radical or anionic) can be manipulated by controlling the protonation state of the 3'-HO group.