A lesson for the maestro of the replication fork: Targeting the protein-binding interface of proliferating cell nuclear antigen for anticancer therapy

Designed peptide binders and nanobodies as PROTAC starting points for targeted degradation of PCNA and BCL6

The efficient degradation of pathogenic proteins, particularly proliferating cell nuclear antigen (PCNA) and B-cell lymphoma 6 protein (BCL6), is crucial for treating various diseases related to cancer. As key biological macromolecules, PCNA plays a critical role in DNA replication and repair, while BCL6 acts as a transcriptional repressor involved in B-cell lymphoma. To enhance the efficiency and specificity of protein degradation, we developed a RS80E-based bioPROTACs system that consists of truncated variants of Ring-B-boxed coiled-coil (RBCC) domains (RS80E) with improved degradation efficiency fused to an AI-driven binder/nanobody targeting specific antigens. Combining state-of-the-art methodologies such as ProteinMPNN, RFdiffusion, AlphaFold3, AlphaFold2, and HADDOCK, we identified binders for PCNA and predicted spatial interrelationships. Employing fragment-based and alanine scanning methods, we designed nanobodies targeting PCNA and BCL6 by combinatorially designing CDR3 and grafting them onto nanobody scaffolds. Significantly, our results demonstrate the utility of bioPROTACs in degrading PCNA and BCL6, thereby activating p53 and promoting apoptosis. This highlights the therapeutic potential of targeting PCNA and BCL6 degradation and lays the groundwork for developing PCNA and BCL6-degrading therapeutics. In summary, our system offers a modular and rapid pathway for exploration other intractable therapeutic targets, and emphasizes the importance of interdisciplinary methods in advancing therapeutic interventions.

Drug-Loaded Microspheres on NIR-Responsive PLA/MXene Scaffolds: Controlled Release and Bone Tissue Regeneration

The resection of bone tumors results in large bone defects with some residual tumor cells, and the treatment of this type of bone defect area often faces a dilemma, namely, the trade-off between bone repair and antitumor after the resection of bone tumors. In order to promote local bone repair, and at the same time inhibit tumor recurrence by continuous and controlled drug administration, we developed a multifunctional NIR-responsive scaffold, whose main components are polylactic acid and MXene, and loaded with PLGA/DOX microspheres, and we hope that the scaffold can take into account both antitumor and bone repair in the bidirectional modulation effect of NIR. The results showed that the scaffold with 1% MXene content had relatively good performance in photothermal therapy (PT) and other aspects, and it could be smoothly increased to 50 °C within 2 min under NIR illumination, and the drug release of microspheres was increased by 10% after illumination compared with that at body temperature. In vivo experiments in animals showed that this scaffold effectively limited the in situ recurrence of tumor cells and lung metastasis and was able to promote osteogenic differentiation under NIR irradiation. Therefore, this scaffold can not only control the release of antitumor drugs but also enhance the antitumor effect through the bidirectional modulation effect of PT and at the same time promote bone formation, which provides a good application solution for the integrated treatment of the bone defect area after bone tumor surgery.

Mechanisms of breast cancer treatment using Gentiana robusta: evidence from comprehensive bioinformatics investigation

This study investigates the potential treatment of breast cancer utilizing Gentiana robusta King ex Hook. f. (QJ) through an integrated approach involving network pharmacology, molecular docking, and molecular dynamics simulation. Building upon prior research on QJ's chemical constituents, we conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the DAVID database. Network interactions and core genes were identified using Cytoscape 3.9.1. Key target genes, including Interleukin-6 (IL-6), tumour suppressor gene P53 (TP53), and epidermal growth factor receptor (EGFR), were selected for molecular docking with QJ's active components, 2'-O-β-D-glucopyranosyl-gentiopicroside and macrophylloside D, employing Schrodinger Maestro 13.5. Molecular dynamics (MD) simulations were performed using the Desmond program. A total of 270 intersection targets of active ingredients and diseases were identified, with three core targets determined through network topology screening. Enrichment analysis highlighted the involvement of QJ in breast cancer treatment, primarily through the hsa05200 cancer signaling pathway and the hsa04066 HIF-1 signaling pathway. Molecular docking and dynamics simulations demonstrated the close interaction of 2'-O-β-D-glucopyranosyl-gentiopicroside (QJ17) and macrophylloside D (QJ25) with IL6, TP53, and EGFR, and other target genes, showcasing a stabilizing effect. In conclusion, this study unveils the effective components and potential mechanisms of 2'-O-β-D-glucopyranosyl-gentiopicroside and macrophylloside D in breast cancer prevention and treatment. The identified components act on target genes such as IL6, TP53, and EGFR, regulating crucial pathways including the cancer signaling and Hypoxia-inducible factor 1 (HIF-1) signaling pathways. These findings provide valuable insights into the therapeutic potential of QJ in breast cancer management. However, further experimental research are needed to validate the computational findings of QJ.

Molecular targets and mechanisms of Sijunzi decoction in the treatment of Parkinson's disease: evidence from network pharmacology, molecular docking, molecular dynamics simulation, and experimental validation

Aim

To explore the molecular mechanism of Sijunzi Decoction (SJZD) in the treatment of Parkinson's disease (PD) through the application of network pharmacology, molecular docking, and molecular dynamics simulations, complemented by experimental verification.

Methods

The BATMAN-TCM, GeneCards, and DisGeNet databases were searched to screen the active components and therapeutic targets of SJZD. Cytoscape (3.7.1) was used to create a network diagram of the components and targets. The STRING platform was used to construct a protein-protein interaction (PPI) network. The Bioconductor database and RX64 (4.0.0) software were used to conduct Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis on the core target genes. The binding sites and binding energies between SJZD active components and the target were analyzed by molecular docking and dynamic simulation. Finally, the therapeutic effect and mechanism of SJZD were verified by Cell Counting Kit-8 (CCK-8) and Western blotting (WB).

Results

This research identified 188 active compounds in SJZD, 1568 drug targets, 2069 PD targets, and 451 intersection targets related to PD. According to network analysis, Adenosine Triphosphate, Tridecanoic Acid, Hexadecanoic Acid, Pentadecanoic Acid, and Adenosine were identified as the core components of SJZD in the treatment of PD. The five targets with the highest Degree values in the PPI network were AKT1, INS, TNF, IL-6, and TP53. The GO and KEGG enrichment analyses, in turn, determined that the administration of SJZD for the treatment of PD may engage processes such as xenobiotic stimulation and biological stimulus response. Furthermore, AGE-RAGE and cAMP signaling pathways related to diabetic complications may be involved. Molecular docking and kinetic simulations showed that IL-6 and AKT1 bind best to Adenosine. Experimental results showed that SJZD significantly reduced 6-OHDA-induced apoptosis of SH⁃SY5Y cells by activating the PI3K/AKT signaling pathway and regulating the expression of apoptosis factors such as Bcl⁃2 and Bax.

Conclusion

SJZD is essential in the processes of apoptosis and neuronal protection, acting through various components that target multiple pathways. Notably, the PI3K/AKT pathway is a verified SJZD-PD target, providing a reference for clinical precision drug use for PD.

Do All Paths Lead to Rome? How Reliable is Umbrella Sampling Along a Single Path?

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

Comparison of interactions between alpha-lactalbumin and three protopanaxadiol ginsenosides: Impacts on the structure and antitumor properties

In this research, interactions between α-lactalbumin (ALA) and three protopanaxadiol ginsenosides [20(S)-Rg3, 20(S)-Rh2, and 20(S)-PPD] were compared to explore the effects of similar ligand on structure and cytotoxicity of ALA. Multi-spectroscopy revealed the binding between ALA and ginsenoside changed the conformation of ALA, which related to different structures and solubility of ligands. Scanning electron microscope illustrated that all ALA-ginsenoside complexes exhibited denser structures via hydrophobic interactions. Additionally, the cytotoxic experiments confirmed that the cytotoxicity of ginsenoside was enhanced after binding with ALA. Molecular docking showed all three ginsenosides were bound to the sulcus depression region of ALA via hydrogen bonding and hydrophobic interaction. Furthermore, molecular dynamics simulation elucidated the precise binding sites and pertinent system properties. Among all three composite systems, 20(S)-Rh2 had optimal binding affinity. These findings enhanced understanding of the synergistic utilization of ALA and ginsenosides as functional ingredients in food, medicine, and cosmetics.

A phenotypic screen with Trypanosoma brucei for discovering small molecules that target the SLiM-binding pocket of proliferating cell nuclear antigen orthologs

Proliferating cell nuclear antigen (PCNA) is a homo-trimeric protein complex that clamps around DNA to tether DNA polymerases to the template during replication and serves as a hub for many other interacting proteins. It regulates DNA metabolic processes and other vital cellar functions through the binding of proteins having short linear motifs (SLiMs) like the PIP-box (PCNA-interacting protein-box) or the APIM (AlkB homolog 2 PCNA-interacting motif) in the hydrophobic pocket where SLiMs bind. However, overproducing TbPCNA or human PCNA (hPCNA) in the pathogenic protist Trypanosoma brucei triggers a dominant-negative phenotype of arrested proliferation. The mechanism for arresting T. brucei proliferation requires the overproduced PCNA orthologs to have functional intact SLiM-binding pocket. Sight-directed mutagenesis studies showed that T. brucei overproducing PCNA variants with disrupted SLiM-binding pockets grew normally. We hypothesized that chemically disrupting the SLiM-binding pocket would restore proliferation in T. brucei, overproducing PCNA orthologs. Testing this hypothesis is the proof-of-concept for a T. brucei-based PCNA screening assay. The assay design is to discover bioactive small molecules that restore proliferation in T. brucei strains that overproduce PCNA orthologs, likely by disrupting interactions in the SLiM-binding pocket. The pilot screen for this assay discovered two hit compounds that linked to predetermined PCNA targets. Compound #1, a known hPCNA inhibitor, had selective bioactivity to hPCNA overproduced in T. brucei, validating the assay. Compound #6 had promiscuous bioactivity for hPCNA and TbPCNA but is the first compound discovered with bioactivity for inhibiting TbPCNA.

Understanding the Role of Activation Loop Mutants in Drug Efficacy for FLT3-ITD

The type III receptor tyrosine kinase FLT3 is a pivotal kinase for hematopoietic progenitor cell regulation, with significant implications in acute myeloid leukemia (AML) through mutations like internal tandem duplication (ITD). This study delves into the structural intricacies of FLT3, the roles of activation loop mutants, and their interaction with tyrosine kinase inhibitors. Coupled with this, the research leverages molecular contrastive learning and protein language modeling to examine interactions between small molecule inhibitors and FLT3 activation loop mutants. Utilizing the ConPLex platform, over 5.7 million unique FLT3 activation loop mutants-small molecule pairs were analyzed. The binding free energies of three inhibitors were assessed, and cellular apoptotic responses were evaluated under drug treatments. Notably, the introduction of the Xepto50 scoring system provides a nuanced metric for drug efficacy. The findings underscore the modulation of molecular interactions and cellular responses by Y842 mutations in FLT3-KD, highlighting the need for tailored therapeutic approaches in FLT3-ITD-related malignancies.

Elucidating the mutational impact in causing Niemann-Pick disease type C: an in silico approach

Niemann-Pick disease type C is a rare autosomal recessive of lysosomal storage disorder characterized by impaired intracellular lipid transport and has a tendency to accumulate the fatty acids and glycosphingolipids in a variety of neurovisceral tissues. This work includes computational tools to deciphere the mutational effect in NPC protein. The study initiated with the collection of 471 missense mutations from various databases, which were then analyzed using computational tools. The mutations (G549V, F703S, Q775P and L1244P) were said to be disease associated, altering the biophysical properties, in highly conserved regions and reduces the stability using several in silico methods and were subjected to molecular docking analysis. To analyze the ligand (Itraconazole: a small molecule of antifungal drug class, which is known to inhibit cholesterol export from lysosomes) activity Molecular docking study was performed for all the complex proteins. The average binding affinity was taken and found to be -10.76 kcal/mol (native) and -11.06 kcal/mol (Q775P was located in transmembrane region IV which impacts the sterol-sensing domain of the NPC1 protein and associated with a severe infantile neurological form). Finally, molecular dynamic simulation was performed in duplicate and trajectories were built for the backbone of the RMSD, RMSF, the number of intramolecular hydrogen bonds, the radius of gyration and the SSE percent for both the complex proteins. This work contributes to understand the effectiveness and may provide an insight on the stability of the drug with the complex variant structures.Communicated by Ramaswamy H. Sarma.

Phylogenetic investigation of Gammaproteobacteria proteins involved in exogenous long-chain fatty acid acquisition and assimilation

Background

The incorporation of exogenous fatty acids into the cell membrane yields structural modifications that directly influence membrane phospholipid composition and indirectly contribute to virulence. FadL and FadD are responsible for importing and activating exogenous fatty acids, while acyltransferases (PlsB, PlsC, PlsX, PlsY) incorporate fatty acids into the cell membrane. Many Gammaproteobacteria species possess multiple homologs of these proteins involved in exogenous fatty acid metabolism, suggesting the evolutionary acquisition and maintenance of this transport pathway.

Methods

This study developed phylogenetic trees based on amino acid and nucleotide sequences of homologs of FadL, FadD, PlsB, PlsC, PlsX, and PlsY via Mr. Bayes and RAxML algorithms. We also explored the operon arrangement of genes encoding for FadL. Additionally, FadL homologs were modeled via SWISS-MODEL, validated and refined by SAVES, Galaxy Refine, and GROMACS, and docked with fatty acids via AutoDock Vina. Resulting affinities were analyzed by 2-way ANOVA test and Tukey's post-hoc test.

Results

Our phylogenetic trees revealed grouping based on operon structure, original homolog blasted from, and order of the homolog, suggesting a more ancestral origin of the multiple homolog phenomena. Our molecular docking simulations indicated a similar binding pattern for the fatty acids between the different FadL homologs.

General significance

Our study is the first to illustrate the phylogeny of these proteins and to investigate the binding of various FadL homologs across orders with fatty acids. This study helps unravel the mystery surrounding these proteins and presents topics for future research.

Exploring the interactive mechanism of acarbose with the amylase SusG in the starch utilization system of the human gut symbiont Bacteroides thetaiotaomicron through molecular modeling

The α-amylase, SusG, is a principal component of the Bacteroides thetaiotaomicron (Bt) starch utilization system (Sus) used to metabolize complex starch molecules in the human gastrointestinal (GI) tract. We previously reported the non-microbicidal growth inhibition of Bt by the acarbose-mediated arrest of the Sus as a potential therapeutic strategy. Herein, we report a computational approach using density functional theory (DFT), molecular docking, and molecular dynamics (MD) simulation to explore the interactive mechanism between acarbose and SusG at the atomic level in an effort to understand how acarbose shuts down the Bt Sus. The docking analysis reveals that acarbose binds orthosterically to SusG with a binding affinity of -8.3 kcal/mol. The MD simulation provides evidence of conformational variability of acarbose at the active site of SusG and also suggests that acarbose interacts with the main catalytic residues via a general acid-base double-displacement catalytic mechanism. These results suggest that small molecule competitive inhibition against the SusG protein could impact the entire Bt Sus and eliminate or reduce the system's ability to metabolize starch. This computational strategy could serve as a potential avenue for structure-based drug design to discover other small molecules capable of inhibiting the Sus of Bt with high potency, thus providing a holistic approach for selective modulation of the GI microbiota.

Anti-breast cancer drugs targeting cell-surface glucose-regulated protein 78: a drug repositioning in silico study

Breast cancer (BC) is prevalent worldwide and is a leading cause of death among women. However, cell-surface glucose-regulated protein 78 (cs-GRP78) is overexpressed in several types of cancer and during pathogen infections. This study examines two well-known BC drugs approved by the FDA as BC treatments to GRP78. The first type consists of inhibitors of cyclin-based kinases 4/6, including abemaciclib, palbociclib, ribociclib, and dinaciclib. In addition, tunicamycin, and doxorubicin, which are among the most effective anticancer drugs for early and late-stage BC, are tested against GRP78. As (-)-epiGallocatechin gallate inhibits GRP78, it is also being evaluated (used as positive control). Thus, using molecular dynamics simulation approaches, this study aims to examine the advantages of targeting GRP78, which represents a promising cancer therapy regime. In light of recent advances in computational drug response prediction models, this study aimed to examine the benefits of GRP78 targeting, which represents a promising cancer therapy regime, by utilizing combined molecular docking and molecular dynamics simulation approaches. The simulated protein (50 ns) was docked with the drugs, then a second round of dynamics simulation was performed for 100 ns. After that, the binding free energies were calculated from 30 to 100 ns for each complex during the simulation period. These findings demonstrate the efficacy of abemaciclib, ribociclib, and tunicamycin in binding to the nucleotide-binding domain of the GRP78, paving the way for elucidating the mode of interactions between these drugs and cancer (and other stressed) cells that overexpress GRP78.Communicated by Ramaswamy H. Sarma.

Target validation and structure-based virtual screening to Discover potential lead molecules against the oncogenic NSD1 histone methyltransferase

The aim of the study was to validate Nuclear receptor-binding SET Domain NSD1 as a cancer drug target followed by the design of lead molecules against NSD1. TCGA clinical data, molecular expression techniques were used to validate the target and structure-based virtual screening was performed to design hits against NSD1. Clinical data analysis suggests the role of NSD1 in metastasis, prognosis and influence on overall survival in various malignancies. Furthermore, the mRNA and protein expression profile of NSD1 was evaluated in various cell lines. NSD1 was exploited as a target protein for in silico design of inhibitors using two major databases including ZINC15 and ChemDiv by structure-based virtual screening approach. Virtual screening was performed using the pharmacophore hypothesis designed with a protein complex S-adenosyl-l-methionine (SAM) as an endogenous ligand. Subsequently, a combined score was used to distinguish the top 10 compounds from the docking screened compounds having high performance in all four scores (docking score, XP, Gscore, PhaseScreenScore, and MMGBSA delta G Bind). Finally, the top three Zinc compounds were subjected to molecular dynamic simulation. The binding MMGBSA data suggests that ZINC000257261703 and ZINC000012405780 can be taken for in vitro and in vivo studies as they have lesser MMGBSA energy towards the cofactor binding site of NSD1 than the sinefungin. Our data validates NSD1 as a cancer drug target and provides promising structures that can be utilized for further lead optimization and rational drug design to open new gateways in the field of cancer therapeutics.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00158-0.

In Silico Characterization and Analysis of Clinically Significant Variants of Lipase-H (LIPH Gene) Protein Associated with Hypotrichosis

Hypotrichosis is an uncommon type of alopecia (hair loss) characterized by coarse scalp hair caused by the reduced or fully terminated activity of the Lipase-H (LIPH) enzyme. LIPH gene mutations contribute to the development of irregular or non-functional proteins. Because several cellular processes, including cell maturation and proliferation, are inhibited when this enzyme is inactive, the hair follicles become structurally unreliable, undeveloped, and immature. This results in brittle hair, as well as altered hair shaft development and structure. Because of these nsSNPs, the protein's structure and/or function may be altered. Given the difficulty in discovering functional SNPs in genes associated with disease, it is possible to assess potential functional SNPs before conducting broader population investigations. As a result, in our in silico analysis, we separated potentially hazardous nsSNPs of the LIPH gene from benign representatives using a variety of sequencing and architecture-based bioinformatics approaches. Using seven prediction algorithms, 9 out of a total of 215 nsSNPs were shown to be the most likely to cause harm. In order to distinguish between potentially harmful and benign nsSNPs of the LIPH gene, in our in silico investigation, we employed a range of sequence- and architecture-based bioinformatics techniques. Three nsSNPs (W108R, C246S, and H248N) were chosen as potentially harmful. The present findings will likely be helpful in future large population-based studies, as well as in drug discovery, particularly in the creation of personalized medicine, since this study provides an initial thorough investigation of the functional nsSNPs of LIPH.

Molecular mechanisms augmenting resistance to current therapies in clinics among cervical cancer patients

Cervical cancer (CC) is the fourth leading cause of cancer death (~ 324,000 deaths annually) among women internationally, with 85% of these deaths reported in developing regions, particularly sub-Saharan Africa and Southeast Asia. Human papillomavirus (HPV) is considered the major driver of CC, and with the availability of the prophylactic vaccine, HPV-associated CC is expected to be eliminated soon. However, female patients with advanced-stage cervical cancer demonstrated a high recurrence rate (50-70%) within two years of completing radiochemotherapy. Currently, 90% of failures in chemotherapy are during the invasion and metastasis of cancers related to drug resistance. Although molecular target therapies have shown promising results in the lab, they have had little success in patients due to the tumor heterogeneity fueling resistance to these therapies and bypass the targeted signaling pathway. The last two decades have seen the emergence of immunotherapy, especially immune checkpoint blockade (ICB) therapies, as an effective treatment against metastatic tumors. Unfortunately, only a small subgroup of patients (< 20%) have benefited from this approach, reflecting disease heterogeneity and manifestation with primary or acquired resistance over time. Thus, understanding the mechanisms driving drug resistance in CC could significantly improve the quality of medical care for cancer patients and steer them to accurate, individualized treatment. The rise of artificial intelligence and machine learning has also been a pivotal factor in cancer drug discovery. With the advancement in such technology, cervical cancer screening and diagnosis are expected to become easier. This review will systematically discuss the different tumor-intrinsic and extrinsic mechanisms CC cells to adapt to resist current treatments and scheme novel strategies to overcome cancer drug resistance.

Computational analysis of protein-ligand interaction by targeting a cell cycle restrainer

BACKGROUND AND OBJECTIVE

The cyclin-dependent kinases 4/6 (CDK4/6) are among the most crucial controllers of the cell cycle, and their abnormal activity may induce uncontrolled cell multiplication, leading to cancers. The FDA currently approved three CDK4/6 inhibitors, however, they are associated with a variety of side effects. Thus it is required to design/develop novel potent and safe CDK4/6 inhibitors.

METHODS

In the present work, we furnished an integrated in-silico approach followed by steered molecular dynamics (SMD) simulations to identify molecules that can be developed into novel CDK4/6 inhibitors.

RESULTS

Out of thirty-two 3-methyleneisoindolin-1-one molecules we selected top three M18, M24, and M32 molecules as potential drug candidates based on their respective interaction energies. According to the robust 250 ns MD simulations and thermodynamic free energy, M24 was the best molecule in comparison to palbociclib. In SMD, M24 required ∼205.587 kJ/mol/nm external pulling force, while palbociclib needed ∼160.97 kJ/mol/nm to dissociate from the binding pocket of the CDK4.

CONCLUSIONS

The high pulling force required for M24 dissociation from the binding site denotes stronger binding with CDK4. Therefore, M24 offers the possibility of a critical starting structure in developing effective CDK4 inhibitors.

Conformational dynamics of A30G α-synuclein that causes familial Parkinson disease

The first gene shown to be responsible for autosomal-dominant Parkinson's disease (PD) is the SNCA gene, which encodes for alpha synuclein (α-Syn). Recently, a novel heterozygous A30G mutation of the SNCA gene associated with familial PD has been reported. However, little research has been done on how the A30G mutation affects the structure of α-Syn. So, using atomistic molecular dynamics (MD) simulation, we demonstrate here the key structural characteristics of A30G α-Syn in the free monomer form and membrane associated state. From the MD trajectory analysis, the structure of A30G α-Syn was noticed to exhibit rapid conformational change, increase in backbone flexibility near the site of mutation and decrease in α-helical propensity. The typical torsion angles in residues (Val26 and Glu28) near the mutation site were observed to deviate significantly in A30G α-Syn. In the case of membrane bound A30G α-Syn, the regions that were submerged in the lipid bilayer (N-helix (3-37) and turn region (38-44)) found to contain higher helical content than the elevated region above the lipid surface. The bending angle in the helix-N and helix-C regions were noticed to be relatively higher in the free form of A30G α-Syn (38.50) than in the membrane bound form (370). The A30G mutation in α-Syn was predicted to have an impact on the stability and function of the protein based on ΔΔG values obtained from the online servers. Our results demonstrate that the A30G mutation in α-Syn altered the protein's α-helical structure and slightly altered the membrane binding.Communicated by Ramaswamy H. Sarma.

Leaf Extract of Perilla frutescens (L.) Britt Promotes Adipocyte Browning via the p38 MAPK Pathway and PI3K-AKT Pathway

The leaf of Perilla frutescens (L.) Britt (PF) has been reported to negatively affect adipocyte formation, inhibit body-fat formation, and lower body weight. However, its effect on adipocyte browning remains unknown. Thus, the mechanism of PF in promoting adipocyte browning was investigated. The ingredients of PF were acquired from the online database and filtered with oral bioavailability and drug-likeness criteria. The browning-related target genes were obtained from the Gene Card database. A Venn diagram was employed to obtain the overlapped genes that may play a part in PF promoting adipocyte browning, and an enrichment was analysis conducted based on these overlapped genes. A total of 17 active ingredients of PF were filtered, which may regulate intracellular receptor-signaling pathways, the activation of protein kinase activity, and other pathways through 56 targets. In vitro validation showed that PF promotes mitochondrial biogenesis and upregulates brite adipocyte-related gene expression. The browning effect of PF can be mediated by the p38 MAPK pathway as well as PI3K-AKT pathway. The study revealed that PF could promote adipocyte browning through multitargets and multipathways. An in vitro study validated that the browning effect of PF can be mediated by both the P38 MAPK pathway and the PI3K-AKT pathway.

Structure restoration and aggregate inhibition of V30M mutant transthyretin protein by potential quinoline molecules

Transthyretin (TTR) is a tetrameric protein found in human plasma and cerebrospinal fluid that functions as a transporter of thyroxine (T4) and retinol. A mutation resulting in the substitution of valine to methionine at position 30 (V30M) is the most common mutation that destabilizes the tetramer structure of TTR protein resulting in a fatal neuropathy known as TTR amyloidosis. The V30M TTR-induced neuropathy can be inhibited through stabilization of the TTR tetramer by the binding of small molecules. We accessed the potential of in-house synthesized quinoline molecules to stabilize the V30M TTR structure and analyzed the impact of protein-ligand interactions through molecular docking, molecular dynamics (MD) simulations, steered MD, and umbrella sampling simulations. This study revealed that the binding of quinoline molecules reverted back the structural changes including the residual flexibility, changes in secondary structural elements, and also restored the alterations in the electrostatic surface potential induced by the V30M mutation. Further, the top-most 4G and 4R molecules were compared with an FDA-approved drug (Tafamidis) and a reference quinoline molecule 14C. Here, we intend to suggest that the quinoline molecules could revert the structural changes, cease tetramer dissociation, prevent abnormal oligomerization and therefore could be developed as an effective therapeutics against TTR amyloidosis.

Insights into the binding mechanism between α-TOH and CYP4F2: A homology modeling, molecular docking, and molecular dynamics simulation study

α-Tocopherol (α-TOH) is a potent antioxidant. The concentrations of α-TOH in plasma are closely related to human health. α-TOH can be regulated by the metabolism of cytochrome P450 4F2 (CYP4F2). However, the atomic-level basis for this regulation process remains elusive. Here, we successfully constructed the structure of CYP4F2 by homology modeling and obtained the α-TOH-CYP4F2 complex models using molecular docking. Three parallel 500 ns molecular dynamics simulations were performed on each complex model to investigate the details of the interaction between α-TOH and CYP4F2. MM-GBSA method combined with principal component analysis shows that 8 key residues establish a hydrophobic cavity stabilizing α-TOH in the pocket of CYP4F2 and S423 forms an important hydrogen bond with α-TOH anchoring α-TOH in the favorable position for ω-hydroxylation. Based on our simulation results and the experimental facts, we designed mutation simulation experiments to clarify the important role of two key residues (S423 and V433) in the binding of α-TOH with CYP4F2. The results show that the mutations directly or indirectly change the binding mode of α-TOH and decrease its binding affinity with CYP4F2, which is unfavorable for ω-hydroxylation. Our results could enrich the information on structure-function relationships of CYP4F2 and provide valuable insights into the regulatory mechanism of CYP4F2 on the metabolism of α-TOH.

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, β-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of β-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (N^α-H---N^δ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

Elucidating activation and deactivation dynamics of VEGFR-2 transmembrane domain with coarse-grained molecular dynamics simulations

The vascular endothelial growth factor receptor 2 (VEGFR-2) is a member of receptor tyrosine kinases (RTKs) and is a dimeric membrane protein that functions as a primary regulator of angiogenesis. As is usual with RTKs, spatial alignment of its transmembrane domain (TMD) is essential toward VEGFR-2 activation. Experimentally, the helix rotations within TMD around their own helical axes are known to participate importantly toward the activation process in VEGFR-2, but the detailed dynamics of the interconversion between the active and inactive TMD forms have not been clearly elucidated at the molecular level. Here, we attempt to elucidate the process by using coarse grained (CG) molecular dynamics (MD) simulations. We observe that inactive dimeric TMD in separation is structurally stable over tens of microseconds, suggesting that TMD itself is passive and does not allow spontaneous signaling of VEGFR-2. By starting from the active conformation, we reveal the mechanism of TMD inactivation through analyzing the CG MD trajectories. We observe that interconversions between a left-handed overlay and a right-handed one are essential for the process of going from an active TMD structure to the inactive form. In addition, our simulations find that the helices can rotate properly when the overlaying structure of the helices interconverts and when the crossing angle of the two helices changes by larger than ~40 degrees. As the activation right after the ligand attachment on VEGFR-2 will take place in the reverse manner of this inactivation process, these structural aspects will also appear importantly for the activation process. The rather large change in helix configuration for activation also explains why VEGFR-2 rarely self-activate and how the activating ligand structurally drive the whole VEGFR-2. This mechanism of TMD activation / inactivation within VEGFR-2 may help in further understanding the overall activation processes of other RTKs.

Effects of the pathological E200K mutation on human prion protein: A computational screening and molecular dynamics approach

The human prion protein gene (PRNP) is mapped to the short arm of chromosome 20 (20pter-12). Prion disease is associated with mutations in the prion protein-encoding gene sequence. Earlier studies found that the mutation G127V in the PRNP increases protein stability. In contrast, the mutation E200K, which has the highest mutation rate in the prion protein, causes Creutzfeldt-Jakob disease (CJD) in humans and induces protein aggregation. We aimed to identify the structural mechanisms of E200k and G127V mutations causing CJD. We used a variety of bioinformatic algorithms, including SIFT, PolyPhen, I-Mutant, PhD-SNP, and SNP& GO, to predict the association of the E200K mutation with prion disease. MD simulation is performed, and graphs for root mean square deviation, root mean square fluctuation, radius of gyration, DSSP, principal component analysis, porcupine, and free energy landscape are generated to confirm and prove the stability of the wild-type and mutant protein structures. The protein is analyzed for aggregation, and the results indicate more fluctuations in the protein structure during the simulation owing to the E200K mutation; however, the G127V mutation makes the protein structure stable against aggregation during the simulation.

Computational Design of Plant-Based Antistress Agents Targeting Nociceptin Receptor

Stress is the body's reaction to the challenges it faces, and it produces a multitude of chemical molecules known as stressors as a result of these reactions. It's also a misalignment of the sympathetic and parasympathetic nervous systems causing changes in a variety of physiological reactions and perhaps leading to stress disorders. The reduction in neurotransmitter & neurohormonal hormones is mainly governed by the nociceptin receptor as G-protein coupled receptor and increased the level of reactive oxygen species. Various synthetic medicines that target nociceptin receptors were utilized to reduce the effects of stress but they come up with a variety of side effects. Because of the widespread utilization and renewed interest in medicinal herbal plants considered to be alternative antistress therapy. Our present work is an approach to decipher the molecular nature of novel herbal leads by targeting nociceptin receptor, under which herbal compounds were screened and validated through in-silico methods. Among screened leads, withanolide-B showed stable association in the active site of the nociceptin receptor as an antistress agent with no side effects. Furthermore, the selected lead was also evaluated for stability by molecular dynamic stimulation as well as for pharmacokinetics and toxicity profile. It has been concluded stable conformation of withanolide-B without presence of any major toxic effects. As a result, the in silico molecular docking technique is a highly successful method for selecting a prospective herbal lead molecule with respect to a specific target, and future research can pave the way for further exploration in the drug development field.

An In silico Approach towards Finding the Cancer-Causing Mutations in Human MET Gene

Mesenchymal-epithelial transition (MET) factor is a proto-oncogene encoding tyrosine kinase receptor with hepatocyte growth factor (HGF) or scatter factor (SF). It is found on the human chromosome number 7 and regulates the diverse cellular mechanisms of the human body. The impact of mutations occurring in the MET gene is demonstrated by their detrimental effects on normal cellular functions. These mutations can change the structure and function of MET leading to different diseases such as lung cancer, neck cancer, colorectal cancer, and many other complex syndromes. Hence, the current study focused on finding deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) and their subsequent impact on the protein's structure and functions, which may contribute to the emergence of cancers. These nsSNPs were first identified utilizing computational tools like SIFT, PROVEAN, PANTHER-PSEP, PolyPhen-2, I-Mutant 2.0, and MUpro. A total of 45359 SNPs of MET gene were accumulated from the database of dbSNP, and among them, 1306 SNPs were identified as non-synonymous or missense variants. Out of all 1306 nsSNPs, 18 were found to be the most deleterious. Moreover, these nsSNPs exhibited substantial effects on structure, binding affinity with ligand, phylogenetic conservation, secondary structure, and post-translational modification sites of MET, which were evaluated using MutPred2, RaptorX, ConSurf, PSIPRED, and MusiteDeep, respectively. Also, these deleterious nsSNPs were accompanied by changes in properties of MET like residue charge, size, and hydrophobicity. These findings along with the docking results are indicating the potency of the identified SNPs to alter the structure and function of the protein, which may lead to the development of cancers. Nonetheless, Genome-wide association study (GWAS) studies and experimental research are required to confirm the analysis of these nsSNPs.

Network Pharmacology, Molecular Docking, and Molecular Dynamic-Based Investigation on the Mechanism of Compound Chrysanthemum in the Treatment of Asthenopia

As a clinical empirical prescription for ophthalmology, compound chrysanthemum has been used gradually and has a good effect on eye fatigue. However, the detailed mechanisms of antiasthenopia have not been studied. In order to clarify the mechanisms of the compound chrysanthemum in the treatment of asthenopia, network pharmacology was combined with experimental study in this paper. A total of 593 genes and 39 active chemicals were identified, and both were considered to be essential to the advancement of asthenopia research. The results of the molecular docking analysis demonstrated a certain affinity between PRKACA, PRKCA, PRKCB, and their related compounds; molecular dynamic simulations assessed the stability of these receptors and ligands. The effects of compound chrysanthemum extract on ciliary muscle were studied in vitro and in vivo. By using the MTT assay, compound chrysanthemum extracts (50, 100, 200, 400, and 800 g·mL-1) showed no effect on the proliferation of rCSMCs for 24 and 48 hours. It raised nitric oxide and decreased Ca²⁺ in ciliary muscle cells isolated from the eyeballs of rats. Besides, compound chrysanthemum extract had a direct relaxing effect on the isolated gastric smooth muscle of rats by reducing the contractile tension. Furthermore, in vivo experiment results showed that, compared to the incandescent lamp-irradiated rats (model group), SD rats treated with compound chrysanthemum extracts (660 mg·kg^-1 and 1320 mg·kg^-1, orally) displayed considerably retracted pupils and increased NO content. It is also found that compound chrysanthemum extract can downregulate the mRNA expression of PKA and PKC in the calcium signaling pathway. Overall, our results suggested that compound chrysanthemum extract may lessen visual fatigue through multiple components, multiple targets, and multiple pathways.

Predicting and Validating the Mechanism of Qingyi II Granules in the Treatment of Acute Pancreatitis by Network Pharmacology

Network pharmacology, reverse molecular docking, and rat acute pancreatitis (AP) models were used to analyze the mechanism of protection by Qingyi II granules. The chemical components of 7 Chinese herbal medicines in Qingyi II granules were searched through the TCMSP (traditional Chinese medicine systems pharmacology database and analysis platform) database. The active ingredients were screened out in the OB (oral bioavailability) and DL (drug likeness) filters as a condition for inclusion. Then, the prediction analysis of potential targets was performed through databases. A GO (gene ontology) enrichment analysis of target proteins related to AP and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway annotation was performed using the DAVID (The Database for Annotation, Visualization, and Integrated Discovery) database. Finally, the “Herbal-Compound-Target” network was constructed using Cytoscape software. The active component structure and target name were uploaded to the Systems Dock database for reverse molecular docking. With octreotide as a positive control, Qingyi II decoction and Qingyi II granules were administered to AP rats at low, medium, and high doses. The pathological changes in the pancreas were observed using HE staining. The levels of Bcl-2, AMS, BAX, IL-2, and CASP3 in plasma were determined by an ELISA kit. Real-time PCR detected the expression of AKT1 and PIK3CA mRNA in the pancreas. The database predicted 94 active components of Qingyi II granules, 76 potential targets, and 64 signaling pathways. Twenty pathways were directly or indirectly associated with acute pancreatitis, including the TNF signaling pathway and the PI3K-AKT signaling pathway. In the reverse molecular docking experiment, the matching scores of the active components and the target were mainly between 6.0 and 7.0, with strong binding activity. Compared to the normal group, the plasma concentrations of BAX, IL-2, Bcl-2, AMS, and CASP3 in the model group were significantly increased ( $P<0.05$ ). Compared with the model group, the low-dose group of Qingyi II granules only significantly reduced IL-2 levels and had no effect on other indicators. The other groups could significantly reduce the levels of AMS, BAX, and CASP3 ( $P<0.05$ ). Compared with the model group, the octreotide group and Qingyi II granules high-dose group significantly increased the Bcl-2 level ( $P<0.05$ ), and there was no statistical difference in other drug-administered groups. Compared with the normal group, the expression of AKT1 and PIK3CA in the pancreas of the model group was significantly higher. Compared to the model group, the expression of PIK3CA was low in all drug-administered groups. In addition to the low-dose group, the other drug-administered groups significantly reduced the expression of AKT1. Qingyi can reduce the levels of AMS, BAX, IL-2, and CASP3 and increase the levels of Bcl-2. This mechanism may be related to the PI3K- AKT signaling pathway.

PHF8 promotes osteogenic differentiation of BMSCs in old rat with osteoporosis by regulating Wnt/β-catenin pathway

Osteoporosis is a progressive bone disorder with a higher incidence in the elderly and has become a major public health concern all over the world. Therefore, it is urgent to investigate the mechanisms underlying the pathogenesis of osteoporosis. In this study, the osteoporosis animal model was established, and then rat bone marrow mesenchymal stem cells (rBMSCs) were cultured. The results showed that PHF8 expression was decreased in osteoporosis rats compared to controls. Overexpression of PHF8 promoted BMSC osteogenic differentiation and the expression of osteogenesis-related genes. In addition, the Wnt/β-catenin signaling pathway in BMSCs was inhibited in osteoporosis rats, which was rescued by overexpression of PHF8. After treatment with the Wnt pathway antagonist, the improved osteogenic differentiation of BMSCs induced by overexpression of PHF8 was blocked. Collectively, our data revealed that the decreased expression of PHF8 in osteoporosis rats suppressed the osteogenic differentiation of BMSCs, which was then restored by PHF8 overexpression. Furthermore, the inhibition of the Wnt/β-catenin signaling pathway in BMSCs suppressed osteogenic differentiation. Thus, these findings indicated that PHF8 plays a role in osteogenic differentiation through the Wnt/β-catenin signaling pathway.

Point mutation consideration in CcO protein of the electron transfer chain by MD simulation

In Acidithiobacillus ferrooxidans, proteins such as CcO are present in the electron transport pathway. They cause ferrous iron oxidation to ferric leading to the electron release. CcO has two copper atoms (CuA, CuB). CuA plays an important role in electron transfer. According to previous studies, the conversion of histidine to methionine in a similar protein increased the redox potential and was directly related to the number of electrons received. Also, the binding of methionine 233 to CuA and CuB in the wild protein structure is the reason for the selection of the H230 M mutation in the CuA site. Then, wild-type and H230 M mutant were simulated in the presence of a bilayer membrane POPC using the gromacs version 5.1.4. The changes performed in the H230 M mutant were evaluated by MD simulations analyzes. CcO and CoxA proteins are the last two proteins in the chain and were docked by the PatchDock server. By H230 M mutation, the connection between CuA and M230 weakens. The M230 moves further away from CuA, resulting become more flexible. Therefore, the Methionine gets closer to E149 of the CoxA leading to the higher stability of the CcO/CoxA complex. The results of RMSF analysis at the mutation point showed a significant increase. This indicates more flexibility in the active site. And leads to an increase in E₀ in the mutation point, an increase in the rate of electron reception, and an improved bioleaching process.

Gαs and Gαq/11 protein coupling bias of two AVPR2 mutants (R68W and V162A) that cause nephrogenic diabetes insipidus

Loss-of-function mutations of the arginine vasopressin receptor 2 gene (AVPR2) cause Nephrogenic diabetes insipidus (NDI). AVPR2 is a kind of G protein coupled receptor (GPCR) and mainly couples with Gαs protein leading to cAMP accumulation in the cell as a secondary messenger. Recent studies showed that some AVPR2 mutations could cause biased Gαq/11 protein coupling rather than Gαs. Investigation into the characterization of biased receptors may give insights into the relationship between the conformational change of the receptor because of the mutation and related downstream signaling. In this study, R68W and V162A were analyzed to whether they show a bias to Gαs or Gαq/11 proteins. Their functionality in terms of cAMP production via Gαs protein coupling was decreased compared to the wild-type receptor. On the other hand, they showed the ability to couple with Gαq/11 protein and make Ca²⁺ mobilization at different levels in the cell. R68W showed bias to coupling with Gαq/11 protein rather than V162A and wild-type receptor. Studies about the Gα protein coupling bias of mutant AVPR2s may broaden our understanding of the relationship between the changed conformation of the receptor and consequently activated signaling pathways, and also may shed light on the development of more effective new therapeutics.

Virtual Screening and Network Pharmacology-Based Study to Explore the Pharmacological Mechanism of Clerodendrum Species for Anticancer Treatment

BACKGROUND

Cancer is a second leading cause of death in the world, killing approximately 3500 per million people each year. Therefore, the drugs with multitarget pharmacology based on biological networks are crucial to investigate the molecular mechanisms of cancer drugs and repurpose the existing drugs to reduce adverse effects. Clerodendrum is a diversified genus with a wide range of economic and pharmacological properties. Limited studies were conducted on the genus's putative anticancer properties and the mechanisms of action based on biological networks remains unknown. This study was aimed to construct the possible compound/target/pathway biological networks for anticancer effect of Clerodendrum sp. using docking weighted network pharmacological approach and to investigate its potential mechanism of action.

METHODS

A total of 194 natural Clerodendrum sp. Compounds were retrieved from public databases and screened using eight molecular descriptors. The cancer-associated gene targets were retrieved from databases and the function of the target genes with related pathways were examined. Cytoscape v3.7.2 was used to build three major networks: compound-target network, target-target pathway network, and compound-target-pathway network.

RESULTS

Our finding indicates that the anticancer activity of Clerodendrum sp. involves 6 compounds, 9 targets, and 63 signaling pathways, resulting in multicompounds, multitargets, and multipathways networks. Additionally, molecular dynamics (MD) simulations were used to estimate the binding affinity of the best hit protein-ligand complexes. Conclusion. This study suggests the potential anticancer activity of Clerodendrum sp. which could further contribute to scavenger novel compounds for the development of new alternative anticancer drugs.

Network pharmacology-based analysis of signaling pathways of an anti-osteoporotic triterpenoid from Acyranthes bidentata Blume root

In Korea folk remedies, Acyranthes bidentata Blume is a functional food plant to treat bone diseases; especially, its roots have been used to alleviate osteoporosis (OP), but its key chemical compound(s) and mechanism of action against osteoporosis have not reported yet. This study suggests that Acyranthes bidentata Blume root (ABBR) has promising compound(s) against OP. We utilized network pharmacology to evaluate the therapeutic value. The chemical compounds from Acyranthes bidentata Blume root (ABBR) were identified by gas chromatography-mass spectrum (GC-MS); their physicochemical properties have been evaluated by SwissADME. Next, the target(s) related to a triterpenoid or OP-related targets were investigated by public databases. The signaling pathways from final targets were visualized, constructed, and analyzed by RPackage. Finally, we performed a molecular docking (MD) to explore key target(s) and compound(s) by employing AutoDockVina tools; the residues of amino acids interacted with ligands were identified by LigPlot + v.22. A total of 24 chemicals were accepted by the Lipinski's rules. We found a sole triterpenoid from ABBR via GC-MS, suggesting that might be a potent ligand to alleviate OP. Thereby, the 42 targets were associated with the triterpenoid; the 19 targets among them were connected to OP-targets (1426). The final 19 targets were related directly to 8 signaling pathways on STRING database. On Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and a key signaling pathway (PPAR signaling pathway), four key targets (PPARA, PPARD, FABP3, and FABP4) and a key compound (Methyl 3β-hydroxyolean-18-en-28-oate) were selected via MD. Collectively, the triterpenoid from ABBR might have potent anti-osteoporotic efficacy by activating PPARA, PPARD, FABP3, and FABP4 on PPAR signaling pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03362-5.

Decoding the Key Functional Combined Components Group and Uncovering the Molecular Mechanism of Longdan Xiegan Decoction in Treating Uveitis

Objective

Longdan Xiegan Decoction (LXD) is a famous herbal formula in China. It has been proved that LXD has been shown to have a significant inhibitory effect on suppresses the inflammatory cells associated with uveitis. However, the key functional combination of component groups and their possible mechanisms remain unclear.

Methods

The community detecting model of the network, the functional response space, and reverse prediction model were utilized to decode the key components group (KCG) and possible mechanism of LXD in treating uveitis. Finally, MTT assay, NO assay and ELISA assay were applied to verify the effectiveness of KCG and the accuracy of our strategy.

Results

In the components-targets-pathogenic genes-disease (CTP) network, a combination of Huffman coding and random walk algorithm was used and eight foundational acting communities (FACs) were discovered with important functional significance. Verification has shown that FACs can represent the corresponding C-T network for treating uveitis. A novel node importance calculation method was designed to construct the functional response space and pick out 349 effective proteins. A total of 54 components were screened and defined as KCG. The pathway enrichment results showed that KCG and their targets enriched signal pathways of IL-17, Toll-like receptor, and T cell receptor played an important role in the pathogenesis of uveitis. Furthermore, experimental verification results showed that important KCG quercetin and sitosterol markedly inhibited the production of nitric oxide and significantly regulated the level of TNF-α and IFN-γ in Lipopolysaccharide-induced RAW264.7 cells.

Discussion

In this research, we decoded the potential mechanism of the multi-components-genes-pathways of LXD’s pharmacological action mode against uveitis based on an integrated pharmacology approach. The results provided a new perspective for the future studies of the anti-uveitis mechanism of traditional Chinese medicine.

Potential of paracetamol for reproductive disruption: Molecular interaction, dynamics and MM-PBSA based In-silico assessment

Paracetamol is generally recommended for pain and fever. However, as per experimental and epidemiological data, widespread and irrational or long-term use of paracetamol may be harmful to human endocrine homeostasis, especially during pregnancy. Some researchers suggest that prenatal exposure to paracetamol might alter foetal development and also enhance the risk of reproductive disorders. An imbalance in the levels of these hormones may play a significant role in the emergence of various diseases, including infertility. Therefore, in this study, the interaction mechanism of paracetamol with reproductive hormone receptors was investigated by molecular docking, molecular dynamics (MD) simulations, and poisson-Boltzmann surface area (MM-PBSA) for assessing paracetamol's potency to disrupt reproductive hormones. The results indicate that paracetamol has the ability to interact with reproductive hormone receptors (estrogen 1XP9; 1QKM with binding energy of -5.61 kcal/mol; -5.77 kcal/mol; androgen 5CJ6 -5.63 kcal/mol; and progesterone 4OAR-5.60 kcal/mol) by hydrogen bonds as well as hydrophobic and van der Waals interactions to maintain its stability. In addition, the results of the MD simulations and MM-PBSA confirm that paracetamol and reproductive receptor complexes are stable. This research provides a molecular and atomic level understanding of how paracetamols disrupt reproductive hormone synthesis. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), Radius of Gyration and hydrogen bonding exhibited that paracetamol mimic at various attribute to bisphenol and native ligand.

Network Toxicology Guided Mechanism Study on the Association between Thyroid Function and Exposures to Polychlorinated Biphenyls Mixture

Polychlorinated biphenyls (PCBs) are persistent and highly toxic pollutants, which can accumulate in organisms and produce toxic effects, especially damaging the function of thyroid hormones. So far, the molecular mechanism of PCBs mixture and their metabolites interfering with thyroid hormones has not been studied thoroughly except for individual compounds. In this study, PubMed, Web of Science, and STITCH databases were used to search PCBs and their corresponding target proteins. The intersection of PCBs and thyroid hormone dysfunction target proteins was obtained from GeneCards. The “compounds-targets-pathways” network was constructed by Cytoscape software. And KEGG and Go analyses were performed for key targets. Finally, molecular docking was used to verify the binding effect. Four major active components, five key targets, and 10 kernel pathways were successfully screened by constructing the network. Functional enrichment analysis showed that the interference was mediated by cancer, proteoglycans, PI3K-Akt, thyroid hormone, and FoxO signaling pathways. The molecular docking results showed that the binding energies were less than -5 kcal·mol-1. PCBs and their metabolites may act on the key targets of MAPK3, MAPK1, RXRA, PIK3R1, and TP53. The toxic effect of sulfated and methyl sulfone PCBs is greater. The method of screening targets based on the simultaneous action of multiple PCBs can provide a reference for other research. The targets were not found in previous metabolite toxicity studies. It also provides a bridge for the toxic effects and experimental research of PCBs and their metabolites in the future.

Screening of Sepsis Biomarkers Based on Bioinformatics Data Analysis

Methods

Gene expression profiles of GSE13904, GSE26378, GSE26440, GSE65682, and GSE69528 were obtained from the National Center for Biotechnology Information (NCBI). The differentially expressed genes (DEGs) were searched using limma software package. Gene Ontology (GO) functional analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and protein-protein interaction (PPI) network analysis were performed to elucidate molecular mechanisms of DEGs and screen hub genes.

Results

A total of 108 DEGs were identified in the study, of which 67 were upregulated and 41 were downregulated. 15 superlative diagnostic biomarkers (CCL5, CCR7, CD2, CD27, CD274, CD3D, GNLY, GZMA, GZMH, GZMK, IL2RB, IL7R, ITK, KLRB1, and PRF1) for sepsis were identified by bioinformatics analysis.

Conclusion

15 hub genes (CCL5, CCR7, CD2, CD27, CD274, CD3D, GNLY, GZMA, GZMH, GZMK, IL2RB, IL7R, ITK, KLRB1, and PRF1) have been elucidated in this study, and these biomarkers may be helpful in the diagnosis and therapy of patients with sepsis.

Network Pharmacology and Molecular Docking Analyses Unveil the Mechanisms of Yiguanjian Decoction against Parkinson’s Disease from Inner/Outer Brain Perspective

Objective

This study aims to explore the pharmacodynamic mechanism of Yiguanjian (YGJ) decoction against Parkinson's disease (PD) through integrating the central nervous (inner brain) and peripheral system (outer brain) relationship spectrum.

Methods

The active components of YGJ were achieved from the TCMSP, TCMID, and TCM@Taiwan databases. The blood-brain barrier (BBB) permeability of the active components along with their corresponding targets was evaluated utilizing the existing website, namely, SwissADME and SwissTargetPrediction. The targets of PD were determined through database retrieval. The interaction network was constructed upon the STRING database, followed by the visualization using Cytoscape software. Then, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses on potential targets. Finally, the molecular docking approach was employed to assess the binding affinity between key components and key targets.

Results

Overall, we identified 79 active components, 128 potential targets of YGJ, and 97 potential targets of YGJ-BBB potentially suitable for the treatment of PD. GO and KEGG analyses showed that the YGJ treatment of PD mainly relied on PI3K-Akt pathway while the YGJ-BBB was mostly involved in endocrine resistance. The molecular docking results displayed high affinity between multiple compounds and targets in accordance with previous observations.

Conclusions

Our study unveiled the potential mechanisms of YGJ against PD from a systemic perspective: (1) for the YGJ, they have potential exerting effects on the peripheral system and inhibiting neuronal apoptosis through regulating the PI3K-Akt pathway; (2) for the YGJ-BBB, they can directly modulate endocrine resistance of the central nervous and holistically enhance body resistance to PD along with YGJ on PI3K-Akt pathway.

Enhancement of Vitamin K2 Efflux in Bacillus subtilis Natto via a Potential Protein Receptor for Increased Yield

Bacillus subtilis is one of the few strains that can secrete synthetic menaquinone-7 (MK-7) to the outside of the cell, and its purpose and mechanism have not been clearly studied. As an amphiphilic protein naturally synthesized by Bacillus subtilis, the BslA protein may be involved in the inversion of extracellular vitamin K2 solubility. The protein structure in UniProt was used to search for the possible binding sites of MK-7, and the analysis of the higher ranking results of the genetic algorithm showed that the ASP166 residue was likely to be the binding site. They could form a stable hydrogen bond connection through ASP166, and approximately 7 proteins formed the conformation of a fixed naphthoquinone ring. We isolated and obtained the BslA protein by Ni-NTA affinity chromatography. Then, MK-7 was modified by BslA in vitro. A series of experiments, such as SEM, XPS, and WCA, showed that MK-7 and BslA proteins can realize self-assembly and transform from fat-soluble to water-soluble complexes. When the bslA protein in Bacillus subtilis natto was overexpressed, its MK-7 synthesis ability was further improved, especially the extracellular MK-7 content, which increased by 16%. This finding suggests that the BslA protein in Bacillus subtilis is likely to be involved in the extracellular secretion of MK-7 as a receptor.

Exploring the Role of Chemical Reactions in the Selectivity of Tyrosine Kinase Inhibitors

A variety of diseases are associated with tyrosine kinase enzymes that activate many proteins via signal transduction cascades. The similar ATP-binding pockets of these tyrosine kinases make it extremely difficult to design selective covalent inhibitors. The present study explores the contribution of the chemical reaction steps to the selectivity of the commercialized inhibitor acalabrutinib over the Bruton's tyrosine kinase (BTK) and the interleukin-2-inducible T-cell kinase (ITK). Ab initio and empirical valence bond (EVB) simulations of the two kinases indicate that the most favorable reaction path involves a water-assisted mechanism of the 2-butynamide reactive group of acalabrutinib. BTK reacts with acalabrutinib with a substantially lower barrier than ITK, according to our calculated free-energy profile and kinetic simulations. Such a difference is due to the microenvironment of the active site, as further supported by a sequence-based analysis of specificity determinants for several commercialized inhibitors. Our study involves a new approach of simulating directly the IC50 and inactivation efficiency k_eff, instead of using the standard formulas. This new strategy is particularly important in studies of covalent inhibitors with a very exothermic bonding step. Overall, our results demonstrate the importance of understanding the chemical reaction steps in designing selective covalent inhibitors for tyrosine kinases.

Therapeutic Mechanism and Key Active Ingredients of Shenfu Injection in Sepsis: A Network Pharmacology and Molecular Docking Approach

At present, although the early treatment of sepsis is advocated, the treatment effect of sepsis is unsatisfactory, and the mortality rate remains high. Shenfu injection (SFI) has been used to treat sepsis with good clinical efficacy. Based on network pharmacology, this study adopted a new research strategy to identify the potential therapeutic targets and key active ingredients of SFI for sepsis from the perspective of the pathophysiology of sepsis. This analysis identified 28 active ingredients of SFI based on UHPLC-QQQ MS, including 18 ginsenosides and 10 aconite alkaloids. 59 targets were associated with the glycocalyx and sepsis pathways. Based on the number of targets related to the pathophysiological process of sepsis, we identified songorine, ginsenoside Rf, ginsenoside Re, and karacoline as the key active ingredients of SFI for the treatment of sepsis. According to the cluster analysis of MCODE and the validation on the GEO dataset, LGALS3, BCHE, AKT1, and IL2 were identified as the core targets. This study further explored the therapeutic mechanism and the key active ingredients of SFI in sepsis and provided candidate compounds for drug development.

Synthesis, Biological Evaluation, and In Silico Studies of Novel Coumarin-Based 4H,5H-pyrano[3,2-c]chromenes as Potent β-Glucuronidase and Carbonic Anhydrase Inhibitors

The search for novel heterocyclic compounds with a natural product skeleton as potent enzyme inhibitors against clinical hits is our prime concern in this study. Here, a simple and facile two-step strategy has been designed to synthesize a series of novel coumarin-based dihydropyranochromenes (12a-12m) in a basic moiety. The synthesized compounds were thus characterized through spectroscopic techniques and screened for inhibition potency against the cytosolic hCA II isoform and β-glucuronidase. Few of these compounds were potent inhibitors of hCA II and β-glucuronidase with varying IC₅₀ values ranging from 4.55 ± 0.22 to 21.77 ± 3.32 μM and 440.1 ± 1.17 to 971.3 ± 0.05 μM, respectively. Among the stream of synthesized compounds, 12e and 12i were the most potent inhibitors of β-glucuronidase, while 12h, 12i, and 12j showed greater potency against hCA II. In silico docking studies illustrated the significance of substituted groups on the pyranochromene skeleton and binding pattern of these highly potent compounds inside enzyme pockets.

Identification of potential inhibitors, conformational dynamics, and mechanistic insights into mutant Kirsten rat sarcoma virus (G13D) driven cancers

The mutations at the hotspot region of K-Ras result in the progression of cancer types. Our study aimed to explore the small molecule inhibitors against the G13D mutant K-Ras model with anti-cancerous activity from food and drug administration (FDA)-approved drug compounds. We implemented several computational strategies such as pharmacophore-based virtual screening, molecular docking, absorption, distribution, metabolism and excretion features, and molecular simulation to ensure the identified hit compounds have potential efficacy against G13D K-Ras. We found that the FDA-approved compounds, namely, azelastine, dihydrocodeine, paroxetine, and tramadol, are potential candidates to inhibit the action of G13D mutant K-Ras. All four compounds exhibited similar binding patterns of sotorasib, and a structural binding mechanism with significant hydrophobic contacts. The descriptor features from the QikProp of all four compounds are within allowable limits compared to sotorasib drug. Consequently, a molecular simulation result emphasized that the dihydrocodeine and tramadol exhibited less fluctuation, minimal basin, significant h-bonds, and potent inhibition against G13D K-Ras. As a result, the current research identifies prospective K-Ras inhibitors that could be further improved with biochemical analysis for precision medicine against K-Ras-driven cancers.

Investigating the Role of Dahuang in Hepatoma Treatment Using Network Pharmacology, Molecular Docking, and Survival Analysis

Hepatoma is one of the most common malignant tumors. The incidence rate is high in developing countries, and China has the most significant number of cases. Dahuang is a classic traditional antitumor drug commonly used in China and has also been applied to treat hepatoma. However, the potential mechanism of Dahuang in treating hepatoma is not clear. Therefore, this study is aimed at elucidating the possible molecular mechanism and key targets of Dahuang using methods of network pharmacology, molecular docking, and survival analysis. Firstly, the active ingredients and key targets of Dahuang were analyzed through public databases, and then the drug-ingredient-target-disease network diagram of Dahuang against hepatoma was constructed. Five main active components and five core targets were determined according to the enrichment degree. Enrichment analysis demonstrated that Dahuang treated hepatoma through the multiple pathways in cancer. Additionally, molecular docking predicted that aloe-emodin and PIK3CG depicted the best binding energy. Survival analysis indicated that a high/ESR1 gene expression had a relatively good prognosis for patients with hepatoma (p < 0.05). In conclusion, the current study results demonstrated that Dahuang could treat hepatoma through a variety of active ingredients, targets, and multiantitumor pathways. Moreover, it effectively improved the prognosis of hepatoma patients. ESR1 is the potential key gene that is beneficial for the survival of hepatoma patients. Also, aloe-emodin and beta-sitosterol are the two main active crucial ingredients for hepatoma treatment. The study also provided some functional bases and references for the development of new drugs, target mining, and experimental animal research of hepatoma in the future.

Comparative molecular dynamics analyses on PIK3CA hotspot mutations with PI3Kα specific inhibitors and ATP

PI3K pathway is heavily emphasized in cancer where PIK3CA, which encodes for the p110α subunit of PI3Kα, presents itself as the second most common mutated gene. A lot of effort has been put in developing PI3K inhibitors, opening promising avenues for the treatment of cancer. Among these, PI3Kα specific inhibitor alpelisib was approved by FDA for breast cancer and other α-isoform specific inhibitors such as inavolisib and serabelisib reached clinical trials. However, the mode of action of these inhibitors on mutated PI3Kα and how they interact with mutant structures has not been fully elucidated yet. In this study, we are revealing the calculated interactions and binding affinities of these inhibitors within the context of PIK3CA hotspot mutations (E542K, E545K and H1047R) by employing molecular dynamics (MD) simulations. We performed principal component analysis to understand the motions of the protein complex during our simulations and also checked the correlated motions of all amino acids. Binding affinity calculations with MM-PBSA confirmed the consistent binding of alpelisib across mutations and revealed relatively higher affinities for inavolisib towards wild-type and H1047R mutant structures in comparison to other inhibitors. On the other hand, E542K mutation significantly impaired the interaction of inavolisib and serabelisib with PI3Kα. We also investigated the structural relationship of the natural ligand ATP with PI3Kα, and interestingly realized a significant reduction in binding affinity for the mutants, with potentially unexpected implications on the mechanisms that render these mutations oncogenic. Moreover, correlated motions of all residues were generally higher for ATP except the H1047R mutation which exhibited a distinguishable reduction. The results presented here could be guiding for pre-clinical and clinical studies of personalized medicine where individual mutations are a strong consideration point.

The soluble N-terminal autoinhibitory module of the A1 domain in von Willebrand factor partially suppresses its catch bond with glycoprotein Ibα in a sandwich complex

von Willebrand factor (VWF) senses and responds to the hemodynamic forces to interact with the circulatory system and platelets in hemostasis and thrombosis. The dark side of this mechanobiology is implicated in atherothrombosis, stroke, and, more recently, the COVID-19 thrombotic symptoms. The force-responsive element controlling VWF activation predominantly resides in the N terminal auto-inhibitory module (N-AIM) flanking its A1 domain. Nevertheless, the detailed mechano-chemistry of soluble VWF N-AIM is poorly understood at the sub-molecular level as it is assumed to be unstructured loops. Using the free molecular dynamics (MD) simulations, we first predicted a hairpin-like structure of the soluble A1 N-AIM derived polypeptide (Lp; sequences Q1238-E1260). Then we combined molecular docking and steered molecular dynamics (SMD) simulations to examine how Lp regulates the A1-GPIbα interaction under tensile forces. Our simulation results indicate that Lp suppresses the catch bond in a sandwich complex of A1-Lp-GPIbα yet contributes an additional catch-bond residue D1249. To experimentally benchmark the binding kinetics for A1-GPIbα in the absence or presence of Lp, we conducted the force spectroscopy-biomembrane force probe (BFP) assays. We found similar suppression on the A1-GPIbα catch bond with soluble Lp in presence. Clinically, as more and more therapeutic candidates targeting the A1-GPIbα axis have entered clinical trials to treat patients with TTP and acute coronary syndrome, our work represents an endeavor further towards an effective anti-thrombotic approach without severe bleeding side effects as most existing drugs suffer.