Exploring the formation pathway and antioxidant properties of the sugar-smoking pigment 5-GGMF

The present study aimed to elucidate the pathway of pigment formation and identify the source of antioxidant activity during sugar smoking. Building upon previous research, this investigation replicated the sucrose cleavage process involved in sugar-smoking through model reactions to obtain distinct model reaction products. The products were analyzed using various techniques such as ultraviolet-visible spectrometry, Fourier-transform infrared spectroscopy, high-performance liquid chromatography, and high-performance liquid chromatography-tandem mass spectrometry. The findings revealed that the pyrolysis of sucrose at 330 °C yielded glucose and fructose, with fructose pyrolysis producing significantly more 5-HMF than glucose. Moreover, the antioxidant capacity of 5-HMF was found to make a substantial contribution. The primary source of 5-HMF was identified as fructose resulting from the cleavage of sucrose at 330 °C, while the primary pathway for the formation of the sugar-smoking pigment 5-GGMF was attributed to the intermolecular dehydration of 5-HMF and glucose at 150 °C.

Understanding cellulose pyrolysis via ab initio deep learning potential field

Comprehensive and dynamic studies of cellulose pyrolysis reaction mechanisms are crucial in designing experiments and processes with enhanced safety, efficiency, and sustainability. The details of the pyrolysis mechanism are not readily available from experiments but can be better described via molecular dynamics (MD) simulations. However, the large size of cellulose molecules challenges accurate ab initio MD simulations, while existing reactive force field parameters lack precision. In this work, precise ab initio deep learning potentials field (DPLF) are developed and applied in MD simulations to facilitate the study of cellulose pyrolysis mechanisms. The formation mechanism and production rate of both valuable and greenhouse products from cellulose at temperatures larger than 1073 K are comprehensively described. This study underscores the critical role of advanced simulation techniques, particularly DLPF, in achieving efficient and accurate understanding of cellulose pyrolysis mechanisms, thus promoting wider industrial applications.

Chemical kinetic behaviors at the chain initiation stage of CH4/H2/air mixture

To intensively investigate chemical kinetic behaviors at the initial stage of CH₄/H₂/air mixture thoroughly, the density functional theory (CAMB3LYP/6-31 G) and a detailed mechanism (GRI-Mech3.0) were used to obtain kinetic and thermodynamic parameters. The reaction paths during the explosion process were analyzed, and the reaction rates of elementary reactions were compared with different ratios of CH₄/H₂/air mixture. The key reactions at the initiation stage of CH₄/H₂/air mixture explosion were determined, and their configurations were optimized. The reaction mechanism, reaction channel and configuration parameters of key reactions were obtained, which was verified by the intrinsic reaction coordinate (IRC) theory. Results show that H₂ addition increases the laminar burning velocity, while it shortens the ignition delay time of H₂/CH₄/air mixture. The addition of hydrogen greatly accelerated the explosion reaction from sample 1 to sample 4. Moreover, CH₄ still plays a key role at the chain initiation stage in H₂/CH₄/air mixture system; the addition of H₂ would not compete with CH₄ for triggering the explosion reaction, nor will it suppress the explosion of CH₄. H₂ could not replace or take precedence over the chain branching reactant (CH₂O) of CH₄ explosion to react with O₂. Besides, H₂ takes precedence over CH₄ in the process of chain transfer after the chain reaction beginning, although CH₄ has a distinct advantage in the chain initiation stage. The present results can provide theoretical guidance for the prevention and control of gas explosion, which may effectively reduce the explosion hazards.

The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning

BACKGROUND

Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest.

METHOD

A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition.

RESULTS

The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space.

CONCLUSIONS

The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements.

GENERAL SIGNIFICANCE

The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.

Theoretical rate constant of methane oxidation from the conventional transition-state theory

The potential energy surface for the first step of the methane oxidation CH₄ + O₂➔CH₃ + HO₂ was studied using the London-Eyring-Polanyi-Sato equation (LEPS) and the conventional transition-state theory (CTST). The calculated activation energy and rate constant values were in good agreement with the experimental and theoretical values reported in the literature using the shock tube technique and coupled cluster method respectively. The rate equation from CTST, although simple, provides good results to study the H-shift between methane and the oxygen molecules.

A minimization principle for transition paths of maximum flux for collective variables

Considered is the construction of transition paths of conformational changes for proteins and other macromolecules, using methods that do not require the generation of dynamics trajectories. Special attention is given to the use of a reduced set of collective variables for describing such paths. A favored way to define transition paths is to seek channels through the transition state having cross sections with a high reactive flux (density of last hitting points of reactive trajectories). Given here is a formula for reactive flux that is independent of the parameterization of "collective variable space." This formula is needed for the principal curve of the reactive flux (as in the revised finite temperature string method) and for the maximum flux transition (MaxFlux) path. Additionally, a resistance functional is derived for narrow tubes, which when minimized yields a MaxFlux path. A strategy for minimization is outlined in the spirit of the string method. Finally, alternative approaches based on determining trajectories of high probability are considered, and it is observed that they yield paths that depend on the parameterization of collective variable space, except in the case of zero temperature, where such a path coincides with a MaxFlux path.

Computational Biochemistry-Enzyme Mechanisms Explored

Understanding enzyme mechanisms is a major task to achieve in order to comprehend how living cells work. Recent advances in biomolecular research provide huge amount of data on enzyme kinetics and structure. The analysis of diverse experimental results and their combination into an overall picture is, however, often challenging. Microscopic details of the enzymatic processes are often anticipated based on several hints from macroscopic experimental data. Computational biochemistry aims at creation of a computational model of an enzyme in order to explain microscopic details of the catalytic process and reproduce or predict macroscopic experimental findings. Results of such computations are in part complementary to experimental data and provide an explanation of a biochemical process at the microscopic level. In order to evaluate the mechanism of an enzyme, a structural model is constructed which can be analyzed by several theoretical approaches. Several simulation methods can and should be combined to get a reliable picture of the process of interest. Furthermore, abstract models of biological systems can be constructed combining computational and experimental data. In this review, we discuss structural computational models of enzymatic systems. We first discuss various models to simulate enzyme catalysis. Furthermore, we review various approaches how to characterize the enzyme mechanism both qualitatively and quantitatively using different modeling approaches.

Toward Determining ATPase Mechanism in ABC Transporters: Development of the Reaction Path-Force Matching QM/MM Method

Adenosine triphosphate (ATP)-binding cassette (ABC) transporters are ubiquitous ATP-dependent membrane proteins involved in translocations of a wide variety of substrates across cellular membranes. To understand the chemomechanical coupling mechanism as well as functional asymmetry in these systems, a quantitative description of how ABC transporters hydrolyze ATP is needed. Complementary to experimental approaches, computer simulations based on combined quantum mechanical and molecular mechanical (QM/MM) potentials have provided new insights into the catalytic mechanism in ABC transporters. Quantitatively reliable determination of the free energy requirement for enzymatic ATP hydrolysis, however, requires substantial statistical sampling on QM/MM potential. A case study shows that brute force sampling of ab initio QM/MM (AI/MM) potential energy surfaces is computationally impractical for enzyme simulations of ABC transporters. On the other hand, existing semiempirical QM/MM (SE/MM) methods, although affordable for free energy sampling, are unreliable for studying ATP hydrolysis. To close this gap, a multiscale QM/MM approach named reaction path-force matching (RP-FM) has been developed. In RP-FM, specific reaction parameters for a selected SE method are optimized against AI reference data along reaction paths by employing the force matching technique. The feasibility of the method is demonstrated for a proton transfer reaction in the gas phase and in solution. The RP-FM method may offer a general tool for simulating complex enzyme systems such as ABC transporters.

Methods for Efficiently and Accurately Computing Quantum Mechanical Free Energies for Enzyme Catalysis

Enzyme activity is inherently linked to free energies of transition states, ligand binding, protonation/deprotonation, etc.; these free energies, and thus enzyme function, can be affected by residue mutations, allosterically induced conformational changes, and much more. Therefore, being able to predict free energies associated with enzymatic processes is critical to understanding and predicting their function. Free energy simulation (FES) has historically been a computational challenge as it requires both the accurate description of inter- and intramolecular interactions and adequate sampling of all relevant conformational degrees of freedom. The hybrid quantum mechanical molecular mechanical (QM/MM) framework is the current tool of choice when accurate computations of macromolecular systems are essential. Unfortunately, robust and efficient approaches that employ the high levels of computational theory needed to accurately describe many reactive processes (ie, ab initio, DFT), while also including explicit solvation effects and accounting for extensive conformational sampling are essentially nonexistent. In this chapter, we will give a brief overview of two recently developed methods that mitigate several major challenges associated with QM/MM FES: the QM non-Boltzmann Bennett's acceptance ratio method and the QM nonequilibrium work method. We will also describe usage of these methods to calculate free energies associated with (1) relative properties and (2) along reaction paths, using simple test cases with relevance to enzymes examples.

Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes

BACKGROUND

Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting.

METHODS

Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes.

RESULTS

Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation.

CONCLUSIONS

The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures.

GENERAL SIGNIFICANCE

Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.