Combined Quantum Research Articles

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an SN2 Reaction in Water.

We present a tutorial to carry out umbrella-sampling free-energy simulations with a combined quantum mechanical and molecular mechanical (QM/MM) potential, which may also be used in a computational or biophysical chemistry curriculum for first-year graduate and undergraduate students. In this article, we choose the Type II SN2 Menshutkin reaction between ammonia and chloromethane to construct the potential of mean force (PMF) for the reaction in aqueous solution. In this exercise, we wish to accomplish three tasks: (1) an understanding of the concept of PMF and the umbrella-sampling free-energy simulation method, (2) the use of a combined QM/MM potential in molecular dynamics simulation of chemical reactions, and (3) an understanding of solvent effects and intermolecular interactions on chemical reactions through comparison with gas-phase results. Analysis of the simulation results allows students to appreciate the difference between the transition state along a PMF from statistical simulations versus the optimized saddle point in the gas phase, the shift of transition state with respect to solvation, and the significance of electronic polarization. Through this tutorial, students will gain the necessary skills to carry out free energy simulations of chemical reactions in solution and in enzymes both in a classroom and in the laboratory.

Gas-liquid interfacial reaction mechanisms of typical small α-dicarbonyls in the neutral and acidic droplets: Implications for secondary organic aerosol formation

Small α-dicarbonyls (SαDs) are well-known as the important precursors of secondary organic aerosol (SOA). Hence, it is imperative to understand the atmospheric chemistry of SαDs to contribute to SOA formation. In this work, we investigated the interfacial chemistry of typical SαDs, including methylglyoxal (MG) and biacetyl (BA) in the neutral and acidic droplets by combined molecular dynamics and quantum chemical calculations. The trans configurations of MG and BA are found to be the favorable configurations at the interfaces and are prone to stay at the gas-liquid interface of the acidic droplet. The C=O group exhibits a preferential uptake orientation towards the interface because the carbonyl-O atom has a strong interaction with interfacial H2O. The uptakes and accumulations of MG and BA at the interfaces are promoted by the acidic condition. Subsequent interfacial hydrations of MG and BA in the acidic droplet are beneficial to yield diols, which can engage in oligomerization in the droplet interior to contribute SOA formation. Our results provide the theoretical insight into the interfacial chemistry of SαDs and their role in SOA formation.

Observation of a V-Shape Superconductivity Evolution on Tungsten-Intercalated 2H-Type Niobium Diselenide.

van der Waals-layered niobium diselenide (NbSe2) intercalated by d-electron transition metals is an ideal test bed for the exploration of their diversiform evolution of ground states. These intercalations are mostly viewed as ordered structures aligned with periodicities of their host materials that enable control of the electronic phases via gradually changing of intercalation ratios. Here, we present the structure and superconductivity in tungsten (W)-intercalated 2H-NbSe2 crystals, which reveals an order to disorder distribution of W atoms with increasing confined intercalating amounts, leading to an approximate V-shape suppression of superconductivity. Aided by density functional theory calculations, we demonstrate that the local magnetic moment around W intercalants induced by the charge redistribution gives rise to the quick superconductivity suppression in 2H-NbSe2 below a certain dilute amount (W% = 0.06). Simultaneously, W intercalants also induce structural aberration due to aggregation effects and inhibit the generation of an ordered structure in 2H-NbSe2, resulting in a recovery of its superconductivity. The alteration of structure and electronic phases in 2H-NbSe2 via intercalation of nonmagnetic transition metals in the van der Waals gap enables the exploration of combined magnetic quantum criticality, superconductivity, and other related electronic correlations.

Observation of Circular Dichroism Induced by Electronic Chirality.

Chiral phases of matter, characterized by a definite handedness, abound in nature, ranging from the crystal structure of quartz to spiraling spin states in helical magnets. In 1T-TiSe_{2} a source of chirality has been proposed that stands apart from these classical examples as it arises from combined electronic charge and quantum orbital fluctuations. This may allow its chirality to be accessed and manipulated without imposing either structural or magnetic handedness. However, direct bulk evidence that broken inversion symmetry and chirality are intrinsic to TiSe_{2} remains elusive. Here, employing resonant elastic x-ray scattering technique, we reveal the presence of circular dichroism, i.e., polarization dependence of the resonant diffraction intensity, up to ∼40% at forbidden Bragg peaks that emerge at the charge and orbital ordering transition. The dichroism varies dramatically with incident energy and azimuthal angle. Comparison to calculated scattering intensities traces its origin to bulk chiral electronic order in TiSe_{2} and establishes resonant elastic x-ray scattering as a sensitive probe to electronic chirality.

Molecular modelling of luciferyl adenylate deprotonation in the active site of Photinus pyralis luciferase

Firefly bioluminescence is a product of chemical reactions that involve luciferin chromophore oxidation in the active site of luciferase proteins. We perform a series of classical molecular dynamic simulations and combined quantum and molecular mechanics (QM/MM) calculations to expose the molecular mechanism of C4 carbon atom deprotonation in luciferyl adenylate molecule. QM/MM calculations confirm that ND-protonated His245 residue is a suitable proton acceptor in the wild-type Photinus pyralis luciferase. Classical molecular dynamic simulations reveal oxygen-binding cavities inside the protein including the one located close to the C4 atom of luciferin. In mutant forms that lack direct interactions with the His245 side chain, a proton wire comprising water molecules stimulates either protonation of the phosphate group of the luciferyl adenylate forming an unstable intermediate or keto–enol tautomerisation of luciferin. The comparison of the ionisation potentials of molecular systems with the ‘deprotonated’ C4 carbon atom revealed that the ionisation energy for the enol tautomeric form is close to the system with the His245 proton acceptor. Thus, the existence of the keto–enol tautomerisation channel might explain bioluminescence in case of the absence of the amino acid proton acceptor.

Identification and molecular mechanism of novel antioxidant peptide from fish sauce: A combined quantum chemistry and molecular simulation

Fish sauce, derived from fermented fish, exhibits a notable antioxidant effect after a six-month fermentation process, and we propose that potential antioxidant peptides were present in the fish sauce. We isolated, purified, and identified potential bioactive antioxidant peptides by using fish sauce fermented for 6 months. Additionally, molecular simulation was employed to investigate the antioxidant action mechanism of these bioactive peptides. The molecular docking results revealed that FS4-1 (MHQLSKK), FS4-2 (VLDNSPER), FS4-3 (MNPPAASIK), FS6-1(VLKQAAAGR), and FS6-2 (SPDVSPRR), could dock with the Keap1 receptor. The primary force (Van der Waals' force and hydrogen bonds) and key sites (GLY509 and ALA510) of Keap1 binding to peptides were determined. The active center was located in the side chain of amino acid Met at positions C7H78 and C7H79. We here identified antioxidant peptides in fish sauce and revealed the antioxidant mechanism through molecular simulations.

Monte Carlo QM/MM simulation studies of the Cannizzaro reaction in ionic liquids for improved biofuel production.

The conversion of biomass to 5-hydroxymethylfurfural (HMF) holds substantial promise as a renewable energy source. Notably, HMF can be transformed into 2,5-bis(hydroxymethyl)furan (BHMF), a crucial reactant in biofuel production, but requires harsh operating conditions, H2, and precious metal catalysts. A recently reported Cannizzaro reaction of HMF to BHMF, characterized by its efficiency, mild conditions, and eco-friendliness, instead employed ionic liquids (ILs) to achieve high yields. In this study, combined quantum mechanical and molecular mechanical (QM/MM) simulations in conjunction with Metropolis Monte Carlo statistical mechanics and free-energy perturbation theory utilized M06-2X/6-31+G(d), PDDG/PM3, and the OPLS-VSIL force field to uncover important solute-solvent interactions present in the HMF to BHMF reaction pathway. The Cannizzaro reaction was examined in water and in five ILs composed of the 1-butyl-3-methylimidazolium [BMIM] cation coupled to hexafluorophosphate, tetrafluoroborate, thiocyanate, chloride, and bromide. Energetic and structural analysis of the rate-determining hydride transfer between HMF and the hydride-donor anion HMFOH- attributed the enhanced reactivity to highly organized solvent interactions featuring (1) hydrogen bonding between the ring protons of [BMIM] and the negatively charged carbonyl oxygen atoms on the transition structure, (2) favorable electrostatic interactions between the IL anions and solute hydroxyl groups, and (3) beneficial π-π stacking interactions between [BMIM] and the two aromatic rings present in HMF and HMFOH-.

Reaction Mechanism for CO Reduction by Mo-Nitrogenase Studied by QM/MM.

We have studied the conversion of two molecules of carbon monoxide to ethylene catalyzed by nitrogenase. We start from a recent crystal structure showing the binding of two carbon monoxide molecules to nitrogenase and employ the combined quantum mechanics and molecular mechanics approach. Our results indicate that the reaction is possible only if S2B dissociates as H2S (i.e., the charge of the FeMo cluster remains the same as in the E0 state, indicating that the Fe ions are formally reduced two steps when CO binds). Eight electrons and protons are needed for the reaction, and our mechanism suggests that the first four bind alternatively to the two carbon atoms. The C-C bond formation takes place already after the first protonation (in the E3 state). The next two protons bind to the same O atom, which then dissociates as water. In the same state (E8), the second C-O bond is cleaved, forming the ethylene product. The last two electrons and protons are used to form a water molecule that can be exchanged by S2B or by two CO molecules to start a new reaction cycle.

Photoelectron Spectroscopy and Computational Study on Microsolvated [B10H10]2- Clusters and Comparisons to Their [B12H12]2- Analogues.

Microhydrated closo-boranes have attracted great interest due to their superchaotropic activity related to the well-known Hofmeister effect and important applications in biomedical and battery fields. In this work, we report a combined negative ion photoelectron spectroscopy and quantum chemical investigation on hydrated closo-decaborate clusters [B10H10]2-·nH2O (n = 1-7) with a direct comparison to their analogues [B12H12]2-·nH2O and free water clusters. A single H2O molecule is found to be sufficient to stabilize the intrinsically unstable [B10H10]2- dianion. The first two water molecules strongly interact with the solute forming B-H···H-O dihydrogen bonds while additional water molecules show substantially reduced binding energies. Unlike [B12H12]2-·nH2O possessing a highly structured water network with the attached H2O molecules arranged in a unified pattern by maximizing B-H···H-O dihydrogen bonding, distinct structural arrangements of the water clusters within [B10H10]2-·nH2O are achieved with the water cluster networks from trimer to heptamer resembling free water clusters. Such a distinct difference arises from the variations in size, symmetry, and charge distributions between these two dianions. The present finding again confirms the structural diversity of hydrogen-bonding networks in microhydrated closo-boranes and enriches our understanding of aqueous borate chemistry.

The combined effects of quantum and strong coupling on the nonlinear collective excitation in quantum dusty plasmas

The quantum hydrodynamic model for electrons and ions and the generalized hydrodynamic model for the strongly coupled dust particles are proposed in the strongly coupled quantum dusty plasma, where the combined quantum effects of quantum diffraction, quantum statistic pressure, as well as electron exchange and correlation effects are all considered in the quantum hydrodynamic model. The shear and bulk viscosity effects are included in the viscoelastic relaxation, which leads to the decay of the dust-ion-acoustic waves. The approximate time-dependent solitary solution is obtained by the momentum conservation law in the presence of viscosity.

Pentaphosphorylation via the Anhydride of Dihydrogen Pentametaphosphate: Access to Nucleoside Hexa- and Heptaphosphates and Study of Their Interaction with RibonucleaseA.

Pentametaphosphate is the little studied cyclic pentamer of the metaphosphate ion, [PO3]5 5-. We show that the doubly protonated form of this pentamer can be selectively dehydrated to provide the anhydride [P5O14]3- (1). This trianion is the well-defined condensed phosphate component of a novel reagent for attachment of a pentaphosphate chain to biomolecules all in one go. Here, we demonstrate by extending adenosine monophosphate (AMP) and uridine monophosphate (UMP) to their corresponding nucleoside hexaphosphates, while adenosine diphosphate (ADP) and uridine diphosphate (UDP) are phosphate chain-extended to the corresponding nucleoside heptaphosphates. Such constructs are of interest for their potential biological function with respect to RNA-processing enzymes. Thus, we go on to investigate in detail the interaction of the polyanionic constructs with ribonuclease A, a model protein containing a polycationic active site and for which X-ray crystal structures are relatively straightforward to obtain. This work presents a combined experimental and quantum chemical approach to understanding the interactions of RNase A with the new nucleoside hexa- and heptaphosphate constructs.

Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study.

Hydrogen bonding is a crucial feature of biomolecules, but its characterization in glycans dissolved in aqueous solutions is challenging due to rapid hydrogen exchange between hydroxyl groups and H2O. In principle, the scalar (J) coupling constant can reveal the relative orientation of the atoms in the molecule. In contrast to J-coupling through H-bonds reported in proteins and nucleic acids, research on J-coupling through H-bonds in glycans dissolved in water is lacking. Here, we use sucrose as a model system for H-bonding studies; its structure, which consists of glucose (Glc) and fructose (Frc), is well-studied, and it is readily available. We apply the in-phase, antiphase-HSQC-TOCSY and quantify previously unreported through H-bond J-values for Frc-OH1-Glc-OH2 in H2O. While earlier reports of Brown and Levy indicate this H-bond as having only a single direction, our reported findings indicate the potential presence of two involving these same atoms, namely, G2OH➔F1O and F1OH➔G2O (where F and G stand for Frc and Glc, respectively). The calculated density functional theory J-values for the G2OH➔F1O agree with the experimental values. Additionally, we detected four other possible H-bonds in sucrose, which require different phi, psi (ϕ, ψ) torsion angles. The ϕ, ψ values are consistent with previous predictions of du Penhoat et al. and Venable et al. Our results will provide new insights into the molecular structure of sucrose and its interactions with proteins.

Putative reaction mechanism of nitrogenase with a half-dissociated S2B ligand.

We have studied whether dissociation of the S2B sulfide ligand from one of its two coordinating Fe ions may affect the later parts of the reaction mechanism of nitrogenase. Such dissociation has been shown to be favourable for the E2-E4 states in the reaction mechanism, but previous studies have assumed that S2B either remains bridging or has fully dissociated from the active-site FeMo cluster. We employ combined quantum mechanical and molecular mechanical (QM/MM) calculations with two density-functional theory methods, r2SCAN and TPSSh. To make dissociation of S2B possible, we have added a proton to this group throughout the reaction. We study the reaction starting from the E4 state with N2H2 bound to the cluster. Our results indicate that half-dissociation of S2B is unfavourable in most steps of the reaction mechanism. We observe favourable half-dissociation of S2B only when NH or NH2 is bound to the cluster, bridging Fe2 and Fe6. However, the former state is most likely not involved in the reaction mechanism and the latter state is only an intermittent intermediate of the E7 state. Therefore, half-dissociation of S2B seems to play only a minor role in the later parts of the reaction mechanism of nitrogenase. Our suggested mechanism with a protonated S2B is alternating (the two N atoms of the substrate is protonated in an alternating manner) and the substrate prefers to bind to Fe2, in contrast to the preferred binding to Fe6 observed when S2B is unprotonated and bridging Fe2 and Fe6.

STRUCTURE AND DYNAMICS OF THE ENZYME-SUBSTRATE COMPLEX OF N-ACETYLASPARTYLGLUTAMATE SYNTHASE ACCORDING TO COMPUTER SIMULATION DATA

N-acetylaspartilglutamate is the most common dipeptide in brain cells, which is synthesized using the enzyme N-acetylaspartilglutamate synthase. Herein we utilize bioinformatics methods to predict the protein structure from the primary sequence of the coding gene, classical molecular dynamics to obtain a stable protein complex with N-acetylaspartate and glutamate ligands within the trajectory, as well as machine learning methods to analyze, describe and select potential reactive and non-reactive conformations of the model system describing the enzyme-substrate complex. Molecular dynamics simulations with combined quantum mechanics / molecular mechanics potentials were performed for a set of selected conformations and the potential reaction mechanism were characterized.

QM/MM Simulations of Afatinib-EGFR Addition: The Role of β-Dimethylaminomethyl Substitution.

Acrylamides are the most commonly used warheads of targeted covalent inhibitors (TCIs) directed at cysteines; however, the reaction mechanisms of acrylamides in proteins remain controversial, particularly for those involving protonated or unreactive cysteines. Using the combined semiempirical quantum mechanics (QM)/molecular mechanics (MM) free energy simulations, we investigated the reaction between afatinib, the first TCI drug for cancer treatment, and Cys797 in the EGFR kinase. Afatinib contains a β-dimethylaminomethyl (β-DMAM) substitution which has been shown to enhance the intrinsic reactivity and potency against EGFR for related inhibitors. Two hypothesized reaction mechanisms were tested. Our data suggest that Cys797 becomes deprotonated in the presence of afatinib, and the reaction proceeds via a classical Michael addition mechanism, with Asp800 stabilizing the ion-pair reactant state β-DMAM+/C797- and the transition state of the nucleophilic attack. Our work elucidates an important structure-activity relationship of acrylamides in proteins.

Hydrostatic pressure effect on excited state properties of room temperature phosphorescence molecules: A QM/MM study

Stimulus-responsive organic room temperature phosphorescence (RTP) materials exhibit variations in their luminescent characteristics (lifetime and efficiency) upon exposure to external stimuli, including force, heat, light and acid-base conditions, the development of stimulus-responsive RTP molecules becomes imperative. However, the inner responsive mechanism is unclear, theoretical investigations to reveal the relationship among hydrostatic pressures, molecular structures and photophysical properties are highly desired. Herein, taking the Se-containing RTP molecule (SeAN) as a model, based on the dispersion corrected density functional theory (DFT-D), the combined quantum mechanics and molecular dynamics (QM/MM) method and thermal vibration correlation function (TVCF) theory, the influences of hydrostatic pressure on molecular structures, transition properties as well as lifetimes and efficiencies of RTP molecule are theoretically studied. Results show that extended lifetime and enhanced efficiency are observed at 2 Gpa compared with molecule at normal pressure, and this is related with the small reorganization energy and large oscillator strength. Moreover, due to the small energy gap (0.34 eV) and remarkable spin-orbit coupling (SOC) constant (8.56 cm-1) between first singlet excited state and triplet state, fast intersystem crossing (ISC) process is determined for molecule at 6 Gpa. Furthermore, the intermolecular interactions are visualized using independent gradient model based on Hirshfeld partition (IGMH) and the changes of molecular packing modes, SOC values, lifetimes and efficiencies with pressures are detected. These results reveal the relationship between molecular structures and RTP properties. Our work provides theoretical insights into the hydrostatic pressure response mechanism and could promote the development new efficient stimulus-responsive molecules.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials.

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

The Unique Role of the Second Coordination Sphere to Unlock and Control Catalysis in Nonheme Fe(II)/2-Oxoglutarate Histone Demethylase KDM2A.

Nonheme Fe(II) and 2-oxoglutarate (2OG)-dependent histone lysine demethylases 2A (KDM2A) catalyze the demethylation of the mono- or dimethylated lysine 36 residue in the histone H3 peptide (H3K36me1/me2), which plays a crucial role in epigenetic regulation and can be involved in many cancers. Although the overall catalytic mechanism of KDMs has been studied, how KDM2 catalysis takes place in contrast to other KDMs remains unknown. Understanding such differences is vital for enzyme redesign and can help in enzyme-selective drug design. Herein, we employed molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) to explore the complete catalytic mechanism of KDM2A, including dioxygen diffusion and binding, dioxygen activation, and substrate oxidation. Our study demonstrates that the catalysis of KDM2A is controlled by the conformational change of the second coordination sphere (SCS), specifically by a change in the orientation of Y222, which unlocks the 2OG rearrangement from off-line to in-line mode. The study demonstrates that the variant Y222A makes the 2OG rearrangement more favorable. Furthermore, the study reveals that it is the size of H3K36me3 that prevents the 2OG rearrangement, thus rendering the enzyme inactivity with trimethylated lysine. Calculations show that the SCS and long-range interacting residues that stabilize the HAT transition state in KDM2A differ from those in KDM4A, KDM7B, and KDM6A, thus providing the basics for the enzyme-selective redesign and modulation of KDM2A without influencing other KDMs.

Theoretical Insights into N-Glycoside Bond Cleavage of 5-Carboxycytosine by Thymine DNA Glycosylase: A QM/MM Study.

Thymine DNA glycosylase (TDG)-mediated excision of 5-formylcytosine and 5-carboxylcytosine (5-caC) is a critical step in active DNA demethylation. Herein, we employed a combined quantum mechanics/molecular mechanics approach to investigate the reaction mechanism of TDG-catalyzed N-glycosidic bond cleavage of 5-caC. The calculated results show that TDG-catalyzed 5-caC excision follows a concerted (SN2) mechanism in which glycosidic bond dissociation is coupled with nucleophile attack. Protonation of the 5-caC anion contributes to the cleavage of the N-glycoside bond, in which the N3-protonated zwitterion and imino tautomers are more favorable than carboxyl-protonated amino tautomers. This is consistent with the experimental data. Furthermore, our results reveal that the configuration rearrangement process of the protonated 5-caC would lower the stability of the N-glycoside bond and substantially reduce the barrier height for the subsequent C1'-N1 bond cleavage. This should be attributed to the smaller electrostatic repulsion between the leaving base and the negative phosphate group as a result of the structural rearrangement.

Theoretical Insights into the Photophysical Properties of 4CzIPN Doped in Different Hosts: A Multiscale Study.

As a typical thermally activated delayed fluorescence (TADF) emitter with green emission, 4CzIPN has attracted much attention recently. Most studies indicated that 4CzIPN doped in different hosts presented different performances; thus, the hosts should have an obvious influence on its photophysical properties. Herein, the influence of four kinds of hosts, including m-CzPym, m-CzTrz, p-CzPym, and p-CzTrz, on the photophysical properties of 4CzIPN is investigated. Molecular dynamics simulations were performed to simulate the host-guest conformations, and the photophysical properties were studied using the combined quantum mechanics/molecular mechanics method coupled with the thermal-vibration correlation function method. It is found that 4CzIPN in doped films has larger transition dipole moments and spin-orbital coupling constants compared to that in nondoped films. Faster radiative decay, intersystem crossing rates, and higher fluorescence efficiency could be obtained in doped films. Our work helps to better understand the photophysical properties of 4CzIPN in doped films and may favor the design of new hosts.