Is It Reliable to Take the Molecular Docking Top Scoring Position as the Best Solution without Considering Available Structural Data?

Exploring fatty acids from royal jelly as a source of histone deacetylase inhibitors: from the hive to applications in human well-being and health

A differential diet with royal jelly (RJ) during early larval development in honeybees shapes the phenotype, which is probably mediated by epigenetic regulation of gene expression. Evidence indicates that small molecules in RJ can modulate gene expression in mammalian cells, such as the fatty acid 10-hydroxy-2-decenoic acid (10-HDA), previously associated with the inhibition of histone deacetylase enzymes (HDACs). Therefore, we combined computational (molecular docking simulations) and experimental approaches for the screening of potential HDAC inhibitors (HDACi) among 32 RJ-derived fatty acids. Biochemical assays and gene expression analyses (Reverse Transcriptase - quantitative Polymerase Chain Reaction) were performed to evaluate the functional effects of the major RJ fatty acids, 10-HDA and 10-HDAA (10-hydroxy-decanoic acid), in two human cancer cell lines (HCT116 and MDA-MB-231). The molecular docking simulations indicate that these fatty acids might interact with class I HDACs, specifically with the catalytic domain of human HDAC2, likewise well-known HDAC inhibitors (HDACi) such as SAHA (suberoylanilide hydroxamic acid) and TSA (Trichostatin A). In addition, the combined treatment with 10-HDA and 10-HDAA inhibits the activity of human nuclear HDACs and leads to a slight increase in the expression of HDAC-coding genes in cancer cells. Our findings indicate that royal jelly fatty acids collectively contribute to HDAC inhibition and that 10-HDA and 10-HDAA are weak HDACi that facilitate the acetylation of lysine residues of chromatin, triggering an increase in gene expression levels in cancer cells.

Computational screening of chemical constituents derived from berry fruits as allosteric caspace-3/-7 inhibitors

With the aim of finding the plant-derived allosteric inhibitors of caspase-3/-7, we conducted computational investigations of bioactive compounds present in various berry fruits. In a molecular docking study, perulactone demonstrated excellent binding affinity scores of -12.1 kcal/mol and -9.1 kcal/mol for caspase 7 and 3, respectively, whereas FDA-approved allosteric inhibitors (DICA and FICA) were found to show lower docking scores (-5.6 and -6.1 kcal/mol) against caspase 7 while (-5.0 and -5.1 kcal/mol) for caspase 3, respectively. MD simulations were used to validate the binding stability of perulactone in the active sites of caspase-7/-3, and the results showed outstanding stability with lower ligand RMSDs of 1.270-3.088 Å and 2.426-9.850 Å against the targeted receptor. Furthermore, we performed MMGBSA free binding energy, where the perulactone values of ΔG Bind were determined to be -63.98 kcal/mol and -66.32 kcal/mol for both receptors (3IBF and 1NME), which are significantly better than the -45.16 kcal/mol and -39.51 kcal/mol for DICA as well as -26.37 kcal/mol and -15.50 kcal/mol for FICA, respectively. The drug resemblance of perulactone was effectively evaluated by ADMET. Thus, our findings indicated that perulactone could be an orally administered therapeutic candidate for regulating apoptosis in a variety of disorders. However, there may be an urgent need to study using in vitro and in vivo experiments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04067-7.

Screening of promising molecules against potential drug targets in Yersinia pestis by integrative pan and subtractive genomics, docking and simulation approach

This study focuses on Yersinia pestis, the bacterium responsible for plague, which posed a severe threat to public health in history. Despite the availability of antibiotics treatment, the emergence of antibiotic resistance in this pathogen has increased challenges of controlling the infections and plague outbreaks. The development of new drug targets and therapies is urgently needed. This research aims to identify novel protein targets from 28 Y. pestis strains by the integrative pan-genomic and subtractive genomics approach. Additionally, it seeks to screen out potential safe and effective alternative therapies against these targets via high-throughput virtual screening. Targets should lack homology to human, gut microbiota, and known human 'anti-targets', while should exhibit essentiality for pathogen's survival and virulence, druggability, antibiotic resistance, and broad spectrum across multiple pathogenic bacteria. We identified two promising targets: the aminotransferase class I/class II domain-containing protein and 3-oxoacyl-[acyl-carrier-protein] synthase 2. These proteins were modeled using AlphaFold2, validated through several structural analyses, and were subjected to molecular docking and ADMET analysis. Molecular dynamics simulations determined the stability of the ligand-target complexes, providing potential therapeutic options against Y. pestis.

Rigidifying Flexible Sites: A Promising Strategy to Improve Thermostability of Lysophospholipase From Pyrococcus abyssi

High thermostability of the enzymes is one of the distinguishing characteristics that increase their industrial utility. In the current research work, rigidifying the flexible amino acid residues of a lysophospholipase (Pa-LPL) from Pyrococcus abyssi was used as a protein engineering approach to improve its thermostability. A truncated variant of Pa-LPL (t-LPL∆12) was constructed by trimming its 12 amino acid residues (50-61) through overlap extension PCR. The truncated enzyme worked optimally at 65°C and pH 6.5 with remarkable thermostability at 65°C-85°C. In comparison to wild-type Pa-LPL, 5.8 and 1.2-fold increase in half-life (t1/2) of t-LPL∆12 was observed at 65 (optimum temperature) and 95°C, respectively. The activity of t-LPL∆12 was stimulated by 1 mM Cu2+ followed by Ca2+, Ni2+, Co2+, and Mg2+. Both substrate docking and experimental results indicated that the truncated enzyme could hydrolyze a variety of p-nitrophenyl esters. Km, Vmax, and Kcat for enzymatic hydrolysis of p-nitrophenyl butyrate were calculated to be 1 ± 0.087 mM, 1456 ± 36.474 U/mg, and 1.397 × 1011 min-1, respectively. In short, broad substrate specificity and thermostability of t-LPL∆12 are some of the distinctive features that make it an ideal candidate for degumming of vegetable oils.

Exploring Triazole-Connected Steroid-Pyrimidine Hybrids: Synthesis, Spectroscopic Characterization, and Biological Assessment

Molecules originating from natural sources are physicochemically and biologically diverse. The conjugation of two active biomolecules has become the foundation for medical and pharmaceutical sciences. An effective synthesis of 11 new steroid-pyrimidine conjugates containing 1,2,3-triazole rings was carried out. The group of 3α-OH bile acids (lithocholic, deoxycholic, cholic) and 3β-OH sterols (cholesterol, cholestanol) were respectively modified to azidoacetates. 2-thiouracil was converted into N(1)S and N(3)S dipropargyl derivatives. Azide-alkyne cycloaddition in the presence of copper(I) of the obtained compounds led to the preparation of 1,2,3-triazole derivatives. Based on a series of spectroscopic (1H NMR, 13C NMR, Fourier-transform infrared (FT-IR)), spectrometric analyses (Electrospray ionization-mass spectrometry (ESI-MS), electron impact-mass spectrometry (EI-MS)), and semiempirical calculations, the structures of all compounds were confirmed. In silico biological tests and molecular docking (for domain 1KZN, 2H94, 5V5Z, 1EZF, 2Q85) were performed for selected compounds. The tests performed indicate the theoretical antimicrobial potential of the obtained ligands.

Discovery of Indole-Thiourea Derivatives as Tyrosinase Inhibitors: Synthesis, Biological Evaluation, Kinetic Studies, and In Silico Analysis

Tyrosinase, a key enzyme in melanin synthesis, represents a crucial therapeutic target for hyperpigmentation disorders due to excessive melanin production. This study aimed to design and evaluate a series of indole-thiourea derivatives by conjugating thiosemicarbazones with strong tyrosinase inhibitory activity to indole. Among these derivatives, compound 4b demonstrated tyrosinase inhibitory activity with an IC50 of 5.9 ± 2.47 μM, outperforming kojic acid (IC50 = 16.4 ± 3.53 μM). Kinetic studies using Lineweaver-Burk plots confirmed competitive inhibition by compound 4b. Its favorable ADMET and drug-likeness properties make compound 4b a promising therapeutic candidate with a reduced risk of toxicity. Molecular docking revealed that the compounds bind strongly to mushroom tyrosinase (mTYR) and human tyrosinase-related protein 1 (TYRP1), with compound 4b showing superior binding energies of -7.0 kcal/mol (mTYR) and -6.5 kcal/mol (TYRP1), surpassing both kojic acid and tropolone. Molecular dynamics simulations demonstrated the stability of the mTYR-4b complex with low RMSD and RMSF and consistent Rg and SASA values. Persistent strong hydrogen bonds with mTYR, along with favorable Gibbs free energy and MM/PBSA calculations (-19.37 kcal/mol), further support stable protein-ligand interactions. Overall, compound 4b demonstrated strong tyrosinase inhibition and favorable pharmacokinetics, highlighting its potential for treating pigmentary disorders.

Novel peptide and hyaluronic acid coated biomimetic liposomes for targeting bacterial infections and sepsis

Sepsis is a life-threatening syndrome resulting from an imbalanced immune response to severe infections. Despite advances in nanomedicines, effective treatments for sepsis are still lacking. Herein, vancomycin free base (VCM)-loaded dual functionalized biomimetic liposomes based on a novel TLR4-targeting peptide (P3) and hyaluronic acid (HA) (HA-P3-Lipo) were developed to enhance sepsis therapy. The nanocarrier revealed appropriate physicochemical parameters, good stability, and biocompatibility. The release of VCM from HA-P3-Lipo was found to be sustained with 76 % VCM released in 48 h. The biomimicry was elucidated by in silico tools and MST and results confirmed strong binding between the system and TLR4. Furthermore, HA-P3-Lipo revealed 2-fold enhanced antibacterial activity against S. aureus, sustained antibacterial activity against MRSA over 72 h and 5-fold better MRSA biofilm inhibition compared to bare VCM. Bacterial-killing kinetics and flow cytometry confirmed the superiority of HA-P3-Lipo in eliminating MRSA faster than VCM. The in vivo potential of the nanocarrier was elucidated in an MRSA-induced sepsis mice model, and the results confirmed the superiority of HA-P3-Lipo compared to free VCM in eliminating bacteria and down-regulating the proinflammatory markers. Therefore, HA-P3-Lipo exhibits potential as a promising novel multi-functional nanosystem against sepsis and could significantly contribute to the transformation of sepsis therapy.

Identification of Hub of the Hub-Genes From Different Individual Studies for Early Diagnosis, Prognosis, and Therapies of Breast Cancer

Breast cancer (BC) is a complex disease, which causes of high mortality rate in women. Early diagnosis and therapeutic improvements may reduce the mortality rate. There were more than 74 individual studies that have suggested BC-causing hub-genes (HubGs) in the literature. However, we observed that their HubG sets are not so consistent with each other. It may be happened due to the regional and environmental variations with the sample units. Therefore, it was required to explore hub of the HubG (hHubG) sets that might be more representative for early diagnosis and therapies of BC in different country regions and their environments. In this study, we selected top-ranked 10 HubGs (CCNB1, CDK1, TOP2A, CCNA2, ESR1, EGFR, JUN, ACTB, TP53, and CCND1) as the hHubG set by the protein-protein interaction network analysis based on all of 74 individual HubG sets. The hHubG set enrichment analysis detected some crucial biological processes, molecular functions, and pathways that are significantly associated with BC progressions. The expression analysis of hHubGs by box plots in different stages of BC progression and BC prediction models indicated that the proposed hHubGs can be considered as the early diagnostic and prognostic biomarkers. Finally, we suggested hHubGs-guided top-ranked 10 candidate drug molecules (SORAFENIB, AMG-900, CHEMBL1765740, ENTRECTINIB, MK-6592, YM201636, masitinib, GSK2126458, TG-02, and PAZOPANIB) by molecular docking analysis for the treatment against BC. We investigated the stability of top-ranked 3 drug-target complexes (SORAFENIB vs ESR1, AMG-900 vs TOP2A, and CHEMBL1765740 vs EGFR) by computing their binding free energies based on 100-ns molecular dynamic (MD) simulation based Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach and found their stable performance. The literature review also supported our findings much more for BC compared with the results of individual studies. Therefore, the findings of this study may be useful resources for early diagnosis, prognosis, and therapies of BC.

Staphylococcus aureus major cell division protein FtsZ assembly is inhibited by silibinin, a natural flavonolignan that also blocked bacterial growth and biofilm formation

The bacterial cell division protein FtsZ has been considered a potential therapeutic target due to its rapid treadmilling that induces cellular wall construction in bacteria. The current study discovered a novel antimicrobial compound, silibinin, a natural flavonolignan and its impact on the recombinant S. aureus FtsZ (SaFtsZ). Silibinin inhibited S. aureus Newman growth in a dose-dependent manner. The IC50 and MIC values for silibinin were 75 μM and 200 μM, respectively. It had no cytotoxicity against HEK293 cells in vitro. Silibinin also enlarged the bacterial cell morphology by ∼40 folds and showed antibiofilm property. It perturbed the S. aureus membrane potential both at IC50 conc. and at MIC conc. Further, it inhibited both the polymerization and GTPase activity of SaFtsZ. It did not inhibit tubulin assembly, a eukaryotic FtsZ homolog. A fluorescence quenching study yielded the Kd value for SaFtsZ-Silibinin interaction and binding stoichiometry 0.857 ± 0.188 μM and 1:1, respectively. Both in silico study and competition assay indicated that silibinin binds at the GTP binding site on SaFtsZ. The Ki value for the silibinin-mediated inhibition of SaFtsZ was 8.8 μM. Therefore, these findings have comprehensively shown the antimicrobial behavior of silibinin on S. aureus Newman cells targeting SaFtsZ.

Identification of mycobacterial Thymidylate kinase inhibitors: a comprehensive pharmacophore, machine learning, molecular docking, and molecular dynamics simulation studies

Thymidylate kinase (TMK) is a pivotal enzyme in Mycobacterium tuberculosis (Mtb), crucial for phosphorylating thymidine monophosphate (dTMP) to thymidine diphosphate (dTDP), thereby playing a critical role in DNA biosynthesis. Dysregulation or inhibition of TMK activity disrupts DNA replication and cell division, making it an attractive target for anti-tuberculosis drug development. In this study, the statistically validated pharmacophore mode was developed from a set of known TMK inhibitors. Further, the robust pharmacophore was considered for screening the Enamine database. The chemical space was reduced through multiple molecular docking approaches, pharmacokinetics, and absolute binding energy estimation. Two different molecular docking algorithms favor the strong binding affinity of the proposed molecules towards TMK. Machine learning-based absolute binding energy also showed the potentiality of the proposed molecules. The binding interactions analysis exposed the strong binding affinity between the proposed molecules and active site amino residues of TMK. Several statistical parameters from all atoms MD simulation explained the stability between proposed molecules and TMK in the dynamic states. The MM-GBSA approach also found a strong binding affinity for each proposed molecule. Therefore, the proposed molecules might be crucial TMK inhibitors for managing Mtb inhibition subjected to in vitro/in vivo validations.

Recombination Events Among SARS-CoV-2 Omicron Subvariants: Impact on Spike Interaction With ACE2 Receptor and Neutralizing Antibodies

The recombination plays a key role in promoting evolution of RNA viruses and emergence of potentially epidemic variants. Some studies investigated the recombination occurrence among SARS-CoV-2, without exploring its impact on virus-host interaction. In the aim to investigate the burden of recombination in terms of frequency and distribution, the occurrence of recombination was first explored in 44 230 Omicron sequences among BQ subvariants and the under investigation "ML" (Multiple Lineages) denoted sequences, using 3seq software. Second, the recombination impact on interaction between the Spike protein and ACE2 receptor as well as neutralizing antibodies (nAbs), was analyzed using docking tools. Recombination was detected in 56.91% and 82.20% of BQ and ML strains, respectively. It took place mainly in spike and ORF1a genes. For BQ recombinant strains, the docking analysis showed that the spike interacted strongly with ACE2 and weakly with nAbs. The mutations S373P, S375F and T376A constitute a residue network that enhances the RBD interaction with ACE2. Thirteen mutations in RBD (S373P, S375F, T376A, D405N, R408S, K417N, N440K, S477N, P494S, Q498R, N501Y, and Y505H) and NTD (Y240H) seem to be implicated in immune evasion of recombinants by altering spike interaction with nAbs. In conclusion, this "in silico" study demonstrated that the recombination mechanism is frequent among Omicron BQ and ML variants. It highlights new key mutations, that potentially implicated in enhancement of spike binding to ACE2 (F376A) and escape from nAbs (RBD: F376A, D405N, R408S, N440K, S477N, P494S, and Y505H; NTD: Y240H). Our findings present considerable insights for the elaboration of effective prophylaxis and therapeutic strategies against future SARS-CoV-2 waves.

Leonurus japonicus Houtt. modulates neuronal apoptosis in intracerebral hemorrhage: Insights from network pharmacology and molecular docking

ETHNOPHARMACOLOGICAL RELEVANCE

Leonurus japonicus Houtt. (Labiatae), commonly known as Chinese motherwort, is a herbaceous flowering plant that is native to Asia. It is widely acknowledged in traditional medicine for its diuretic, hypoglycemic, antiepileptic properties and neuroprotection. Currently, Leonurus japonicus (Leo) is included in the Pharmacopoeia of the People's Republic of China. Traditional Chinese Medicine (TCM) recognizes Leo for its myriad pharmacological attributes, but its efficacy against ICH-induced neuronal apoptosis is unclear.

AIMS OF THE STUDY

This study aimed to identify the potential targets and regulatory mechanisms of Leo in alleviating neuronal apoptosis after ICH.

MATERIALS AND METHODS

The study employed network pharmacology, UPLC-Q-TOF-MS technique, molecular docking, pharmacodynamic studies, western blotting, and immunofluorescence techniques to explore its potential mechanisms.

RESULTS

Leo was found to assist hematoma absorption, thus improving the neurological outlook in an ICH mouse model. Importantly, molecular docking highlighted JAK as Leo's potential therapeutic target in ICH scenarios. Further experimental evidence demonstrated that Leo adjusts JAK1 and STAT1 phosphorylation, curbing Bax while augmenting Bcl-2 expression.

CONCLUSION

Leo showcases potential in mitigating neuronal apoptosis post-ICH, predominantly via the JAK/STAT mechanism.

Unveiling the characteristics of D4 and R4 aptamers for their future use in prostate cancer clinical practice

The DNA and RNA aptamers D4 and R4, respectively, emerged from the modification of PC-3 cell-binding aptamer A4. Our objective was to characterize the aptamers in silico and in vitro and begin to identify their target molecules. We represented their structures using computational algorithms; evaluated their binding to several prostate cell lines and their effects on the viability and migration of these cells; and determined their dissociation constant by flow cytometry. We analyzed circulating prostate tumor cells from patients using D4, R4, anti-CD133 and anti-CD44. Finally, the target proteins of both aptamers were precipitated and identified by mass spectrometry to simulate their in silico docking. The aptamers presented similar structures and bound to prostate tumor cells without modifying the cellular parameters studied, but with different affinities. The ligand cells for both aptamers were CD44+, indicating that they could identify cells in the mesenchymal stage of the metastatic process. The possible target proteins NXPE1, ADAM30, and MUC6 need to be further studied to better understand their interaction with the aptamers. These results support the development of new assays to determine the clinical applications of D4 and R4 aptamers in prostate cancer.

Computational Assessment of the Phytochemicals of Panax ginseng C.A. Meyer Against Dopamine Receptor D1 for Early Huntington's Disease Prophylactics

A herb, Panax ginseng C.A. Meyer has been used traditionally for the treatment of various diseases. In this work, its chemical components have been explored by computational methods for the possibility of therapeutic potential against early Huntington's disease. The molecular docking calculations against dopamine receptor D1 (PDB ID: 7X2F) involved in pathogenesis of early Huntington's disease gave the binding affinities (kcal/mol) of schizandrin (-10.530), ergosterol (-10.124), protopanaxadiol (-9.650), panaxydol (-9.399), diphenhydramine (-9.358), and panasenoside (-9.358). The values for native ligand (-7.748) and some selected drugs, Nefazodone (-9.880), Risperidone (-9.752), and Haloperidol (-9.712) were higher revealing weaker interactions. The stability assessment of top protein-ligand adducts in terms of various geometrical and thermodynamical parameters extracted from 200 ns molecular dynamics simulations pointed to schizandrin, protopanaxadiol, and panasenoside as hit molecules. The minimal translational and rotational motion of the docked ligands at orthosteric pocket of the receptor at near physiological conditions hinted at the probability of it restricting or inhibiting over-activation of DRD1. The sustained thermodynamic spontaneity of complex formation reaction augmented the inferences derived from spatial results. The phytochemicals from Panax ginseng could be used in the prophylactics of early Huntington's disease and recommendation is made for further evaluation by experimental work.

Analysis of the anti-Alzheimer potential of bioactive compounds from Citrus hystrix DC. peel, leaf, and essential oil by network pharmacology

Alzheimer's disease (AD) is the most known neurodegenerative disease, and its prevalence is predicted to increase significantly. Discovering novel drugs and treatments for AD is urgently needed. Drugs from natural products have been preferred lately due to their high potential and low toxicity. Citrus hystrix DC. (kaffir lime; KL) is one such herbal plant that is found abundantly in Southeast Asia with many biological activities. In this study, the potential of bioactive compounds from KL peel, leaf, and essential oil as anti-AD agents was explored using network pharmacology. First, the compounds were identified with KNApSAcK database and related literature. Subsequently, the targets of each corresponding compound were determined with SEA Search Server and Swiss Target Prediction, while the proteins associated with AD were identified using OMIM, GenCLiP3, and DisGeNET. Furthermore, a protein-protein interaction network and a compound-target interaction network were constructed to identify the most crucial proteins and compounds in the network by employing Cytoscape v3.9.1. The study continued with pathway enrichment analysis using STRING v1.7.1, molecular docking with PyRx and SwissDock, and molecular dynamics simulation with YASARA for further confirmation. Our results showed that almost all the secondary metabolites of KL targeted AD-associated genes, with oxypeucedanin and citrusoside A showing the highest anti-AD potential and targeting essential genes, EGFR and MAPK14, respectively. These targets were associated with inflammatory and oxidative stress pathways, indicating the potential mechanism of KL in attenuating AD clinical manifestation.

Optimization, gene cloning, expression, and molecular docking insights for enhanced cellulase enzyme production by Bacillus amyloliquefaciens strain elh1

BACKGROUND

In this study, we isolated a cellulase-producing bacterium, Bacillus amyloliquefaciens strain elh, from rice peel. We employed two optimization methods to enhance the yield of cellulase. Firstly, we utilized a one-variable-at-a-time (OVAT) approach to evaluate the impact of individual physical and chemical parameters. Subsequently, we employed response surface methodology (RSM) to investigate the interactions among these factors. We heterologously expressed the cellulase encoding gene using a cloning vectorin E. coli DH5α. Moreover, we conducted in silico molecular docking analysis to analyze the interaction between cellulase and carboxymethyl cellulose as a substrate.

RESULTS

The bacterial isolate eh1 exhibited an initial cellulase activity of 0.141 ± 0.077 U/ml when cultured in a specific medium, namely Basic Liquid Media (BLM), with rice peel as a substrate. This strain was identified as Bacillus amyloliquefaciens strain elh1 through 16S rRNA sequencing, assigned the accession number OR920278 in GenBank. The optimal incubation time was found to be 72 h of fermentation. Urea was identified as the most suitable nitrogen source, and dextrose as the optimal sugar, resulting in a production increase to 5.04 ± 0.120 U/ml. The peak activity of cellulase reached 14.04 ± 0.42 U/ml utilizing statistical optimization using Response Surface Methodology (RSM). This process comprised an initial screening utilizing the Plackett-Burman design and further refinement employing the BOX -Behnken Design. The gene responsible for cellulase production, egl, was effectively cloned and expressed in E. coli DH5α. The transformed cells exhibited a cellulase activity of 22.3 ± 0.24 U/ml. The egl gene sequence was deposited in GenBank with the accession number PP194445. In silico molecular docking revealed that the two hydroxyl groups of carboxymethyl cellulose bind to the residues of Glu169 inside the binding pocket of the CMCase. This interaction forms two hydrogen bonds, with an affinity score of -5.71.

CONCLUSIONS

Optimization of cultural conditions significantly enhances the yield of cellulase enzyme when compared to unoptimized culturing conditions. Additionally, heterologous expression of egl gene showed that the recombinant form of the cellulase is active and that a valid expression system can contribute to a better yield of the enzyme.

Phytochemicals and enzymes inhibitory potentials of leaves and rootbarks of Sarcocephallus latifolius (smith): In vitro and in silico investigations

In this study, bioactive compounds were isolated and characterized from the leaves and root-barks extracts of S.latifolius, with subsequent in vitro experimental investigations for antihyperglycemic potentials on α-amylase and α-glucosidase enzymes. Thirteen bioactive compounds were identified, including 10-Hydroxystrictosamide (2) and Quinovic acid-3-O-α-L-rhamnosyl-28-O-β-d-glucopyranosyl ester (8), using chromatographic, nuclear magnetic resonance spectroscopy (NMR), and mass spectrometry (MS) techniques. Experimental assays revealed that compounds 1-4 exhibited potent inhibition of α-amylase and α-glucosidase, with compound 2 demonstrating the most potent α-amylase inhibition (IC50 value of 0.52 ± 0.003 μg/mL). Compound 8 showed a lower IC50 value (0.098 ± 0.016 μg/mL) against α-glucosidase compared to compound 2 and acarbose. Synergistic effects among the compounds could enhance their inhibitory actions on the enzymes, positioning them as potential anti-hyperglycemia agents. Compound 2 displayed the highest binding affinity (-7.970 kcal/mol) when docked against α-amylase (PDB ID: 2QV4), comparable to acarbose (-8.515 kcal/mol). It also ranked among the top ligands against α-glucosidase (PDB ID 3TOP), although compound 13 and acarbose had higher docking scores. All compounds exhibited ideal ADMET properties, suggesting good bioavailability and low toxicity. In conclusion, the isolated compounds exhibit promising antihyperglycemic potential and favourable safety profiles for further exploration.

Building polarization into protein-inhibitor binding dynamics in rational drug design for rheumatoid arthritis

Standard force field-based simulations to accomplish structure-based evaluations of lead molecules is a powerful tool. Combining protein fragmentation into tractable sub-systems with continuum solvation method is envisaged to enable quantum mechanics-based electronic structure calculations of macromolecules in their realistic environment. This along with incorporation of many-body polarization effect in molecular dynamics simulations may augment an accurate description of electrostatics of protein-inhibitor systems for effective drug design. Rheumatoid arthritis (RA) is a complex autoimmune disorder plagued by the ceiling effect of current targeted therapies, encouraging identification of new druggable targets and corresponding drug design to tackle the refractory form of disease. In this study, polarization-inclusive force field approach has been used to model protein solvation and ligand binding for 'Mitogen-activated protein kinase' (MAP3K8), a regulatory node of notable pharmacological relevance in RA synovial biology. For MAP3K8 inhibitors belonging to different scaffold series, the calculations illustrated differential electrostatic contribution to their relative binding affinities and successfully explained examples from available structure-activity relationship studies. Results from this study exemplified i) the advantage of this approach in reliably ranking inhibitors having close nanomolar range activities for the same target; and ii) its prospective application in lead molecule identification aiding drug discovery efforts in RA.Communicated by Ramaswamy H. Sarma.

1,4-Diol Hq (TBHQ) vs 1,4-dithiol (TBDT); simulation of safe antioxidant with a lower carcinogenic activity

OBJECTIVES

tert-Butylhydroquinone (TBHQ) is an antioxidant and preservative used in unsaturated vegetable oils and processed foods. However, when consumed in higher doses daily, it may pose a threat to public health by potentially increasing the risk of cancer, as it has an affinity with both the aryl hydrocarbon receptor (AhR) and the estrogen receptor alpha (ERα).

METHODS

This study aimed to examine the impact of substituting the 1,4-diol of TBHQ with 1,4-dithiol, referred to as TBDT, on the carcinogenic and antioxidant systems using computational methods. The binding affinity of TBHQ and TBDT to the two carcinogenic receptors, AhR and ERα, as well as to the antioxidant receptor Keap1 alone and in connection with Nrf2 (Nrf2-Keap1) was investigated through docking analysis.

RESULTS

The results indicated a decrease in TBDT's binding strength to ERα and AhR when assessed using Molegro Virtual Docker (P-value: 0.0001 and 0.00001, respectively), AutoDock Vina (P-value: 0.0001 and 0.0001), and the online server Fast DRH (P-value: 0.0001 and 0.0001). However, TBDT's binding affinity to Keap1 was predicted to be significantly stronger than TBHQ's by both MVD and AutoDock Vina (P-value: 0.0001 and 0.04), while its binding to Nrf2-Keap1 assessed to be stronger only by MVD (P-value: 0.0001).

CONCLUSION

These findings suggest that TBDT not only exhibits higher antioxidant activity as a better ligand for the antioxidant system but also shows lower affinity with the AhR and ERα receptors. Therefore, TBDT can be considered a safer compound than TBHQ.

Identification of potential biogenic chalcones against antibiotic resistant efflux pump (AcrB) via computational study

The cases of bacterial multidrug resistance are increasing every year and becoming a serious concern for human health. Multidrug efflux pumps are key players in the formation of antibiotic resistance, which transfer out a broad spectrum of drugs from the cell and convey resistance to the host. Efflux pumps have significantly reduced the efficacy of the previously available antibiotic armory, thereby increasing the frequency of therapeutic failures. In gram-negative bacteria, the AcrAB-TolC efflux pump is the principal transporter of the substrate and plays a major role in the formation of antibiotic resistance. In the current work, advanced computer-aided drug discovery approaches were utilized to find hit molecules from the library of biogenic chalcones against the bacterial AcrB efflux pump. The results of the performed computational studies via molecular docking, drug-likeness prediction, pharmacokinetic profiling, pharmacophore mapping, density functional theory, and molecular dynamics simulation study provided ZINC000004695648, ZINC000014762506, ZINC000014762510, ZINC000095099506, and ZINC000085510993 as stable hit molecules against the AcrB efflux pumps. Identified hits could successfully act against AcrB efflux pumps after optimization as lead molecules.Communicated by Ramaswamy H. Sarma.

cFLIP - An interacting partner and a novel substrate for pro-apoptotic serine protease HtrA2

Background

HtrA2, a pro-apoptotic protease, plays a crucial role in apoptosis by cleaving inhibitory and anti-apoptotic proteins by translocating from mitochondria to the cytosol. Prior studies in ischemic cells have indicated that cytosolic HtrA2 triggers cFLIP degradation, plausibly through direct interaction. In this study, we have characterized the cFLIP protein, validated its interaction with HtrA2, and demonstrated that cFLIP is also a substrate of HtrA2.

Methods

We have identified the probable cleavage sites of cFLIP through gel-based assays and mass spectrometric analysis of the cleaved fragments.

Results

Our findings shed light on a key protein-protein interaction involving pro-apoptotic HtrA2, confirming cFLIP as its interacting partner and substrate.

Conclusion

Understanding the nuances of HtrA2's interaction with cFLIP (a decoy protein of the initiator procaspase-8 in the extrinsic apoptotic pathway) and deciphering the cFLIP's mode of cleavage, would provide an excellent alternative to modulate the pathway for therapeutic benefits toward diseases like ischemia and cancer.

Exploring the mechanism of action of spirooxindoles as a class of CDK2 inhibitors: a structure-based computational approach

Cyclin-dependent kinase 2 (CDK2) regulates cell cycle checkpoints in the synthesis and mitosis phases and plays a pivotal role in cancerous cell proliferation. The activation of CDK2, influenced by various protein signaling pathways, initiates the phosphorylation process. Due to its crucial role in carcinogenesis, CDK2 is a druggable hotspot target to suppress cancer cell proliferation. In this context, several studies have identified spirooxindoles as an effective class of CDK2 inhibitors. In the present study, three spirooxindoles (SOI1, SOI2, and SOI3) were studied to understand their inhibitory mechanism against CDK2 through a structure-based approach. Molecular docking and molecular dynamics (MD) simulations were performed to explore their interactions with CDK2 at the molecular level. The calculated binding free energy for the spirooxindole-based CDK2 inhibitors aligned well with experimental results regarding CDK2 inhibition. Energy decomposition (ED) analysis identified key binding residues, including I10, G11, T14, R36, F82, K89, L134, P155, T158, Y159, and T160, in the CDK2 active site and T-loop phosphorylation. Molecular mechanics (MM) energy was identified as the primary contributor to stabilizing inhibitor binding in the CDK2 protein structure. Furthermore, the analysis of binding affinity revealed that the inhibitor SOI1 binds more strongly to CDK2 compared to the other inhibitors under investigation. It demonstrated a robust interaction with the crucial residue T160 in the T-loop phosphorylation site, responsible for kinase activation. These insights into the inhibitory mechanism are anticipated to contribute to the development of potential CDK2 inhibitors using the spirooxindole scaffold.

In silico screening of a series of 1,6-disubstituted 1H-pyrazolo[3,4-d]pyrimidines as potential selective inhibitors of the Janus kinase 3

Rheumatoid arthritis is a common chronic disabling inflammatory disease that is characterized by inflammation of the synovial membrane and leads to discomfort. In the current study, twenty-seven 1,6-disubstituted 1H-pyrazolo[3,4-d]pyrimidines were tested as potential selective inhibitors of the tyrosine-protein kinase JAK3 using a number of molecular modeling methods. The activity of the screened derivatives was statistically quantified using multiple linear regression and artificial neural networks. To assess the quality, robustness, and predictability of the generated models, the leave-one-out cross-validation method was applied with favorable results (Q2 = 0.75) and Y-randomization. In addition, the evaluation of the predictive ability of the established model was confirmed by means of an external validation using a composite test set and an applicability domain approach. The covalent docking indicated that the tested 1H-pyrazolo[3,4-d]pyrimidines containing the acrylic aldehyde moiety had irreversible interaction with the residue Cys909 in the active sites of the tyrosine-protein kinase JAK3 by Michael addition. The molecular dynamics for three selected derivatives (compounds 9, 12, and 18) were used to verify the covalent docking by determining the stability of hydrogen bonding interactions with active sites, which are needed to stop tyrosine-protein kinase JAK3. The results obtained showed that the tested compounds containing acrylic aldehyde moiety had favorable binding free energies, indicating a strong affinity for the JAK3 enzyme. Overall, this current study suggests that the tested compounds containing the acrylic aldehyde moiety have the potential to act as anti-JAK3 inhibitors. They could be explored further to be used as treatment options for rheumatoid arthritis.Communicated by Ramaswamy H. Sarma.

Site-specific mutagenesis screening in KRASG12D mutant library to uncover resistance mechanisms to KRASG12D inhibitors

KRAS plays a crucial role in regulating cell survival and proliferation and is one of the most commonly mutated oncogenes in human cancers. The novel KRASG12D inhibitor, MRTX1133, demonstrates promising antitumor efficacy in vitro and in vivo. However, the development of acquired resistance in treated patients presents a considerable challenge to sustained therapeutic effectiveness. In response to this challenge, we conducted site-specific mutagenesis screening to identify potential secondary mutations that could induce resistance to MRTX1133. We screened a range of KRASG12D variants harboring potential secondary mutations, and 44 representative variants were selected for in-depth validation of the pooled screening outcomes. We identified eight variants (G12D with V9E, V9W, V9Q, G13P, T58Y, R68G, Y96W, and Q99L) that exhibited substantial resistance, with V9W showing notable resistance, and downstream signaling analyses and structural modeling were conducted. We observed that secondary mutations in KRASG12D can lead to acquired resistance to MRTX1133 and BI-2865, a novel pan-KRAS inhibitor, in human cancer cell lines. This evidence is critical for devising new strategies to counteract resistance mechanisms and, ultimately, enhance treatment outcomes in patients with KRASG12D-mutant cancers.

Carbohydrate polymer degradation derivatives as possible natural mannanase inhibitors

The role of mannanases is diverse and they are used in many industrial applications, in animal feed, in the food industry and in healthcare. They are also applied in biomass processing, because they play an important role in the breakdown of hemicellulose. Among the mannanase inhibitors, heavy metal ions and general enzyme inhibitors are mainly mentioned. Unfortunately, almost no data are available on carbohydrate-based natural inhibitors of mannanases. According to the literature, carbohydrates do not play an important role in the inhibition of mannanases, so neither do oligosaccharides. This is in contrast to the action and inhibition of other O-glycosyl hydrolases. My hypothesis is that mannanases, like other polysaccharide-degrading enzymes, work in the same way and can be inhibited by oligosaccharides. Evidence from docking and modeling results supports and makes probable the hypothesis that oligosaccharides can inhibit the activity of mannanases, similar to the inhibition of other O-glycosyl hydrolases. Among natural carbohydrate oligomers, several potential mannanase inhibitors have been identified and characterized. In addition to expensive research, it is very important to use research based on cheaper modeling to explore the processes. The results obtained are novel and forward-looking, enabling in-depth and targeted research to be carried out.

Unraveling DEHP influence on hemoglobin S polymerization in sickle cell disease: Ex vivo, in vitro and in silico analysis

Sickle cell disease (SCD) is a hereditary hemoglobinopathy, caused by a mutation at position 6 of the β-globin chain and patients are frequently exposed to several blood transfusions in order to maintain physiological function. Transfusion blood bags are composed of PVC and phthalates (as DEHP) are often introduced to the material in order to confer malleability. In this sense, DEHP can easily elute to the blood and cause harmful effects. This study aimed to unravel DEHP effect on SCD patient's hemoglobin function. We found that HbS polymerization using whole erythrocytes is decreased by DEHP in ex vivo experiments and this effect might be mediated by the DEHP-VAL6 interaction, evaluated in silico. Isolated HbS exhibited less polymerization at low DEHP concentrations and increased polymerization rate at higher concentration. When analyzing the propensity to aggregate, HbS is more inclined to aggregate when compared to HbA due to the residue 6 mutation. Circular dichroism showed characteristic hemoglobin peaks for oxygenated HbS that are lost when oxygen is sequestered, and DEHP at higher concentration mildly recovers a peak close to the second hemoglobin one. Finally, by transmission electron microscopy we demonstrated that high DEHP concentration increased polymer formation with a more organized structure. These findings show for the first-time the beneficial effect of low-dose DEHP on HbS polymerization.

Role of CSF1R 550th-tryptophan in kusunokinin and CSF1R inhibitor binding and ligand-induced structural effect

Binding affinity is an important factor in drug design to improve drug-target selectivity and specificity. In this study, in silico techniques based on molecular docking followed by molecular dynamics (MD) simulations were utilized to identify the key residue(s) for CSF1R binding affinity among 14 pan-tyrosine kinase inhibitors and 15 CSF1R-specific inhibitors. We found tryptophan at position 550 (W550) on the CSF1R binding site interacted with the inhibitors' aromatic ring in a π-π way that made the ligands better at binding. Upon W550-Alanine substitution (W550A), the binding affinity of trans-(-)-kusunokinin and imatinib to CSF1R was significantly decreased. However, in terms of structural features, W550 did not significantly affect overall CSF1R structure, but provided destabilizing effect upon mutation. The W550A also did not either cause ligand to change its binding site or conformational changes due to ligand binding. As a result of our findings, the π-π interaction with W550's aromatic ring could be still the choice for increasing binding affinity to CSF1R. Nevertheless, our study showed that the increasing binding to W550 of the design ligand may not ensure CSF1R specificity and inhibition since W550-ligand bound state did not induce significantly conformational change into inactive state.

G. G, K. P. C, Natarajamurthy S, K. N, Pradeep S, Shivamallu C, Elossaily GM, Achar RR, Silina E, Stupin V, Manturova N, A. Shati A, Y. Alfaifi M, I. Elbehairi SE, Kestur Nagaraj A, Mahadevamurthy M, Kollur SP. Computational exploration of Picrasma quassioides compounds as CviR-mediated quorum sensing inhibitors against Chromobacterium violaceum

Chromobacterium violaceum an opportunistic human pathogenic bacterium, exhibits resistance to conventional antibiotics by exploiting its quorum sensing mechanism to regulate virulence factor expression. In light of this, disrupting the quorum sensing mechanism presents a promising avenue for treating infections caused by this pathogen. The study focused on using the cytoplasmic quorum sensing receptor CviR from C. violaceum as a model target to identify novel quorum sensing inhibitors from P. quassioides through in silico computational approaches. Molecular docking analyses unveiled that several phytochemicals derived from Picrasma quassioides exhibit the potential to inhibit quorum sensing by binding to CviR protein. Notably, the compounds such as Quassidine I (- 8.8 kcal/mol), Quassidine J (- 8.8 kcal/mol), Kumudine B (- 9.1 kcal/mol) and Picrasamide A (- 8.9 kcal/mol) exhibited high docking scores, indicating strong binding affinity to the CviR protein. The native ligand C6-HSL (N-hexanoyl-L-homoserine lactone) as a positive control/co-crystal inhibitor also demonstrated a significant binding energy of-7.7 kcal/mol. The molecular dynamics simulation for 200 ns showed the thermodynamic stability and binding affinity refinement of the top-ranked CviR inhibitor (Kumudine B) with its stable binding and minor fluctuations compared to positive control (C6-HSL). Pharmacokinetic predictions indicated that Kumudine B possesses favourable drug-like properties, which suggest its potential as a drug candidate. The study highlight Kumudine B as a potential agent for inhibiting the CviR protein in C. violaceum. The comprehensive evaluation of Kumudine B provides valuable insights into its pharmacological profiles, facilitating its assessment for diverse therapeutic applications and guiding future research activities, particularly as antibacterial agents for clinical drug development.

Novel C3-Methylene-Bridged Indole Derivatives with and without Substituents at N1: The Influence of Substituents on Their Hemolytic, Cytoprotective, and Antimicrobial Activity

Alkaloids are natural compounds useful as scaffolds for discovering new bioactive molecules. This study utilized alkaloid gramine to synthesize two groups of C3-substituted indole derivatives, which were either functionalized at N1 or not. The compounds were characterized by spectroscopic methods. The protective effects of the new compounds against in vitro oxidative hemolysis induced by standard oxidant 2,2'-azobis(2-amidinopropane dihydro chloride (AAPH) on human erythrocytes as a cell model were investigated. Additionally, the compounds were screened for antimicrobial activity. The results indicated that most of the indole derivatives devoid of the N1 substitution exhibited strong cytoprotective properties. The docking studies supported the affinities of selected indole-based ligands as potential antioxidants. Furthermore, the derivatives obtained exhibited potent fungicidal properties. The structures of the eight derivatives possessing indole moiety bridged to the imidazole-, benzimidazole-, thiazole-, benzothiazole-, and 5-methylbenzothiazoline-2-thiones were determined by X-ray diffraction. The C=S bond lengths in the thioamide fragment pointed to the involvement of zwitterionic structures of varying contribution. The predominance of zwitterionic mesomers may explain the lack of cytoprotective properties, while steric effects, which limit multiple the hydrogen-bond acceptor properties of a thione sulfur, seem to be responsible for the high hemolytic activity.

Exploration of isoxazole analogs: Synthesis, COX inhibition, anticancer screening, 3D multicellular tumor spheroids, and molecular modeling

In this study, a new series of Isoxazole-carboxamide derivatives were synthesized and characterized via HRMS, 1H-, 13CAPT-NMR, and MicroED. The findings revealed that nearly all of the synthesized derivatives exhibited potent inhibitory activities against both COX enzymes, with IC50 values ranging from 4.1 nM to 3.87 μM. Specifically, MYM1 demonstrated the highest efficacy among the compounds tested against the COX-1, displaying an IC50 value of 4.1 nM. The results showed that 5 compounds possess high COX-2 isozyme inhibitory effects with IC50 value in range 0.24-1.30 μM with COX-2 selectivity indexes (2.51-6.13), among these compounds MYM4 has the lowest IC50 value against COX-2, with selectivity index around 4. Intriguingly, this compound displayed significant antiproliferative effects against CaCo-2, Hep3B, and HeLa cancer cell lines, with IC50 values of 10.22, 4.84, and 1.57 μM, respectively, which was nearly comparable to that of doxorubicin. Compound MYM4 showed low cytotoxic activities on normal cell lines LX-2 and Hek293t with IC50 values 20.01 and 216.97 μM respectively, with safer values than doxorubicin. Furthermore, compound MYM4 was able to induce the apoptosis, suppress the colonization of both HeLa and HepG2 cells. Additionally, the induction of Reactive oxygen species (ROS) production could be the mechanism underlying the apoptotic effect and the cytotoxic activity of the compound. In the 3D multicellular tumor spheroid model, results revealed that MYM4 compound hampered the spheroid formation capacity of Hep3B and HeLa cancer cells. Moreover, the molecular docking of MYM4 compound revealed a high affinity for the COX2 enzyme, with energy scores (S) -7.45 kcal/mol, which were comparable to celecoxib (S) -8.40 kcal/mol. Collectively, these findings position MYM4 as a promising pharmacological candidate as COX inhibitor and anticancer agent.

Searching possible SARS-CoV-2 main protease inhibitors in constituents from herbal medicines using in silico studies

The largest threat to civilization since the Second World War is the spread of the new coronavirus disease (COVID-19). Therefore, there is an urgent need for innovative therapeutic medicines to treat COVID-19. Reusing bio-actives is a workable and efficient strategy in the battle against new epidemics because the process of developing new drugs is time-consuming. This research aimed to identify which herbal remedies had the highest affinity for the receptor and assess a variety of them for potential targets to suppress the SARS-CoV-2 Mpro. The use of AutoDock Vina for structure-based virtual screening was done first due to the importance of protein interactions in the development of drugs. Molecular docking was used in the comparative study to assess 89 different chemicals from medicinal herbs. To anticipate their effectiveness against the primary protease of SARS-CoV-2, more analysis was done on the ADMET profile, drug-likeness, and Lipinski's rule of five. The next step involved three replicas of 100 ns-long molecular dynamics simulations on the potential candidates, which were preceded by calculations of the binding free energy of MM-GBSA. The outcomes showed that Achyrodimer A, Cinchonain Ib, Symphonone F, and Lupeol acetate all performed well and had the highest 6LU7 binding affinities. Using RMSD, RMSF, and protein-ligand interactions, the stability of the protein-ligand complex was assessed. The studies indicate that bioactive substances obtained from herbal medicines may function as a COVID-19 therapeutic agent, necessitating additional wet lab research to confirm their therapeutic potential, efficacy, and pharmacological capacity against the condition.Communicated by Ramaswamy H. Sarma.

PROTACable Is an Integrative Computational Pipeline of 3-D Modeling and Deep Learning To Automate the De Novo Design of PROTACs

Proteolysis-targeting chimeras (PROTACs) that engage two biological targets at once are a promising technology in degrading clinically relevant protein targets. Since factors that influence the biological activities of PROTACs are more complex than those of a small molecule drug, we explored a combination of computational chemistry and deep learning strategies to forecast PROTAC activity and enable automated design. A new method named PROTACable was developed for the de novo design of PROTACs, which includes a robust 3-D modeling workflow to model PROTAC ternary complexes using a library of E3 ligase and linker and an SE(3)-equivariant graph transformer network to predict the activity of newly designed PROTACs. PROTACable is available at https://github.com/giaguaro/PROTACable/.

Molecular interaction and MD-simulations: investigation of Sizofiran as a promising anti-cancer agent targeting eIF4E in colorectal cancer

CRC has a major global health impact due to high mortality rates. CRC shows high expression of eukaryotic translation initiation factor (eIF4E) protein, the rapid development of lung, bladder, colon, prostate, breast, head, and neck cancer is attributed to the dysregulation of eIF4E making an important target for treatment. Targeting eIF4E-mediated translation is a promising anti-cancer strategy. Many organic compounds that inhibit eIF4E are being studied clinically. The compound Sizofiran has emerged as a promising eIF4E inhibitor candidate, but its exact mechanism of action is unclear. In an effort to close this discrepancy by clarifying the mechanism of the interactions between phytochemical substances and eIF4E, molecular docking and dynamics studies were conducted. Molecular docking studies found Sizofiran (- 12.513 kcal/mol) has the most affinity eIF4E binding energy out of 93 phytochemicals, 5 current drugs, and 4 known inhibitors. This positions it as a top eIF4E inhibitor candidate. An alignment of eIF4E protein sequences from multiple pathogens revealed that the glutamate103 interacting residues are evolutionarily conserved across the different eIF4E proteins. Further insights from 100 ns of MD simulations supported Sizofiran having superior stability and eIF4E inhibition compared to reference compounds. Designed Sizofiran-related compounds showed better activity than the current drugs such as Camptosar, Sorafenib, Regorafenib, Doxorubicin, and Kenpaullone, indicating strong potential to suppress CRC progression by targeting eIF4E. This research aims to significantly aid development of improved eIF4E-targeting drugs for cancer treatment.

Graphical abstract

Showing the Graphical abstract of the complete study.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00206-3.

Deep Generative Models in De Novo Drug Molecule Generation

The discovery of new drugs has important implications for human health. Traditional methods for drug discovery rely on experiments to optimize the structure of lead molecules, which are time-consuming and high-cost. Recently, artificial intelligence has exhibited promising and efficient performance for drug-like molecule generation. In particular, deep generative models achieve great success in de novo generation of drug-like molecules with desired properties, showing massive potential for novel drug discovery. In this study, we review the recent progress of molecule generation using deep generative models, mainly focusing on molecule representations, public databases, data processing tools, and advanced artificial intelligence based molecule generation frameworks. In particular, we present a comprehensive comparison of state-of-the-art deep generative models for molecule generation and a summary of commonly used molecular design strategies. We identify research gaps and challenges of molecule generation such as the need for better databases, missing 3D information in molecular representation, and the lack of high-precision evaluation metrics. We suggest future directions for molecular generation and drug discovery.

In silico investigation of falcipain-2 inhibition by hybrid benzimidazole-thiosemicarbazone antiplasmodial agents: A molecular docking, molecular dynamics simulation, and kinetics study

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 n M) and favorable binding free energies (-30.32 to -17.17 kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18 kcal/mol in water and 28.42 kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

A new paradigm for applying deep learning to protein-ligand interaction prediction

Protein-ligand interaction prediction presents a significant challenge in drug design. Numerous machine learning and deep learning (DL) models have been developed to accurately identify docking poses of ligands and active compounds against specific targets. However, current models often suffer from inadequate accuracy or lack practical physical significance in their scoring systems. In this research paper, we introduce IGModel, a novel approach that utilizes the geometric information of protein-ligand complexes as input for predicting the root mean square deviation of docking poses and the binding strength (pKd, the negative value of the logarithm of binding affinity) within the same prediction framework. This ensures that the output scores carry intuitive meaning. We extensively evaluate the performance of IGModel on various docking power test sets, including the CASF-2016 benchmark, PDBbind-CrossDocked-Core and DISCO set, consistently achieving state-of-the-art accuracies. Furthermore, we assess IGModel's generalizability and robustness by evaluating it on unbiased test sets and sets containing target structures generated by AlphaFold2. The exceptional performance of IGModel on these sets demonstrates its efficacy. Additionally, we visualize the latent space of protein-ligand interactions encoded by IGModel and conduct interpretability analysis, providing valuable insights. This study presents a novel framework for DL-based prediction of protein-ligand interactions, contributing to the advancement of this field. The IGModel is available at GitHub repository https://github.com/zchwang/IGModel.

Novel 2-Aminobenzothiazole Derivatives: Docking, Synthesis, and Biological Evaluation as Anticancer Agents

Sixteen novel 2-aminobenzothiazole compounds with different amines or substituted piperazine moieties were designed, synthesized, and tested using various methods. Potential interactions were assessed by docking new compounds in the adenosine triphosphate (ATP) binding domain of the PI3Kγ enzyme (PDB code: 7JWE) by nucleophilic substitution or solvent-free/neat fusion for docked compound synthesis. Final 2-aminobenzothiazole compounds were characterized by direct probe gas chromatography-mass spectrometry (GC-MS), proton (1H-NMR), carbon-13 (13C-NMR), and attenuated total reflectance-infrared Fourier transform infrared (ATR FT-IR). The synthesized compounds were investigated for anticancer activities on lung cancer (A549) and breast cancer (MCF-7) cell lines. The compounds' PI3Kγ inhibition was evaluated at a 100 μM concentration. 4-Nitroaniline and piperazine-4-nitroaniline combination in OMS5 and OMS14 reduced lung and breast cancer cell line growth. IC50 values for OMS5 and OMS14, the strongest compounds, ranged from 22.13 to 61.03 μM. OMS1 and OMS2 inhibited PI3Kγ at the highest rates (47 and 48%, respectively) at a 100 μM concentration. Results show that the PI3Kγ enzyme suppression is not the main mechanism behind these OMS5 and OMS14 anticancer effects. CDK2, Akt, mTOR, and p42/44 MAPK are affected. EGF receptor suppression matters. AKT1, AKT3, CDK1/cyclin B, PDK1 direct, PIK3CA E542 K/PIK3R1 (p110 α/p85 α), PIK3CD/PIK3R1 (p110 δ/p85 α), and PKN inhibition were measured to evaluate the possible mechanism of compound OMS14. PIK3CD/PIK3R1 (p110 δ/p85 α) is the most, with 65% inhibition, suggesting a possible mechanism of anticancer properties. Furthermore, the NCI 60-cell line inhibition demonstrates promising broad anticancer inhibition against numerous cancer cell lines of OMS5 and OMS14, which could be good lead compounds for future development.

Gene Cloning, Heterologous Expression, and In Silico Analysis of Chitinase B from Serratia marcescens for Biocontrol of Spodoptera frugiperda Larvae Infesting Maize Crops

Spodoptera frugiperda, the fall armyworm (FAW), is a highly invasive polyphagous insect pest that is considered a source of severe economic losses to agricultural production. Currently, the majority of chemical insecticides pose tremendous threats to humans and animals besides insect resistance. Thus, there is an urgent need to develop new pest management strategies with more specificity, efficiency, and sustainability. Chitin-degrading enzymes, including chitinases, are promising agents which may contribute to FAW control. Chitinase-producing microorganisms are reported normally in bacteria and fungi. In the present study, Serratia marcescens was successfully isolated and identified from the larvae of Spodoptera frugiperda. The bacterial strain NRC408 displayed the highest chitinase enzyme activity of 250 units per milligram of protein. Subsequently, the chitinase gene was cloned and heterologously expressed in E. coli BL21 (DE3). Recombinant chitinase B was overproduced to 2.5-fold, driven by the T7 expression system. Recombinant chitinase B was evaluated for its efficacy as an insecticidal bioagent against S. frugiperda larvae, which induced significant alteration in subsequent developmental stages and conspicuous malformations. Additionally, our study highlights that in silico analyses of the anticipated protein encoded by the chitinase gene (ChiB) offered improved predictions for enzyme binding and catalytic activity. The effectiveness of (ChiB) against S. frugiperda was evaluated in laboratory and controlled field conditions. The results indicated significant mortality, disturbed development, different induced malformations, and a reduction in larval populations. Thus, the current study consequently recommends chitinase B for the first time to control FAW.

Preliminary insight on diarylpentanoids as potential antimalarials: In silico, in vitro pLDH and in vivo zebrafish toxicity assessment

Malaria remains a major public health problem worldwide, including in Southeast Asia. Chemotherapeutic agents such as chloroquine (CQ) are effective, but problems with drug resistance and toxicity have necessitated a continuous search for new effective antimalarial agents. Here we report on a virtual screening of ∼300 diarylpentanoids and derivatives, in search of potential Plasmodium falciparum lactate dehydrogenase (PfLDH) inhibitors with acceptable drug-like properties. Several molecules with binding affinities comparable to CQ were chosen for in vitro validation of antimalarial efficacy. Among them, MS33A, MS33C and MS34C are the most promising against CQ-sensitive (3D7) with EC50 values of 1.6, 2.5 and 3.1 μM, respectively. Meanwhile, MS87 (EC50 of 1.85 μM) shown the most active against the CQ-resistant Gombak A strain, and MS33A and MS33C the most effective P. knowlesi inhibitors (EC50 of 3.6 and 5.1 μM, respectively). The in vitro cytotoxicity of selected diarylpentanoids (MS33A, MS33C, MS34C and MS87) was tested on Vero mammalian cells to evaluate parasite selectivity (SI), showing moderate to low cytotoxicity (CC50 > 82 μM). In addition, MS87 exhibited a high SI and the lowest resistance index (RI), suggesting that MS87 may exert effective parasite inhibition with low resistance potential in the CQ-resistant P. falciparum strain. Furthermore, the in vivo toxicity of the molecules on early embryonic development, the cardiovascular system, heart rate, motor activity and apoptosis were assessed in a zebrafish animal model. The overall results indicate the preliminary potential of diarylpentanoids, which need further investigation for their development as new antimalarial agents.

Identifying potential monkeypox virus inhibitors: an in silico study targeting the A42R protein

Monkeypox (now Mpox), a zoonotic disease caused by the monkeypox virus (MPXV) is an emerging threat to global health. In the time span of only six months, from May to October 2022, the number of MPXV cases breached 80,000 and many of the outbreaks occurred in locations that had never previously reported MPXV. Currently there are no FDA-approved MPXV-specific vaccines or treatments, therefore, finding drugs to combat MPXV is of utmost importance. The A42R profilin-like protein of the MPXV is involved in cell development and motility making it a critical drug target. A42R protein is highly conserved across orthopoxviruses, thus A42R inhibitors may work for other family members. This study sought to identify potential A42R inhibitors for MPXV treatment using computational approaches. The energy minimized 3D structure of the A42R profilin-like protein (PDB ID: 4QWO) underwent virtual screening using a library of 36,366 compounds from Traditional Chinese Medicine (TCM), AfroDb, and PubChem databases as well as known inhibitor tecovirimat via AutoDock Vina. A total of seven compounds comprising PubChem CID: 11371962, ZINC000000899909, ZINC000001632866, ZINC000015151344, ZINC000013378519, ZINC000000086470, and ZINC000095486204, predicted to have favorable binding were shortlisted. Molecular docking suggested that all seven proposed compounds have higher binding affinities to A42R (-7.2 to -8.3 kcal/mol) than tecovirimat (-6.7 kcal/mol). This was corroborated by MM/PBSA calculations, with tecovirimat demonstrating the highest binding free energy of -68.694 kJ/mol (lowest binding affinity) compared to the seven shortlisted compounds that ranged from -73.252 to -97.140 kJ/mol. Furthermore, the 7 compounds in complex with A42R demonstrated higher stability than the A42R-tecovirimat complex when subjected to 100 ns molecular dynamics simulations. The protein-ligand interaction maps generated using LigPlot+ suggested that residues Met1, Glu3, Trp4, Ile7, Arg127, Val128, Thr131, and Asn133 are important for binding. These seven compounds were adequately profiled to be potential antivirals via PASS predictions and structural similarity searches. All seven potential lead compounds were scored Pa > Pi for antiviral activity while ZINC000001632866 and ZINC000015151344 were predicted as poxvirus inhibitors with Pa values of 0.315 and 0.215, and Pi values of 0.052 and 0.136, respectively. Further experimental validations of the identified lead compounds are required to corroborate their predicted activity. These seven identified compounds represent solid footing for development of antivirals against MPXV and other orthopoxviruses.

Evaluating the Therapeutic Potential of Curcumin and Synthetic Derivatives: A Computational Approach to Anti-Obesity Treatments

Natural compounds such as curcumin, a polyphenolic compound derived from the rhizome of turmeric, have gathered remarkable scientific interest due to their diverse metabolic benefits including anti-obesity potential. However, curcumin faces challenges stemming from its unfavorable pharmacokinetic profile. To address this issue, synthetic curcumin derivatives aimed at enhancing the biological efficacy of curcumin have previously been developed. In silico modelling techniques have gained significant recognition in screening synthetic compounds as drug candidates. Therefore, the primary objective of this study was to assess the pharmacokinetic and pharmacodynamic characteristics of three synthetic derivatives of curcumin. This evaluation was conducted in comparison to curcumin, with a specific emphasis on examining their impact on adipogenesis, inflammation, and lipid metabolism as potential therapeutic targets of obesity mechanisms. In this study, predictive toxicity screening confirmed the safety of curcumin, with the curcumin derivatives demonstrating a safe profile based on their LD50 values. The synthetic curcumin derivative 1A8 exhibited inactivity across all selected toxicity endpoints. Furthermore, these compounds were deemed viable candidate drugs as they adhered to Lipinski's rules and exhibited favorable metabolic profiles. Molecular docking studies revealed that both curcumin and its synthetic derivatives exhibited favorable binding scores, whilst molecular dynamic simulations showed stable binding with peroxisome proliferator-activated receptor gamma (PPARγ), csyclooxygenase-2 (COX2), and fatty acid synthase (FAS) proteins. The binding free energy calculations indicated that curcumin displayed potential as a strong regulator of PPARγ (-60.2 ± 0.4 kcal/mol) and FAS (-37.9 ± 0.3 kcal/mol), whereas 1A8 demonstrated robust binding affinity with COX2 (-64.9 ± 0.2 kcal/mol). In conclusion, the results from this study suggest that the three synthetic curcumin derivatives have similar molecular interactions to curcumin with selected biological targets. However, in vitro and in vivo experimental studies are recommended to validate these findings.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule's binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R2 = 0.8521, Q2LOO = 0.8438, and Q2LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: -8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies.Communicated by Ramaswamy H. Sarma.

Identification of natural xanthine oxidase inhibitors: Virtual screening, anti-xanthine oxidase activity, and interaction mechanism

Xanthine oxidase (XO) is a crucial target for hyperuricemia treatment(s). Naturally occurred XO inhibitors with minimal toxicity and high efficacy have attracted researchers' attention. With the goal of quickly identifying natural XO inhibitors, an integrated computational screening strategy was constructed by molecular docking and calculating the free energy of binding. Twenty-seven hits were achieved from a database containing 19,377 natural molecules. This includes fourteen known XO inhibitors and four firstly-reported inhibitors (isolicoflavonol, 5,7-dihydroxycoumarin, parvifolol D and clauszoline M, IC50 < 40 μM). Iolicoflavonol (hit 8, IC50 = 8.45 ± 0.68 μM) and 5,7-dihydroxycoumarin (hit 25, IC50 = 10.91 ± 0.71 μM) displayed the great potency as mixed-type inhibitors. Docking study and molecular dynamics simulation revealed that both hits could interact with XO's primarily active site residues ARG880, MOS1328, and ASN768 of XO. Fluorescence spectroscopy studies showed that hit 8 bound to the active cavity region of XO, causing changes in XO's conformation and hydrophobicity. Hits 8 and 25 exhibit favorable Absorption, Distribution, Metabolism, and Excretion (ADME) properties. Additionally, no cytotoxicity against human liver cells was observed at their median inhibition concentrations against XO. Therefore, the present study offers isolicoflavonol and 5,7-dihydroxycoumarin with the potential to be disease-modifying agents for hyperuricemia.

Structure-based virtual screening and molecular dynamics approaches to identify new inhibitors of Staphylococcus aureus sortase A

Staphylococcus aureus is a prevalent Gram-positive bacteria leading cause of a wide range of human pathologies. Moreover, antibiotic résistance of pathogenesis bacteria is one of the worldwide health problems. In Gram-positive bacteria, the enzyme of SrtA, is responsible for the anchoring of surface-exposed proteins to the cell wall peptidoglycan. Because of its critical role in Gram-positive bacterial pathogenesis, SrtA is an attractive target for anti-virulence during drug development. To date, some SrtA inhibitors have been discovered most of them being derived from flavonoid compounds, like Myricetin. In order to provide potential hit molecules against SrtA for clinical use, we obtained a total of 293 compounds by performing in silico shape-based screening of compound libraries against Myristin as a reference structure. Employing molecular docking and scoring functions, the top 3 compounds Apigenin, Efloxate, and Compound 8261032 were screened by comparing their docking scores with Myricetin. Furthermore, MD simulations and MM-PBSA binding energy calculation studies revealed that only Compound 8261032 strongly binds to the catalytic core of the SrtA enzyme than Myricetin, and stable behavior was consistently observed in the docking complex. Compound 8261032 showed a good number of hydrogen bonds with SrtA and higher MM-PBSA binding energy when compared to all three molecules. Also, it makes strength interactions with Arg139 and His62, which are critical for SrtA biological activity. This study showed that the development of this inhibitor could be a fundamental strategy against resistant bacteria, but further studies in vitro are needed to confirm this claim.Communicated by Ramaswamy H. Sarma.

Assessment of polystyrene nano plastics effect on human salivary α-amylase structural alteration: Insights from an in vitro and in silico study

The study found that the enzyme activity of human salivary α-amylase (α-AHS) was competitively inhibited by nanoplastic polystyrene (PS-NPs), with a half-inhibitory concentration (IC50) of 92 μg/mL, while the maximum reaction rate (Vmax) remained unchanged at 909 μg/mL•min. An increase in the concentration of PS-NPs led to a quenching of α-AHS fluorescence with a slight red shift, indicating a static mechanism. The binding constant (Ka) and quenching constant (Kq) were calculated to be 2.92 × 1011 M-1 and 1.078 × 1019 M-1• S-1 respectively, with a hill coefficient (n) close to one and an apparent binding equilibrium constant (KA) of 1.54 × 1011 M-1. Molecular docking results suggested that the interaction between α-AHS and PS-NPs involved π-anion interactions between the active site Asp197, Asp300 residues, and van der Waals force interactions affecting the Tyr, Trp, and other residues. Fourier transform infrared (FT-IR) and circular dichroism (CD) analyses revealed conformational changes in α-AHS, including a loss of secondary structure α-helix and β-sheet. The study concludes that the interaction between α-AHS and PS-NPs leads to structural and functional changes in α-AHS, potentially impacting human health. This research provides a foundation for further toxicological analysis of MPs/NPs in the human digestive system.

Recombinant production and characterization of a metal ion-independent Lysophospholipase from a hyperthermophilic archaeon Pyrococcus abyssi DSM25543

Genome sequence of Pyrococcus abyssi DSM25543 contains a coding sequence (PAB_RS01410) for α/β hydrolase (WP_010867387.1). Structural analysis revealed the presence of a consensus motif GXSXG and a highly conserved catalytic triad in the amino acid sequence of α/β hydrolase that were characteristic features of lysophospholipases. A putative lysophospholipase from P. abyssi with its potential applications in oil degumming and starch processing was heterologously produced in E. coli Rosetta (DE3) pLysS in soluble form followed by its purification and characterization. The recombinant enzyme was found to be active at temperature of 40-90 °C and pH 5.5-7.0. However, the enzyme exhibited its optimum activity at 65 °C and pH 6.5. None of the metal ions (Mn2+, Mg2+, Ni2+, Cu2+, Fe2+, Co2+, Zn2+ and Ca2+) being tested had stimulatory effect on lysophospholipase activity. Km and Vmax for hydrolysis of 4-nitrophenyl butyrate were calculated to be 1 ± 0.089 mM and 1637 ± 24.434 U/mg, respectively. It is the first report on the soluble production and characterization of recombinant lysophospholipase from P. abyssi which exhibits its lipolytic activity in the absence of divalent metal ions. Broad substrate specificity, activity and stability at elevated temperatures make recombinant lysophospholipase an ideal candidate for potential industrial applications.

How do physicochemical properties contribute to inhibitory activity of promising peptides against Zika Virus NS3 protease?

CONTEXT AND RESULTS

Flavivirus diseases' cycles, especially Dengue and Yellow Fever, can be observed all over Brazilian territory, representing a great health concern. Additionally, there are no drugs available in therapy. In this scenario, in silico methodologies were applied to obtain physicochemical properties, as well as to better understand the ligand-biological target interaction mode of 20 previously reported NS2B/NS3 protease inhibitors of Dengue virus. Since catalytic site of flavivirus hold similarities, such as the same catalytic triad (His51, Asp75 e Ser135), the ability of this series of molecules to fit in Zika NS3 domains can be achieved. We performed an exploratory data analysis, using statistical methodologies, such as PCA (Principal Component Analysis) and HCA (Hierarchical Component Analysis), to assist the comprehension of how physicochemical properties impact the interaction observed by the docking studies, as well as to build a correlation between the respective ranked characteristics. Based on these previous studies, peptides were selected for the dynamics simulations, which were useful to better understand the ligand-protein interactions. Information relating to, for instance, energy, ΔG, average number of hydrogen bonds and distance from Ser135 (one of the main amino acids in the catalytic pocket) were discussed. In this sense, peptides 15 (considering ΔG value and Hbond number), 7 (ΔG and energy) and 1, 6, 7 and 15 (the proximity to Ser135 throughout the dynamics simulation) were highlighted as promising. Those interesting results could contribute to future studies regarding Zika virus drug design, since this infection represents a great concern in neglected populations.

METHODS

The models were constructed in the ChemDraw software. The ligand parametrization was performed in the CHEM3D 17.0, UCSF Chimera. Docking simulations were carried out in the GOLD software, after the redocking validation. We used ASP as the function score. Additionally, for dynamics simulations we applied GROMACS software, exploring, mainly, free binding energy calculations. Exploratory analysis was carried out in Minitab 17.3.1 statistical software. Prior to the exploratory analysis, data of quantum chemical properties of the peptides were collected in Microsoft Excel spreadsheet and organized to obtain Hierarchical Cluster Analysis (HCA) and Principal Component Analysis (PCA).

Secondary metabolites of Trichoderma spp. as EGFR tyrosine kinase inhibitors: Evaluation of anticancer efficacy through computational approach

The present study explores the epidermal growth factor receptor (EGFR) tyrosine kinase inhibition efficacy of secondary metabolites in Trichoderma spp. through molecular docking, molecular dynamics (MD) simulation and MM-PBSA approach. The result of molecular docking confirmed that out of 200 metabolites screened, three metabolites such as Harzianelactone A, Pretrichodermamide G and Aspochalasin M, potentially bound with the active binding site of EGFR tyrosine kinase domain(PDB ID: 1M17) with a threshold docking score of ≤- 9.0 kcal/mol when compared with the standard EGFR inhibitor (Erlotinib). The MD simulation was run to investigate the potential for stable complex formation in EGFR tyrosine kinase domain-unbound/lead metabolite (Aspochalasin M)-bound/standard inhibitor (Erlotinib)-bound complex. The MD simulation analysis at 100 ns revealed that Aspochalasin M formed the stable complex with EGFR. Besides, the in silico predication of pharmacokinetic properties further confirmed that Aspochalasin M qualified the drug-likeness rules with no harmful side effects (viz., hERG toxicity, hepatotoxicity and skin sensitization), non-mutagenicity and favourable logBB value. Moreover, the BOILED-Egg model predicted that Aspochalasin M showed a higher gastrointestinal absorption with improved bioavailability when administered orally and removed from the central nervous system (CNS). The results of the computational studies concluded that Aspochalasin M possessed significant efficacy in binding EGFR's active sites compared to the known standard inhibitor (Erlotinib). Therefore, Aspochalasin M can be used as a possible anticancer drug candidate and further in vitro and in vivo experimental validation of Aspochalasin M of Trichoderma spp. are required to determine its anticancer potential.

Revealing Propolis Potential Activity on Inhibiting Estrogen Receptor and Heat Shock Protein 90 Overexpressed in Breast Cancer by Bioinformatics Approaches

Breast cancer is the most commonly diagnosed cancer globally, with the highest incidence of breast cancer occurring in Asian countries including Indonesia. Among the types of breast cancer, the estrogen receptor (ER)-positive subtype which is prominent with estrogen receptor alpha (ERα) and heat shock protein 90 (HSP90) overexpression genes becomes the most prevalent than the others, approximately 75% of all breast cancer cases. ERα and HSP90 play a role in breast cancer activities including breast tumor growth, invasion, and metastasis mechanism. Propolis, a natural bee product, has been explored for its anticancer activity. However, there is lack of studies that evaluated the potential inhibitor from propolis compounds to the ERα and HSP90 proteins. Therefore, this article focuses on examining the correlation between ERα and HSP90's role in breast cancer and investigating the potential of 93 unique propolis compositions in inhibiting these genes in breast cancer using in silico approaches. This study revealed the positive correlation between ERα and HSP90 genes in breast cancer disease development. Furthermore, we also found novel potential bioactive compounds of propolis against breast cancer through binding with ERα and HSP90; they were 3',4',7-trihydroxyisoflavone and baicalein-7-O-β-D glucopyranoside, respectively. Further research on these compounds is needed to elucidate deeper mechanisms and activity in the real biological system to develop new breast cancer drug treatments.