Combined Quantum Mechanics/Molecular Mechanics Study on the Reversible Isomerization of Glucose and Fructose Catalyzed by Pyrococcus furiosus Phosphoglucose Isomerase

Short-Term Exposure to 2-Amino-1-methyl-6-phenylimidazo[4,5-b]pyridine Induces Colonic Energy Metabolism Disorders in Rats

2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) is one of the most abundant heterocyclic aromatic amines generated in thermally processed meat products, and the toxicities of its short-term exposure in the intestines remain unclear. This study aimed to elucidate the short-term PhIP toxicity in colons through administering PhIP orally to rats for 4 weeks. The results indicated that short-term PhIP exposure induced colonic oxidative stress, a significant decrease of serum triglyceride, and a disrupted colonic gene expression pattern associated with mitochondrial electron transport chain and energy metabolism. Thirteen energy metabolites, including lactate and d-erythrose-4-phosphate, showed significant changes under short-term PhIP effects. Energy metabolism pathway analysis revealed that PhIP-induced colonic energy metabolism disorders are characterized by inhibited glycolysis and enhanced tricarboxylic acid cycle. Further investigation found that PhIP altered the energy metabolic phenotype of colon epithelial cells to increase aerobic respiration. In summary, our study provides new insights into the colon toxicity induced by a short-term PhIP exposure.

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer

Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.

Isomerase and epimerase: overview and practical application in production of functional sugars

The biosynthesis of functional sugars has gained significant attention due to their potential health benefits and increasing demand in the food industry. Enzymatic synthesis has emerged as a promising approach, offering high catalytic efficiency, chemoselectivity, and stereoselectivity. However, challenges such as poor thermostability, low catalytic efficiency, and food safety concerns have limited the commercial production of functional sugars. Protein engineering, including directed evolution and rational design, has shown promise in overcoming these barriers and improving biocatalysts for large-scale production. Furthermore, enzyme immobilization has proven effective in reducing costs and facilitating the production of functional sugars. To ensure food safety, the use of food-grade expression systems has been explored. However, downstream technologies, including separation, purification, and crystallization, still pose challenges in terms of efficiency and cost-effectiveness. Addressing these challenges is crucial to optimize the overall production process. Despite the obstacles, the future outlook for functional sugars is promising, driven by increasing awareness of their health benefits and continuous technological advancements. With further research and technological breakthroughs, industrial-scale production of functional sugars through biosynthesis will become a reality, leading to their widespread incorporation in various industries and products.

Biochemical and Structural Characterisation of a Novel D-Lyxose Isomerase From the Hyperthermophilic Archaeon Thermofilum sp

A novel D-lyxose isomerase has been identified within the genome of a hyperthermophilic archaeon belonging to the Thermofilum species. The enzyme has been cloned and over-expressed in Escherichia coli and biochemically characterised. This enzyme differs from other enzymes of this class in that it is highly specific for the substrate D-lyxose, showing less than 2% activity towards mannose and other substrates reported for lyxose isomerases. This is the most thermoactive and thermostable lyxose isomerase reported to date, showing activity above 95°C and retaining 60% of its activity after 60 min incubation at 80°C. This lyxose isomerase is stable in the presence of 50% (v/v) of solvents ethanol, methanol, acetonitrile and DMSO. The crystal structure of the enzyme has been resolved to 1.4–1.7 A. resolution in the ligand-free form and in complexes with both of the slowly reacting sugar substrates mannose and fructose. This thermophilic lyxose isomerase is stabilised by a disulfide bond between the two monomers of the dimeric enzyme and increased hydrophobicity at the dimer interface. These overall properties of high substrate specificity, thermostability and solvent tolerance make this lyxose isomerase enzyme a good candidate for potential industrial applications.

Distant Non-Obvious Mutations Influence the Activity of a Hyperthermophilic Pyrococcus furiosus Phosphoglucose Isomerase

The cupin-type phosphoglucose isomerase (PfPGI) from the hyperthermophilic archaeon Pyrococcus furiosus catalyzes the reversible isomerization of glucose-6-phosphate to fructose-6-phosphate. We investigated PfPGI using protein-engineering bioinformatics tools to select functionally-important residues based on correlated mutation analyses. A pair of amino acids in the periphery of PfPGI was found to be the dominant co-evolving mutation. The position of these selected residues was found to be non-obvious to conventional protein engineering methods. We designed a small smart library of variants by substituting the co-evolved pair and screened their biochemical activity, which revealed their functional relevance. Four mutants were further selected from the library for purification, measurement of their specific activity, crystal structure determination, and metal cofactor coordination analysis. Though the mutant structures and metal cofactor coordination were strikingly similar, variations in their activity correlated with their fine-tuned dynamics and solvent access regulation. Alternative, small smart libraries for enzyme optimization are suggested by our approach, which is able to identify non-obvious yet beneficial mutations.

Computational Analysis of Host-Pathogen Protein Interactions between Humans and Different Strains of Enterohemorrhagic Escherichia coli

Serotype O157:H7, an enterohemorrhagic Escherichia coli (EHEC), is known to cause gastrointestinal and systemic illnesses ranging from diarrhea and hemorrhagic colitis to potentially fatal hemolytic uremic syndrome. Specific genetic factors like ompA, nsrR, and LEE genes are known to play roles in EHEC pathogenesis. However, these factors are not specific to EHEC and their presence in several non-pathogenic strains indicates that additional factors are involved in pathogenicity. We propose a comprehensive effort to screen for such potential genetic elements, through investigation of biomolecular interactions between E. coli and their host. In this work, an in silico investigation of the protein–protein interactions (PPIs) between human cells and four EHEC strains (viz., EDL933, Sakai, EC4115, and TW14359) was performed in order to understand the virulence and host-colonization strategies of these strains. Potential host–pathogen interactions (HPIs) between human cells and the “non-pathogenic” E. coli strain MG1655 were also probed to evaluate whether and how the variations in the genomes could translate into altered virulence and host-colonization capabilities of the studied bacterial strains. Results indicate that a small subset of HPIs are unique to the studied pathogens and can be implicated in virulence. This subset of interactions involved E. coli proteins like YhdW, ChuT, EivG, and HlyA. These proteins have previously been reported to be involved in bacterial virulence. In addition, clear differences in lineage and clade-specific HPI profiles could be identified. Furthermore, available gene expression profiles of the HPI-proteins were utilized to estimate the proportion of proteins which may be involved in interactions. We hypothesized that a cumulative score of the ratios of bound:unbound proteins (involved in HPIs) would indicate the extent of colonization. Thus, we designed the Host Colonization Index (HCI) measure to determine the host colonization potential of the E. coli strains. Pathogenic strains of E. coli were observed to have higher HCIs as compared to a non-pathogenic laboratory strain. However, no significant differences among the HCIs of the two pathogenic groups were observed. Overall, our findings are expected to provide additional insights into EHEC pathogenesis and are likely to aid in designing alternate preventive and therapeutic strategies.

Understanding the conformation transition in the activation pathway of β2 adrenergic receptor via a targeted molecular dynamics simulation

The conformation transition in the activation pathway of β2 adrenergic receptor was explored mainly using a target molecular dynamics simulation.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

A QM/MM MD study of the pH-dependent ring-opening catalysis and lid motif flexibility in glucosamine 6-phosphate deaminase

QM/MM MD and MM MD simulations reveal pH-dependent proton-shuttle ring-opening mechanisms of GlcN6P and dynamical behavior of the lid motif inSmuNagB.

A full picture of enzymatic catalysis by hydroxynitrile lyases from Hevea brasiliensis: protonation dependent reaction steps and residue-gated movement of the substrate and the product

Hydroxynitrile lyases (HNLs) defend plants from herbivores and microbial attack by releasing cyanide from hydroxynitriles.

Effects of water content on the tetrahedral intermediate of chymotrypsin - trifluoromethylketone in polar and non-polar media: observations from molecular dynamics simulation

The work uses MD simulation to study effects of five water contents (3 %, 10 %, 20 %, 50 %, 100 % v/v) on the tetrahedral intermediate of chymotrypsin - trifluoromethyl ketone in polar acetonitrile and non-polar hexane media. The water content induced changes in the structure of the intermediate, solvent distribution and H-bonding are analyzed in the two organic media. Our results show that the changes in overall structure of the protein almost display a clear correlation with the water content in hexane media while to some extent U-shaped/bell-shaped dependence on the water content is observed in acetonitrile media with a minimum/maximum at 10-20 % water content. In contrast, the water content change in the two organic solvents does not play an observable role in the stability of catalytic hydrogen-bond network, which still exhibits high stability in all hydration levels, different from observations on the free enzyme system [Zhu L, Yang W, Meng YY, Xiao X, Guo Y, Pu X, Li M (2012) J Phys Chem B 116(10):3292-3304]. In low hydration levels, most water molecules mainly distribute near the protein surface and an increase in the water content could not fully exclude the organic solvent from the protein surface. However, the acetonitrile solvent displays a stronger ability to strip off water molecules from the protein than the hexane. In a summary, the difference in the calculated properties between the two organic solvents is almost significant in low water content (<10 %) and become to be small with increasing water content. In addition, some structural properties at 10 ~ 20 % v/v hydration zone, to large extent, approach to those in aqueous solution.

Can hydridic-to-protonic hydrogen bonds catalyze hydride transfers in biological systems?

Catalysis of hydride transfer by hydridic-to-protonic hydrogen (HHH) bonding in α-hydroxy carbonyl isomerization reactions was examined computationally in the lithium salts of 7-substituted endo-3-hydroxybicyclo[2.2.1]hept-5-en-2-ones. The barrier for intramolecular hydride transfer in the parent system was calculated to be 17.2 kcal/mol. Traditional proton donors, such as OH and NH(3)(+), stabilized the metal cation-bridged transition state by 1.4 and 3.3 kcal/mol, respectively. Moreover, among the conformers of the OH systems, the one in which the proton donor is able to interact with the migrating hydride (H(m)) has an activation barrier lower by 1.3 and 1.7 kcal/mol than the other possible OH conformers. By contrast, the presence of an electronegative group such as F, which disfavors the migration electronically by opposing development of hydridic charge, destabilizes the hydride migration by 1.5 kcal/mol relative to the epimeric exo system. In both ground and transition states the H(m)···H distance decreased with increasing acidity of the proton donor, reaching a minimum of 1.58 Å at the transition state for NH(3)(+). Both Mulliken and NPA charges show enhancement of negative character of the migrating hydride in the cases in which HHH bonding is possible.

Structure-based annotation of a novel sugar isomerase from the pathogenic E. coli O157:H7

Prokaryotes can use a variety of sugars as carbon sources in order to provide a selective survival advantage. The gene z5688 found in the pathogenic Escherichia coli O157:H7 encodes a "hypothetical" protein of unknown function. Sequence analysis identified the gene product as a putative member of the cupin superfamily of proteins, but no other functional information was known. We have determined the crystal structure of the Z5688 protein at 1.6 A resolution and identified the protein as a novel E. coli sugar isomerase (EcSI) through overall fold analysis and secondary-structure matching. Extensive substrate screening revealed that EcSI is capable of acting on d-lyxose and d-mannose. The complex structure of EcSI with fructose allowed the identification of key active-site residues, and mutagenesis confirmed their importance. The structure of EcSI also suggested a novel mechanism for substrate binding and product release in a cupin sugar isomerase. Supplementation of a nonpathogenic E. coli strain with EcSI enabled cell growth on the rare pentose d-lyxose.

DL-FIND: an open-source geometry optimizer for atomistic simulations

Geometry optimization, including searching for transition states, accounts for most of the CPU time spent in quantum chemistry, computational surface science, and solid-state physics, and also plays an important role in simulations employing classical force fields. We have implemented a geometry optimizer, called DL-FIND, to be included in atomistic simulation codes. It can optimize structures in Cartesian coordinates, redundant internal coordinates, hybrid-delocalized internal coordinates, and also functions of more variables independent of atomic structures. The implementation of the optimization algorithms is independent of the coordinate transformation used. Steepest descent, conjugate gradient, quasi-Newton, and L-BFGS algorithms as well as damped molecular dynamics are available as minimization methods. The partitioned rational function optimization algorithm, a modified version of the dimer method and the nudged elastic band approach provide capabilities for transition-state search. Penalty function, gradient projection, and Lagrange-Newton methods are implemented for conical intersection optimizations. Various stochastic search methods, including a genetic algorithm, are available for global or local minimization and can be run as parallel algorithms. The code is released under the open-source GNU LGPL license. Some selected applications of DL-FIND are surveyed.

Determination of the structure form of the fourth ligand of zinc in Acutolysin A using combined quantum mechanical and molecular mechanical simulation

Acutolysin A, which is isolated from the snake venom of Agkistrodon acutus, is a member of the SVMPs subfamily of the metzincin family, and it is a snake venom zinc metalloproteinase possessing only one catalytic domain. The catalytic zinc ion, in the active site, is coordinated in a tetrahedral manner with three imidazole nitrogen atoms of histidine and one oxygen atom. It is uncertain whether this oxygen atom is a water molecule or a hydroxide ion just from the three-dimensional X-ray crystal structure. The identity of the fourth ligand of zinc is theoretically determined for the first time by performing both combined quantum mechanical and molecular mechanical (QM/MM) simulation and high-level quantum mechanical calculations. All of the results obtained indicate that the fourth ligand in the active site of the reported X-ray crystal structure is a water molecule rather than a hydroxide anion. On the basis of these theoretical results, we note that the experimental observed pH dependence of the proteolytic and hemorrhagic activity of Acutolysin A can be attributed to the deprotonation of the zinc-bound water to yield a better nucleophile, the hydroxide ion. Structural analyses revealed structural details useful for the understanding of acutolysin catalytic mechanism.