New Insights into the Catalytic Mechanism of Aldose Reductase: A QM/MM Study.
Ontology highlight
ABSTRACT:
SUBMITTER: Dreanic MP
PROVIDER: S-EPMC6044956 | biostudies-literature | 2017 Sep
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the ?-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and ?-aminoadipate-?-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, ?-aminoadipate-?-semialdehyde and NADPH in a manner and polar environment conducive to reaction.
| S-EPMC6264447 | biostudies-literature
Project description:SARS-CoV-2 Mpro, also known as the main protease or 3C-like protease, is a key enzyme involved in the replication process of the virus that is causing the COVID-19 pandemic. It is also the most promising antiviral drug target targeting SARS-CoV-2 virus. In this work, the catalytic mechanism of Mpro was studied using the full model of the enzyme and a computational QM/MM methodology with a 69/72-atoms QM region treated at DLPNO-CCSD(T)/CBS//B3LYP/6-31G(d,p):AMBER level and including the catalytic important oxyanion-hole residues. The transition state of each step was fully characterized and described together with the related reactants and products. The rate-limiting step of the catalytic process is the hydrolysis of the thioester-enzyme adduct, and the calculated barrier closely agrees with the available kinetic data. The calculated Gibbs free energy profile, together with the full atomistic detail of the structures involved in catalysis, can now serve as valuable models for the rational drug design of transition state analogs as new inhibitors targeting the SARS-CoV-2 virus.
| S-EPMC8224256 | biostudies-literature
Project description:Catalysis by S-adenosylmethionine synthetase has been investigated by quantum mechanical/molecular mechanical calculations, exploiting structures of the active crystalline enzyme. The transition state energy of +19.1 kcal/mol computed for a nucleophilic attack of the methionyl sulfur on carbon-5' of the nucleotide was indistinguishable from the experimental (solution) value when the QM residues were an uncharged histidine that hydrogen bonds to the leaving oxygen-5' and an aspartate that chelates a Mg2+ ion, and was similar (+18.8 kcal/mol) when the QM region also included the active site arginine and lysines. The computed energy difference between reactant and product was also consistent with their equimolar abundance in co-crystals. The calculated geometrical changes support catalysis of a S(N)2 reaction through hydrogen bonding of the liberated oxygen-5' to the histidine, charge neutralization by the two Mg2+ ions, and stabilization of the product sulfonium cation through a close, non-bonded, contact between the sulfur and the ribose oxygen-4'.
| S-EPMC2788016 | biostudies-literature
Project description:Sulfite oxidase is a mononuclear molybdenum enzyme that oxidises sulfite to sulfate in many organisms, including man. Three different reaction mechanisms have been suggested, based on experimental and computational studies. Here, we study all three with combined quantum mechanical (QM) and molecular mechanical (QM/MM) methods, including calculations with large basis sets, very large QM regions (803?atoms) and QM/MM free-energy perturbations. Our results show that the enzyme is set up to follow a mechanism in which the sulfur atom of the sulfite substrate reacts directly with the equatorial oxo ligand of the Mo ion, forming a Mo-bound sulfate product, which dissociates in the second step. The first step is rate limiting, with a barrier of 39-49?kJ/mol. The low barrier is obtained by an intricate hydrogen-bond network around the substrate, which is preserved during the reaction. This network favours the deprotonated substrate and disfavours the other two reaction mechanisms. We have studied the reaction with both an oxidised and a reduced form of the molybdopterin ligand and quantum-refinement calculations indicate that it is in the normal reduced tetrahydro form in this protein.
| S-EPMC5856855 | biostudies-literature
Project description:Sortases are key virulence factors in Gram-positive bacteria. These enzymes embed surface proteins in the cell wall through a transpeptidation reaction that involves recognizing a penta-peptide "sorting signal" in a target protein, cleaving it, and covalently attaching it to a second substrate that is later inserted into the cell wall. Although well studied, several aspects of the mechanism by which sortases perform these functions remains unclear. In particular, experiments have revealed two potential sorting signal binding motifs: a "Threonine-Out" (Thr-Out) structure in which the catalytically critical threonine residues protrudes into solution, and a "Threonine-In" (Thr-In) configuration in which this residue inserts into the binding site. To determine which of these is the biologically relevant state, we have performed a series of conventional and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations of the Staphylococcus aureus sortase A (SrtA) enzyme bound to a sorting signal substrate. Through the use of multi-dimensional metadynamics, our simulations were able to both map the acylation mechanism of SrtA in the Thr-In and Thr-Out states, as well as determine the free energy minima and barriers along these reactions. Results indicate that in both states the catalytic mechanisms are similar, however the free energy barriers are lower in the Thr-In configuration, suggesting that Thr-In is the catalytically relevant state. This has important implications for advancing our understanding of the mechanisms of sortase enzymes, as well we for future structure based drug design efforts aimed at inhibiting sortase function in vivo.
| S-EPMC4917473 | biostudies-literature
Project description:Multidrug-resistant organisms contain antibiotic-modifying enzymes that facilitate resistance to a variety of antimicrobial compounds. Particularly, the fosfomycin (FOF) drug can be structurally modified by several FOF-modifying enzymes before it reaches the biological target. Among them, FosB is an enzyme that utilizes l-cysteine or bacillithiol in the presence of a divalent metal to open the epoxide ring of FOF and, consequently, inactivate the drug. Here, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of the reaction involving FosB and FOF. The calculated free-energy profiles show that the cost to open the epoxide ring of FOF at the C2 atom is ∼3.0 kcal/mol higher than that at the C1 atom. Besides, our QM/MM MD results revealed the critical role of conformation change of Cys9 and Asn50 to release the drug from the active site. Overall, the present study provides insights into the mechanism of FOF-resistant proteins.
| S-EPMC8154160 | biostudies-literature
Project description:The environmental problems derived from the generalized plastic consumption and disposal could find a friendly solution in enzymatic biodegradation. Recently, two hydrolases from Ideonella sakaiensis 201-F6 and the metagenome-derived leaf-branch compost cutinase (LCC), more specially the improved ICCG variant, have revealed degradation activity toward poly ethylene terephthalate (PET). In the present study, the reaction mechanism of this polymer breakage is studied at an atomic level by multiscale QM/MM molecular dynamics simulations, using semiempirical and DFT Hamiltonians to describe the QM region. The obtained free energy surfaces confirmed a characteristic four-step path for both systems, with activation energies in agreement with the experimental observations. Structural analysis of the evolution of the active site along the reaction progress and the study of electrostatic effects generated by the proteins reveal the similarity in the behavior of the active site of these two enzymes. The origin of the apparent better performance of the LCC-ICCG protein over PETase must be due to its capabilities of working at higher temperature and its intrinsic relationship with the crystallinity grade of the polymer. Our results may be useful for the development of more efficient enzymes in the biodegradation of PET.
| S-EPMC8576761 | biostudies-literature
Project description:Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1-2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.
| S-EPMC2919153 | biostudies-literature
Project description:Lipoyl synthase (LipA) catalyses the final step of the biosynthesis of the lipoyl cofactor by insertion of two sulfur atoms at the C6 and C8 atoms of the protein-bound octanoyl substrate. In this reaction, two [4Fe4S] clusters and two molecules of S-adenosyl-L-methionine are used. One of the two FeS clusters is responsible for the generation of a powerful oxidant, the 5'-deoxyadenosyl radical (5'-dA•). The other (the auxiliary cluster) is the source of both sulfur atoms that are inserted into the substrate. In this paper, the spin state of the FeS clusters and the reaction mechanism is investigated by the combined quantum mechanical and molecular mechanics approach. The calculations show that the ground state of the two FeS clusters, both in the [4Fe4S]2+ oxidation state, is a singlet state with antiferromagnetically coupled high-spin Fe ions and that there is quite a large variation of the energies of the various broken-symmetry states, up to 40 kJ/mol. For the two S-insertion reactions, the highest energy barrier is found for the hydrogen-atom abstraction from the octanoyl substrate by 5'-dA•. The formation of 5'-dA• is very facile for LipA, with an energy barrier of 6 kJ/mol for the first S-insertion reaction and without any barrier for the second S-insertion reaction. In addition, the first S ion attack on the C6 radical of octanoyl was found to take place directly by the transfer of the H6 from the substrate to 5'-dA•, whereas for the second S-insertion reaction, a C8 radical intermediate was formed with a rate-limiting barrier of 71 kJ/mol.
| S-EPMC5816104 | biostudies-literature
Project description:Using a structure generated by induced fit modeling of the protein-ligand complex, the reaction path for hydrogen atom abstraction in P450 BM3 is studied by means of mixed QM/MM methods to determine the structures and energetics along the reaction path. The IFD structure is suitable for hydrogen atom abstraction at the ?-1 position. The electronic structures obtained are similar to those observed in P450 cam. We show that the barrier for the hydrogen abstraction step from QM/MM modeling is 13.3 kcal/mol in quartet and 15.6 kcal/mol in doublet. Although there is some strain energy present in the ligand, the activation barrier is not dramatically affected. A crystal water molecule, HOH502, plays a role as catalyst and decreases the activation barrier by about 2 kcal/mol and reaction energy by about 3-4 kcal/mol. In order to achieve reactive chemistry at the remaining experimentally observed positions in the hydrocarbon tail of the ligand, other structures would have to be utilized as a starting point for the reaction. Finally, the present results still leave open the question of whether DFT methods provide an accurate computation of the barrier height in the P450 hydrogen atom abstraction reaction.
| S-EPMC2757081 | biostudies-literature