Is It Reliable to Use Common Molecular Docking Methods for Comparing the Binding Affinities of Enantiomer Pairs for Their Protein Target?

Astaxanthin targets IL-6 and alleviates the LPS-induced adverse inflammatory response of macrophages

Numerous natural compounds are recognized for their anti-inflammatory properties attributed to antioxidant effects and the modulation of key inflammatory factors. Among them, astaxanthin (AST), a potent carotenoid antioxidant, remains relatively underexplored regarding its anti-inflammatory mechanisms and specific molecular targets. In this study, human monocytic leukemia cell-derived macrophages (THP-1) were selected as experimental cells, and lipopolysaccharides (LPS) served as inflammatory stimuli. Upon LPS treatment, the oxidative stress was significantly increased, accompanied by remarkable cellular damage. Moreover, LPSs escalated the expression of inflammation-related molecules. Our results demonstrate that AST intervention could effectively alleviate LPS-induced oxidative stress, facilitate cellular repair, and significantly attenuate inflammation. Further exploration of the anti-inflammatory mechanism revealed AST could substantially inhibit NF-κB translocation and activation, and mitigate inflammatory factor production by hindering NF-κB through the antioxidant mechanism. We further confirmed that AST exhibited protective effects against cell damage and reduced the injury from inflammatory cytokines by activating p53 and inhibiting STAT3. In addition, utilizing network pharmacology and in silico calculations based on molecular docking, molecular dynamics simulation, we identified interleukin-6 (IL-6) as a prominent core target of AST anti-inflammation, which was further validated by the RNA interference experiment. This IL-6 binding capacity actually enabled AST to curb the positive feedback loop of inflammatory factors, averting the onset of possible inflammatory storms. Therefore, this study offers a new possibility for the application and development of astaxanthin as a popular dietary supplement of anti-inflammatory or immunomodulatory function.

Novel hydrazone compounds with broad-spectrum antiplasmodial activity and synergistic interactions with antimalarial drugs

The development of novel antiplasmodial compounds with broad-spectrum activity against different stages of Plasmodium parasites is crucial to prevent malaria disease and parasite transmission. This study evaluated the antiplasmodial activity of seven novel hydrazone compounds (referred to as CB compounds: CB-27, CB-41, CB-50, CB-53, CB-58, CB-59, and CB-61) against multiple stages of Plasmodium parasites. All CB compounds inhibited blood stage proliferation of drug-resistant or sensitive strains of Plasmodium falciparum in the low micromolar to nanomolar range. Interestingly, CB-41 exhibited prophylactic activity against hypnozoites and liver schizonts in Plasmodium cynomolgi, a primate model for Plasmodium vivax. Four CB compounds (CB-27, CB-41, CB-53, and CB-61) inhibited P. falciparum oocyst formation in mosquitoes, and five CB compounds (CB-27, CB-41, CB-53, CB-58, and CB-61) hindered the in vitro development of Plasmodium berghei ookinetes. The CB compounds did not inhibit the activation of P. berghei female and male gametocytes in vitro. Isobologram assays demonstrated synergistic interactions between CB-61 and the FDA-approved antimalarial drugs, clindamycin and halofantrine. Testing of six CB compounds showed no inhibition of Plasmodium glutathione S-transferase as a putative target and no cytotoxicity in HepG2 liver cells. CB compounds are promising candidates for further development as antimalarial drugs against multidrug-resistant parasites, which could also prevent malaria transmission.

In silico exploration of phenolics as modulators of penicillin binding protein (PBP) 2× of Streptococcus pneumoniae

Infections caused by multidrug-resistant Streptococcus pneumoniae remain the leading cause of pneumonia-related deaths in children < 5 years globally, and mutations in penicillin-binding protein (PBP) 2 × have been identified as the major cause of resistance in the organism to beta-lactams. Thus, the development of new modulators with enhanced binding of PBP2x is highly encouraged. In this study, phenolics, due to their reported antibacterial activities, were screened against the active site of PBP2x using structure-based pharmacophore and molecular docking techniques, and the ability of the top-hit phenolics to inhibit the active and allosteric sites of PBP2x was refined through 120 ns molecular dynamic simulation. Except for gallocatechin gallate and lysidicichin, respectively, at the active and allosteric sites of PBP2x, the top-hit phenolics had higher negative binding free energy (ΔGbind) than amoxicillin [active site (- 19.23 kcal/mol), allosteric site (- 33.75 kcal/mol)]. Although silicristin had the best broad-spectrum effects at the active (- 38.41 kcal/mol) and allosteric (- 50.54 kcal/mol) sites of PBP2x, the high thermodynamic entropy (4.90 Å) of the resulting complex might suggest the need for its possible structural refinement for enhanced potency. Interestingly, silicristin had a predicted synthetic feasibility score of < 5 and quantum calculations using the DFT B3LYP/6-31G+ (dp) revealed that silicristin is less stable and more reactive than amoxicillin. These findings point to the possible benefits of the top-hit phenolics, and most especially silicristin, in the direct and synergistic treatment of infections caused by S. pneumoniae. Accordingly, silicristin is currently the subject of further confirmatory in vitro research.

Compromise in Docking Power of Liganded Crystal Structures of Mpro SARS-CoV-2 Surpasses 90% Success Rate

Herein, we present the capacity of three different molecular docking programs (AutoDock, AutoDock Vina, and PLANTS) to identify and reproduce the binding modes of ligands present in 247 covalent and 169 noncovalent complex crystal structures of the severe acute respiratory syndrome coronavirus 2 main protease (Mpro). The compromise in docking power is evaluated with respect to their ability to generate poses similar to the crystal structure binding mode (heavy atoms' root-mean-square deviation < 2 Å) and their ability to recognize the native binding mode with an included compensation for the scoring function error. Noncovalently bound inhibitors are best modeled by AutoDock Vina (90.6% success rate in the active site), while the most relevant results for covalently bound inhibitors are produced by PLANTS (93.0%). AutoDock shows acceptable performance for both types of ligands, 81.1 and 76.4% for noncovalent and covalent complexes, respectively. All three programs manifest worse performance when reproducing surface-bound ligands. Comparison with other works illustrates the importance of crystal structure processing (12% of noncovalent and 26% of covalent ligands had to be manually corrected), proper sampling protocol settings, and inclusion of root-mean-square deviation (RMSD)/scoring function error compensations in crystal structure pose identification. Results are analyzed with respect to a clustering scheme of the noncovalently bound ligands and the chemical reaction type of the covalent ligand bound to the Cys145 residue. A comparison of screening power based on the docking scores of noncovalent ligands from the crystal structures with a "Directory of Useful Decoys, Enhanced" set of known decoys (6562 compounds) and ZINC15 in vivo subset (60,394 compounds) is provided. Ligand and protein input files are provided for future benchmarking purposes.

Waste to Medicine: Evidence from Computational Studies on the Modulatory Role of Corn Silk on the Therapeutic Targets Implicated in Type 2 Diabetes Mellitus

Type 2 diabetes mellitus (T2DM) is characterized by insulin resistance and/or defective insulin production in the human body. Although the antidiabetic action of corn silk (CS) is well-established, the understanding of the mechanism of action (MoA) behind this potential is lacking. Hence, this study aimed to elucidate the MoA in different samples (raw and three extracts: aqueous, hydro-ethanolic, and ethanolic) as a therapeutic agent for the management of T2DM using metabolomic profiling and computational techniques. Ultra-performance liquid chromatography-mass spectrometry (UP-LCMS), in silico techniques, and density functional theory were used for compound identification and to predict the MoA. A total of 110 out of the 128 identified secondary metabolites passed the Lipinski's rule of five. The Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analysis revealed the cAMP pathway as the hub signaling pathway, in which ADORA1, HCAR2, and GABBR1 were identified as the key target genes implicated in the pathway. Since gallicynoic acid (-48.74 kcal/mol), dodecanedioc acid (-34.53 kcal/mol), and tetradecanedioc acid (-36.80 kcal/mol) interacted well with ADORA1, HCAR2, and GABBR1, respectively, and are thermodynamically stable in their formed compatible complexes, according to the post-molecular dynamics simulation results, they are suggested as potential drug candidates for T2DM therapy via the maintenance of normal glucose homeostasis and pancreatic β-cell function.

Cheminformatics bioprospection of selected medicinal plants metabolites against trypsin cleaved VP4 (spike protein) of rotavirus A

Rotaviruses have continued to be the primary cause of acute dehydrating diarrhoea in children under five years of age despite the global introduction of four World Health Organization (WHO) prequalified oral vaccines in over 106 countries. Currently, no medication is approved by the Food and Drug Administration (FDA) specifically for treating rotavirus A-induced diarrhoea. Consequently, it is important to focus on developing prophylactic and curative therapeutics to combat rotaviral infections. For the first time, this study computationally screened and identified metabolites from Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum as potential novel inhibitors with broad-spectrum activity against VP5* and VP8* (spike protein) of rotavirus A (RVA). The initial top 20 metabolites identified through molecular docking were further filtered using drug-likeness and pharmacokinetics parameters. The molecular properties of the resulting top-ranked compounds were predicted by conducting density functional theory (DFT) calculations, while molecular dynamics (MD) simulation revealed their thermodynamic compatibility with a significant affinity towards VP8* than VP5*. Except for ellagic acid (-11.78 kcal/mol), the lead compounds had higher binding free energy than the reference standard (VP5* (-11.81 kcal/mol), VP8* (-14.12 kcal/mol)) with 2SG (-20.98 kcal/mol) and apigenin-4'-glucoside (-23.56 kcal/mol) having the highest affinity towards VP5* and VP8*, respectively. Of all the top-ranked compounds, better broad-spectrum affinities for both VP5* and VP8* than tizoxanide were observed in 2SG (VP5* (-20.98 kcal/mol), VP8* (-20.13 kcal/mol)) and sericetin (VP5* (-20.46 kcal/mol), VP8* (-18.31 kcal/mol)). While the identified leads could be regarded as potential modulators of the investigated therapeutic targets for effective management of rotaviral infection, additional in vitro and in vivo evaluation is strongly recommended, and efforts are on-going in this regard.Communicated by Ramaswamy H. Sarma.

Exploration of CviR-mediated quorum sensing inhibitors from Cladosporium spp. against Chromobacterium violaceum through computational studies

An opportunistic human pathogenic bacterium, Chromobacterium violaceum resists the potency of most antibiotics by exploiting the quorum sensing system within their community to control virulence factor expression. Therefore, blocking the quorum sensing mechanism could help to treat several infectious caused by this organism. The quorum sensing receptor (CviR) of C. violaceum was used as a model target in the current investigation to identify potentially novel quorum sensing inhibitors from Cladosporium spp. through in silico computational approaches. The molecular docking results confirmed the anti-quorum sensing potential of bioactive compounds from Cladosporium spp. through binding to CviR with varying docking scores between - 5.2 and - 9.5 kcal/mol. Relative to the positive control [Azithromycin (- 7.4 kcal/mol)], the top six metabolites of Cladosporium spp. had higher docking scores and were generally greater than - 8.5 kcal/mol. The thermodynamic stability and binding affinity refinement of top-ranked CviR inhibitors were further studied through a 160 ns molecular dynamic (MD) simulation. The Post-MD simulation analysis confirmed the top-ranked compounds' affinity, stability, and biomolecular interactions with CviR at 50 ns, 100 ns, and 160 ns with Coniochaetone K of the Cladosporium spp. having the highest binding free energy (- 30.87 kcal/mol) and best interactions (two consistent hydrogen bond contact) following the 160 ns simulation. The predicted pharmacokinetics properties of top selected compounds point to their drug likeliness, potentiating their chance as a possible drug candidate. Overall, the top-ranked compounds from Cladosporium spp., especially Coniochaetone K, could be identified as potential C. violaceum CviR inhibitors. The development of these compounds as broad-spectrum antibacterial medicines is thus possible in the future following the completion of further preclinical and clinical research.

Identification and validation of a small molecule targeting ROR1 for the treatment of triple negative breast cancer

Introduction: Breast cancer is the most common cancer in women, with roughly 10-15% of new cases classified as triple-negative breast cancer (TNBC). Traditional chemotherapies are often toxic to normal cells. Therefore, it is important to discover new anticancer compounds that target TNBC while causing minimal damage to normal cells. Receptor tyrosine kinase-like Orphan Receptor 1 (ROR1) is an oncofetal protein overexpressed in numerous human malignancies, including TNBC. This study investigated potential small molecules targeting ROR1. Methodology: Using AutoDock Vina and Glide, we screened 70,000 chemicals for our investigation. We obtained 10 representative compounds via consensus voting, deleting structural alerts, and clustering. After manual assessment, compounds 2 and 4 were chosen for MD simulation and cell viability experiment. Compound 4 showed promising results in the viability assay, which led us to move further with the apoptosis assay and immunoblotting. Results: Compound 4 (CID1261330) had docking scores of -6.635 and -10.8. It fits into the pocket and shows interactions with GLU64, ASP174, and PHE93. Its RMSD fluctuates around 0.20 nm and forms two stable H-bonds indicating compound 4 stability. It inhibits cell proliferation in MDA-MB-231, HCC1937, and HCC1395 cell lines, with IC50 values of approximately 2 μM to 10 μM, respectively. Compound 4 did not kill non-malignant epithelial breast cells MCF-10A (IC50 > 27 μM). These results were confirmed by the significant number of apoptotic cells in MDA-MB-231 cells (47.6%) but not in MCF-10A cells (7.3%). Immunoblot analysis provided additional support in the same direction. Discussion: These findings collectively suggest that compound 4 has the potential to effectively eliminate TNBC cells while causing minimal harm to normal breast cells. The promising outcomes of this study lay the groundwork for further testing of compound 4 in other malignancies characterized by ROR1 upregulation, serving as a proof-of-concept for its broader applicability.

Molecular Modeling Unveils the Effective Interaction of B-RAF Inhibitors with Rare B-RAF Insertion Variants

The Food and Drug Administration (FDA) has approved MAPK inhibitors as a treatment for melanoma patients carrying a mutation in codon V600 of the BRAF gene exclusively. However, BRAF mutations outside the V600 codon may occur in a small percentage of melanomas. Although these rare variants may cause B-RAF activation, their predictive response to B-RAF inhibitor treatments is still poorly understood. We exploited an integrated approach for mutation detection, tumor evolution tracking, and assessment of response to treatment in a metastatic melanoma patient carrying the rare p.T599dup B-RAF mutation. He was addressed to Dabrafenib/Trametinib targeted therapy, showing an initial dramatic response. In parallel, in-silico ligand-based homology modeling was set up and performed on this and an additional B-RAF rare variant (p.A598_T599insV) to unveil and justify the success of the B-RAF inhibitory activity of Dabrafenib, showing that it could adeptly bind both these variants in a similar manner to how it binds and inhibits the V600E mutant. These findings open up the possibility of broadening the spectrum of BRAF inhibitor-sensitive mutations beyond mutations at codon V600, suggesting that B-RAF V600 WT melanomas should undergo more specific investigations before ruling out the possibility of targeted therapy.

Computational insight into stability-enhanced systems of anthocyanin with protein/peptide

Anthocyanins, which belong to the flavonoid group, are commonly found in the organs of plants native to South and Central America. However, these pigments are unstable under conditions of varying pH, heat, etc., which limits their potential applications. One method for preserving the stability of anthocyanins is through encapsulation using proteins or peptides. Nevertheless, the complex and diverse structure of these molecules, as well as the limitation of experimental technologies, have hindered a comprehensive understanding of the encapsulation processes and the mechanisms by which stability is enhanced. To address these challenges, computational methods, such as molecular docking and molecular dynamics simulation have been used to study the binding affinity and dynamics of interactions between proteins/peptides and anthocyanins. This review summarizes the mechanisms of interaction between these systems, based on computational approaches, and highlights the role of proteins and peptides in the stability enhancement of anthocyanins. It also discusses the current limitations of these methods and suggests possible solutions.

Insight into the Potential Mechanisms of Endocrine Disruption by Dietary Phytoestrogens in the Context of the Etiopathogenesis of Endometriosis

Phytoestrogens (PEs) are estrogen-like nonsteroidal compounds derived from plants (e.g., nuts, seeds, fruits, and vegetables) and fungi that are structurally similar to 17β-estradiol. PEs bind to all types of estrogen receptors, including ERα and ERβ receptors, nuclear receptors, and a membrane-bound estrogen receptor known as the G protein-coupled estrogen receptor (GPER). As endocrine-disrupting chemicals (EDCs) with pro- or antiestrogenic properties, PEs can potentially disrupt the hormonal regulation of homeostasis, resulting in developmental and reproductive abnormalities. However, a lack of PEs in the diet does not result in the development of deficiency symptoms. To properly assess the benefits and risks associated with the use of a PE-rich diet, it is necessary to distinguish between endocrine disruption (endocrine-mediated adverse effects) and nonspecific effects on the endocrine system. Endometriosis is an estrogen-dependent disease of unknown etiopathogenesis, in which tissue similar to the lining of the uterus (the endometrium) grows outside of the uterus with subsequent complications being manifested as a result of local inflammatory reactions. Endometriosis affects 10-15% of women of reproductive age and is associated with chronic pelvic pain, dysmenorrhea, dyspareunia, and infertility. In this review, the endocrine-disruptive actions of PEs are reviewed in the context of endometriosis to determine whether a PE-rich diet has a positive or negative effect on the risk and course of endometriosis.

Metabolites profiling and cheminformatics bioprospection of selected medicinal plants against the main protease and RNA-dependent RNA polymerase of SARS-CoV-2

Despite the existence of some vaccines, SARS-CoV-2 (S-2) infections persist for various reasons relating to vaccine reluctance, rapid mutation rate, and an absence of specific treatments targeted to the infection. Due to their availability, low cost and low toxicity, research into potentially repurposing phytometabolites as therapeutic alternatives has gained attention. Therefore, this study explored the antiviral potential of metabolites of some medicinal plants [Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum (Sesame plant)] identified using liquid chromatography-mass spectrometry (LCMS) as possible inhibitory agents against the S-2 main protease (S-2 MP) and RNA-dependent RNA polymerase (RP) using computational approaches. Molecular docking was used to identify the compounds with the best affinities for the selected therapeutics targets. Afterwards, compounds with poor physicochemical characteristics, pharmacokinetics, and drug-likeness were screened out. The top-ranked compounds were further subjected to a 120-ns molecular dynamics (MD) simulation. Only quercetin 3-O-rhamnoside (-48.77 kcal/mol) had higher binding free energy than the reference standard (zafirlukast) (-44.99 kcal/mol) against S-2 MP. Conversely, all the top-ranked compounds (ellagic acid hexoside, spiraeoside, apigenin-4'-glucoside and chrysoeriol 7-glucuronide) except gnetin L (-24.24 kcal/mol) had higher binding free energy (-55.19 kcal/mol, -52.75 kcal/mol, -47.22 kcal/mol and -43.35 kcal/mol) respectively, against S-2 RP relative to the reference standard (-34.79 kcal/mol). The MD simulations study further revealed that the investigated inhibitors are thermodynamically stable and form structurally compatible complexes that impede the regular operation of the respective S-2 therapeutic targets. Although, these S-2 therapeutic candidates are promising, further in vitro and in vivo evaluation is required and highly recommended.Communicated by Ramaswamy H. Sarma.

Discovery of Selective Dopamine Receptor Ligands Derived from (-)-Stepholidine via C-3 Alkoxylation and C-3/C-9 Dialkoxylation

We evaluated C-3 alkoxylated and C-3/C-9 dialkoxylated (-)-stepholidine analogues to probe the tolerance at the C-3 and C-9 positions of the tetrahydroprotoberberine (THPB) template toward affinity for dopamine receptors. A C-9 ethoxyl substituent appears optimal for D1R affinity since high D1R affinities were observed for compounds that contain an ethyl group at C-9, with larger C-9 substituents tending to decrease D1R affinity. A number of novel ligands were identified, such as compounds 12a and 12b, with nanomolar affinities for D1R and no affinity for either D2R or D3R, with compound 12a being identified as a D1R antagonist for both G-protein- and β-arrestin-based signaling. Compound 23b was identified as the most potent and selective D3R ligand containing a THPB template to date and functions as an antagonist for both G-protein- and β-arrestin-based signaling. Molecular docking and molecular dynamics studies validated the D1R and D3R affinity and selectivity of 12a, 12b, and 23b.

Updates on drug designing approach through computational strategies: a review

The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.

Marine-Derived Bioactive Peptides Self-Assembled Multifunctional Materials: Antioxidant and Wound Healing

Peptide self-assembling materials have received significant attention from researchers in recent years, emerging as a popular field in biological, environmental, medical, and other new materials studies. In this study, we utilized controllable enzymatic hydrolysis technology (animal proteases) to obtain supramolecular peptide self-assembling materials (CAPs) from the Pacific oyster (Crassostrea gigas). We conducted physicochemical analyses to explore the pro-healing mechanisms of CAPs on skin wounds in both in vitro and in vivo experiments through a topical application. The results demonstrated that CAPs exhibit a pH-responsive behavior for self-assembly and consist of peptides ranging from 550 to 2300 Da in molecular weight, with peptide chain lengths of mainly 11-16 amino acids. In vitro experiments indicated that CAPs display a procoagulant effect, free radical scavenging activity, and promote the proliferation of HaCaTs (112.74% and 127.61%). Moreover, our in vivo experiments demonstrated that CAPs possess the ability to mitigate inflammation, boost fibroblast proliferation, and promote revascularization, which accelerates the epithelialization process. Consequently, a balanced collagen I/III ratio in the repaired tissue and the promotion of hair follicle regeneration were observed. With these remarkable findings, CAPs can be regarded as a natural and secure treatment option with high efficacy for skin wound healing. The potential of CAPs to be further developed for traceless skin wound healing is an exciting area for future research and development.

SVSBI: sequence-based virtual screening of biomolecular interactions

Virtual screening (VS) is a critical technique in understanding biomolecular interactions, particularly in drug design and discovery. However, the accuracy of current VS models heavily relies on three-dimensional (3D) structures obtained through molecular docking, which is often unreliable due to the low accuracy. To address this issue, we introduce a sequence-based virtual screening (SVS) as another generation of VS models that utilize advanced natural language processing (NLP) algorithms and optimized deep K-embedding strategies to encode biomolecular interactions without relying on 3D structure-based docking. We demonstrate that SVS outperforms state-of-the-art performance for four regression datasets involving protein-ligand binding, protein-protein, protein-nucleic acid binding, and ligand inhibition of protein-protein interactions and five classification datasets for protein-protein interactions in five biological species. SVS has the potential to transform current practices in drug discovery and protein engineering.

Cheminformatics-Based Study Identifies Potential Ebola VP40 Inhibitors

The Ebola virus (EBOV) is still highly infectious and causes severe hemorrhagic fevers in primates. However, there are no regulatorily approved drugs against the Ebola virus disease (EVD). The highly virulent and lethal nature of EVD highlights the need to develop therapeutic agents. Viral protein 40 kDa (VP40), the most abundantly expressed protein during infection, coordinates the assembly, budding, and release of viral particles into the host cell. It also regulates viral transcription and RNA replication. This study sought to identify small molecules that could potentially inhibit the VP40 protein by targeting the N-terminal domain using an in silico approach. The statistical quality of AutoDock Vina’s capacity to discriminate between inhibitors and decoys was determined, and an area under the curve of the receiver operating characteristic (AUC-ROC) curve of 0.791 was obtained. A total of 29,519 natural-product-derived compounds from Chinese and African sources as well as 2738 approved drugs were successfully screened against VP40. Using a threshold of −8 kcal/mol, a total of 7, 11, 163, and 30 compounds from the AfroDb, Northern African Natural Products Database (NANPDB), traditional Chinese medicine (TCM), and approved drugs libraries, respectively, were obtained after molecular docking. A biological activity prediction of the lead compounds suggested their potential antiviral properties. In addition, random-forest- and support-vector-machine-based algorithms predicted the compounds to be anti-Ebola with IC50 values in the micromolar range (less than 25 μM). A total of 42 natural-product-derived compounds were identified as potential EBOV inhibitors with desirable ADMET profiles, comprising 1, 2, and 39 compounds from NANPDB (2-hydroxyseneganolide), AfroDb (ZINC000034518176 and ZINC000095485942), and TCM, respectively. A total of 23 approved drugs, including doramectin, glecaprevir, velpatasvir, ledipasvir, avermectin B1, nafarelin acetate, danoprevir, eltrombopag, lanatoside C, and glycyrrhizin, among others, were also predicted to have potential anti-EBOV activity and can be further explored so that they may be repurposed for EVD treatment. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area calculations corroborated the stability and good binding affinities of the complexes (−46.97 to −118.9 kJ/mol). The potential lead compounds may have the potential to be developed as anti-EBOV drugs after experimental testing.

Chemical Compositions and Experimental and Computational Modeling of the Anticancer Effects of Cnidocyte Venoms of Jellyfish Cassiopea andromeda and Catostylus mosaicus on Human Adenocarcinoma A549 Cells

Nowadays, major attention is being paid to curing different types of cancers and is focused on natural resources, including oceans and marine environments. Jellyfish are marine animals with the ability to utilize their venom in order to both feed and defend. Prior studies have displayed the anticancer capabilities of various jellyfish. Hence, we examined the anticancer features of the venom of Cassiopea andromeda and Catostylus mosaicus in an in vitro situation against the human pulmonary adenocarcinoma (A549) cancer cell line. The MTT assay demonstrated that both mentioned venoms have anti-tumoral ability in a dose-dependent manner. Western blot analysis proved that both venoms can increase some pro-apoptotic factors and reduce some anti-apoptotic molecules that lead to the inducing of apoptosis in A549 cells. GC/MS analysis demonstrated some compounds with biological effects, including anti-inflammatory, antioxidant and anti-cancer activities. Molecular docking and molecular dynamic showed the best position of each biologically active component on the different death receptors, which are involved in the process of apoptosis in A549 cells. Ultimately, this study has proven that both venoms of C. andromeda and C. mosaicus have the capability to suppress A549 cells in an in vitro condition and they might be utilized in order to design and develop brand new anticancer agents in the near future.

Design, synthesis, molecular docking studies and biological evaluation of thiazole carboxamide derivatives as COX inhibitors

BACKGROUND

Nonsteroidal anti-inflammatory drugs (NSAIDs) have been the most commonly used class of medications worldwide for the last three decades.

OBJECTIVES

This study aimed to design and synthesize a novel series of methoxyphenyl thiazole carboxamide derivatives and evaluate their cyclooxygenase (COX) suppressant and cytotoxic properties.

METHODS

The synthesized compounds were characterized using ¹H, ¹³C-NMR, IR, and HRMS spectrum analysis and were evaluated for their selectivity towards COX-1 and COX-2 using an in vitro COX inhibition assay kit. Besides, their cytotoxicity was evaluated using the Sulforhodamine B (SRB) assay. Moreover, molecular docking studies were conducted to identify the possible binding patterns of these compounds within both COX-1 and COX-2 isozymes, utilizing human X-ray crystal structures. The density functional theory (DFT) analysis was used to evaluate compound chemical reactivity, which was determined by calculating the frontier orbital energy of both HOMO and LUMO orbitals, as well as the HOMO-LUMO energy gap. Finally, the QiKProp module was used for ADME-T analysis.

RESULTS

The results revealed that all synthesized molecules have potent inhibitory activities against COX enzymes. The percentage of inhibitory activities at 5 µM concentration against the COX2 enzyme was in the range of 53.9-81.5%, while the percentage against the COX-1 enzyme was 14.7-74.8%. That means almost all of our compounds have selective inhibition activities against the COX-2 enzyme, and the most selective compound was 2f, with selectivity ratio (SR) value of 3.67 at 5 µM concentration, which has a bulky group of trimethoxy on the phenyl ring that could not bind well with the COX-1 enzyme. Compound 2h was the most potent, with an inhibitory activity percentage at 5 µM concentration of 81.5 and 58.2% against COX-2 and COX-1, respectively. The cytotoxicity of these compounds was evaluated against three cancer cell lines: Huh7, MCF-7, and HCT116, and negligible or very weak activities were observed for all of these compounds except compound 2f, which showed moderate activities with IC₅₀ values of 17.47 and 14.57 µM against Huh7 and HCT116 cancer cell lines, respectively. Analysis of the molecular docking suggests 2d, 2e, 2f, and 2i molecules were bound to COX-2 isozyme favorably over COX-1 enzyme, and their interaction behaviors within COX-1 and COX-2 isozymes were comparable to celecoxib, as an ideal selective COX-2 drug, which explained their high potency and COX-2 selectivity. The molecular docking scores and expected affinity using the MM-GBSA approach were consistent with the recorded biological activity. The calculated global reactivity descriptors, such as HOMO and LUMO energies and the HOMO-LUMO gaps, confirmed the key structural features required to achieve favorable binding interactions and thus improve affinity. The in silico ADME-T studies asserted the druggability of molecules and have the potential to become lead molecules in the drug discovery process.

CONCLUSION

In general, the series of the synthesized compounds had a strong effect on both enzymes (COX-1 and COX-2) and the trimethoxy compound 2f was more selective than the other compounds.

Computational Chemistry for the Identification of Lead Compounds for Radiotracer Development

The use of computer-aided drug design (CADD) for the identification of lead compounds in radiotracer development is steadily increasing. Traditional CADD methods, such as structure-based and ligand-based virtual screening and optimization, have been successfully utilized in many drug discovery programs and are highlighted throughout this review. First, we discuss the use of virtual screening for hit identification at the beginning of drug discovery programs. This is followed by an analysis of how the hits derived from virtual screening can be filtered and culled to highly probable candidates to test in in vitro assays. We then illustrate how CADD can be used to optimize the potency of experimentally validated hit compounds from virtual screening for use in positron emission tomography (PET). Finally, we conclude with a survey of the newest techniques in CADD employing machine learning (ML).

Integrated in silico and experimental discovery of trimeric peptide ligands targeting Butyrylcholinesterase

Butyrylcholinesterase (BChE) is recognized as a high value biotherapeutic in the treatment of Alzheimer's disease and drug addiction. This study presents the rational design and screening of an in-silico library of trimeric peptides against BChE and the experimental characterization of peptide ligands for purification. The selected peptides consistently afforded high BChE recovery (> 90 %) and purity, yielding up to a 1000-fold purification factor. This study revealed a marked anti-correlated conformational movement governed by the ionic strength and pH of the aqueous environment, which ultimately controls BChE binding and release during chromatographic purification; and highlighted the role of residues within and allosteric to the catalytic triad of BChE in determining biorecognition, thus providing useful guidance for ligand design and affinity maturation.

Discovery, Validation, and Target Prediction of Antibacterial and Antidiabetic Components of Archidendron clypearia Based on a Combination of Multiple Analytical Methods

Archidendron clypearia (A. clypearia), a Fabaceae family member, is widely used as an anti-inflammatory herbal medicine; however, its antibacterial and antidiabetic properties have not been extensively investigated. This study aimed to systematically analyze the antibacterial and antidiabetic components of A. clypearia by utilizing a combination of analytical methods. First, ten different polarity extracts were analyzed through ultra-performance liquid chromatography (UPLC), and their antibacterial and antidiabetic activities were evaluated. Then the spectrum-effect relationship between the biological activity and UPLC chromatograms was analyzed by partial least squares regression and gray relational analysis, followed by corresponding validation using isolated components. Finally, network pharmacology and molecular docking were implemented to predict the main antibacterial target components of A. clypearia and the enzyme inhibition active sites of α-amylase and α-glucosidase. P15, P16, and P20 were found to be the antibacterial and antidiabetic active components. The inhibitory effect of 7-O-galloyltricetiflavan (P15) on six bacterial species may be mediated through the lipid and atherosclerosis pathway, prostate cancer, adherens junctions, and targets such as SRC, MAPK1, and AKT1. The molecular docking results revealed that 7-O-galloyltricetiflavan and 7,4'-di-O-galloyltricetiflavan (P16/P20) can bind to α-amylase and α-glucosidase pockets with binding energies lower than -6 kcal/mol. Our study provides guidance for the development of antibacterial and antidiabetic products based on A. clypearia and can be used as a reference for the evaluation of bioactivity of other herbs.

Computational design of new tacrine analogs: an in silico prediction of their cholinesterase inhibitory, antioxidant, and hepatotoxic activities

Tacrine, the first drug approved for the treatment of Alzheimer's disease (AD), is a non-competitive cholinesterase inhibitor withdrawn due to its acute hepatotoxicity. However, new non-hepatotoxic forms of tacrine have been actively researched. Moreover, several recent reports have shown that oxidative stress is the cause of damage and plays a role in the pathogenesis of several neurodegenerative diseases including AD. The aim of the present study is the design of new easily synthesized tacrine analogs with less hepatotoxicity and potent antioxidant activity. In this context, a library of 34 novel tacrine analogs bearing an antioxidant fragment was designed and evaluated for its hepatotoxicity as well as anticholinesterase and antioxidant activities using computational methods. As a result, six new tacrine analogs have been proposed as potential inhibitors of cholinesterase with antioxidant activity and low or no hepatotoxicity. Furthermore, ADME calculations suggest that these compounds are promising oral drug candidates. Communicated by Ramaswamy H. Sarma.

Synthesis of New Ester Derivatives of Salicylic Acid and Evaluation of Their COX Inhibitory Potential

Salicylic acid is an NSAID with serious side effects on the GIS. The side effects of salicylic acid on the GIS are slightly reduced by acetylating salicylic acid. 12 new ester analogs of salicylic acid were synthesized with high yields in this study. The chemical structures of the synthesized compounds were characterized by ¹ H-NMR, ¹³ C-NMR, and HRMS spectra. The inhibitory potential of the compounds was evaluated on COXs by in vitro and in silico studies. The COX2 inhibitory activity of the most potent inhibitor MEST1 (IC₅₀ : 0.048 μM) was found to be much higher than the COX2 inhibitory activity of aspirin (IC₅₀ : 2.60 μM). In docking studies, the strongest inhibitor among the compounds synthesized was predicted to be MEST1, with the lowest binding energy. Docking studies revealed that MEST1 extends from the hydrophobic channel to the top of the cyclooxygenase active site, forming various interactions with residues in the binding pocket.

Pyridyl Methylsulfinyl Benzimidazole Derivatives as Promising Agents against Giardia lamblia and Trichomonas vaginalis

Protozoan parasites, such as Giardia lamblia and Trichomonas vaginalis, cause the most prevalent infections in humans in developing countries and provoke significant morbidity and mortality in endemic countries. Despite its side-effects, metronidazole is still the drug of choice as a giardiacidal and trichomonacidal tissue-active agent. However, the emergence of metronidazole resistance and its evolved strategies of parasites to evade innate host defenses have hindered the identification and development of new therapeutic strategies against these parasites. Here, we tested five synthesized benzimidazole derivatives as possible drugs for treating giardiasis and trichomoniasis, probing the bifunctional enzyme glucose 6-phosphate dehydrogenase::6-phosphogluconolactone from G. lamblia (GlG6PD::6PGL) and T. vaginalis (TvG6PD::6PGL) as a drug target. The investigated benzimidazole derivatives were H-B2M1, H-B2M2, H₂N-BZM6, O₂N-BZM7, and O₂N-BZM9. The recombinant enzymes were used in inhibition assays, and in silico computational predictions and spectroscopic studies were applied to follow the structural alteration of the enzymes and identify the possible mechanism of inhibition. We identified two potent benzimidazole compounds (O₂N-BZM7 and O₂N-BZM9), which are capable of inhibiting both protozoan G6PD::6PGL enzymes and in vitro assays with these parasites, showing that these compounds also affect their viability. These results demonstrate that other therapeutic targets of the compounds are the enzymes GlG6PD::6PGL and TvG6PD::6PGL, which contribute to their antiparasitic effect and their possible use in antigiardial and trichomonacidal therapies.

Deffini: A family-specific deep neural network model for structure-based virtual screening

Deep learning-based virtual screening methods have been shown to significantly improve the accuracy of traditional docking-based virtual screening methods. In this paper, we developed Deffini, a structure-based virtual screening neural network model. During training, Deffini learns protein-ligand docking poses to distinguish actives and decoys and then to predict whether a new ligand will bind to the protein target. Deffini outperformed Smina with an average AUC ROC of 0.92 and AUC PRC of 0.44 in 3-fold cross-validation on the benchmark dataset DUD-E. However, when tested on the maximum unbiased validation (MUV) dataset, Deffini achieved poor results with an average AUC ROC of 0.517. We used the family-specific training approach to train the model to improve the model performance and concluded that family-specific models performed better than the pan-family models. To explore the limits of the predictive power of the family-specific models, we constructed Kernie, a new protein kinase dataset consisting of 358 kinases. Deffini trained with the Kernie dataset outperformed all recent benchmarks on the MUV kinases, with an average AUC ROC of 0.745, which highlights the importance of quality datasets in improving the performance of deep neural network models and the importance of using family-specific models.

Structural insights into the inhibition of the nsP2 protease from Chikungunya virus by molecular modeling approaches

Chikungunya virus (CHIKV) is the etiological agent of the Chikungunya fever which has spread worldwide. Clinically, this disease may lead to prolonged incapacitating joint pain that can compromise remarkably the patients' quality of life. However, there are no licensed vaccines or specific drugs to fight this infection yet, making the search for novel therapies an imperative need. In this scenario, the CHIKV nsP2 protease emerged as an attractive therapeutic target once this protein plays a pivotal role in viral replication and pathogenesis. Hence, we investigated the structural basis for the inhibition of this enzyme by using molecular docking and dynamics simulations. Compounds with inhibitory activities against CHIKV nsP2 protease determined experimentally were selected from the literature. Docking studies with a set of stereoisomers showed that trans isomers, but not cis ones, bound close to the catalytic dyad which may explain isomerism requirements to the enzyme's inhibition. Further, binding mode analyses of other known inhibitors revealed highly conserved contacts between inhibitors and enzyme residues like N1011, C1013, A1046, Y1079, N1082, W1084, L1205, and M1242. Molecular dynamics simulations reinforced the importance of some of these interactions and pointed to nonpolar interactions as the main forces for inhibitors' binding. Finally, we observed that true inhibitors exhibited lower structural fluctuation, higher ligand efficiency and did not induce significant changes in protein correlated motions. Collectively, our findings might allow discerning true inhibitors from false ones and can guide drug development efforts targeting the nsP2 protease to fight CHIKV infections in the future.

Cheminformatics Identification of Phenolics as Modulators of Penicillin-Binding Protein 2a of Staphylococcus aureus: A Structure–Activity-Relationship-Based Study

The acquisition of penicillin-binding protein (PBP) 2a in resistant strains of Staphylococcus aureus allows for the continuous production of cell walls even after the inactivation of intrinsic PBPs. Thus, the discovery of novel therapeutics with enhanced modulatory activity on PBP2a is crucial, and plant secondary metabolites, such as phenolics, have found relevance in this regard. In this study, using computational techniques, phenolics were screened against the active site of PBP2a, and the ability of the lead phenolics to modulate PBP2a’s active and allosteric sites was studied. The top-five phenolics (leads) identified through structure–activity-based screening, pharmacokinetics and synthetic feasibility evaluations were subjected to molecular dynamics simulations. Except for propan-2-one at the active site, the leads had a higher binding free energy at both the active and allosteric sites of PBP2a than amoxicillin. The leads, while promoting the thermodynamic stability of PBP2a, showed a more promising affinity at the allosteric site than the active site, with silicristin (−25.61 kcal/mol) and epicatechin gallate (−47.65 kcal/mol) having the best affinity at the active and allosteric sites, respectively. Interestingly, the modulation of Tyr446, the active site gatekeeper residue in PBP2a, was noted to correlate with the affinity of the leads at the allosteric site. Overall, these observations point to the leads’ ability to inhibit PBP2a, either directly or through allosteric modulation with conventional drugs. Further confirmatory in vitro studies on the leads are underway.

Anticholinesterase Inhibition, Drug-Likeness Assessment, and Molecular Docking Evaluation of Milk Protein-Derived Opioid Peptides for the Control of Alzheimer’s Disease

The drug-likeness and pharmacokinetic properties of 23 dairy-protein-derived opioid peptides were studied using SwissADME and ADMETlab in silico tools. All the opioid peptides had poor drug-like properties based on violations of Lipinski’s rule-of-five. Moreover, prediction of their pharmacokinetic properties showed that the peptides had poor intestinal absorption and bioavailability. Following this, two well-known opioid peptides (βb-casomorphin-5, βb-casomorphin-7) from A1 bovine milk and caffeine (positive control) were selected for in silico molecular docking and in vitro inhibition study with two cholinesterase enzyme receptors important for the pathogenesis of Alzheimer’s disease. Both peptides showed higher binding free energies and inhibitory activities to butyrylcholinesterase (BChE) than caffeine, but in vitro binding energy values were lower than those from the docking model. Moreover, the two casomorphins had lower inhibitory properties against acetylcholinesterase (AChE) than caffeine, although the docking model predicted the opposite. At 1 mg/mL concentrations, βb-casomorphin-5 and βb-casomorphin-7 showed promising results in inhibiting both cholinesterases (i.e., respectively 34% and 43% inhibition of AChE, and 67% and 81% inhibition of BChE). These dairy-derived opioid peptides have the potential to treat Alzheimer’s disease via cholinesterase inhibition. However, appropriate derivatization may be required to improve their poor predicted intestinal absorption and bioavailability.

Macroalgae Bioactive Compounds for the Potential Antiviral of SARS-CoV-2: An In Silico Study

Coronavirus disease (COVID-19), which was due to novel coronavirus was detected in December 2019 in Wuhan, China for the first time and spread rapidly became a global pandemic. This study aimed to predict the potential of macroalgae compounds as SARS-CoV-2 antiviral by inhibiting of ACE2 receptor through in silico approach. Twenty-seven macroalgae compounds were obtained from PubChem (NCBI, USA), while target protein ACE2 receptor was collected from Protein Data Bank (PDB). Then the initial screening study drug-likeness conducted by Lipinski rule of five web server and prediction of bioactive probability carried out by PASS (Prediction of activity spectra for biologically active substances) Online web server. After those compounds were approved by Lipinski’s rule of five and PASS online prediction web server, the blind docking simulation was performed using PyRx 0.8 software to show binding energy value. Molecular interaction analysis was done using BIOVIA Discovery Studio 2016 v16.1.0 and PyMOL v2.4.1 software. There are six macroalgae compounds approved by Lipinski’s rule of five and PASS Online Analysis. The result is that macroalgae compound siphonaxanthin among 27 macroalgae compound showed strong binding energy to bind ACE2 receptor with -8.8 kcal/mol. This study also used the SARS-CoV-2 drugs as positive control: remdesivir, molnupiravir, baricitinib, lopinavir, oseltamivir, and favipiravir. The result shows that siphonaxanthin has lowest binding energy than the common SARS-CoV-2 drug. Macroalgae compounds are predicted to have potential as SARS-CoV-2 antiviral. Thus, extension studies need to investigate by in vitro and in vivo analysis for confirmation the siphonaxanthin’s inhibitory activity in combat SARS-CoV-2.

In silico mutational analysis of ACE2 to check the susceptibility of lung cancer patients towards COVID-19

Being the second major cause of death worldwide, lung cancer poses a significant threat to the health of patients. This worsened during the era of pandemic since lung cancer is found to be more prone to SARS-CoV-2 infection. Many recent studies imply a high frequency of COVID-19 infection associated severe outcome. However, molecular studies are still lacking in this respect. Hence the current study is designed to investigate the binding affinities of ACE2 lung cancer mutants with the viral spike protein to find the susceptibility of respective mutants carrying patients in catching the virus. Quite interestingly, our study found lesser binding affinities of all the selected mutants thus implying that these cancer patients might be less affected by the virus than others. These results are opposed to the recent studies’ propositions and open new avenues for more in-depth studies.

Study on the Mechanism of Üstikuddus Sherbiti in Ischemic Cerebrovascular Diseases: Based on Network Pharmacology

This paper aims to study the potential biological mechanism of Üstikuddus Sherbiti (ÜS) in the treatment of ischemic cerebrovascular diseases (ICVD) by the network pharmacology method. Traditional Chinese Medicine Systems Pharmacology (TCMSP) database was used to obtain effective constituents of ÜS by screening eligible oral utilization, drug similarity, and blood-brain barrier permeability threshold. By drug target prediction and stroke treatment target mining, 2 target data sets were analyzed to find intersection targets and the corresponding constituents were used as active constituents. An active constituent target network and an effective constituent target network were constructed by using Cytoscape 3.7.2 software. Degree parameters of the effective constituent target network were analyzed to find important effective constituents and targets. Through protein-protein interaction (PPI) analysis/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, potential signaling pathways of ÜS in ischemic stroke were found out. AutoDock was used for molecular docking verification. A total of 90 active constituents of ÜS were screened out. There were 10 active constituents against ICVD, including quercetin, luteolin, kaempferol, and naringenin, and 10 important targets for anticerebral ischemia, namely, PIK3CA, APP, PIK3R1, MAPK1, MAPK3, AKT1, PRKCD, Fyn, RAC1, and NF-κB1. Based on the protein interaction network, the important targets of ÜS were significantly enriched in PI3K-Akt signaling pathway, neuroactive ligand-receptor interaction pathway, Ras signaling pathway, etc. ÜS in ICVD has characteristics like multiple targets, multiple approaches, and multiple pathways. Results of molecular docking showed that the active components in ICVD had a good binding ability with the key targets. Its main biological mechanism may be related to the PI3K-Akt and Ras-MAPK centered signaling pathway. Our study demonstrated that ÜS exerted the effect of treating ICVD by regulating multiple targets and multiple channels with multiple components through the method of network pharmacology and molecular docking.

Molecular docking-based interactions in QSAR studies on Mycobacterium tuberculosis ATP synthase inhibitors

Tuberculosis (TB) is a global threat with a large burden across the continents in terms of mortality, morbidity, and financial losses. The disease has evolved into multi-drug-resistant (MDR-TB) and extensively drug-resistant (XDR-TB) tuberculosis owing to numerous factors ranging from patients' non-compliance to demographical implications. There have been very few new drugs for resistant TB. Resistance has already been reported even for the newly introduced drug bedaquiline. An attempt has been made to integrate both structure-based and QSAR drug design techniques (QSAR-SBDD) for the identification of novel leads. The docking scores normally do not correlate with the activity. Hence, the docking results have been analysed in terms of the number of interactions rather than docking scores. The parameters derived from interactions have been used in developing the QSAR models. The best model shows a good correlation (r = 0.908) between the activity and interaction parameter 'C' describing the sum of all the interactions with each amino acid residue. This model also predicts external dataset with a good correlation (r_ext = 0.851) and can be used for the identification of novel chemical entities (NCEs) and repurposed drugs for TB therapeutics.

Selectivity of enzymes involved in the formation of opposite enantiomeric series of p-menthane monoterpenoids in peppermint and Japanese catnip

Peppermint (Mentha x piperita L.) and Japanese catnip (Schizonepeta tenuifolia (Benth.) Briq.) accumulate p-menthane monoterpenoids with identical functionalization patterns but opposite stereochemistry. In the present study, we investigate the enantioselectivity of multiple enzymes involved in monoterpenoid biosynthesis in these species. Based on kinetic assays, mint limonene synthase, limonene 3-hydroxylase, isopiperitenol dehydrogenase, isopiperitenone reductase, and menthone reductase exhibited significant enantioselectivity toward intermediates of the pathway that proceeds through (-)-4S-limonene. Limonene synthase, isopiperitenol dehydrogenase and isopiperitenone reductase of Japanese catnip preferred intermediates of the pathway that involves (+)-4R-limonene, whereas limonene 3-hydroxylase was not enantioselective, and the activities of pulegone reductase and menthone reductase were too low to acquire meaningful kinetic data. Molecular modeling studies with docked ligands generally supported the experimental data obtained with peppermint enzymes, indicating that the preferred enantiomer was aligned well with the requisite cofactor and amino acid residues implicated in catalysis. A striking example for enantioselectivity was peppermint (-)-menthone reductase, which binds (-)-menthone with exquisite affinity but was predicted to bind (+)-menthone in a non-productive orientation that positions its carbonyl functional group at considerable distance to the NADPH cofactor. The work presented here lays the groundwork for structure-function studies aimed at unraveling how enantioselectivity evolved in closely related species of the Lamiaceae and beyond.

Astaxanthin-Mediated Bacterial Lethality: Evidence from Oxidative Stress Contribution and Molecular Dynamics Simulation

The involvement of cellular oxidative stress in antibacterial therapy has remained a topical issue over the years. In this study, the contribution of oxidative stress to astaxanthin-mediated bacterial lethality was evaluated in silico and in vitro. For the in vitro analysis, the minimum inhibitory concentration (MIC) of astaxanthin was lower than that of novobiocin against Staphylococcus aureus but generally higher than those of the reference antibiotics against other test organisms. The level of superoxide anion of the tested organisms increased significantly following treatment with astaxanthin when compared with DMSO-treated cells. This increase compared favorably with those observed with the reference antibiotics and was consistent with a decrease in the concentration of glutathione (GSH) and corresponding significant increase in ADP/ATP ratio. These observations are suggestive of probable involvement of oxidative stress in antibacterial capability of astaxanthin and in agreement with the results of the in silico evaluations, where the free energy scores of astaxanthins' complexes with topoisomerase IV ParC and ParE were higher than those of the reference antibiotics. These observations were consistent with the structural stability and compactness of the complexes as astaxanthin was observed to be more stable against topoisomerase IV ParC and ParE than DNA Gyrase A and B. Put together, findings from this study underscored the nature and mechanism of antibacterial action of astaxanthin that could suggest practical approaches in enhancing our current knowledge of antibacterial arsenal and aid in the novel development of alternative natural topo2A inhibitor.

Molecular Docking Analysis of Ficus religiosa Active Compound with Anti-Inflammatory Activity by Targeting Tumour Necrosis Factor Alpha and Vascular Endothelial Growth Factor Receptor in Diabetic Wound Healing

BACKGROUND: Diabetes mellitus contributes to the delayed healing of wounds causes disturbance of inflammatory cytokine. Tumour necrosis factor alpha (TNF-alpha) and Vascular Endothelial Growth Factor Receptor (VEGFR) both have a role in the persistent inflammation associated with diabetic wounds. Ficus religiosa has developed a reputation as a traditional wound healer among some java people in Indonesia. AIM: Our study aims to discover the molecular interaction between the active constituents of F. religiosa with TNF-alpha and VEGFR. MATERIALS AND METHODS: This research was conducted in computerized molecular docking using Protein database, Pymol, Discovery studio, and Pyrex software. A thorough literature search was conducted to identify the potential compound and molecular target for diabetic wounds. Analysis of its anti-inflammatory properties was also carried out using a passonline webserver. Pharmacokinetic analysis was performed using the Lipinski Rule of Five websites and the PreADMET website. RESULTS: Each of the study’s active compounds has a good pharmacokinetic profile. The predictions of the compound’s structure indicate that it has a strong anti-inflammatory impact. Lupenyl acetate and Lanosterol bind more strongly to the TNF-alpha than the natural ligand, but Piperine binds more strongly to VEGFR. CONCLUSIONS: Lupenyl acetate, Lanosterol, and Piperine compounds have anti-inflammatory effects through inhibition of TNF-alpha and VEGFR. In addition, this compound has potential to become a drug because it has good pharmacokinetics. Future studies are required to determine the effectiveness and toxicity of Lupenyl acetate, Lanosterol, and Piperine as potential treatment in diabetic wounds.

Truncated human angiotensin converting enzyme 2; a potential inhibitor of SARS-CoV-2 spike glycoprotein and potent COVID-19 therapeutic agent

The current pandemic of Covid-19 caused by SARS-CoV-2 is continued to spread globally and no potential drug or vaccine against it is available. Spike (S) glycoprotein is the structural protein of SARS-CoV-2 located on the envelope surface, involve in interaction with angiotensin converting enzyme 2 (ACE2), a cell surface receptor, followed by entry into the host cell. Thereby, blocking the S glycoprotein through potential inhibitor may interfere its interaction with ACE2 and impede its entry into the host cell. Here, we present a truncated version of human ACE2 (tACE2), comprising the N terminus region of the intact ACE2 from amino acid position 21-119, involved in binding with receptor binding domain (RBD) of SARS-CoV-2. We analyzed the in-silico potential of tACE2 to compete with intact ACE2 for binding with RBD. The protein-protein docking and molecular dynamic simulation showed that tACE2 has higher binding affinity for RBD and form more stabilized complex with RBD than the intact ACE2. Furthermore, prediction of tACE2 soluble expression in E. coli makes it a suitable candidate to be targeted for Covid-19 therapeutics. This is the first MD simulation based findings to provide a high affinity protein inhibitor for SARS-CoV-2 S glycoprotein, an important target for drug designing against this unprecedented challenge.Communicated by Ramaswamy H. Sarma.

Studies on the recombinant production and anticancer activity of thermostable L- asparaginase I from Pyrococcus abyssi

L-Asparaginase catalysing the breakdown of L-Asparagine to L-Aspartate and ammonia is an enzyme of therapeutic importance in the treatment of cancer, especially the lymphomas and leukaemia. The present study describes the recombinant production, properties and anticancer potential of enzyme from a hyperthermophilic archaeon Pyrococcus abyssi. There are two genes coding for asparaginase in the genome of this organism. A 918 bp gene encoding 305 amino acids was PCR amplified and cloned in BL21 (DE3) strain of E. coli using pET28a (+) plasmid. The production of recombinant enzyme was induced under 0.5mM IPTG, purified by selective heat denaturation and ion exchange chromatography. Purified enzyme was analyzed for kinetics, in silico structure and anticancer properties. The recombinant enzyme has shown a molecular weight of 33 kDa, specific activity of 1175 U/mg, KM value 2.05mM, optimum temperature and pH 80°C and 8 respectively. No detectable enzyme activity found when L-Glutamine was used as the substrate. In silico studies have shown that the enzyme exists as a homodimer having Arg11, Ala87, Thr110, His112, Gln142, Leu172, and Lys232 being the putative active site residues. The free energy change calculated by molecular docking studies of enzyme and substrate was found as ∆G - 4.5 kJ/mole indicating the affinity of enzyme with the substrate. IC50 values of 5U/mL to 7.5U/mL were determined for FB, caco2 cells and HepG2 cells. A calculated amount of enzyme (5U/mL) exhibited 78% to 55% growth inhibition of caco2 and HepG2 cells. In conclusion, the recombinant enzyme produced and characterized in the present study offers a good candidate for the treatment of cancer. The procedures adopted in the present study can be prolonged for in vivo studies.

Assessment of anti-inflammatory-like, antioxidant activities and molecular docking of three alkynyl-substituted 3-ylidene-dihydrobenzo[d]isothiazole 1,1-dioxide derivatives

The presence of enyne and benzoisothiazole functions in the molecular architecture of compounds 1, 2 and 3 were expected to provide biochemical activities. In the present work, we first examined the molecular surface contact of three alkynyl-substituted 3-ylidenedihydrobenzo[d] isothiazole 1,1-dioxides. The analysis of the Hirshfeld surfaces reveals that only compound 3 exhibited a well-defined red spots, indicating intermolecular interactions identified as S-O⋯H, C-H⋯O and C-O⋯H contacts. Comparative fingerprint histograms of the three compounds show that close pair interactions are dominated by C-H⋯H-C contact. By UV-visible analysis, compound 1 showed the most intense absorbances at 407 and 441 nm, respectively. The radical scavenging activity explored in the DPPH test, shows that only 1 exhibited low anti-radical activity. Furthermore, cellular antioxidant capacity of benzoisothiazoles 1-3 was investigated with PMA-activated HL-60 cells using chemiluminescence and fluorescence techniques in the presence of L-012 and Amplex Red probe, respectively. Results highlight that compound 1 exhibited moderate anti-ROS capacity while compounds 2 and 3 enhanced ROS production. The cytotoxicity test performed on HL-60 cells, using the MTS assay, confirmed the lack of toxicity of the tested benzoisothiazole 1 compared to 2 and 3 which show low cytotoxicity (≤30%). Anti-catalytic activity was evaluated by following the inhibitory potential of the benzoisothiazoles on MPO activity and depicted benzoisothiazoles-MPO interactions by docking. Both SIEFED and docking studies demonstrated an anti-catalytic activity of the tested benzoisothiazoles towards MPO with the best activity for compound 2.

Computer-assisted evaluation of plant-derived β-secretase inhibitors in Alzheimer’s disease

Background

Alzheimer’s disease (AD) is a progressive neurodegenerative age-related dementia that results in memory loss of elderly people. Many hypotheses have been formally articulated till now to decipher the pathogenesis of this disease. According to the compelling amyloidogenic hypothesis, β-secretase is a key regulatory enzyme in AD development and is therefore considered as one of the major targets for the development of drugs to treat AD. In this study, 40 plant-derived phytocompounds, proven to have β-secretase inhibitory activity in different laboratory experiments, were evaluated using computational approaches in order to identify the best possible β-secretase inhibitor(s).

Results

Amentoflavone (IFD score: − 7.842 Kcal/mol), Bilobetin (IFD score: − 7.417 Kcal/mol), and Ellagic acid (IFD score: − 6.923 Kcal/mol) showed highest β-secretase inhibitory activities with high binding affinity among all the selected phytocompounds and interacted with key amino acids, i.e., Asp32, Tyr71, and Asp228 in the catalytic site of β-secretase. Moreover, these three molecules exhibited promising results in different drug potential assessment experiments and displayed signs of correlation with significant pharmacological and biological activities.

Conclusion

Amentoflavone, Biolbetin, and Ellagic acid could be investigated further in developing β-secretase-dependent drug for the effective treatment of AD. However, additional in vivo and in vitro experiments might be required to strengthen the findings of this experiment.

An apta-aggregation based machine learning assay for rapid quantification of lysozyme through texture parameters

A novel assay technique that involves quantification of lysozyme (Lys) through machine learning is put forward here. This article reports the tendency of the well- documented Ellington group anti-Lys aptamer, to produce aggregates when exposed to Lys. This property of apta-aggregation has been exploited here to develop an assay that quantifies the Lys using texture and area parameters from a photograph of the elliptical aggregate mass through machine learning. Two assay sets were made for the experimental procedure: one with high Lys concentration between 25–100 mM and another with low concentration between 1–20 mM. The high concentration set had a sample volume of 10 μl while the low concentration set had a higher sample volume of 100 μl, in order to obtain the statistical texture values reliably from the aggregate mass. The platform exhibited an experimental limit of detection of 1 mM and a response time of less than 10 seconds. Further, two potential operating modes for the aptamer were hypothesized for this aggregation property and the more accurate mode among the two was ascertained through bioinformatics studies.

7-Amine-spiro[chromeno[4,3-b]quinoline-6,1'-cycloalkanes]: Synthesis and cholinesterase inhibitory activity of structurally modified tacrines

Five new examples of 9,10-chloro(bromo)-7-amine-spiro[chromeno[4,3-b]quinoline-6,1'-cycloalkanes] - in which cycloalkanes = cyclopentane, cyclohexane, and cycloheptane - were synthesized at yields of 42-56%, using a sequential one-pot two-step cyclocondensation reaction of three different scaffolds of 2-aminobenzonitriles and the respective spiro[chroman-2,1'-cycloalkan]-4-ones, and using AlCl₃ as the catalyst in a solvent-free method. Subsequently, the five new spirochromeno-quinolines and nine quinolines previously published by us (14 modified tacrine scaffolds) were subjected to AChE and BChE inhibitory activity evaluation. The molecule containing a spirocyclopentane derivative had the highest AChE and BChE inhibitory activity (IC₅₀ = 3.60 and 4.40 μM, respectively), and in general, the non-halogenated compounds were better inhibitors of AChE and BChE than the halogenated molecules. However, the inhibitory potency of compounds 3a-n was weaker than that of tacrine. By molecular docking simulations, it was found that the size of the spirocarbocyclic moieties is inversely proportional to the inhibitory activity of the cholinesterases, probably because an increase in the size of the spirocyclic component sterically hindered the interaction of tacrine derivatives with the active site of tested cholinesterases. The findings obtained here may help in the design and development of new anticholinesterase drugs.

The latest automated docking technologies for novel drug discovery

INTRODUCTION

Molecular docking has been consolidated as one of the most important methods in the molecular modeling field. It has been recognized as a prominent tool in the study of protein-ligand complexes, to describe intermolecular interactions, to accurately predict poses of multiple ligands, to discover novel promising bioactive compounds. Molecular docking methods have evolved in terms of their accuracy and reliability; but there are pending issues to solve for improving the connection between the docking results and the experimental evidence.

AREAS COVERED

In this article, the author reviews very recent innovative molecular docking applications with special emphasis on reverse docking, treatment of protein flexibility, the use of experimental data to guide the selection of docking poses, the application of Quantum mechanics(QM) in docking, and covalent docking.

EXPERT OPINION

There are several issues being worked on in recent years that will lead to important breakthroughs in molecular docking methods in the near future These developments are related to more efficient exploration of large datasets and receptor conformations, advances in electronic description, and the use of structural information for guiding the selection of results.

Crystal structure of AmpC BER and molecular docking lead to the discovery of broad inhibition activities of halisulfates against β-lactamases

AmpC BER is an extended-spectrum (ES) class C β-lactamase with a two-amino-acid insertion in the H10 helix region located at the boundary of the active site compared with its narrow spectrum progenitor. The crystal structure of the wild-type AmpC BER revealed that the insertion widens the active site by restructuring the flexible H10 helix region, which is the structural basis for its ES activity. Besides, two sulfates originated from the crystallization solution were observed in the active site. The presence of sulfate-binding subsites, together with the recognition of ring-structured chemical scaffolds by AmpC BER, led us to perform in silico molecular docking experiments with halisulfates, natural products isolated from marine sponge. Inspired by the snug fit of halisulfates within the active site, we demonstrated that halisulfate 3 and 5 significantly inhibit ES class C β-lactamases. Especially, halisulfate 5 is comparable to avibactam in terms of inhibition efficiency; it inhibits the nitrocefin-hydrolyzing activity of AmpC BER with a K_i value of 5.87 μM in a competitive manner. Furthermore, halisulfate 5 displayed moderate and weak inhibition activities against class A and class B/D enzymes, respectively. The treatment of β-lactamase inhibitors (BLIs) in combination with β-lactam antibiotics is a working strategy to cope with infections by pathogens producing ES β-lactamases. Considering the emergence and dissemination of enzymes insensitive to clinically-used BLIs, the broad inhibition spectrum and structural difference of halisulfates would be used to develop novel BLIs that can escape the bacterial resistance mechanism mediated by β-lactamases.

Chirality-Dependent Amino Acid Modulation of RNA Folding

The preponderance of a specific d- or l-chirality in fats, sugars, amino acids, nucleic acids, and so on is ubiquitous in nature, yet the biological origin of such chiral dominance (i.e., with one enantiomer overwhelmingly present) remains an open question. One plausible proposal for the predominance of l-chirality in amino acids could be through evolutionary templating of chiral RNA-folding via chaperone activity. To help evaluate this possibility, single molecule fluorescence experiments have been performed that measure the chiral dependence of chaperone folding dynamics for the simple tetraloop-tetraloop receptor (TL-TLR) tertiary binding motif in the presence of a series of chiral amino acids. Specifically, d- vs l-arginine is found to accelerate the unfolding of this RNA motif in a chirally selective fashion, with temperature-dependent studies of the kinetics performed to extract free energy, enthalpy, and entropy landscapes for the underlying thermodynamics. Furthermore, all-atom molecular dynamics (MD) simulations are pursued to provide additional physical insight into this chiral sensitivity, which reveal enantiomer-specific sampling of nucleic acid surfaces by d- vs l-arginine and support a putative mechanism for chirally specific denaturation of RNA tertiary structure by arginine but not other amino acids.

In vitro anti-leishmanial activity of Prunus armeniaca fractions on Leishmania tropica and molecular docking studies

Prunus armeniaca (L.) is a member of the Rosaceae, subfamily Prunoideae, shows anticancer, antitubercular, antimutagenic, antimicrobial, antioxidant, and cardioprotective activities. Here we fractionated the leaves extract of this highly medicinally important plant for antileishmanial activity. In the current study, the leaves extract was fractionated and characterized using column and thin layer chromatography by n-hexane, ethyl acetate, and methanol solvents. Twelve fractions were isolated and subjected for evaluation of their cytotoxicity and in vitro antileishmanial activity against promastigotes and amastigotes of Leishmania tropica. Among all fractions used, the fraction (F7) exhibited the strongest antileishmanial activity. The bioactive fraction was further characterized by spectroscopy (FTIR, UV-Vis), and GC-MS analysis. The in silico docking was carried out to find the active site of PTR1. All derived fractions exhibited toxicity in the safety range IC₅₀ > 100 μg/ml. The fraction (F7) showed significantly the highest antipromastigotes activity with IC₅₀11.48 ± 0.82 μg/ml and antiamastigotes activity with IC₅₀ 21.03 ± 0.98 μg/ml compared with control i.e. 11.60 ± 0.70 and 22.03 ± 1.02 μg/ml respectively. The UV-Vis spectroscopic analysis revealed the presence of six absorption peaks and the FTIR spectrum revealed the presence of alkane, aldehyde, carboxylic acid, thiols, alkynes, and carbonyls compounds The GC-MS chromatogram exhibited the presence of nine compounds: (a) benzeneethanol, alpha, beta dimethyl, (b)carbazic acid, 3-(1 propylbutylidene)-, ethyl ester, (c)1, 2-benzenedicarboxylic acid, diisooctyl ester, (d)benzeneethanamine a-methyl, (e)2aminononadecane, (f)2-heptanamine-5-methyl, (g)cyclobutanol, (h)cyclopropyl carbine, and (i)nitric acid, nonyl ester. Among all compounds, the 1, 2-benzenedicarboxylic acid, diisooctyl ester bound well to the PTR1 receptor. Fraction (F7) showed acceptable results with no cytotoxicity. However, in vivo studies are required in the future.