Browse
Submit Data
Databases
API
Help

Dataset Information

31 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

QM Cluster or QM/MM in Computational Enzymology: The Test Case of LigW-Decarboxylase.

ABSTRACT: The catalytic mechanism of the decarboxylation of 5-carboxyvanillate by LigW producing vanillic acid has been studied by using QM cluster and hybrid QM/MM methodologies. In the QM cluster model, the environment of a small QM model is treated with a bulky potential while two QM/MM models studies include partial and full protein with and without explicitly treated water solvent. The studied reaction involves two sequential steps: the protonation of the carbon of the 5-carboxy-vanillate substrate and the decarboxylation of the intermediate from which results deprotonated vanillic acid as product. The structures and energetics obtained by using three structural models and two density functionals are quite consistent to each other. This indicates that the small QM cluster model of the presently considered enzymatic reaction is appropriate enough and the reaction is mainly influenced by the active site.

SUBMITTER: Prejano M

PROVIDER: S-EPMC6031855 | biostudies-literature | 2018

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

QM Cluster or QM/MM in Computational Enzymology: The Test Case of LigW-Decarboxylase.

Prejanò Mario M Marino Tiziana T Russo Nino N

Frontiers in chemistry 20180628

The catalytic mechanism of the decarboxylation of 5-carboxyvanillate by LigW producing vanillic acid has been studied by using QM cluster and hybrid QM/MM methodologies. In the QM cluster model, the environment of a small QM model is treated with a bulky potential while two QM/MM models studies include partial and full protein with and without explicitly treated water solvent. The studied reaction involves two sequential steps: the protonation of the carbon of the 5-carboxy-vanillate substrate a ...[more]

PMID: 30003076

Similar Datasets

Comparative Computational Approach To Study Enzyme Reactions Using QM and QM-MM Methods.

Project description:Choline oxidase catalyzes oxidation of choline into glycine betaine through a two-step reaction pathway employing flavin as the cofactor. On the light of kinetic studies, it is proposed that a hydride ion is transferred from α-carbon of choline/hydrated-betaine aldehyde to the N5 position of flavin in the rate-determining step, which is preceded by deprotonation of hydroxyl group of choline/hydrated-betaine aldehyde to one of the possible basic side chains. Using the crystal structure of glycine betaine-choline oxidase complex, we formulated two computational systems to study the hydride-transfer mechanism including main active-site amino acid side chains, flavin cofactor, and choline as a model system. The first system used pure density functional theory calculations, whereas the second approach used a hybrid ONIOM approach consisting of density functional and molecular mechanics calculations. We were able to formulate in silico model active sites to study the hydride-transfer steps by utilizing noncovalent chemical interactions between choline/betaine aldehyde and active-site amino acid chains using an atomistic approach. We evaluated and compared the geometries and energetics of hydride-transfer process using two different systems. We highlighted chemical interactions and studied the effect of protonation state of an active-site histidine base on the energetics of transfer. Furthermore, we evaluated energetics of the second hydride-transfer process as well as hydration of betaine aldehyde.

| S-EPMC6643517 | biostudies-literature

N2H2 binding to the nitrogenase FeMo cluster studied by QM/MM methods.

Project description:We have made a systematic combined quantum mechanical and molecular mechanical (QM/MM) investigation of possible structures of the N2 bound state of nitrogenase. We assume that N2 is immediately protonated to a N2H2 state, thereby avoiding the problem of determining the position of the protons in the cluster. We have systematically studied both end-on and side-on structures, as well as both HNNH and NNH2 states. Our results indicate that the binding of N2H2 is determined more by interactions and steric clashes with the surrounding protein than by the intrinsic preferences of the ligand and the cluster. The best binding mode with both the TPSS and B3LYP density-functional theory methods has trans-HNNH terminally bound to Fe2. It is stabilised by stacking of the substrate with His-195 and Ser-278. However, several other structures come rather close in energy (within 3-35 kJ/mol) at least in some calculations: The corresponding cis-HNNH structure terminally bound to Fe2 is second best with B3LYP. A structure with HNNH2 terminally bound to Fe6 is second most stable with TPSS (where the third proton is transferred to the substrate from the homocitrate ligand). Structures with trans-HNNH, bound to Fe4 or Fe6, or cis-HNNH bound to Fe6 are also rather stable. Finally, with the TPSS functional, a structure with cis-HNNH side-on binding to the Fe3-Fe4-Fe5-Fe7 face of the cluster is also rather low in energy, but all side-on structures are strongly disfavoured by the B3LYP method.

| S-EPMC7186253 | biostudies-literature

A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase.

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

How Large Should the QM Region Be in QM/MM Calculations? The Case of Catechol O-Methyltransferase.

Project description:Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.

| S-EPMC5108028 | biostudies-literature

QM/MM calculations with deMon2k.

Project description:The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.

| S-EPMC6272552 | biostudies-literature

Immersive virtual reality in computational chemistry: Applications to the analysis of QM and MM data.

Project description:The role of Virtual Reality (VR) tools in molecular sciences is analyzed in this contribution through the presentation of the Caffeine software to the quantum chemistry community. Caffeine, developed at Scuola Normale Superiore, is specifically tailored for molecular representation and data visualization with VR systems, such as VR theaters and helmets. Usefulness and advantages that can be gained by exploiting VR are here reported, considering few examples specifically selected to illustrate different level of theory and molecular representation.

| S-EPMC5101850 | biostudies-literature

QM/MM metadynamics study of the direct decarboxylation mechanism for orotidine-5'-monophosphate decarboxylase using two different QM regions: acceleration too small to explain rate of enzyme catalysis.

Project description:Despite decades of study, the mechanism by which orotidine-5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine monophosphate remains unresolved. A computational investigation of the direct decarboxylation mechanism has been performed using mixed quantum mechanical/molecular mechanical (QM/MM) dynamics simulations. The study was performed with the program CP2K that integrates classical dynamics and ab initio dynamics based on the Born-Oppenheimer approach. Two different QM regions were explored. The free energy barriers for direct decarboxylation of orotidine-5'-monophosphate (OMP) in solution and in the enzyme (using the larger QM region) were determined with the metadynamics method to be 40 and 33 kcal/mol, respectively. The calculated change in activation free energy (DeltaDeltaG++) on going from solution to the enzyme is therefore -7 kcal/mol, far less than the experimental change of -23 kcal/ mol (for k(cat.)/k(uncat.): Radzicka, A.; Wolfenden, R., Science 1995, 267, 90-92). These results do not support the direct decarboxylation mechanism that has been proposed for the enzyme. However, in the context of QM/MM calculations, it was found that the size of the QM region has a dramatic effect on the calculated reaction barrier.

| S-EPMC3163499 | biostudies-literature

CHARMM-GUI QM/MM Interfacer for a Quantum Mechanical and Molecular Mechanical (QM/MM) Simulation Setup: 1. Semiempirical Methods.

Project description:Quantum mechanical (QM) treatments, when combined with molecular mechanical (MM) force fields, can effectively handle enzyme-catalyzed reactions without significantly increasing the computational cost. In this context, we present CHARMM-GUI QM/MM Interfacer, a web-based cyberinfrastructure designed to streamline the preparation of various QM/MM simulation inputs with ligand modification. The development of QM/MM Interfacer has been achieved through integration with existing CHARMM-GUI modules, such as PDB Reader and Manipulator, Solution Builder, and Membrane Builder. In addition, new functionalities have been developed to facilitate the one-stop preparation of QM/MM systems and enable interactive and intuitive ligand modifications and QM atom selections. QM/MM Interfacer offers support for a range of semiempirical QM methods, including AM1(+/d), PM3(+/PDDG), MNDO(+/d, +/PDDG), PM6, RM1, and SCC-DFTB, tailored for both AMBER and CHARMM. A nontrivial setup related to ligand modification, link-atom insertion, and charge distribution is automatized through intuitive user interfaces. To illustrate the robustness of QM/MM Interfacer, we conducted QM/MM simulations of three enzyme-substrate systems: dihydrofolate reductase, insulin receptor kinase, and oligosaccharyltransferase. In addition, we have created three tutorial videos about building these systems, which can be found at https://www.charmm-gui.org/demo/qmi. QM/MM Interfacer is expected to be a valuable and accessible web-based tool that simplifies and accelerates the setup process for hybrid QM/MM simulations.

| S-EPMC11209942 | biostudies-literature

QM/MM Benchmarking of Cyanobacteriochrome Slr1393g3 Absorption Spectra.

Project description:Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that are promising candidates for biotechnological applications. Computational studies can contribute to an understanding at a molecular level of their wide spectral tuning and diversity. In this contribution, we benchmark methods to model a 110 nm shift in the UV/Vis absorption spectrum from a red- to a green-absorbing form of the cyanobacteriochrome Slr1393g3. Based on an assessment of semiempirical methods to describe the chromophore geometries of both forms in vacuo, we find that DFTB2+D leads to structures that are the closest to the reference method. The benchmark of the excited state calculations is based on snapshots from quantum mechanics/molecular mechanics molecular dynamics simulations. In our case, the methods RI-ADC(2) and sTD-DFT based on CAM-B3LYP ground state calculations perform the best, whereas no functional can be recommended to simulate the absorption spectra of both forms with time-dependent density functional theory. Furthermore, the difference in absorption for the lowest energy absorption maxima of both forms can already be modelled with optimized structures, but sampling is required to improve the shape of the absorption bands of both forms, in particular for the second band. This benchmark study can guide further computational studies, as it assesses essential components of a protocol to model the spectral tuning of both cyanobacteriochromes and the related phytochromes.

| S-EPMC6540152 | biostudies-literature

QM/MM X-ray refinement of zinc metalloenzymes.

Project description:Zinc metalloenzymes play an important role in biology. However, due to the limitation of molecular force field energy restraints used in X-ray refinement at medium or low resolutions, the precise geometry of the zinc coordination environment can be difficult to distinguish from ambiguous electron density maps. Due to the difficulties involved in defining accurate force fields for metal ions, the QM/MM (quantum-mechanical/molecular-mechanical) method provides an attractive and more general alternative for the study and refinement of metalloprotein active sites. Herein we present three examples that indicate that QM/MM based refinement yields a superior description of the crystal structure based on R and R(free) values and on the inspection of the zinc coordination environment. It is concluded that QM/MM refinement is an useful general tool for the improvement of the metal coordination sphere in metalloenzyme active sites.

| S-EPMC2843799 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data