Spin-projected QM/MM Free Energy Simulations for Oxidation Reaction of Guanine in B-DNA by Singlet Oxygen.
Ontology highlight
ABSTRACT:
SUBMITTER: Saito T
PROVIDER: S-EPMC8048875 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A projected hybrid orbital (PHO) method was described to model the covalent boundary in a hybrid quantum mechanical and molecular mechanical (QM/MM) system. The PHO approach can be used in ab initio wave function theory and in density functional theory with any basis set without introducing system-dependent parameters. In this method, a secondary basis set on the boundary atom is introduced to formulate a set of hybrid atomic orbtials. The primary basis set on the boundary atom used for the QM subsystem is projected onto the secondary basis to yield a representation that provides a good approximation to the electron-withdrawing power of the primary basis set to balance electronic interactions between QM and MM subsystems. The PHO method has been tested on a range of molecules and properties. Comparison with results obtained from QM calculations on the entire system shows that the present PHO method is a robust and balanced QM/MM scheme that preserves the structural and electronic properties of the QM region.
| S-EPMC4306490 | biostudies-literature
Project description:Indoleamine 2,3-dioxygenase (IDO) and tryptophan dioxygenase (TDO) are two heme proteins that catalyze the oxidation reaction of tryptophan (Trp) to N-formylkynurenine (NFK). Human IDO (hIDO) has recently been recognized as a potent anticancer drug target, a fact that triggered intense research on the reaction and inhibition mechanisms of hIDO. Our recent studies revealed that the dioxygenase reaction catalyzed by hIDO and TDO is initiated by addition of the ferric iron-bound superoxide to the C(2)?C(3) bond of Trp to form a ferryl and Trp-epoxide intermediate, via a 2-indolenylperoxo radical transition state. The data demonstrate that the two atoms of dioxygen are inserted into the substrate in a stepwise fashion, challenging the paradigm of heme-based dioxygenase chemistry. In the current study, we used QM/MM methods to decipher the mechanism by which the second ferryl oxygen is inserted into the Trp-epoxide to form the NFK product in hIDO. Our results show that the most energetically favored pathway involves proton transfer from Trp-NH(3)(+) to the epoxide oxygen, triggering epoxide ring opening and a concerted nucleophilic attack of the ferryl oxygen to the C(2) of Trp that leads to a metastable reaction intermediate. This intermediate subsequently converts to NFK, following C(2)-C(3) bond cleavage and the associated back proton transfer from the oxygen to the amino group of Trp. A comparative study with Xantomonas campestris TDO (xcTDO) indicates that the reaction follows a similar pathway, although subtle differences distinguishing the two enzyme reactions are evident. The results underscore the importance of the NH(3)(+) group of Trp in the two-step ferryl-based mechanism of hIDO and xcTDO, by acting as an acid catalyst to facilitate the epoxide ring-opening reaction and ferryl oxygen addition to the indole ring.
| S-EPMC3304497 | biostudies-literature
Project description:Human arginase I (HARGI) is a metalloprotein highly expressed in the liver cytosol and catalyzes the hydrolysis of l-arginine to form l-ornithine and urea. Understanding the reaction mechanism would be highly helpful to design new inhibitor molecules for HARGI as it is a target for heart- and blood-related diseases. In this study, we explored the hydrolysis reaction mechanism of HARGI with antiferromagnetic and ferromagnetic coupling between two Mn(II) ions at the catalytic site by employing molecular dynamics simulations coupled with quantum mechanics and molecular mechanics (QM/MM). The spin states, high-spin ferromagnetic couple (S Mn1 = 5/2, S Mn2 = 5/2), low-spin ferromagnetic couple (S Mn1 = 1/2, S Mn2 = 1/2), high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2), and low-spin antiferromagnetic couple (S Mn1 = 1/2, S Mn2 = -1/2) are considered, and the calculated energetics for the complex of the substrate and HARGI are compared. The results show that the high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) is more stable than other spin states. The low-spin ferromagnetic and antiferromagnetic coupled states are highly unstable compared with the corresponding high-spin states. The high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) is stabilized by 0.39 kcal/mol compared with the ferromagnetic couple (S Mn1 = 5/2, S Mn2 = 5/2). The reaction mechanism is independent of spin states; however, the energetics of transition states and intermediates are more stable in the case of the high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) than the corresponding ferromagnetic state. It is evident that the calculated coupling constants are higher for antiferromagnetic states and, interestingly, superexchange coupling is found to occur between Mn(II) ions via hydroxide ions in a reactant. The hydroxide ion enhances the coupling interaction and initiates the catalytic reaction. It is also noted that the first intermediate structure where there is no superexchange coupling is similar to the known inhibitor 2(S)-amino-6-boronohexanoic acid.
| S-EPMC9475637 | biostudies-literature
Project description:Among all bioluminescent organisms, the firefly is the most famous, with a high luminescent efficiency of 41%, which is widely used in the fields of biotechnology, biomedicine and so on. The entire bioluminescence (BL) process involves a series of complicated in-vivo chemical reactions. The BL is initiated by the enzymatic oxidation of luciferin (LH2). However, the mechanism of the efficient spin-forbidden oxygenation is far from being totally understood. Via MD simulation and QM/MM calculations, this article describes the complete process of oxygenation in real protein. The oxygenation of luciferin is initiated by a single electron transfer from the trivalent anionic LH2 (L3-) to O2 to form 1[L•2-…O2•-]; the entire reaction is carried out along the ground-state potential energy surface to produce the dioxetanone (FDO-) via three transition states and two intermediates. The low energy barriers of the oxygenation reaction and biradical annihilation involved in the reaction explain this spin-forbidden reaction with high efficiency. This study is helpful for understanding the BL initiation of fireflies and the other oxygen-dependent bioluminescent organisms.
| S-EPMC8307557 | 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:Bioluminescence of a variety of marine organisms, mostly cnidarians and ctenophores, is carried out by Ca2+-dependent photoproteins. The mechanism of light emission operates via the same reaction in both animal families. Despite numerous studies on the ctenophore photoprotein family, the detailed catalytic mechanism and arrangement of amino acid residues surrounding the chromophore in this family are a mystery. Here, we report the crystal structure of Cd2+-loaded apo-mnemiopsin1, a member of the ctenophore family, at 2.15 Å resolution and used quantum mechanics/molecular mechanics (QM/MM) to investigate its reaction mechanism. The simulations suggested that an Asp-156-Arg-39-Tyr-202 triad creates a hydrogen-bonded network to facilitate the transfer of a proton from the 2-hydroperoxy group of the chromophore coelenterazine to bulk solvent. We identified a water molecule in the coelenteramide-binding cavity that forms a hydrogen bond with the amide nitrogen atom of coelenteramide, which, in turn, is hydrogen-bonded via another water molecule to Tyr-131. This observation supports the hypothesis that the function of the coelenteramide-bound water molecule is to catalyze the 2-hydroperoxycoelenterazine decarboxylation reaction by protonation of a dioxetanone anion, thereby triggering the bioluminescence reaction in the ctenophore photoprotein family.
| S-EPMC6322872 | biostudies-literature
Project description:Singlet oxygen (1O2) plays a pivotal role in numerous catalytic oxidation processes utilized in water purification and chemical synthesis. The spin-trapping method based on electron paramagnetic resonance (EPR) analysis is commonly employed for 1O2 detection. However, it is often limited to time-independent acquisition. Recent studies have raised questions about the reliability of the 1O2 trapper, 2,2,6,6-tetramethylpiperidine (TEMP), in various systems. In this study, we introduce a comprehensive, kinetic examination to monitor the spin-trapping process in EPR analysis. The EPR intensity of the trapping product was used as a quantitative measurement to evaluate the concentration of 1O2 in aqueous systems. This in situ kinetic study was successfully applied to a classical photocatalytic system with exceptional accuracy. Furthermore, we demonstrated the feasibility of our approach in more intricate 1O2-driven catalytic oxidation processes for water decontamination and elucidated the molecular mechanism of direct TEMP oxidation. This method can avoid the false-positive results associated with the conventional 2D 1O2 detection techniques, and provide insights into the reaction mechanisms in 1O2-dominated catalytic oxidation processes. This work underscores the necessity of kinetic studies for spin-trapping EPR analysis, presenting an avenue for a comprehensive exploration of the mechanisms governing catalytic oxidation processes.
| S-EPMC10372693 | biostudies-literature
Project description:We present a new development in quantum mechanics/molecular mechanics (QM/MM) methods by replacing conventional MM models with data-driven many-body (MB) representations rigorously derived from high-level QM calculations. The new QM/MM approach builds on top of mutually polarizable QM/MM schemes developed for polarizable force fields with inducible dipoles and uses permutationally invariant polynomials to effectively account for quantum-mechanical contributions (e.g., exchange-repulsion and charge transfer and penetration) that are difficult to describe by classical expressions adopted by conventional MM models. Using the many-body MB-pol and MB-DFT potential energy functions for water, which include explicit two-body and three-body terms fitted to reproduce the corresponding CCSD(T) and PBE0 two-body and three-body energies for water, we demonstrate a smooth energetic transition as molecules are transferred between QM and MM regions, without the need of a transition layer. By effectively elevating the accuracy of both the MM region and the QM/MM interface to that of the QM region, the new QM/MB-MM approach achieves an accuracy comparable to that obtained with a fully QM treatment of the entire system.
| S-EPMC8131112 | biostudies-literature
Project description:The use of the most accurate (i.e., QM or QM/MM) levels of theory for free energy simulations (FES) is typically not possible. Primarily, this is because the computational cost associated with the extensive configurational sampling needed for converging FES is prohibitive. To ensure the feasibility of QM-based FES, the "indirect" approach is generally taken, necessitating a free energy calculation between the MM and QM/MM potential energy surfaces. Ideally, this step is performed with standard free energy perturbation (Zwanzig's equation) as it only requires simulations be carried out at the low level of theory; however, work from several groups over the past few years has conclusively shown that Zwanzig's equation is ill-suited to this task. As such, many approximations have arisen to mitigate difficulties with Zwanzig's equation. One particularly popular notion is that the convergence of Zwanzig's equation can be improved by using interaction energy differences instead of total energy differences. Although problematic numerical fluctuations (a major problem when using Zwanzig's equation) are indeed reduced, our results and analysis demonstrate that this "interaction energy approximation" (IEA) is theoretically incorrect, and the implicit approximation invoked is spurious at best. Herein, we demonstrate this via solvation free energy calculations using IEA from two different low levels of theory to the same target high level. Results from this proof-of-concept consistently yield the wrong results, deviating by ∼1.5 kcal/mol from the rigorously obtained value.
| S-EPMC6745021 | biostudies-literature
Project description:Taking long-range electrostatic effects into account in classical and hybrid quantum mechanics-molecular mechanics (QM/MM) simulations is necessary for an accurate description of the system under study. We have recently developed a method, termed long-range electrostatic corrections (LREC), for monopolar QM/MM calculations. Here, we present an extension of LREC for multipolar/polarizable QM/MM simulations within the LICHEM software package. Reaction barriers and QM-MM interaction energies converge with a LREC cutoff between 20 and 25 Å, in agreement with our previous results. Additionally, the LREC approach for the QM-MM interactions can be smoothly combined with standard shifting or Ewald summation methods in the MM calculations. We recommend the use of QM(LREC)/MM(PME), where the QM region is treated with LREC and the MM region is treated with particle mesh Ewald (PME) summation. This combination is an excellent compromise between simplicity, speed, and accuracy for large QM/MM simulations.
| S-EPMC5373107 | biostudies-literature