Unraveling the atomic mechanisms underlying glyphosate insensitivity in EPSPS: implications of distal mutations.

5-enolpyruvyl shikimate-3-phosphate synthase (EPSPS), as an indispensable enzyme in the shikimate pathway, is the specific target of grasser killer glyphosate (GPJ). GPJ is a competitive inhibitor of phosphoenolpyruvate (PEP), which is the natural substrate of EPSPS. A novel Ls-EPSPS gene variant discovered from Liliaceae, named ELs-EPSPS, includes five distal mutations, E112V, D142N, T351S, D425G, and R496G, endowing high GPJ insensitivity. However, the implicit molecular mechanism of the enhanced tolerance/insensitivity of GPJ in ELs-EPSPS is not fully understood. Herein, we try to interpret the hidden molecular mechanism using computational methods. Computational results reveal the enhanced flexibility of apo EPSPS upon mutations. The enhanced affinity of the initial binding substrate shikimate-3-phosphate (S3P), and the higher probability of second ligands PEP/GPJ entering the pocket are observed in the ELs-EPSPS-S3P system. Docking and MD results further confirmed the decreased GPJ-induced EPSPS inhibition upon mutations. And, the alterations of K98 and R179 side-chain orientations upon mutations are detrimental to GPJ binding at the active site. Additionally, the oscillation of side chain K98, in charge of PEP location, improves the proximity effect for substrates in the dual-substrate systems upon mutations. Our results clarify that the enhanced GPJ tolerance of EPSPS is achieved from decreased competitive inhibition of GPJ at the atomic perspective, and this finding further contributes to the cultivation of EPSPS genes with higher GPJ tolerance/insensitivity and a mighty renovation for developing glyphosate-resistant crops.Communicated by Ramaswamy H. Sarma.

Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models.

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 µs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

Characterizing the Molecular Mechanism of the Lethal C423D Mutation in FgMyoI: A Molecular Perspective.

The lethal mutation C423D in Fusarium graminearum myosin I (FgMyoI) occurs close to the binding pocket of the allosteric inhibitor phenamacril and causes severe inhibition on mycelial growth of F. graminearum strain PH-1. Here, based on extensive Gaussian accelerated molecular dynamics simulations and wet experiments, we elucidate the underlying molecular mechanism of the abnormal functioning of the FgMyoIC423D mutant at the atomistic level. Our results suggest that the damaging mutation C423D exhibits a synergistic allosteric inhibition mechanism similar to but more robust than that of phenamacril, including effects on the active site and actin binding. Unlike phenamacril-induced closure of Switch2, the mutation results in unfolding of the N-terminal relay helix with a partially opened Switch2 and blocks the structural rearrangement of the relay/SH1 helices, impairing the proper initiation of the recovery stroke. Due to the significant influence of C423D mutation on the function of FgMyoI, designing covalent inhibitors targeting this site holds tremendous potential.

Unraveling the mechanism of ceftaroline-induced allosteric regulation in penicillin-binding protein 2a: insights for novel antibiotic development against methicillin-resistant Staphylococcus aureus.

Methicillin-resistant Staphylococcus aureus (MRSA) acquires high-level resistance against ß-lactam antibiotics by expressing penicillin-binding protein 2a (PBP2a). PBP2a is a cell wall-synthesizing protein whose closed active site exhibits a reduced binding affinity toward ß-lactam antibiotics. Ceftaroline (CFT), a fifth-generation cephalosporin, can effectively inhibit the PBP2a activity by binding to an allosteric site to trigger the active site opening, allowing a second CFT to access the active site. However, the essential mechanism behind the allosteric behavior of PBP2a remains unclear. Herein, computational simulations are employed to elucidate how CFT allosterically regulates the conformation and dynamics of the active site of PBP2a. While CFT stabilizes the allosteric domain surrounding it, it simultaneously enhances the dynamics of the catalytic domain. Specifically, the study successfully captured the opening process of the active pocket in the allosteric CFT-bound systems and discovered that CFT alters the potential signal-propagating pathways from the allosteric site to the active site. These findings reveal the implied mechanism of the CFT-mediated allostery in PBP2a and provide new insights into dual-site drug design or combination therapy against MRSA targeting PBP2a.

Surface-independent CO2 and CO reduction on two-dimensional kagome metal KV3Sb5.

Two-dimensional kagome metals possess rich band structure characteristics, including Dirac points, flat bands, and van Hove singularities, because of their special geometric structures. Furthermore, kagome metals AV3Sb5 (A = K, Rb, and Cs) have garnered significant attention due to their nontrivial topological electronic structures. In this study, we theoretically demonstrate that the KV3Sb5 (001) surface is conducive to CO2 and CO reduction. The thermodynamic stability and electrochemical states of various surface types are investigated. The reaction paths reveal that the product is identical on different surfaces, and the free energy profiles exhibit low onset potentials. This paper elucidates the effect of two-dimensional topological kagome metals on CO2 and CO reduction.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches.

BACKGROUND: Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear. RESULTS: This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin. CONCLUSION: Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

Computational Insights into the Allosteric Modulation of a Phthalate-Degrading Hydrolase by Distal Mutations.

Phthalate esters (PAEs) are a ubiquitous kind of environmental endocrine that disrupt chemicals, causing environmental and health issues. EstJ6 is an effective phthalate-degrading hydrolase, and its mutant with a combination of three non-conservative distal mutations has an improved activity against PAEs with unknown molecular mechanisms. Herein, we attempt to fill the significant gap between distal mutations and the activity of this enzyme using computational approaches. We found that mutations resulted in a redistribution of the enzyme's preexisting conformational states and dynamic changes of key functional regions, especially the lid over the active site. The outward motion of the lid upon the mutations made it easier for substrates or products to enter or exit. Additionally, a stronger substrate binding affinity and conformational rearrangements of catalytic reaction-associated residues in the mutant, accompanied by the strengthened communication within the protein, could synergistically contribute to the elevated catalytic efficiency. Finally, an attempt was made to improve the thermostability of EstJ6 upon introducing a distal disulfide bond between residues A23 and A29, and the simulation results were as expected. Together, our work explored the allosteric effects caused by distal mutations, which could provide insights into the rational design of esterases for industrial applications in the future.

GPCR Allostery: A View from Computational Biology.

G protein-coupled receptors (GPCRs) represent a large superfamily of cell-surface proteins that mediate cell signaling and regulate virtually various aspects of physiological and pathological processes, therefore serving as a rich source of drug targets. As intrinsically allosteric proteins, numerous functions of GPCRs are regulated via allostery, whereby allosteric modulators binding at a distal site regulate the function of the typical orthosteric site. However, only a few GPCR allosteric ligands have been presently approved as drugs due to the high dynamic structures of GPCRs. Fortunately, the rapid development of computational biology sheds light on understanding the mechanism of GPCR allosteric ligands, which is critical for the discovery of new therapeutic agents. Here, we present a comprehensive overview of the currently available resources and approaches in computational biology related to G protein-coupled receptor allostery and their conformational dynamics. In addition, current limitations and major challenges in the field are also discussed accordingly.

Mechanism of enhanced sensitivity of mutated ß-adrenergic-like octopamine receptor to amitraz in honeybee Apis mellifera: An insight from MD simulations.

BACKGROUND: Amitraz is one of the critical acaricides/insecticides for effective control of pest infestation of Varroa destructor mite, a devastating parasite of Apis mellifera, because of its low toxicity to honeybees. Previous assays verified that a typical G protein-coupled receptor, ß-adrenergic-like octopamine receptor (Octß2R), is the unique target of amitraz, but the honeybee Octß2R resists to amitraz. However, the underlying molecular mechanism of the enhanced sensitivity or toxicity of amitraz to mutated honeybee Octß2RE208V/I335T/I350V is not fully understood. Here, molecular dynamics simulations are employed to explore the implied mechanism of the enhanced sensitivity to amitraz in mutant honeybee Octß2R. RESULTS: We found that amitraz binding stabilized the structure of Octß2R, particularly the intracellular loop 3 associated with the Octß2R signaling. Then, it was further demonstrated that both mutations and ligand binding resulted in a more rigid and compact amitraz binding site, as well as the outward movement of the transmembrane helix 6, which was a prerequisite for G protein coupling and activation. Moreover, mutations were found to promote the binding between Octß2R and amitraz. Finally, community analysis illuminated that mutations and amitraz strengthened the residue-residue communication within the transmembrane domain, which might facilitate the allosteric signal propagation and activation of Octß2R. CONCLUSION: Our results unveiled structural determinants of improved sensitivity in the Octß2R-amitraz complex and may contribute to further structure-based drug design for safer and less toxic selective insecticides. © 2022 Society of Chemical Industry.

A G-quadruplex/hemin structure-undamaged method to inhibit peroxidase-mimic DNAzyme activity for biosensing development.

Damaging the structure of the G-quadruplex (G4) to prevent the formation of the G4/hemin complex is presently the only available method to inhibit the activity of the peroxidase-mimic DNAzyme. In this study, a unique intramolecular inhibitory effect of the adjacent base-pair (InE(N:N)), by installing a rationally adjacent base-pair of the G4 core sequence, is proposed for the inhibition of the DNAzyme activity, which eliminates the need to damage the entire G4 structure. Various base pairs show different abilities to inhibit DNAzyme activity. The adjacent adenine: thymine pair possesses the best inhibitory efficiency (17 times). Through detailed investigations of the InE(N:N), it was revealed that the adjacent adenine: thymine pair downregulated the formation of compound I in the catalytic process, thus inhibiting the G4 DNAzyme activity. The mechanism of inhibition indicated that the carbonyl group on the hexatomic ring of the complementary base played an important role. To further reflect the advantages of the proposed strategy, two InE(N:N)-based biosensors were developed for DNA analysis and Uracil-DNA glycosylase (UDG) detection. Compared with existing DNAzyme-based methods, the application of InE(N: N) facilitates the real-time assay and simplifies the design difficulty. Therefore, InE(N:N) provides new insights into the regulation of the DNAzyme activity and offers an efficient approach for the future application of DNAzyme.

Uncovering the Molecular Basis for the Better Gefitinib Sensitivity of EGFR with Complex Mutations over Single Rare Mutation: Insights from Molecular Simulations.

Epidermal growth factor receptor (EGFR) is an intensively focused target for anti-tumor compounds used in non-small cell lung cancer (NSCLC) therapy. Compared to the classical activating mutations, there are still many uncommon EGFR mutations associated with poorer responses to EGFR inhibitors. A detailed understanding of the molecular basis for multiple EGFR mutants exhibiting diverse responses to inhibitors is of critical importance for related drug development. Herein, we explored the molecular determinants contributing to the distinct responses of EGFR with a single rare mutation (G719S) or combined mutations (G719S/L858R and G719S/l861Q) to Gefitinib (IRE). Our results indicated that interactions, formed within the tetrad of residues S768 (in the αC-helix), D770 (in the αC-ß4 loop), Y827 (in the αE-helix), and R831 (in the catalytic loop), play an important role in the stability of αC-helix and the maintenance of K745-E762 salt bridge in the absence of IRE, which are weakened in the EGFRG719S system and enhanced in the EGFRG719S/L858R system upon IRE binding. Besides, the introduced hydrogen bonds by the co-occurring mutation partner also contribute to the stability of αC-helix. The work done for inhibitor dissociation suggests that IRE exhibits a stronger binding affinity to EGFRG719S/L858R mutant. Together, these findings provide a deeper understanding of minor mutations, which is essential for drug development targeting EGFR with less common mutations.

Investigating the Permeation Mechanism of Typical Phthalic Acid Esters (PAEs) and Membrane Response Using Molecular Dynamics Simulations.

Phthalic acid esters (PAEs) are typical environmental endocrine disrupters, interfering with the endocrine system of organisms at very low concentrations. The plasma membrane is the first barrier for organic pollutants to enter the organism, so membrane permeability is a key factor affecting their biological toxicity. In this study, based on computational approaches, we investigated the permeation and intramembrane aggregation of typical PAEs (dimethyl phthalate, DMP; dibutyl phthalate, DBP; di-2-ethyl hexyl phthalate, DEHP), as well as their effects on membrane properties, and related molecular mechanisms were uncovered. Our results suggested that PAEs could enter the membrane spontaneously, preferring the headgroup-acyl chain interface of the bilayer, and the longer the side chain (DEHP > DBP > DMP), the deeper the insertion. Compared with the shortest DMP, DEHP apparently increased membrane thickness, order, and rigidity, which might be due to its stronger hydrophobicity. Potential of means force (PMF) analysis revealed the presence of an energy barrier located at the water-membrane interface, with a maximum value of 2.14 kcal mol−1 obtained in the DEHP-system. Therefore, the difficulty of membrane insertion is also positively correlated with the side-chain length or hydrophobicity of PAE molecules. These findings will inspire our understanding of structure-activity relationship between PAEs and their effects on membrane properties, and provide a scientific basis for the formulation of environmental pollution standards and the prevention and control of small molecule pollutants.

Critical Extracellular Ca2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations.

The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca2+. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca2+ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca2+ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca2+ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca2+-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca2+ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.

Environmental and sensitization variations among asthma and/or rhinitis patients between 2008 and 2018 in China.

BACKGROUND: Little is known about the changes in allergen sensitization in China secondary to the environmental variations over the past decade. We aimed at investigating the variations in sensitization among asthma and/or rhinitis patients in China between 2008 and 2018. METHODS: This study analyzed cross-sectional data from national surveys conducted in China in 2008 and 2018. After finishing the questionnaire, participants underwent serum specific IgE measurements. A total of 2322 and 2798 patients were enrolled in 2008 and 2018, respectively. The significance of differences in sensitization rates among four regions of China were assessed. Correlation analysis was used to identify the associations of sensitization with climate change and planting of Artemisia desertorum between the two surveys. RESULTS: Compared with 2008, the general sensitization rate to mites significantly increased in 2018, which ranked highest among all tested allergens. Sensitization to pollens, especially Artemisia vulgaris, showed the greatest increase in the north. The annual mean temperature, rainfall and relative humidity in all four regions, and the Artemisia desertorum coverage in the northeastern area, increased significantly in 2018 as compared with 2008. From 2008 to 2018, an increase in Dermatophagoides pteronyssinus sensitization was significantly associated with an increase in relative humidity (r = 0.54, p = 0.037). The increase in A. vulgaris sensitization was significantly associated with the increase in the A. desertorum planting area (r = 0.67, p = 0.006) and with a decrease in rainfall (r = -0.59, p = 0.021). CONCLUSIONS: House dust mites remain the most important allergen in Chinese individuals with asthma and/or rhinitis. Pollen sensitization dramatically increased in northern China. Increases in sensitization to dust mites and Artemisia were related to the increases in humidity and planting area of A. desertorum.

Improved Thermoelectric Performance of P-type SnTe through Synergistic Engineering of Electronic and Phonon Transports.

SnTe has been regarded as a potential alternative to PbTe in thermoelectrics because of its environmentally friendly features. However, it is a challenge to optimize its thermoelectric (TE) performance as it has an inherent high hole concentration (nH∼2 × 1020 cm-3) and low mobility (µH∼18 cm2 V-1 s-1) at room temperature (RT), arising from a high intrinsic Sn vacancy concentration and large energy separation between its light and heavy valence bands. Therefore, its TE figure of merit is only 0.38 at ∼900 K. Herein, both the electronic and phonon transports of SnTe were engineered by alloying species Ag0.5Bi0.5Se and ZnO in succession, thus increasing the Seebeck coefficient and, at the same time, reducing the thermal conductivity. As a result, the TE performance improves significantly with the peak ZT value of ∼1.2 at ∼870 K for the sample (SnGe0.03Te)0.9(Ag0.5Bi0.5Se)0.1 + 1.0 wt % ZnO. This result proves that synergistic engineering of the electronic and phonon transports in SnTe is a good approach to improve its TE performance.

Molecular Mechanism of Ca2+ in the Allosteric Regulation of Human Parathyroid Hormone Receptor-1.

Parathyroid hormone (PTH) is an endogenous ligand that activates the PTH type 1 receptor (PTH1R) signaling. Ca2+, a common second messenger, acts as an allosteric regulator for prolonging the activation of PTH1R. However, a clear picture of the underlying allosteric mechanism is still missing. Herein, extensive molecular dynamics (MD) simulations are performed for PTH1R-PTH complexes with and without Ca2+ ions, allowing us to delineate the molecular details of calcium-induced allostery. Our results indicate that acidic residues in the extracellular loop 1 (ECL1) (D251, E252, E254, and E258-E260) and PTH (E19 and E22) serve as key determinants for local Ca2+-coupling structures and rigidity of ECL1. Moreover, the binding of Ca2+ induces conformational changes of transmembrane domain 6/7 (TM6/7) that are related to PTH1R activation and strengthens the residue-residue communication within PTH and TMD allosterically. Moreover, our results demonstrate that the presence of Ca2+ ions potentiates the interaction between PTH and PTH1R via steered molecular dynamics (SMD) simulations, while the point mutation in the PTH (PTHR25C) weakens the binding of PTH and PTH1R. These results support that Ca2+ ions might further prolong the residence time of PTH on PTH1R and facilitate the positive allostery of PTH1R. Together, the present work provides new insights into the allosteric regulation mechanism of GPCRs induced by ions and related drug design targeting the PTH1R allosteric pathway.

Understanding the Allosteric Modulation of PTH1R by a Negative Allosteric Modulator.

The parathyroid hormone type 1 receptor (PTH1R) acts as a canonical class B G protein-coupled receptor, regulating crucial functions including calcium homeostasis and bone formation. The identification and development of PTH1R non-peptide allosteric modulators have obtained widespread attention. It has been found that a negative allosteric modulator (NAM) could inhibit the activation of PTH1R, but the implied mechanism remains unclear. Herein, extensive molecular dynamics simulations together with multiple analytical approaches are utilized to unravel the mechanism of PTH1R allosteric inhibition. The results suggest that the binding of NAM destabilizes the structure of the PTH1R-PTH-spep/qpep (the C terminus of Gs/Gq proteins) complexes. Moreover, the presence of NAM weakens the binding of PTH/peps (spep and qpep) and PTH1R. The intra- and inter-molecular couplings are also weakened in PTH1R upon NAM binding. Interestingly, compared with our previous study of the positive allosteric effects induced by extracellular Ca2+, the enhanced correlation between the PTH and G-protein binding sites is significantly reduced by the replacement of this negative allosteric regulator. Our findings might contribute to the development of new therapeutic agents for diseases caused by the abnormal activation of PTH1R.

[A review of consensus statements, practice resources, standards and guidelines for clinical applications of next-generation sequencing technologies in the United States].

The use of whole exome sequencing (WES) for the detection of disease-causing variants of genetic diseases and for non-invasive prenatal screening (NIPS) of fetal aneuploidies are two major clinical applications of next generation sequencing (NGS). This article has summarized the official documents developed and updated by the American College of Medical Genetics and Genomics (ACMG) on governing WES and NIPS. These include the development of expert consensus policies and position statements on an ongoing basis to guide clinical application of NGS technology and variant analysis, establish evidence-based practical resources, as well as standards and guidelines to govern diagnosis and screening. These ACMG documents are valuable references to Chinese geneticists, but direct adoption of these standards and guidelines may not be practical due to the differences in disease-associated variant frequencies in Chinese population, socioeconomic status, and medical practice between the two countries. It is hoped that this review could facilitate the development of NGS and NIPS standards and guidelines that are consistent with international standards and concordant with medical genetics practice in China to provide high-quality, efficient and safe clinical services for patients and their families with genetic diseases.

Conformational dynamics is critical for the allosteric inhibition of cGAS upon acetyl-mimic mutations.

Detection of cytosolic dsDNA by cyclic GMP-AMP synthase (cGAS) is critical for the immune system to sense and fight against infection, but chronic activation of cGAS by self-DNA leads to autoimmune diseases without effective treatment yet. It was found that acetylation on either Lys384, Lys394, or Lys414 could inhibit the catalytic production of cGAMP by cGAS, and further suppressed self-DNA-induced autoimmunity. However, the implied mechanism remains unclear. Here, extensive molecular dynamics simulations combined with multiple analytical approaches were employed to uncover the allosteric inhibition mechanisms by using the K-to-Q mutations to mimic acetylation. Results suggested that the exterior loops contributed most to the conformational dynamics of cGAS, and two concerted intrinsic motions were observed: the inward/outward or twisting movement for the outer appendage of lobe 1 and the open/closed swing of the active-site loops. Mutations slightly affected the binding of dsDNA and cGAMP. The shift of the conformational sampling of the active-site loops or residues around cGAMP upon mutation might potentially explain the inhibition of cGAS activity. Moreover, the intra- and inter-molecular coupling was weakened upon mutations more or less but via distinct pathways. Hence, conformational dynamics play a vital role in the allosteric inhibition of cGAS upon the studied acetyl-mimic mutations. As the studied acetyl-mimic mutations are located at either the inter-lobe or inter-molecular interfaces, hence except for acetylation, our findings might help the development of new therapeutics against autoimmune diseases due to abnormal cGAS activation by designing inter-lobe or intermolecular allosteric inhibitors.

Integration of multicomponent characterization, untargeted metabolomics and mass spectrometry imaging to unveil the holistic chemical transformations and key markers associated with wine steaming of Ligustri Lucidi Fructus.

Processing of traditional Chinese medicine (TCM) can enhance the efficacy and/or reduce the toxicity. Currently available approaches regarding TCM processing generally focus on a few markers, rendering a one-sided strategy that fail to unveil the involved global chemical transformation. We herein present a strategy, by integrating enhanced multicomponent characterization, untargeted metabolomics, and mass spectrometry imaging (MSI), to visualize the chemical transformation and identify the markers associated with the wine steaming of Ligustri Lucidi Fructus (LLF), as a case. An ultra-high-performance liquid chromatography/quadrupole-Orbitrap mass spectrometry-based polarity-switching (between the negative and positive modes), precursor ions list-including data-dependent acquisition approach was developed, which enabled the simultaneous targeted/untargeted characterization of 158 components from LLF via one injection analysis. Holistic, continuous, and time-dependent chemical variation trajectory, among different processing time (0-12 h) for LLF, was depicted by principle component analysis. Pattern recognition chemometrics could unveil 20 markers, among which the peak area ratios of eight components to oleuropein aglycone, used as an internal standard, were diagnostic to identify the processed (both the commercial and in-house prepared) from the raw LLF. Four markers (10-hydroxyoleoside dimethylester, 8-demethyl-7-ketoliganin, elenolic acid, and salidroside) showed an increasing trend, while another four (neonuezhenide/isomer, verbascoside/isomer, luteoline, and nuzhenal A) decreased in LLF after processing. MSI visualized the spatial distribution in the fruit and indicated consistent variation trends for four major markers deduced by the untargeted metabolomics approach. This integral strategy, in contrast to the conventional approaches, gives more convincing data supporting the processing mechanism investigations of TCM from a macroscopic perspective.