Uncovering the pharmacological mechanisms of Patchouli essential oil for treating ulcerative colitis.

ETHNOPHARMACOLOGICAL RELEVANCE: Pogostemonis Herba has long been used in traditional Chinese medicine to treat inflammatory disorders. Patchouli essential oil (PEO) is the primary component of Pogostemonis Herba, and it has been suggested to offer curative potential when applied to treat ulcerative colitis (UC). However, the pharmacological mechanisms of PEO for treating UC remain to be clarified. AIM OF THE STUDY: To elucidate the pharmacological mechanisms of PEO for treating UC. METHODS AND RESULTS: In the present study, transcriptomic and network pharmacology approaches were combined to clarify the mechanisms of PEO for treating UC. Our results reveal that rectal PEO administration in UC model mice significantly alleviated symptoms of UC. In addition, PEO effectively suppressed colonic inflammation and oxidative stress. Mechanistically, PEO can ameliorate UC mice by modulating gut microbiota, inhibiting inflammatory targets (OPTC, PTN, IFIT3, EGFR, and TLR4), and inhibiting the PI3K-AKT pathway. Next, the 11 potential bioactive components that play a role in PEO's anti-UC mechanism were identified, and the therapeutic efficacy of the pogostone (a bioactive component) in UC mice was partially validated. CONCLUSION: This study highlights the mechanisms through which PEO can treat UC, providing a rigorous scientific foundation for future efforts to develop and apply PEO for treating UC.

Unveiling the Dynamic Mechanism of SARS-CoV-2 Entry Host Cells at the Single-Particle Level.

Understanding the dynamic features of severe acute respiratory coronavirus 2 (SARS-CoV-2) binding to the cell membrane and entry cells is crucial for comprehending viral pathogenesis and transmission and facilitating the development of effective drugs against COVID-19. Herein, we employed atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) to study the binding dynamics between the virus and cell membrane. Our findings revealed that the Omicron variant of SARS-CoV-2 virus-like particles (VLPs) exhibited a slightly stronger affinity for the angiotensin-converting enzyme-2 (ACE2) receptor compared with the Delta variant and was significantly higher than the wild-type (WT). Using a real-time force-tracing technique, we quantified the dynamic parameters for a single SARS-CoV-2 VLP entry into cells, showing that approximately 200 ms and 60 pN are required. The parameters aligned with the analysis obtained from coarse-grained molecular dynamics (CGMD) simulations. Additionally, the Omicron variant invades cells at a higher entry cell speed, smaller force, and higher probability. Furthermore, single-particle fluorescence tracking visually demonstrated clathrin-dependent endocytosis for SARS-CoV-2 entry into A549 cells. The dynamic features of endocytosis provide valuable insights into the SARS-CoV-2 entry mechanism and possible intervention strategies targeting the viral infection process.

Temperature-Dependent Oxidation Behavior of Silicon Carbide Surface: Reactive Molecular Dynamics Simulations.

The high-temperature oxidation mechanism of silicon carbide (SiC) is crucial for designing thermal protection systems in aircraft. This study explored the oxidative chemical reaction processes and mechanism on SiC surface and interfaces within the temperature range of 300-2300 K by using reactive molecular dynamics simulation. The oxygen impact on the silica (SiOx) growth of the SiC surface was analyzed, which shows a progressive increase of silica thickness. The simulation results indicated that the oxidation process of SiC was a typical passive oxidation mechanism. With the environmental temperature rising and oxygen impact, the increase of oxidation thickness on the SiC surface undergoes three oxidizing reaction processes: little chemical adsorption of oxygen molecules on the initial surface, rapid oxidation of silicon and carbon, and dramatic oxidation of the interface between the oxidation layer and SiC. Additionally, this work studied the mechanism of oxidation thickness growth and chemical diffusion of oxygen. The oxidation rate is weakened according to the oxygen atom diffusion barrier effect of silica repulsion. Moreover, the kinetic parameters were statistically calculated by fitting the growth of Si-O bonds and their reaction rate constants. Subsequently, the activation energy and pre-exponential factors were derived by using the Arrhenius equation to model the chemical reaction kinetics of the thermal oxidation process. The chemical reaction behaviors of the two stages could be concluded as follows: (i) in stage I, the initial oxidation is reaction rate limiting; (ii) in stage II, SiC oxidation is limited by both the oxidation reaction rate and the oxygen diffusion coefficient of the oxidation layer. The activation energy of stage II increased compared with stage I due to the oxygen atoms diffusion barrier between the oxidation layers. This study on the oxidation and ablation mechanism of the SiC surface at the atomic scale would provide insight into understanding thermal oxidation behavior and the design of ceramic materials.

Improving the activity of an inulosucrase by rational engineering for the efficient biosynthesis of low-molecular-weight inulin.

Inulin, a widely recognized prebiotic, has diverse applications across various industrial sectors. Although inulin is primarily produced through plant extraction, there is growing interest in enzymatic synthesis as an alternative. The enzymatic production of inulin from sucrose, which yields polymers with degrees of polymerization similar to those of plant-derived inulin, shows potential as a viable replacement for traditional extraction methods. In this study, an inulosucrase from Neobacillus bataviensis was identified, demonstrating a non-processive mechanism specifically tailored for synthesizing inulin with polymerization degrees ranging from 3 to approximately 40. The enzyme exhibited optimal activity at pH 6.5 and 55 °C, efficiently producing inulin with a yield of 50.6%. Ca2+ can improve the activity and thermostability of this enzyme. To enhance catalytic total activity, site-directed and truncated mutagenesis techniques were applied, resulting in the identification of a mutant, T149S, displaying a significant 57% increase in catalytic total activity. Molecular dynamics simulations unveiled that the heightened flexibility observed in three surface regions positively influenced enzymatic activity. This study not only contributes to the theoretical foundation for inulosucrase engineering but also presents a potential avenue for the production of inulin.

Interface Stability and Reaction Mechanisms of Li3YCl5Br with High-Voltage Cathodes and Li Metal Anode: Insights from Ab Initio Simulations.

Recent advancements in battery technology emphasize the critical role of solid electrolytes in enhancing the performance and safety of next-generation batteries. In this study, we investigate the interface stability and reaction mechanisms of Li3YCl5Br, a promising halide-based solid electrolyte, in contact with high-voltage Ni-Mn-Co (NMC) cathodes and a Li metal anode using ab initio molecular dynamics simulations. Our findings reveal that Li3YCl5Br reacts with charged NMC cathodes. This reaction involves changes in the oxidation states of Br- anions in Li3YCl5Br and d-element cations in NMC, as well as the diffusion of Li ions from the solid electrolyte to the cathode to maintain charge balance. The reaction is confined to the interface, suggesting bulk stability. Conversely, the Li/Li3YCl5Br interface exhibits significant instability, with a chemical reaction that results in substantial structural changes and the formation of LiCl and LiBr at the solid electrolyte surface and metallic Y at the Li anode surface. These insights provide valuable information for optimizing interfacial design, aiming at improving the performance and reliability of all-solid-state batteries using halide solid electrolytes.

Computational insights into the inhibitory mechanism of type 2 diabetes mellitus by bioactive components of Oryza sativa L. indica (black rice).

Background: Type 2 diabetes mellitus is a metabolic disease categorized by hyperglycemia, resistance to insulin, and ß-cell dysfunction. Around the globe, approximately 422 million people have diabetes, out of which 1.5 million die annually. In spite of innovative advancements in the treatment of diabetes, no biological drug has been known to successfully cure and avert its progression. Thereupon, natural drugs derived from plants are emerging as a novel therapeutic strategy to combat diseases like diabetes. Objective: The current study aims to investigate the antidiabetic potential of natural compounds of Oryza sativa L. indica (black rice) in disease treatment. Methods: Antioxidant activity and alpha amylase assays were performed to evaluate the therapeutic potential of the extract of Oryza sativa L. indica. Gas chromatography-mass spectrometry (GC-MS) was used for identification of constituents from the ethanol extract. ADMET profiling (absorption, distribution, metabolism, excretion, and toxicity), network pharmacology, and molecular dynamics simulation were employed in order to uncover the active ingredients and their therapeutic targets in O. sativa L. indica against type 2 diabetes mellitus. Results: GC-MS of the plant extract provided a list of 184 compounds. Lipinski filter and toxicity parameters screened out 18 compounds. The topological parameters of the protein-protein interaction (PPI) were used to shortlist the nine key proteins (STAT3, HSP90AA1, AKT1, SRC, ESR1, MAPK1, NFKB1, EP300, and CREBBP) in the type 2 diabetes mellitus pathways. Later, molecular docking analysis and simulations showed that C14 (1H-purine-8-propanoic acid, .alpha.-amino-2, 3, 6, 7-tetrahydro-1,3,7-trimethyl-2,6-dioxo-) and C18 (cyclohexane-carboxamide, N-furfuryl) bind with AKT1 and ESR1 with a binding energy of 8.1, 6.9, 7.3, and 7.2 kcal/mol, respectively. RMSD (root-mean-square deviation) and RMSF (root-mean-square fluctuation) values for AKT1 and ESR1 have shown very little fluctuation, indicating that proteins were stabilized after ligand docking. Conclusion: This study suggests therapeutic drug candidates against AKT1 and ESR1 to treat type 2 diabetes mellitus. However, further wet-lab analysis is required to discover the best remedy for type 2 diabetes mellitus.

Investigation of the Inhibitory Activity of ß-Arbutin and Its Analogues on Tyrosinase Based on Molecular Docking and Enzyme Inhibition Kinetics.

ß-Arbutin, a natural glucoside hydroquinone derivative known for its skin-whitening properties through tyrosinase inhibition in melanin synthesis, may pose potential risks of allergy and carcinogenicity due to the release of hydroquinone during use. This study explores the inhibitory effects of phenyl-ß-D-pyranoglucoside (compound 1), 4-methoxyphenyl-ß-D-pyranoglucoside (compound 2), 4-hydroxymethylphenyl-ß-D-pyranoglucoside (compound 3), and ß-arbutin (compound 4) on tyrosinase using enzyme kinetics, molecular docking, and molecular dynamics simulations. Results show compounds 1, 3, and 4 exhibit competitive inhibition, while compound 2 shows mixed inhibition. Docking analysis reveals phenyl rings of all compounds interact with the enzyme's active site, with compound 3 forming a metal bond with copper ions. MD simulations indicate high stability for compounds 2, 3, and 4, with compound 3 showing the lowest RMSD and compact Rg, suggesting stronger binding. Compound 1 is less stable and less inhibitory. These insights are valuable for designing effective tyrosinase inhibitors.

Stick, Slide, or Bounce: Charge Density Controls Nanoparticle Diffusion.

The diffusion and interaction dynamics of charged nanoparticles (NPs) within charged polymer networks are crucial for understanding various biological and biomedical applications. Using a combination of coarse-grained molecular dynamics simulations and experimental diffusion studies, we investigate the effects of the NP size, relative surface charge density (ζ), and concentration on the NP permeation length and time. We propose a scaling law for the relative diffusion of NPs with respect to concentration and ζ, highlighting how these factors influence the NP movement within the network. The analyses reveal that concentration and ζ significantly affect NP permeation length and time, with ζ being critical, as critical as concentration. This finding is corroborated by controlled release experiments. Further, we categorize NP dynamics into sticking, sliding, and bouncing regimes, demonstrating how variations in ζ, concentration, and NP size control these behaviors. Through normalized attachment time (NAT) analyses, we elucidate the roles of electrostatic interactions, steric hindrance, and hydrodynamic forces in governing NP dynamics. These insights provide guidance for optimizing NP design in targeted drug delivery and advanced material applications, enhancing our understanding of NP behavior in complex environments.

CDR L3 Loop Rearrangement Switches Multispecific SPE-7 IgE Antibody From Hapten to Protein Binding.

The monoclonal IgE antibody SPE-7 was originally raised against a 2,4-dinitrophenyl (DNP) target. Through its ability to adopt multiple conformations, the antibody is capable of binding to a diverse range of small haptens and large proteins. The present study examines a dataset of experimentally determined crystal structures of the SPE-7 antibody to gain insight into the mechanisms that contribute to its multispecificity. With the emergence of more and more therapeutic antibodies against a huge repertoire of different targets, our research could be of great interest for future drug development. We are able to discriminate between the different paratope-binding states in the conformational ensembles obtained by enhanced sampling molecular dynamics simulations, and to calculate their transition timescales and state probabilities. Furthermore, we describe the key residues responsible for discriminating between the different binding capacities and identify a tryptophan in a central position of the CDR L3 loop as the residue of greatest interest. The overall dynamics of the paratope appear to be mainly influenced by the CDR L3 and CDR L1 loops.

Computer-Aided Site-Specific PEGylation of PET Hydrolases for Enhanced PET Degradation.

The enormous accumulation of poly(ethylene terephthalate) (PET) waste has posed a serious threat to the environment and human health, and biodegradation with PET hydrolase (PETase) can be a possible solution. Herein, we propose site-specifically modifying PETase with amphiphilic polymers to improve the enzyme performance at ambient temperature. For this purpose, we devise a computer-aided strategy to prioritize the conjugation site, and polyethylene glycol (PEG) preparations of 0.55 to 10 kDa are site-specifically conjugated to PETase. The most active conjugate PETase-PEG 5k (PETase-5K) shows an increase of melting temperature (3.88 °C) and significantly improves PET degradation performance (3.5- and 3.1-fold increases at 30 and 40 °C, respectively). Experimental investigation and molecular dynamics simulations reveal that the site-specific PEGylation increases the hydrophobic solvent-accessible surface area and the binding capability to the PET surface, thickens the hydration layer, increases the intramolecular hydrogen bonding, reduces the interactions between water and the conjugated enzyme surface, and rigidifies the enzyme structure via hydrogen bonding and hydrophobic interactions between the polymer and the enzyme, thus leading to improved enzymatic performance of PETase-5K. We further validate the versatility of the site-specific PEGylation in one of the most evolved variants of PETase, FAST-PETase, by 1.8-fold improvement in PET degradation at 30 °C. The presented computer-aided site-specific conjugation strategy has opened a new avenue to enhancing PETase performance at ambient temperature, and the contribution of PEGylation to PETase unraveled in this work laid a foundation for the rational engineering of PET hydrolases.

Assessing AF2's ability to predict structural ensembles of proteins.

Recent breakthroughs in protein structure prediction have enhanced the precision and speed at which protein configurations can be determined. Additionally, molecular dynamics (MD) simulations serve as a crucial tool for capturing the conformational space of proteins, providing valuable insights into their structural fluctuations. However, the scope of MD simulations is often limited by the accessible timescales and the computational resources available, posing challenges to comprehensively exploring protein behaviors. Recently emerging approaches have focused on expanding the capability of AlphaFold2 (AF2) to predict conformational substates of protein. Here, we benchmark the performance of various workflows that have adapted AF2 for ensemble prediction and compare the obtained structures with ensembles obtained from MD simulations and NMR. We provide an overview of the levels of performance and accessible timescales that can currently be achieved with machine learning (ML) based ensemble generation. Significant minima of the free energy surfaces remain undetected.

On the solvation properties of menthol-thymol mixtures. A Molecular Dynamics Investigation.

Using classical molecular dynamics, we have investigated the solvation of catechol, resorcinol, hydroquinone and 1,4-benzoquinone at infinite dilution, in a series of menthol - thymol mixtures in which the molar fraction of thymol (xTHY) has been increased by steps of 0.1, from 0 (pure menthol) to 1 (pure thymol). The evolution of the solvation shell around the solutes reveals that when xTHY is increased, the average number of hydrogen bonds (HB) where the solute acts as HB acceptor (HBA) and the solvent as HB donor (HBD) increases, while the amount of HB, in which the solute acts as HBD and the solvent as HBA, decreases. Overall, the total number of HBs between the different benzenediols and the solvent decreases with an increase of xTHY, while for benzoquinone the total number of HB increases. This points to the fact that "acidic" or HBD molecules are better solvated in mixtures with high menthol proportion, while "basic" or HBA molecules, are better solvated in thymol rich mixtures. The results reported herein follow the same trends as experimentally reported Kamlet-Taft parameters and present insights on how the composition of these "deep eutectic" mixtures maybe tweaked in order to optimize their solvation properties.

An integrative characterisation of proline cis and trans conformers in a disordered peptide.

Intrinsically disordered proteins (IDPs) often contain proline residues, which undergo cis/trans isomerisation. While molecular dynamics (MD) simulations have the potential to fully characterise the proline cis and trans sub-ensembles, they are limited by the slow timescales of isomerisation and force field inaccuracies. Nuclear magnetic resonance (NMR) spectroscopy can report on ensemble-averaged observables for both the cis-proline and trans-proline states, but a full atomistic characterisation of these conformers is challenging. Given the importance of proline cis/trans isomerisation for influencing the conformational sampling of disordered proteins, we employed a combination of all-atom MD simulations with enhanced sampling (metadynamics), NMR, and small-angle X-ray scattering (SAXS) to characterise the two sub-ensembles of the ORF6 C-terminal region (ORF6CTR) from SARS-CoV-2 corresponding to the proline-57 (P57) cis and trans states. We performed MD simulations in three distinct force fields: AMBER03ws, AMBER99SB-disp, and CHARMM36m, which are all optimised for disordered proteins. Each simulation was run for an accumulated time of 180-220 µs until convergence was reached, as assessed by blocking analysis. A good agreement between the cis-P57 populations predicted from metadynamic simulations in AMBER03ws was observed with populations obtained from experimental NMR data. Moreover, we observed good agreement between the radius of gyration predicted from the metadynamic simulations in AMBER03ws and that measured using SAXS. Our findings suggest that both the cis-P57 and trans-P57 conformations of ORF6CTR are extremely dynamic and that interdisciplinary approaches combining both multi-scale computations and experiments offer avenues to explore highly dynamic states that cannot be reliably characterised by either approach in isolation.

Exploring the Origins of Association of Poly(acrylic acid) Polyelectrolyte with Lysozyme in Aqueous Environment through Molecular Simulations and Experiments.

This study provides a detailed picture of how a protein (lysozyme) complexes with a poly(acrylic acid) polyelectrolyte (PAA) in water at the atomic level using a combination of all-atom molecular dynamics simulations and experiments. The effect of PAA and temperature on the protein's structure is explored. The simulations reveal that a lysozyme's structure is relatively stable except from local conformational changes induced by the presence of PAA and temperature increase. The effect of a specific thermal treatment on the complexation process is investigated, revealing both structural and energetic changes. Certain types of secondary structures (i.e., α-helix) are found to undergo a partially irreversible shift upon thermal treatment, which aligns qualitatively with experimental observations. This uncovers the origins of thermally induced aggregation of lysozyme with PAA and points to new PAA/lysozyme bonds that are formed and potentially enhance the stability in the complexes. As the temperature changes, distinct amino acids are found to exhibit the closest proximity to PAA, resulting into different PAA/lysozyme interactions; consequently, a different complexation pathway is followed. Energy calculations reveal the dominant role of electrostatic interactions. This detailed information can be useful for designing new biopolymer/protein materials and understanding protein function under immobilization of polyelectrolytes and upon mild denaturation processes.

Revealing the pH-dependent conformational changes in sol g 2.1 protein and potential ligands binding.

Sol g 2, a major protein found in the venom of the tropical fire ant (Solenopsis geminata), is well-known for its ability to bind various hydrophobic molecules. In this study, we investigate the binding activity of recombinant Sol g 2.1 protein (rSol g 2.1) with potential molecules, including (E)-ß-Farnesene, α-Caryophyllene, and 1-Octen-3-ol at different pH levels (pH 7.4 and 5.5) using fluorescence competitive binding assays (FCBA). Our results revealed that Sol g 2.1 protein has higher affinity binding with these ligands at neutral pH. Relevance to molecular docking and molecular dynamics simulations were utilized to provide insights into the stability and conformational dynamics of Sol g 2.1 and its ligand complexes. After simulation, we found that Sol g 2.1 protein has higher affinity binding with these ligands as well as high structural stability at pH 7.4 than at an acidic pH level, indicating by RMSD, RMSF, Rg, SASA, and principal component analysis (PCA). Additionally, the Sol g 2.1 protein complexes at pH 7.4 showed significantly lower binding free energy (∆Gbind) and higher total residue contributions, particularly from key non-polar amino acids such as Trp36, Met40, Cys62, and Ile104, compared to the lower pH environment. These explain why they exhibited higher binding affinity than the lower pH. Therefore, we suggested that Sol g 2.1 protein is a pH-responsive carrier protein. These findings also expand our understanding of protein-ligand interactions and offer potential avenues for the development of innovative drug delivery strategies targeting Sol g 2.1 protein.

Stereo-selectivity of enantiomeric inhibitors to ubiquitin-specific protease 7 (USP7) dissected by molecular docking, molecular dynamics simulations, and binding free energy calculations.

The ubiquitin-specific protease 7 (USP7), as a member of deubiquitination enzymes, represents an attractive therapeutic target for various cancers, including prostate cancer and liver cancer. The change of the inhibitor stereocenter from the S to R stereochemistry (S-ALM → R-ALM34) markedly improved USP7 inhibitory activity. However, the molecular mechanism for the stereo-selectivity of enantiomeric inhibitors to USP7 is still unclear. In this work, molecular docking, molecular dynamics (MD) simulations, molecular mechanics/Generalized-Born surface area (MM/GBSA) calculations, and free energy landscapes were performed to address this mystery. MD simulations revealed that S-ALM34 showed a high degree of conformational flexibility compared to the R-ALM34 counterpart, and S-ALM34 binding led to the enhanced intradomain motions of USP7, especially the BL1 and BL2 loops and the two helices α4 and α5. MM/GBSA calculations showed that the binding strength of R-ALM34 to USP7 was stronger than that of S-ALM34 by - 4.99 kcal/mol, a similar trend observed by experimental data. MM/GBSA free energy decomposition was further performed to differentiate the ligand-residue spectrum. These analyses not only identified the hotspot residues interacting with R-ALM34, but also revealed that the hydrophobic interactions from F409, K420, H456, and Y514 play the major determinants in the binding of R-ALM34 to USP7. This result is anticipated to shed light on energetic basis and conformational dynamics information to aid in the design of more potent and selective inhibitors targeting USP7.

Pathogenic G6PD variants: Different clinical pictures arise from different missense mutations in the same codon.

G6PD deficiency results from mutations in the X-linked G6PD gene. More than 200 variants are associated with enzyme deficiency: each one of them may either cause predisposition to haemolytic anaemia triggered by exogenous agents (class B variants), or may cause a chronic haemolytic disorder (class A variants). Genotype-phenotype correlations are subtle. We report a rare G6PD variant, discovered in a baby presenting with severe jaundice and haemolytic anaemia since birth: the mutation of this class A variant was found to be p.(Arg454Pro). Two variants affecting the same codon were already known: G6PD Union, p.(Arg454Cys), and G6PD Andalus, p.(Arg454His). Both these class B variants and our class A variant exhibit severe G6PD deficiency. By molecular dynamics simulations, we performed a comparative analysis of the three mutants and of the wild-type G6PD. We found that the tetrameric structure of the enzyme is not perturbed in any of the variants; instead, loss of the positively charged Arg residue causes marked variant-specific rearrangement of hydrogen bonds, and it influences interactions with the substrates G6P and NADP. These findings explain severe deficiency of enzyme activity and may account for p.(Arg454Pro) expressing a more severe clinical phenotype.

Cheminformatics-enhanced discovery of therapeutic agents targeting isocitrate lyase in Mycobacterium tuberculosis infections.

Tuberculosis (TB) is a global health challenge; therefore, there is an urgent requirement to develop a novel and more effective anti-TB therapeutic. This study targeted the isocitrate lyase (ICL) protein due to its pivotal role in the pathogenicity of Mycobacterium tuberculosis (Mtb). Virtual screening of 8752 bioactive compounds used an ML-based QSAR model and molecular docking. ADMET testing was performed on the top three hits to identify the compound most closely mimicking a drug molecule. The top hits, 648 and 2785758, showed high binding affinity towards ICL with -7.3 and -7 kcal/mol, comparable to the control. These molecules also showed strong binding with the residue Asp108, which plays a vital role in ICL activity. Molecular dynamics simulations showed stability for 648 and 2785758, comparable to the control compound used in this study. It was found that 648 bound to the protein maintained the RMSD constant and consistent at 0.3 nm for a complete 100 ns simulation. 2785758 showed a comparable RMSD trend to the control. Both 648 and 2785758 showed high RMSF for critical residue Asp108. Further, PCA and FEL confirmed the formation of a stable complex. MM/GBSA estimations of binding free energy indicated that compounds 648 had an elevated level of stability (ΔGTOTAL = -28.11 kcal/mol) and 2785758 (ΔGTOTAL = -21.05 kcal/mol). This study suggests that compounds 648 and 2785758 can potentially affect the activity of ICL, leading to its inactivation and ultimately preventing the progression of tuberculosis.Communicated by Ramaswamy H. Sarma.

Network Pharmacology and Molecular Dynamics Simulations Reveal the Mechanism of Total Alkaloid Components in Anisodus tanguticus (Maxim.) Pascher in Treating Inflammation and Pain.

This study aimed to elucidate the mechanism that total alkaloids in Anisodus tanguticus (AT)(Maxim.) Pascher play anti-inflammatory and analgesic effects. In this paper, the anti-inflammatory effect in the total alkaloids of AT was confirmed via lipopolysaccharide (LPS)-induced inflammation model in RAW 264.7 cells and the main components of AT were immediately analyzed by UPLC/MS. Disease targets were obtained in GeneCards and DisGeNET. Targets of major compounds were searched in ETCM, TCMSP and other databases. The protein-protein interaction (PPI) network was constructed using STRING database, and Cytoscape was used for core targets screening. GO and KEGG enrichment analysis were performed using Daivid database. Sailvina was used for molecular docking. Molecular dynamics simulation analysis was performed using the Amber 20 program. The results showed that the main components in AT were anisodamine, atropine, fabiatrin, scopolamine, scopoletin and scopolin, possibly exerting anti-inflammatory and analgesic effects through pathways such as EGFR tyrosine kinase inhibitor resistance and IL-17 signaling pathway. Fabiatrin and scopolin could be potential drugs with good anti-inflammatory and analgesic effects.

Molecular Interactions of the Pioneer Transcription Factor GATA3 With DNA.

The GATA3 transcription factor is a pioneer transcription factor that is critical in the development, proliferation, and maintenance of several immune cell types. Identifying the detailed conformational dynamics and interactions of this transcription factor, as well as its clinically important population variants will allow us to unravel its mode of action. In this study, we analyze the molecular interactions of the GATA3 transcription factor bound to dsDNA as well as three clinically important population variants by atomistic molecular dynamics simulations. We identify the effect of the variants on the DNA conformational dynamics and delineate the differences compared to the wildtype transcription factor that could be related to impaired function. We highlight the structural plasticity in the binding of the GATA3 transcription factor and identify important DNA-protein contacts. Although the DNA-protein contacts are persistent and appear to be stable, they exhibit nanosecond timescale fluctuations and several binding/unbinding events. Further, we identify differential DNA binding in the three variants and show that the N-terminal binding is reduced in two of the variants. Our results indicate that reduced minor groove width and DNA diameter are important hallmarks for the binding of GATA3. Our work is an important step towards understanding the functional dynamics of the GATA3 protein and its clinically significant population variants.