Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Functional characterization of novel mutations in the conserved region of EPSPS for herbicide resistance in pigeonpea: structure-based coherent design

Prospectively, agroecosystems for the growth of crops provide the potential fertile, productive, and tropical environment which attracts infestation by weedy plant species that compete with the primary crop plants. Infestation by weed is a major biotic stress factor faced by pigeonpea that hampers the productivity of the crop. In the modern era with the development of chemicals the problem of weed infestation is dealt with armours called herbicides. The most widely utilized, post-emergent, broad-spectrum herbicide has an essential active ingredient called glyphosate. Glyphosate mechanistically inhibits a chloroplastic enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) by competitively interacting with the PEP binding site which hinders the shikimate pathway and the production of essential aromatic amino acids (Phe, Tyr, Trp) and other secondary metabolites in plants. Moreover, herbicide spray for weed management is lethal to both the primary crop and the weeds. Therefore, it is critical to develop herbicide-resistant crops for field purposes to reduce the associated yield and economic losses. In this study, the in-silico analysis drove the selection and validation of the point mutations in the conserved region of the EPSPS gene, which confers efficient herbicide resistance to mutated-CcEPSPS enzyme along with the retention of the normal enzyme function. An optimized in-silico validation of the target mutation before the development of the genome-edited resistant plant lines is a prerequisite for testing their efficacy as a proof of concept. We validated the combination of GATIPS mutation for its no-cost effect at the enzyme level via molecular dynamic (MD) simulation.Communicated by Ramaswamy H. Sarma.

Transformer-Based Molecular Generative Model for Antiviral Drug Design

Since the Simplified Molecular Input Line Entry System (SMILES) is oriented to the atomic-level representation of molecules and is not friendly in terms of human readability and editable, however, IUPAC is the closest to natural language and is very friendly in terms of human-oriented readability and performing molecular editing, we can manipulate IUPAC to generate corresponding new molecules and produce programming-friendly molecular forms of SMILES. In addition, antiviral drug design, especially analogue-based drug design, is also more appropriate to edit and design directly from the functional group level of IUPAC than from the atomic level of SMILES, since designing analogues involves altering the R group only, which is closer to the knowledge-based molecular design of a chemist. Herein, we present a novel data-driven self-supervised pretraining generative model called "TransAntivirus" to make select-and-replace edits and convert organic molecules into the desired properties for design of antiviral candidate analogues. The results indicated that TransAntivirus is significantly superior to the control models in terms of novelty, validity, uniqueness, and diversity. TransAntivirus showed excellent performance in the design and optimization of nucleoside and non-nucleoside analogues by chemical space analysis and property prediction analysis. Furthermore, to validate the applicability of TransAntivirus in the design of antiviral drugs, we conducted two case studies on the design of nucleoside analogues and non-nucleoside analogues and screened four candidate lead compounds against anticoronavirus disease (COVID-19). Finally, we recommend this framework for accelerating antiviral drug discovery.

An investigation of binding interactions of tumor-targeted peptide conjugated polyphenols with the kinase domain of ephrin B4 and B2 receptors

Recent studies have shown that Ephrin receptors may be upregulated in several types of cancers including breast, ovarian and endometrial cancers, making them a target for drug design. In this work, we have utilized a target-hopping approach to design new natural product-peptide conjugates and examined their interactions with the kinase-binding domain of EphB4 and EphB2 receptors. The peptide sequences were generated through point mutations of the known EphB4 antagonist peptide TNYLFSPNGPIA. Their anticancer properties and secondary structures were analyzed computationally. Conjugates of most optimum of peptides were then designed by binding the N-terminal of the peptides with the free carboxyl group of the polyphenols sinapate, gallate and coumarate, which are known for their inherent anticancer properties. To investigate if these conjugates have a potential to bind to the kinase domain, we carried out docking studies and MMGBSA free energy calculations of the trajectories based on the molecular dynamics simulations, with both the apo and the ATP bound kinase domains of both receptors. In most cases binding interactions occurred within the catalytic loop region, while in some cases the conjugates were found to spread out across the N-lobe and the DFG motif region. The conjugates were further tested for prediction of pharmacokinetic properties using ADME studies. Our results indicated that the conjugates were lipophilic and MDCK permeable with no CYP interactions. These findings provide an insight into the molecular interactions of these peptides and conjugates with the kinase domain of the EphB4 and EphB2 receptor. As a proof of concept, we synthesized and carried out SPR analysis with two of the conjugates (gallate-TNYLFSPNGPIA and sinapate-TNYLFSPNGPIA). Results indicated that the conjugates showed higher binding with the EphB4 receptor and minimal binding to EphB2 receptor. Sinapate-TNYLFSPNGPIA showed inhibitory activity against EphB4. These studies reveal that some of the conjugates may be developed for further investigation into in vitro and in vivo studies and potential development as therapeutics.

Proteomics-based vaccine targets annotation and design of multi-epitope vaccine against antibiotic-resistant Streptococcus gallolyticus

Streptococcus gallolyticus is a non-motile, gram-positive bacterium that causes infective endocarditis. S. gallolyticus has developed resistance to existing antibiotics, and no vaccine is currently available. Therefore, it is essential to develop an effective S. gallolyticus vaccine. Core proteomics was used in this study together with subtractive proteomics and reverse vaccinology approach to find antigenic proteins that could be utilized for the design of the S. gallolyticus multi-epitope vaccine. The pipeline identified two antigenic proteins as potential vaccine targets: penicillin-binding protein and the ATP synthase subunit. T and B cell epitopes from the specific proteins were forecasted employing several immunoinformatics and bioinformatics resources. A vaccine (360 amino acids) was created using a combination of seven cytotoxic T cell lymphocyte (CTL), three helper T cell lymphocyte (HTL), and five linear B cell lymphocyte (LBL) epitopes. To increase immune responses, the vaccine was paired with a cholera enterotoxin subunit B (CTB) adjuvant. The developed vaccine was highly antigenic, non-allergenic, and stable for human use. The vaccine's binding affinity and molecular interactions with the human immunological receptor TLR4 were studied using molecular mechanics/generalized Born surface area (MMGBSA), molecular docking, and molecular dynamic (MD) simulation analyses. Escherichia coli (strain K12) plasmid vector pET-28a ( +) was used to examine the ability of the vaccine to be expressed. According to the outcomes of these computer experiments, the vaccine is quite promising in terms of developing a protective immunity against diseases. However, in vitro and animal research are required to validate our findings.

Synthesis, Anticancer Activity, and In Silico Modeling of Alkylsulfonyl Benzimidazole Derivatives: Unveiling Potent Bcl-2 Inhibitors for Breast Cancer

A series of alkylsulfonyl 1H-benzo[d]imidazole derivatives were synthesized and evaluated for anticancer activity against human breast cancer cells, MCF-7 in vitro. The cytotoxic potential was determined using the xCELLigence real-time cell analysis, and expression levels of genes related to microtubule organization, tumor suppression, apoptosis, cell cycle, and proliferation were examined by quantitative real-time polymerase chain reaction. Molecular docking against Bcl-2 was carried out using AutoDock Vina, while ADME studies were performed to predict the physicochemical and drug-likeness properties of the synthesized compounds. The results revealed that compounds 23 and 27 were the most potent cytotoxic derivatives against MCF-7 cells. Gene expression analysis showed that BCL-2 was the most prominent gene studied. Treatment of MCF-7 cells with compounds 23 and 27 resulted in significant downregulation of the BCL-2 gene, with fold changes of 128 and 256, respectively. Docking analysis predicted a strong interaction between the compounds and the target protein. Interestingly, all of the compounds exhibit a higher binding affinity toward Bcl-2 than the standard drug (compound 27 vina score = -9.6 kcal/mol, vincristine = -6.7 kcal/mol). Molecular dynamics simulations of compounds 23 and 27 showed a permanent stabilization in the binding site of Bcl-2 for 200 ns. Based on Lipinski and Veber's filters, all synthesized compounds displayed drug-like characteristics. These findings suggest that compounds 23 and 27 were the most promising cytotoxic compounds and downregulated the expression of the BCL-2 gene. These derivatives could be further explored as potential candidates for the treatment of breast cancer.

Association of epilepsy and neurological impairments with homozygous recessive missense mutations found in the genes responsible for ganglioside biosynthesis (ST3GAL5) and calcium voltage-gated channels (CACNA1H) - insights through molecular dynamic simulations

With over 2.2 million cases, the incidence rate of epilepsy in Pakistan is far higher than the rest of the world due primarily to the frequent, traditionally imposed cousin marriages. In the present study, comprehensive whole exome sequencing (WES) analyses of a three-generation family with four affected members presenting 'unexplained' childhood absence epilepsy (CAE), seizures and dementia, was performed in a quest to identify heritable, epilepsy-causal gene variants to better aid in carrier screening and genetic counselling. The WES data was generated, analyzed, and validated through Sanger's sequencing, molecular dynamic simulation (MDS) analysis, and molecular mechanics with generalized Born and surface area solvation (MM/GBSA) studies. Two homozygous recessive, missense mutations in ST3GAL5 (c.311A > G, p. His104Arg) and CACNA1H (c.6230G > A, p. Arg2077His) genes, earlier regarded as benign or of uncertain significance, have been identified as a potential etiology. Comparative MDS and free binding energy calculations revealed substantial structural perturbations in mutant forms of ST3GAL5 leading to decreased binding and reduced catalytic activity of the p.His104Arg and two other functional variants (p.Val74Glu and p.Arg288Ter) when compared with wild type. Our findings reinforce that WES analyses may uncover 'hidden', heritable variants and together with MDS and MM/GBSA may provide plausible clues to answer the unexplained causes of epilepsy for an effective management and better patient outcome. Further, revisit of epilepsy-associated mutational landscape in population context is imperative as the variants with 'benign' tags may turn out to be 'non-benign', when exist in combination with other benign.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening and molecular docking approaches to identify potential inhibitors against KIF2C to combat glioma

Glioma, a kind of malignant brain tumor, is extremely lethal. Kinesin family member 2C (KIF2C) was found to have an aberrant expression in several cancer types, including lung cancer and glioma. KIF2C may therefore be a useful therapeutic target for the treatment of glioma. In the current study, new drug candidates that may function as KIF2C enzyme inhibitors were discovered. MTi OpenScreen was used to carry out the structure-based virtual screening of an inbuilt drug library containing 150,000 compounds. These compounds belong to different classes, such as natural product-based compounds (NP-lib), purchasable approved drugs (Drugs-lib), and food constituents compound collection (FOOD-lib). Based on their binding affinities, a total of 84 compounds were further pushed to calculate ADMET properties. The compounds (16) meeting the ADMET cutoff ranges were then further docked to the receptor to find their plausible binding modes using the Glide tool's standard precision (SP) technique. The docking results were examined using the Glide gscore, and the best binding compounds (Rimacalib and Sarizotan) were chosen to test their stability with KIF2C protein through molecular dynamics (MD) simulation. Similarly, Principal Component Analysis and cross-correlation matrix were also examined. The MM/GBSA binding free energies showed a considerable energy contribution in the binding of hits with the KIF2C. Collectively, these findings strongly suggest the potential of the lead compounds to inhibit the biological function of KIF2C, emphasizing the need for further investigation in this area.Communicated by Ramaswamy H. Sarma.

Identification of potent anti-Candida metabolites produced by the soft coral associated Streptomyces sp. HC14 using chemoinformatics

Candida albicans is the most common pathogen responsible for both spontaneous and recurrent candidiasis. The available treatment of Candida infections has several adverse effects, and the development of new drugs is critical. The current study looked at the synthesis of anti-Candida metabolites by Streptomyces sp. HC14 recovered from a soft coral. Using the Plackett Burman design, the medium composition was formulated to maximize production. Using GC-MS, the compounds have been identified, and a cheminformatics approach has been used to identify the potential source of activity. The compounds that showed high potential for activity were identified as pyrrolo[1,2-a]pyrazine-1,4-dione, hexahydro-3-(phenylmethyl)-3 and di-n-octyl based on their docking score against the cytochrome monooxygenase (CYP51) enzyme in Candida albicans. As a result of their discovery, fewer molecules need to be chemically synthesized, and fermentation optimization maximizes their synthesis, providing a strong foundation for the development of novel anti-Candida albicans agents.

Assessment of the inhibitory potential of anabolic steroids towards human AKR1D1 by computational methods and in vitro evaluation

Exposure to xenobiotics can adversely affect biochemical reactions, including hepatic bile acid synthesis. Bile acids are essential for dissolving lipophilic compounds in the hydrophilic environment of the gastrointestinal tract. The critical micellar concentration of bile acids depends on the Δ4-reduction stereochemistry, with the 3-oxo-5β-steroid-Δ4-dehydrogenase (AKR1D1) introducing the cis ring A/B conformation. Loss-of-function mutations in AKR1D1 cause hepatic cholestasis, which, if left untreated can progress into steatosis and liver cirrhosis. Furthermore, AKR1D1 is involved in clearing steroids with an A-ring Δ4-double bond. Here, we tested whether anabolic-androgenic steroids (AAS), often taken off-label at high doses, might inhibit AKR1D1, thereby potentially causing hepatotoxicity. A computational molecular model was established and used for virtual screening of the DrugBank database consisting of 2740 molecules, yielding mainly steroidal hits. Fourteen AAS were selected for in vitro evaluation, as such compounds can reach high hepatic concentrations in an abuse situation. Nandrolone, clostebol, methasterone, drostanolone, and methenolone inhibited to various extent the AKR1D1-mediated reduction of testosterone. Molecular modeling suggests that 9 out of 14 investigated AAS are competitive inhibitors. Moreover quantum mechanical calculations show that nadrolone and clostebol are substrates of AKR1D1 with different activation energy barriers for the hydrogen transfer from cofactor to the C5 position affecting their turnover. In this multidisciplinary approach, we established a molecular model of AKR1D1, identified several AAS as inhibitors, and described their binding mode. This approach may be applied to study other classes of inhibitors including non-steroidal compounds.

Antioxidant evaluation and computational prediction of prospective drug-like compounds from polyphenolic-rich extract of Hibiscus cannabinus L. seed as antidiabetic and neuroprotective targets: assessment through in vitro and in silico studies

BACKGROUND

Reports have implicated diabetes mellitus (DM) and Alzheimer's disease (AD) as some of the global persistent health challenges with no lasting solutions, despite of significant inputs of modern-day pharmaceutical firms. This study therefore, aimed to appraise the in vitro antioxidant potential, enzymes inhibitory activities, and as well carry out in silico study on bioactive compounds from polyphenolic-rich extract of Hibiscus cannabinus seed (PEHc).

METHODS

In vitro antioxidant assays were performed on PEHc using standard methods while the identification of phytoconstituents was carried out with high performance liquid chromatography (HPLC). For the in silico molecular docking using Schrodinger's Grid-based ligand docking with energetics software, seven target proteins were retrieved from the database ( https://www.rcsb.org/ ).

RESULTS

HPLC technique identified twelve chemical compounds in PEHc, while antioxidant quantification revealed higher total phenolic contents (243.5 ± 0.71 mg GAE/g) than total flavonoid contents (54.06 ± 0.09 mg QE/g) with a significant (p < 0.05) inhibition of ABTS (IC₅₀ = 218.30 ± 0.87 µg/ml) and 1, 1-diphenyl-2-picrylhydrazyl free radicals (IC₅₀ = 227.79 ± 0.74 µg/ml). In a similar manner, the extract demonstrated a significant (p < 0.05) inhibitory activity against α-amylase (IC₅₀ = 256.88 ± 6.15 µg/ml) and α-glucosidase (IC₅₀ = 183.19 ± 0.23 µg/ml) as well as acetylcholinesterase (IC₅₀ = 262.95 ± 1.47 µg/ml) and butyrylcholinesterase (IC₅₀ = 189.97 ± 0.82 µg/ml), respectively. Furthermore, In silico study showed that hibiscetin (a lead) revealed a very strong binding affinity energies for DPP-4, (PDB ID: 1RWQ) and α-amylase (PDB ID: 1SMD), gamma-tocopherol ( for peptide-1 receptor; PDB ID: 3C59, AChE; PDB ID: 4EY7 and BChE; PDB ID: 7B04), cianidanol for α-glucosidase; PDB ID: 7KBJ and kaempferol for Poly [ADP-ribose] polymerase 1 (PARP-1); PDB ID: 6BHV, respectively. More so, ADMET scores revealed drug-like potentials of the lead compounds identified in PEHc.

CONCLUSION

As a result, the findings of this study point to potential drug-able compounds in PEHc that could be useful for the management of DM and AD.

HR-LCMS and evaluation of anti-diabetic activity of Hemidesmus indicus (anantmool): Kinetic study, and molecular modelling approach

This study delved into the exploration of novel antidiabetic medications acquired from natural resources, utilizing the Ayurvedic Rasayana herb Hemidesmus indicus through cutting-edge chemoprofiling and molecular modelling techniques. The methanolic extract of Hemidesmus indicus root exhibited the highest extractive yield (24.70 ± 0.08 %) and contained substantial levels of total phenolic and flavonoid content as 154.15 ± 1.24 mg Gallic Acid Equivalent/g extract and 70.61 ± 0.35 Quercetin Equivalent/g extract respectively. Invitro study revealed the potent inhibitory potential of methanolic extract of the herb against essential carbohydrate hydrolytic enzymes α-amylase (IC₅₀ = 4.19 ± 0.04 mg/ml) and α-glucosidase (IC₅₀ = 5.78 ± 0.10 mg/ml). Further, the enzyme kinetic study demonstrated the competitive mode of inhibition of both enzymes. HR-LCMS analysis identified the major phytoconstituents present in the extracts, including Solanocapsine, Cyclovirobuxine C, Lucidine B, Zygadenine, Aspidospermidine, silychristin, 3beta-3-Hydroxy-18-lupen-21-one, Manglupenone, and 19-Noretiocholanolone. Molecular docking, molecular dynamic simulation, and MM/GBSA analysis have proved stable, rigid, compact, and folded form of complexes during the entire 100 ns simulation, illustrating Zygadenine, Solanocapsine, and Cyclovirobuxine C as the superior inhibitors of α-A protein, while Zygadenine, Plumieride, and Phlegmarine exhibited greater inhibitory behaviour towards α-G protein than the FDA-approved drug acarbose. Collectively, our findings indicate that the Hemidesmus indicus could be a promising source of α-A and α-G inhibitors, potentially serving as a lead in order to develop medications for type-2 diabetes.

QSAR-based virtual screening of traditional Chinese medicine for the identification of mitotic kinesin Eg5 inhibitors

Cell division is a crucial process for the growth and development of all living organisms. Unfortunately, uncontrolled cell division and growth is a hallmark of cancer, leading to the formation of tumors. The Human Eg5 protein, also known as the mitotic kinesin Eg5, plays a vital role in the regulation of cell division and its dysfunction has been linked to cancer development. This study aimed to identify new inhibitors of the Human Eg5 protein. Over 2000 Traditional Chinese Medicine (TCM) compounds were screened through a combination of virtual and structure-based screening methods. The top five compounds (Compounds 1-5) showed improved binding affinity to Human Eg5 compared to the standard drug Monastrol, as demonstrated by docking and MMGBSA scores, as well as interactions with key amino acids GLY 116 and GLY 118. The potential absorption and bioactivity of these compounds were also predicted through ADMET properties and a QSAR model, respectively, and showed improved results compared to the standard. Further quantum mechanics docking confirmed the better binding affinity of the lead compound, Compound 1. Our findings highlight Compound 1-5 as promising hits for inhibiting Human Eg5 and the need for experimental validation of their potential in treating cancer.

Molecular dynamics, MMGBSA, and docking studies of natural products conjugated to tumor-targeted peptide for targeting BRAF V600E and MERTK receptors

Recent studies have revealed that MERTK and BRAF V600E receptors have been found to be over-expressed in several types of cancers including melanoma, making these receptors targets for drug design. In this study, we have designed novel peptide conjugates with the natural products vanillic acid, thiazole-2-carboxylic acid, cinnamic acid, theanine, and protocatechuic acid. Each of these compounds was conjugated with the tumor targeting peptide sequence TAASGVRSMH, known to bind to NG2 and target tumor neovasculature. We examined their binding affinities and stability with MERTK and BRAF V600E receptors using molecular docking and molecular dynamics studies. Compared to the neat compounds, the peptide conjugates displayed higher binding affinity toward both receptors. In the case of MERTK, the most stable complexes were formed with di-theaninate-peptide, vanillate-peptide, and thiazole-2-amido peptide conjugates and binding occurred in the hinge region. Additionally, it was discovered that the peptide alone also had high binding ability and stability with the MERTK receptor. In the case of BRAF V600E, the peptide conjugates of protocatechuate, vanillate and thiazole-2-amido peptide conjugates showed the formation of the most stable complexes and binding occurred in the ATP binding cleft. Further analysis revealed that the number of hydrogen bonds and hydrophobic interactions played a critical role in enhanced stability of the complexes. Docking studies also revealed that binding affinities for NG2 were similar to MERTK and higher for BRAF V600E. MMGBSA studies of the trajectories revealed that the protocatechuate-peptide conjugate showed the highest binding energy with BRAF V600E while the peptide-TAASGVRSMH showed the highest binding energy with MERTK. ADME studies revealed that each of the compounds showed medium to high permeability toward MDCK cells and were not hERG blockers. Furthermore, the conjugates were not CYP inhibitors or substrates, but they were found to be Pgp substrates. Our results indicated that the protocatechuate-TAASGVRSMH, thiazole-2-amido-TAASGVRSMH, and vanillate-TAASGVRSMH conjugates may be furthered developed for in vitro and in vivo studies as novel tumor targeting compounds for tumor cells over-expressing BRAF V600E, while di-theaninate-amido-TAASGVRSMH and thiazole-2-amido-TAASGVRSMH conjugates may be developed for targeting MERTK receptors. These studies provide insight into the molecular interactions of natural product-peptide conjugates and their potential for binding to and targeting MERTK and BRAF V600E receptors in developing new therapeutics for targeting cancer.

Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods

Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.

Drug repurposing combined with MM/PBSA based validation strategies towards MEK inhibitors screening

Emergence of oncogenic mutations in the MAPK pathway gaining more impact in the recent years. Importantly, MEK is a core element of this pathway as it is easy to inhibit and is a gatekeeper of multiple malignancies. Therefore, we performed in-silico strategy to screen repurposed candidate for MEK protein using a library of 11,808 compounds from different clusters in the DrugBank database. Glide docking, Prime-MM/GBSA and QikProp analysis were implemented to retrieve the hits with high precision. The stability of the binding mode and binding affinity of the resultant hit were explored using molecular dynamic simulations and MM/PBSA approach. The results highlight that Nebivolol (DB04861) not only achieved a stable conformation in the MEK binding pocket but also displayed highest binding affinity than the other molecules investigated in our study. Taken together, we hypothesized that Nebivolol is an excellent candidate for the inhibition of MEK in NSCLC patients in future.Communicated by Ramaswamy H. Sarma.

Design of short peptides and peptide amphiphiles as collagen mimics and an investigation of their interactions with collagen using molecular dynamics simulations and docking studies

Short peptide sequences and bolaamphiphiles derived from natural proteins are gaining importance due to their ability to form unique nanoscale architectures for a variety of biological applications. In this work, we have designed six short peptides (triplet or monomeric forms) and two peptide bolaamphiphiles that either incorporate the bioactive collagen motif (Gly-X-Y) or sequences where Gly, Pro, or hydroxyproline (Hyp) are replaced by Ala or His. For the bolaamphiphiles, a malate moiety was used as the aliphatic linker for connecting His with Hyp to create collagen mimics. Stability of the assemblies was assessed through molecular dynamics simulations and results indicated that (Pro-Ala-His)₃ and (Ala-His-Hyp)₃ formed the most stable structures, while the amphiphiles and the monomers showed some disintegration over the course of the 200 ns simulation, though most regained structural integrity and formed fibrillar structures, and micelles by the end of the simulation, likely due to the formation of more thermodynamically stable conformations. Multiple replica simulations (REMD) were also conducted where the sequences were simulated at different temperatures. Our results showed excellent convergence in most cases compared to constant temperature molecular dynamics simulation. Furthermore, molecular docking and MD simulations of the sequences bound to collagen triple helix structure revealed that several of the sequences had a high binding affinity and formed stable complexes, particularly (Pro-Ala-His)₃ and (Ala-His-Hyp)₃. Thus, we have designed new hybrid-peptide-based sequences which may be developed for potential applications as biomaterials for tissue engineering or drug delivery.

Exploration of natural product database for the identification of potent inhibitor against IDH2 mutational variants for glioma therapy

Mutation in isocitrate dehydrogenase 2 (mIDH2) is an oncogenic driver prevalently reported in various cancer types including gliomas. To date, enasidenib is the only FDA-approved drug widely used as a mIDH2 (R140Q) inhibitor. However, dose-limiting toxicity and modest brain penetrating capability restrict its use as a plausible mIDH2 inhibitor. Furthermore, secondary site mutations (Q316E and I319M) were identified in patients with enasidenib treatments resulting in acquired therapeutic resistance. Hence, in the present investigation, we aimed to identify novel and potent drug-like compounds to overcome the existing drawbacks using an integrated in-silico strategy. A sum of 1574 natural compounds from the naturally occurring plant-based anti-cancerous compound activity target (NPACT) database was proclaimed and subjected to molecular docking. The binding affinities of the resultant natural compounds were rescored using MM-GBSA scoring functions. The resultant lead molecules were subjected to anticancer activity prediction using the machine-learning model. Furthermore, the toxicity and drug-likeliness of the lead compounds were investigated using ADMET properties. Eventually, the integrated in silico approach resulted in a lead molecule, namely squalene (NPACT00954) against mIDH2 protein. The screened compound was subjected to mutational analysis accomplishing second-site mutations. Interestingly, squalene exhibited appreciable binding affinity alongside good brain penetrating potential than enasidenib. Indeed, the reproducibility and significance of our results are examined by running 3 replicas of 100-ns simulations per system using the random initial velocities of the atoms generated by Maxwell distribution at a given temperature. Thus, we hypothesize from our results that further optimization of squalene could be beneficial for the treatment and management of glioma in the near future.

A combination strategy of structure-based virtual screening, MM-GBSA, cross docking, molecular dynamics and metadynamics simulations used to investigate natural compounds as potent and specific inhibitors of tumor linked human carbonic anhydrase IX

Cancer remains a serious health concern representing one of the leading causes of deaths worldwide. The enzyme human carbonic anhydrase IX (hCA IX) is found to be over-expressed in many cancer types and its selective inhibition over its cytosolic off-target isoform, human carbonic anhydrase II (hCA II), represents a potential area of research in the development of novel anticancer compounds. This work is concerned with the use of various in silico tools for the identification of natural product based molecules that can selectively inhibit hCA IX over hCA II. MM-GBSA assisted structure-based virtual screening against hCA IX was performed for nearly 225,000 natural products imported from the ZINC15 database. The obtained hits were checked for their potency by considering acetazolamide, the bound inhibitor of hCA IX, as the reference molecule, and 121 molecules were identified as potent hCA IX inhibitors. After ensuring their potency, cross-docking, followed by MM-GBSA calculations of the hits with hCA II, was performed, and their selectivity was assessed by considering the hCA IX selective compound SLC-0111 as the reference molecule, and 50 natural products were identified as potent as well as selective hCA IX inhibitors. Molecules with the quinoline scaffold showed the highest selectivity, and their selectivity was attributed to the strong electrostatic interactions of the zinc binding group (ZBG) with the active site Zn(II) ion. Furthermore, the stability of the binding modes of the top hCA IX selective hits was ensured by performing molecular dynamics (MD) simulations, which clearly proved that one of the short-listed molecules is truly selective, as it does not interact with the active site Zn(II) ion of hCA II, but interacts strongly with this ion in hCA IX. Bonding pose metadynamics studies revealed that the ligand moves to a more stable binding site from the one predicted by the docking studies and shows stronger interaction with the protein and Zn(II) at this binding site. The ligand is not likely to have issues with bioavailability. As a result, this ligand can be taken for bioassay testing and subsequently used as a feasible therapeutic treatment for a variety of cancer types. Communicated by Ramaswamy H. Sarma.

Design, Synthesis, Molecular Docking Analysis and Biological Evaluations of 4-[(Quinolin-4-yl)amino]benzamide Derivatives as Novel Anti-Influenza Virus Agents

In this study, a series of 4-[(quinolin-4-yl)amino]benzamide derivatives as the novel anti-influenza agents were designed and synthesized. Cytotoxicity assay, cytopathic effect assay and plaque inhibition assay were performed to evaluate the anti-influenza virus A/WSN/33 (H1N1) activity of the target compounds. The target compound G07 demonstrated significant anti-influenza virus A/WSN/33 (H1N1) activity both in cytopathic effect assay (EC₅₀ = 11.38 ± 1.89 µM) and plaque inhibition assay (IC₅₀ = 0.23 ± 0.15 µM). G07 also exhibited significant anti-influenza virus activities against other three different influenza virus strains A/PR/8 (H1N1), A/HK/68 (H3N2) and influenza B virus. According to the result of ribonucleoprotein reconstitution assay, G07 could interact well with ribonucleoprotein with an inhibition rate of 80.65% at 100 µM. Furthermore, G07 exhibited significant activity target PA−PB1 subunit of RNA polymerase according to the PA−PB1 inhibitory activity prediction by the best pharmacophore Hypo1. In addition, G07 was well drug-likeness based on the results of Lipinski’s rule and ADMET prediction. All the results proved that 4-[(quinolin-4-yl)amino]benzamide derivatives could generate potential candidates in discovery of anti-influenza virus agents.

Repurposing of Four Drugs as Anti-SARS-CoV-2 Agents and Their Interactions with Protein Targets

Although there are existing vaccines against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), new COVID-19 cases are increasing due to low immunization coverage and the emergence of new variants. For this reason, new drugs to treat and prevent severe COVID-19 are needed. Here, we provide four different FDA-approved drugs against SARS-CoV-2 proteins involved in the entry and replication process, aiming to identify potential drugs to treat COVID-19. We use the main protease (Mpro), the spike glycoprotein (S protein), and RNA-dependent RNA polymerase (RdRp) as protein targets for anti- SARS-CoV-2 drugs. In our constructed database, we selected different drugs against each target (Mpro, S protein, and RdRp) based on their common interactions with relevant residues involved in viral entry at the host cell and replication. Furthermore, their stability inside the binding pocket, as well as their predicted binding-free energy, allow us to provide new insight into the possible drug repurposing of viomycin (interacting with Mpro) due to its interactions with key residues, such as Asn 143, Glu 166, and Gln 189 at the same time as hesperidin (interacting with the S protein) is interacting with residues Tyr 449, Ser 494, and Thr 500, keeping inside the predicted binding pocket, as well as interacting with residues in different variants of concern. Finally, we also suggest nystatin and elvitegravir (interacting with RdRp) as possible drugs due to their stability within the predicted pocket along the simulation and their interaction with key residues, such as Asp 760, Asp 761, and Asp 618. Altogether our results provide new knowledge about the possible mechanism of the inhibition of viomycin, hesperidin, elvitegravir, and nystatin to inhibit the viral life cycle of SARS-CoV-2 and some of its variants of concern (VOC). Additionally, some iodide-based contrast agents were also found to bind the S protein strongly, i.e., iohexol (−58.99 Kcal/mol), iotrolan (−76.19 Kcal/mol), and ioxilan (−62.37 Kcal/mol). Despite the information we report here as the possible strong interaction between these contrast agents and the SARS-CoV-2′s S protein, Mpro, and RdRp, we believe that further investigation, including chemical modifications in their structures, are needed for COVID-19 treatment.

Molecular dynamics simulations, docking and MMGBSA studies of newly designed peptide-conjugated glucosyloxy stilbene derivatives with tumor cell receptors

In this work, for the first time, we designed derivatives of beta-D-glucosyloxy-3-hydroxy-trans-stiblene-2-carboxylic acid (GHS), by conjugating GHS with tumor targeting peptides RPARPAR and GGKRPAR to target over-expressed receptors in tumor cells. The sequences RPARPAR and GGKRPAR are known to target the neuropilin1 (NRP1) receptor due to the C-terminal Arg domain; however, their effectiveness has never been examined with other commonly over-expressed receptors in tumor cells, particularly of chronic lymphocytic leukemia that include integrin α1β1 and CD22. By conjugating these peptides with GHS, which is known for its inherent anti-cancer properties, the goal is to further enhance tumor cell targeting by developing compounds that can target multiple receptors. The physicochemical properties of the conjugates and individual peptides were analyzed using Turbomole and COSMOthermX20 in order to determine their hydrogen bond accepting and donating capabilities. The web server POCASA was used in order to determine the surface cavities and binding pockets of the three receptors. To explore the binding affinities, we conducted molecular docking studies with the peptides and the conjugates with each of the receptors. After molecular docking, the complexes were analyzed using Protein-Ligand Interaction Profiler to determine the types of interactions involved. Molecular dynamics simulation studies were conducted to explore the stability of the receptor-ligand complexes. Our results indicated that in most cases the conjugates showed higher binding and stability with the receptors. Additionally, highly stable complexes of conjugates were obtained with CD22, NRP1 and in most cases with the integrin α1β1 receptor as well. The binding energies were calculated for each of the receptor ligand complexes through trajectory analysis using MMGBSA studies. SwissADME studies revealed that the compounds showed low GI absorption and were not found to be CYP inhibitors and had bioavailability score that would allow them to be considered as potential drug candidates. Overall, our results for the first time show that the designed conjugates can target multiple over-expressed receptors in tumor cells and may be potentially developed as future therapeutics for targeting tumor cells.

In Silico Studies of Tumor Targeted Peptide-Conjugated Natural Products for Targeting Over-Expressed Receptors in Breast Cancer Cells Using Molecular Docking, Molecular Dynamics and MMGBSA Calculations

In this work, in silico studies were carried out for the design of diterpene and polyphenol-peptide conjugates to potentially target over-expressed breast tumor cell receptors. Four point mutations were induced into the known tumor-targeting peptide sequence YHWYGYTPQN at positions 1, 2, 8 and 10, resulting in four mutated peptides. Each peptide was separately conjugated with either chlorogenate, carnosate, gallate, or rosmarinate given their known anti-tumor activities, creating dual targeting compounds. Molecular docking studies were conducted with the epidermal growth factor receptor (EGFR), to which the original peptide sequence is known to bind, as well as the estrogen receptor (ERα) and peroxisome proliferator-activated receptor (PPARα) using both Autodock Vina and FireDock. Based on docking results, peptide conjugates and peptides were selected and subjected to molecular dynamics simulations. MMGBSA calculations were used to further probe the binding energies. ADME studies revealed that the compounds were not CYP substrates, though most were Pgp substrates. Additionally, most of the peptides and conjugates showed MDCK permeability. Our results indicated that several of the peptide conjugates enhanced binding interactions with the receptors and resulted in stable receptor-ligand complexes; Furthermore, they may successfully target ERα and PPARα in addition to EGFR and may be further explored for synthesis and biological studies for therapeutic applications.

Discovery of Bis-sydnone styryl ketone as a selective cyclooxygenase-2 inhibitor

Background

Various literature sources have documented a wide spectrum of therapeutic properties of sydnones including anti-inflammatory, anticancer, antimicrobial activities. Phenyl styryl ketones and their derivatives as members of the chalcone family have also been reported as significant bioactive molecules. The current study was initiated to evaluate the anti-inflammatory activity of sydnone-based compounds including some novel bis-sydnone styryl ketone hybrids.

Results

Twenty-five sydnone-containing compounds were successfully synthesized. Compounds 46-48 and 56-58 were reported as new sydnone derivatives. Whereas, compounds 61-63 were synthesized as novel molecules containing two sydnone rings linked via α,β-unsaturated ketone. The structures of the synthesized compounds were confirmed by FTIR, 1H NMR, 13C NMR and ToF-MS analyses. The in vitro COX inhibition assay showed varied activity. Compounds 47, 51, 58 and 63 showed the most potent COX inhibitory effects at a concentration of 200 μM. The selectivity index revealed that 63 was the best selective COX-2 inhibitor. Acetylation of the sydnone ring at C-4 was fruitful for the COX inhibitory effects. Docking analysis showed that COX-2 selectivity was due to a favourable positive charged interaction occurring between the sydnone ring of 63 and Arg513 of COX-2. Compound 51 was hydrogen bonded to Arg513. On the other hand, the low inhibitory effect of 63 against COX-1 was due to an unfavourable polar interaction with His513 in the binding pocket of COX-1.

Conclusions

The compounds were successfully synthesized and characterized. Compound 63 had a common architecture and pharmacophoric features with known selective COX-2 inhibitors (the coxib family) which make it a suitable candidate for the designing of selective and safe NSAID.

Antibody mix-and-read assays based on fluorescence intensity probes

Antibodies and Fc fusion proteins are a rapidly growing class of pharmaceuticals. Cell culture and purification process development and operation require frequent measurement of product concentrations, commonly by complex enzyme-linked immunosorbent assay and high-performance liquid chromatography methods. Here we report a fast (<30 s), and simple antibody Fc assay based on mix-and-read reporting by fluorescence emission. A soluble fluorescein-labeled Fc-affinity reporter produced by standard peptide synthesis is mixed with an Fc-containing sample to produce an immediate shift in both fluorescence polarization and intensity, compatible with on- and at-line measurements and microbioreactor monitoring. We observed significant shifts in fluorescence intensity in Chinese hamster ovary cell culture fluid spiked with IgG and detected an adalimumab biosimilar down to 100 ng/mL (10^-4 g/L), despite the interferents in the complex sample matrix. Neither the fluorescence polarization nor the fluorescence intensity assay is significantly affected by the addition of clarified lysate of 2 million CHO-k1 cells/mL, suggesting applicability even to cultures of low viability. Biochemical and molecular docking approaches suggest that the fluorescence intensity enhancement is caused by changes in the fluorophore's local microenvironment upon binding to IgG Fc, especially by interactions with Fc His433.Abbreviations: CCF: Cell Culture Fluid; CHO: Chinese Hamster Ovary cells; ELISA: Enzyme Linked Immunosorbent Assay; Fc: Fragment Crystallizable of antibody; HPLC: High-Performance Liquid Chromatography; HPβCD: hydroxypropyl-β-cyclodextrin; IgG: ImmunoglobulinG; mAb: Monoclonal Antibody; PBS: Phosphate-Buffered Saline; PDB: Protein Data Bank; SpA: Staphylococcal protein A; SpG: Staphylococcal protein G.

Multifaceted targeting strategies in cancer against the human notch 3 protein: a computational study

Notch receptors play a significant role in the development and the regulation of cell-fate in several multicellular organisms. For normal differentiation, genomes are essential as their regular roles and play a role in cancer is dysregulated. Notch 3 has been shown to play a major role in lung cancer function and therefore, inhibition of notch 3 protein activation represents a clear plan for cancer treatment. This study accomplished a combined structure- and ligand-based pharmacophore hypothesis to explore novel notch 3 inhibitors. The analysis identified common lead molecule ZINC000013449462 that showed better XP GScore and binding energy score than the reference inhibitor DAPT. The identified lead compound that passed all the druggable characteristics exhibited stable binding. Furthermore, the lead molecule can also form hydrogen and salt bridge interactions with binding site residues Asp1621 and Arg1465 residues, respectively of the active pockets of notch 3 protein. In essence, the inhibitory activity of the hit was validated across 109 NSCLC cell lines by employing a deep neural network algorithm. Our study proposes that ZINC000013449462 would be a possible prototype molecule towards the notch 3 target and further examined by clinical studies to combat NSCLC.

Synthesis of New 3-Arylaminophthalides and 3-Indolyl-phthalides using Ammonium Chloride, Evaluation of their Anti-Mycobacterial Potential and Docking Study

OBJECTIVE

The study aims at the derivatization of "Phthalides" and synthesizes 3- arylaminophthalides & 3-indolyl-phthalides compounds, and evaluates their anti-tubercular and antioxidant activities. The study has also intended to employ the in silico methods for the identification of possible drug targets in Mycobacterium and evaluate the binding affinities of synthesized compounds.

METHODS

This report briefly explains the synthesis of phthalide derivatives using ammonium chloride. The synthesized compounds were characterized using spectral analysis. Resazurin Microtiter Assay (REMA) plate method was used to demonstrate the anti-mycobacterial activity of the synthesized compounds. An in-silico pharmacophore probing approach was used for target identification in Mycobacterium. The structural level interaction between the identified putative drug target and synthesized phthalides was studied using Lamarckian genetic algorithm-based software.

RESULTS AND DISCUSSION

In the present study, we report an effective, environmentally benign scheme for the synthesis of phthalide derivatives. Compounds 5c and 5d from the current series appear to possess good anti-mycobacterial activity. dCTP: deaminasedUTPase was identified as a putative drug target in Mycobacterium. The docking results clearly showed the interactive involvement of conserved residues of dCTP with the synthesized phthalide compounds.

CONCLUSION

On the eve of evolving anti-TB drug resistance, the data on anti-tubercular and allied activities of the compounds in the present study demonstrates the enormous significance of these newly synthesized derivatives as possible candidate leads in the development of novel anti-tubercular agents. The docking results from the current report provide a structural rationale for the promising anti-tubercular activity demonstrated by 3-arylaminophthalides and 3-indolyl-phthalides compounds.

Drug repurposing for ligand-induced rearrangement of Sirt2 active site-based inhibitors via molecular modeling and quantum mechanics calculations

Sirtuin 2 (Sirt2) nicotinamide adenine dinucleotide-dependent deacetylase enzyme has been reported to alter diverse biological functions in the cells and onset of diseases, including cancer, aging, and neurodegenerative diseases, which implicate the regulation of Sirt2 function as a potential drug target. Available Sirt2 inhibitors or modulators exhibit insufficient specificity and potency, and even partially contradictory Sirt2 effects were described for the available inhibitors. Herein, we applied computational screening and evaluation of FDA-approved drugs for highly selective modulation of Sirt2 activity via a unique inhibitory mechanism as reported earlier for SirReal2 inhibitor. Application of stringent molecular docking results in the identification of 48 FDA-approved drugs as selective putative inhibitors of Sirt2, but only top 10 drugs with docking scores > - 11 kcal/mol were considered in reference to SirReal2 inhibitor for computational analysis. The molecular dynamics simulations and post-simulation analysis of Sirt2-drug complexes revealed substantial stability for Fluphenazine and Nintedanib with Sirt2. Additionally, developed 3D-QSAR-models also support the inhibitory potential of drugs, which exclusively revealed highest activities for Nintedanib (pIC50 ≥ 5.90 µM). Conclusively, screened FDA-approved drugs were advocated as promising agents for Sirt2 inhibition and required in vitro investigation for Sirt2 targeted drug development.

Pharmacophore model, docking, QSAR, and molecular dynamics simulation studies of substituted cyclic imides and herbal medicines as COX-2 inhibitors

Cyclooxygenase-2 (COX-2) enzyme inhibitors have not eliminated the necessity for developed drugs not only in the nonsteroidal anti-inflammatory drug (NSAIDs) area, but also in other therapeutic applications including prevention of cancer and Alzheimer's disease. A series of novel substituted cyclic imides have been reported as selective COX-2 inhibitors. To understand the structural features responsible for their activity, a 3D validated pharmacophore and quantitative structure−activity relationship (QSAR) model have been developed. The values of enrichment factor (EF), goodness of hit score (GH), area under the ROC curve (AUC), sensitivity, and specificity refer to the good ability of the pharmacophore model to identify active compounds. Multiple linear regression (MLR) produced statistically significant QSAR model with (R²_training = 0.763, R²_test = 0.96) and predictability (Q²_training = 0.66, Q²_test = 0.84). Then, using the pharmacophore and QSAR models, eight authenticated botanicals in two herbal medicines and the ZINC compounds database, were virtually screened for ligands to COX-2. The retrieved hits which also obey lipinski's rule of five (RO5) were docked in the COX-2 3D structure to investigate their binding mode and affinity. Finally, based on the docking results, nine molecules were prioritized as promising hits that could be used as leads to discover novel COX-2 inhibitors. COX-2 inhibition of most of these hits has not been reported previously. Ten-nanosecond molecular dynamics simulation (10-ns MD) was performed on the initial structure COX-2 complex with ZINC000113253375 and ZINC000043170560 resulted from the docking. Our utilization of the 3D pharmacophore model, QSAR, molecular docking, and molecular dynamics simulation trials can be a potent strategy to successfully predict activity, efficiently design drugs, and screen large numbers of new compounds as active drug candidates.

New Benzoxazole Derivatives as Antiprotozoal Agents: In Silico Studies, Synthesis, and Biological Evaluation

Background. Benzoxazole derivatives have different biological activities. In pursuit of designing novel chemical entities with antiprotozoal and antimicrobial activities, benzoxazolyl aniline was utilized as a privileged scaffold of a series of (3-benzoxazole-2-yl) phenylamine derivatives, 3-benzoxazoloyl acetamide, and butyramide derivatives. Methods. These novel analogs were synthesized in straightforward simple chemistry without any quantitative chromatographic separations in reasonable yields. The biological evaluation of all target compounds as potential antimalarial, antileishmanial, antitrypanosomal, and antimicrobial agents was performed by various well-established cell-based methods. Results. Compounds 6d and 5a showed promising biological screening data. The amidation of 3-benzoxazolyl aniline 1 with the chloroacetyl functional group resulted in a good antimalarial activity and showed moderate inhibitory activities against leishmanial and trypanosomal spp. Moreover, chloroacetyl functionalization of benzoxazolyl aniline serves as a good early goal for constructing and synthesizing new antimicrobial and antiprotozoal agents. The molecular docking study rationalizes the relative inhibitory activity of compound 5a as an antimalarial agent with the deregulation of PfPNP activity which has emerged as a major mechanism of these targets.

Targeting the Autophagy Specific Lipid Kinase VPS34 for Cancer Treatment: An Integrative Repurposing Strategy

The impact of autophagy on cancer treatment and its corresponding responsiveness has galvanized the scientific community to develop novel inhibitors for cancer treatment. Importantly, the discovery of inhibitors that targets the early phase of autophagy was identified as a beneficial choice. Despite the number of research in recent years, screening of the DrugBank repository (9591 molecules) for the Vacuolar protein sorting 34 (VPS34) has not been reported earlier. Therefore, the present study was designed to identify potential VPS34 antagonists using integrated pharmacophore strategies. Primarily, an energy-based pharmacophore and receptor cavity-based analysis yielded five (DHRRR) and seven featured (AADDHRR) pharmacophore hypotheses respectively, which were utilized for the database screening process. The glide score, the binding free energy, pharmacokinetics and pharmacodynamics properties were examined to narrow down the screened compounds. This analysis yielded a hit molecule, DB03916 that exhibited a better docking score, higher binding affinity and better drug-like properties in contrast to the reference compound that suffers from a toxicity property. Importantly, the result was validated using a 50 ns molecular dynamics simulation study. Overall, we conclude that the identified hit molecule DB03916 is believed to serve as a prospective antagonist against VPS34 for cancer treatment.

Identification of a repurposed drug as an inhibitor of Spike protein of human coronavirus SARS-CoV-2 by computational methods

Severe acute respiratory syndrome coronavirus (SARS-CoV-2) is an emerging new viral pathogen that causes severe respiratory disease. SARS-CoV-2 is responsible for the outbreak of COVID-19 pandemic worldwide. As there are no confirmed antiviral drugs or vaccines currently available for the treatment of COVID-19, discovering potent inhibitors or vaccines are urgently required for the benefit of humanity. The glycosylated Spike protein (S-protein) directly interacts with human angiotensin-converting enzyme 2 (ACE2) receptor through the receptor-binding domain (RBD) of S-protein. As the S-protein is exposed to the surface and is essential for entry into the host, the S-protein can be considered as a first-line therapeutic target for antiviral therapy and vaccine development. In silico screening, docking, and molecular dynamics simulation studies were performed to identify repurposing drugs using DrugBank and PubChem library against the RBD of S-protein. The study identified a laxative drug, Bisoxatin (DB09219), which is used for the treatment of constipation and preparation of the colon for surgical procedures. It binds nicely at the S-protein–ACE2 interface by making substantial π-π interactions with Tyr505 in the ‘Site 1’ hook region of RBD and hydrophilic interactions with Glu406, Ser494, and Thr500. Bisoxatin consistently binds to the protein throughout the 100 ns simulation. Taken together, we propose that the discovered molecule, Bisoxatin may be a promising repurposable drug molecule to develop new chemical libraries for inhibiting SARS-CoV-2 entry into the host.

TRAF2 and NCK-Interacting Kinase Inhibitors for Colorectal Cancer: In Vitro and Theoretical Validations

TRAF2 and NCK-interacting kinase (TNIK) is a critical factor in colorectal cancer (CRC) proliferation mediated by Wnt signaling. We attempted to identify efficient TNIK inhibitors using computational high-throughput virtual screening (HTVS) from various drug banks and databases. By performing/on performing e-pharmacophore screening and molecular docking methods, from ∼700 000 molecules, compounds LC_222150, LC_112060, and LC_64796 were identified as potential leads, through molecular dynamics (MD) simulations and density functional theory (DFT). These top 3 structures were commercially procured, and their inhibitory activity was assessed in vitro. Significant TNIK inhibition was observed, with an average IC₅₀ of 18.33 ± 0.75 nM. In terms of anticancer activity, the observed average relative % activity (RPA) of 90.28 ± 1.04 for these compounds compared well with doxorubicin (86.75 ± 1.45) as a standard. Compounds LC_222150, LC_112060, and LC_64796, therefore, warrant further evaluation in vivo to assess their CRC therapeutic effects.

Scaffold morphing of arbidol (umifenovir) in search of multi-targeting therapy halting the interaction of SARS-CoV-2 with ACE2 and other proteases involved in COVID-19

The rapid emergence of novel coronavirus, SARS-coronavirus 2 (SARS-CoV-2), originated from Wuhan, China, imposed a global health emergency. Angiotensin-converting enzyme 2 (ACE2) receptor serves as an entry point for this deadly virus while the proteases like furin, transmembrane protease serine 2 (TMPRSS2) and 3 chymotrypsin-like protease (3CLpro) are involved in the further processing and replication of SARS-CoV-2. The interaction of SP with ACE2 and these proteases results in the SARS-CoV-2 invasion and fast epidemic spread. The small molecular inhibitors are reported to limit the interaction of SP with ACE2 and other proteases. Arbidol, a membrane fusion inhibitor approved for influenza virus is currently undergoing clinical trials against COVID-19. In this context, we report some analogues of arbidol designed by scaffold morphing and structure-based designing approaches with a superior therapeutic profile. The representative compounds A_BR4, A_BR9, A_BR18, A_BR22 and A_BR28 restricted the interaction of SARS-CoV-2 SP with ACE2 and host proteases furin and TMPRSS2. For 3CLPro, Compounds A_BR5, A_BR6, A_BR9 and A_BR18 exhibited high binding affinity, docking score and key residue interactions. Overall, A_BR18 and A_BR28 demonstrated multi-targeting potential against all the targets. Among these top-scoring molecules A_BR9, A_BR18, A_BR22 and A_BR28 were predicted to confer favorable ADME properties.

Nano-ellagic acid: inhibitory actions on aldose reductase and α-glucosidase in secondary complications of diabetes, strengthened by in silico docking studies

Increased blood sugar levels in prolonged diabetes lead to secondary complications such as retinopathy, neuropathy, and nephropathy, which gradually end in death. Synthesis of nano-phytomedicines from active phytoconstituents for novel emerging applications in the field of pharmaceuticals is of huge interest among researchers. In the present investigation, encapsulated ellagic acid (NEA) was synthesized at four different concentrations (0.2%, 0.3%, 0.4%, 0.5%) using ZnO nanoparticles as encapsulating agent. The surface morphology (fiber-like structures) of the nanoparticles were determined by scanning electron microscopy (SEM) and particle size (161-297 nm) and zeta potential (- 54.9-38.4 mV) were determined by dynamic light scattering technique. Further, the α-glucosidase and aldose reductase enzymes were significantly inhibited by the 0.4% of NEA compared to the other concentrations which strengthened our studies in overcoming the secondary complications of diabetes. The interaction analysis between ellagic acid and insulin receptor found Hit 1 among 10 executed ∆G score and energy of - 5.76, - 4.63 kcal/mol and formed polar bond with Arg 113 with - 1.175 Å distance. The residues Arg115, Lys116, Phe118, Ile115, Arg1131, Arg1155, Ile1157, Lys1165 and Phe1186 were found in ligand-protein interactions. ADME/T analysis of hit 1 was within the acceptable range without any toxic functional groups, providing a framework for developing novel therapeutics.

In silico molecular docking studies and MM/GBSA analysis of coumarin-carbonodithioate hybrid derivatives divulge the anticancer potential against breast cancer

Background

There are many biomarkers associated with breast cancer. Higher expression of PIK3CA (Phosphoinositide 3-kinase Cα), in its upregulated form, is associated with Hr+ and Her2− breast cancer; therefore, many drugs were synthesized against this protein to treat breast cancer patients. FDA recently approved that the drug alpelisib also inhibits PI3KCα (PDB ID-5DXT) in BC patients with Hr+ and Her2−. In present study, we have exploited fourteen coumarin-carbonodithioate derivatives and alpelisib against this protein along with eighteen others which are responsible for causing BC through computational analysis. We have used Schrödinger Maestro 11.2 version for our in silico docking study, and to calculate relative binding energies of ligands, we used prime MM-GBSA module.

Result

Docking study revealed that among all fourteen compounds, 2f, 2a, 2d, and 2e showed the highest G score than the alpelisib and coumarin against PI3KCα with − 9.3, − 9.0, − 9.0 and − 9.1 kcal/mol respectively, along with individual G score of alpelisib (− 8.9) and coumarin (− 7.9). Prime MM-GBSA analysis gave the relative binding energies of alpelisib, 2f, and 2e with − 19.94864535, − 18.63076296 and − 13.07341286 kcal/mol sequentially.

Conclusion

This study provides an insight into the coumarin-carbonodithioate derivatives that could act as inhibitors of PI3KCα like alpelisib. Further prime MM-GBSA study revealed ligand binding energies and ligands strain energies.

Modeling and simulation study to identify threonine synthase as possible drug target in Leishmania major

Leishmaniasis is one of the most neglected tropical diseases that demand immediate attention to the identification of new drug targets and effective drug candidates. The present study demonstrates the possibility of using threonine synthase (TS) as a putative drug target in leishmaniasis disease management. We report the construction of an effective homology model of the enzyme that appears to be structurally as well as functionally well conserved. The 200 nanosecond molecular dynamics data on TS with and without pyridoxal phosphate (PLP) shed light on mechanistic details of PLP-induced conformational changes. Moreover, we address some important structural and dynamic interactions in the PLP binding region of TS that are in good agreement with previously speculated crystallographic estimations. Additionally, after screening more than 44,000 compounds, we propose 10 putative inhibitor candidates for TS based on virtual screening data and refined Molecular Mechanics Generalized Born Surface Area calculations. We expect that structural and functional dynamics data disclosed in this study will help initiate experimental endeavors toward establishing TS as an effective antileishmanial drug target.

Designing novel inhibitors against cyclopropane mycolic acid synthase 3 (PcaA): targeting dormant state of Mycobacterium tuberculosis

Mycobacterium tuberculosis can sustain inside the host in dormant (non-replicating) state for years. It suppresses the host immune system by residing in the host alveolar macrophage, resulting in the development of latent tuberculosis. Despite many antibiotics available for the treatment of tuberculosis, the major hurdle in complete elimination is the ability of the bacilli to undergo dormancy and develop resistance against the existing drugs. Cyclopropanation of mycolic acids present in the cell wall of mycobacteria is required for its persistence and virulence. Cyclopropane synthases such as PcaA, CmaA1 and CmaA2, introduce site-specific modifications in mycolic acids. PcaA expression levels are high during dormancy and the gene mutants fails to persist, showing reduced survival in host macrophage. Hence, PcaA appears as a potential target to develop inhibitors against the dormant bacilli. In this study, we have identified compounds with maximum binding affinity against PcaA by in-silico virtual screening of anti-tuberculosis compounds and their structural analogues. In-silico docking followed molecular dynamic simulations and free energy calculations of the compounds with highest docking score in their respective libraries. This study reports novel inhibitors that can act as better anti-tuberculosis compounds targeting PcaA.Communicated by Ramaswamy H. Sarma.

Pharmacophore-enabled virtual screening, molecular docking and molecular dynamics studies for identification of potent and selective histone deacetylase 8 inhibitors

Histone deacetylases (HDACs) play important roles in various biological processes, but are also notorious for their over-expression in numerous cancers and neurological disorders. Therefore, the development of isoform selective HDAC inhibitors is crucial in order to prevent any side effects of pan inhibition. This work focuses on identifying novel inhibitors for the selective inhibition of HDAC8, an isoform implicated in fatal diseases like T-cell lymphoma, colon cancer and childhood neuroblastoma. Virtual screening of the 'In-trials' subset of ZINC database has been carried out with the help of two pharmacophore models signifying potent and selective HDAC8 inhibition. A detailed molecular docking strategy, followed by molecular dynamics simulations and post-scoring with MM-GBSA calculations, has led to the identification of six promising molecules that have excellent binding with the HDAC8 active site. In order to establish the selectivity profile of these molecules, their binding to off-target HDAC isoforms has also been evaluated. Substitution analyses of the proposed inhibitors suggest that aromatic substituents that access the adjacent hydrophobic pocket of the HDAC8 active site have the potential to further enhance the HDAC8 selectivity.

Estimation of the ligand-binding free energy of checkpoint kinase 1 via non-equilibrium MD simulations

Checkpoint kinase 1 (CHK1) is a serine/threonine-protein kinase that is involved in cell cycle regulation in eukaryotes. Inhibition of CHK1 is thus considered as a promising approach in cancer therapy. In this study, the fast pulling of ligand (FPL) process was applied to predict the relative binding affinities of CHK1 inhibitors using non-equilibrium molecular dynamics (MD) simulations. The work of external harmonic forces to pull the ligand out of the binding cavity strongly correlated with the experimental binding affinity of CHK1 inhibitors with the correlation coefficient of R = -0.88 and an overall root mean square error (RMSE) of 0.99 kcal/mol. The data indicate that the FPL method is highly accurate in predicting the relative binding free energies of CHK1 inhibitors with an affordable CPU time. A new set of molecules were designed based on the molecular modeling of interactions between the known inhibitor and CHK1 as inhibitory candidates. Molecular docking and FPL results exhibited that the binding affinities of developed ligands were similar to the known inhibitor in interaction with the catalytic site of CHK1, producing very potential CHK1 inhibitors of that the inhibitory activities should be further evaluated in vitro.

Design, synthesis, anti-tobacco mosaic viral and molecule docking simulations of urea/thiourea derivatives of 2-(piperazine-1-yl)-pyrimidine and 1-(4-Fluoro/4-Chloro phenyl)-piperazine and 1-(4-Chloro phenyl)-piperazine - A study

The objectives of the present work are to design, syhthesize and introduce novel urea/thiourea derivatives of 2-(piperazine-1-yl)-pyrimidine and 1-(4-Fluoro/4-Chloro phenyl)-piperazine molecules as tobacco mosaic virus (TMV) inhibitors. A series of urea/thiourea derivatives containing pyrimidine and piperazine moieties were synthesized, characterized using Fourier-transform infrared (FTIR) mass spectra, nuclear magnetic resonance (NMR) spectroscopy, elemental analysis and evaluated their sustainability using biological experiments. The anti-viral bioassay of the title compounds showed an antiviral activity against TMV. The compounds synthesized, 9j, 6g and 3d, showed highly-potential curative, protective, and inhibitory activities against TMV at 500 mg/mL formulation. All these compounds were allowed to quantum-polarized-ligand (quantum mechanical and molecular mechanical (QM/MM)) docking experiments. The compounds 9j, 6g and 3d structurally exhibited identical higher affinity towards TMV-Helicase and TMV-Coat proteins. The docking interactions proposed had two stage inhibition of TMV virus by binding to coat protein and helicase for inhibition of RNA replication. The long-range molecular dynamics (150 ns) simulations has revealed more consistency by 9j, 6g and 3d. The present study outcomes good binding propensity for active-tunnel of TMV-Hel enzyme, by these thiourea, urea derivatives, 9j, 6g and 3d, to suggest that the designed and synthesized were ideal for proposing as selective novel inhibitors to target for TMV.

Promising inhibitors of main protease of novel corona virus to prevent the spread of COVID-19 using docking and molecular dynamics simulation

Coronavirus disease-2019 (COVID-19) is a global health emergency and the matter of serious concern, which has been declared a pandemic by WHO. Till date, no potential medicine/ drug is available to cure the infected persons from SARS-CoV-2. This deadly virus is named as novel 2019-nCoV coronavirus and caused coronavirus disease, that is, COVID-19. The first case of SARS-CoV-2 infection in human was confirmed in the Wuhan city of the China. COVID-19 is an infectious disease and spread from man to man as well as surface to man . In the present work, in silico approach was followed to find potential molecule to control this infection. Authors have screened more than one million molecules available in the ZINC database and taken the best two compounds based on binding energy score. These lead molecules were further studied through docking against the main protease of SARS-CoV-2. Then, molecular dynamics simulations of the main protease with and without screened compounds were performed at room temperature to determine the thermodynamic parameters to understand the inhibition. Further, molecular dynamics simulations at different temperatures were performed to understand the efficiency of the inhibition of the main protease in the presence of the screened compounds. Change in energy for the formation of the complexes between the main protease of novel coronavirus and ZINC20601870 as well ZINC00793735 at room temperature was determined on applying MM-GBSA calculations. Docking and molecular dynamics simulations showed their antiviral potential and may inhibit viral replication experimentally.

Communicated by Ramaswamy H. Sarma

Potential Inhibitors for Novel Coronavirus Protease Identified by Virtual Screening of 606 Million Compounds

The rapid outbreak of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China followed by its spread around the world poses a serious global concern for public health. To this date, no specific drugs or vaccines are available to treat SARS-CoV-2 despite its close relation to the SARS-CoV virus that caused a similar epidemic in 2003. Thus, there remains an urgent need for the identification and development of specific antiviral therapeutics against SARS-CoV-2. To conquer viral infections, the inhibition of proteases essential for proteolytic processing of viral polyproteins is a conventional therapeutic strategy. In order to find novel inhibitors, we computationally screened a compound library of over 606 million compounds for binding at the recently solved crystal structure of the main protease (Mpro) of SARS-CoV-2. A screening of such a vast chemical space for SARS-CoV-2 Mpro inhibitors has not been reported before. After shape screening, two docking protocols were applied followed by the determination of molecular descriptors relevant for pharmacokinetics to narrow down the number of initial hits. Next, molecular dynamics simulations were conducted to validate the stability of docked binding modes and comprehensively quantify ligand binding energies. After evaluation of potential off-target binding, we report a list of 12 purchasable compounds, with binding affinity to the target protease that is predicted to be more favorable than that of the cocrystallized peptidomimetic compound. In order to quickly advise ongoing therapeutic intervention for patients, we evaluated approved antiviral drugs and other protease inhibitors to provide a list of nine compounds for drug repurposing. Furthermore, we identified the natural compounds (-)-taxifolin and rhamnetin as potential inhibitors of Mpro. Rhamnetin is already commercially available in pharmacies.

Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures

The current outbreak of a novel coronavirus, named as SARS-CoV-2 causing COVID-19 occurred in 2019, is in dire need of finding potential therapeutic agents. Recently, ongoing viral epidemic due to coronavirus (SARS-CoV-2) primarily affected mainland China that now threatened to spread to populations in most countries of the world. In spite of this, there is currently no antiviral drug/ vaccine available against coronavirus infection, COVID-19. In the present study, computer-aided drug design-based screening to find out promising inhibitors against the coronavirus (SARS-CoV-2) leads to infection, COVID-19. The lead therapeutic molecule was investigated through docking and molecular dynamics simulations. In this, binding affinity of noscapines(23B)-protease of SARS-CoV-2 complex was evaluated through MD simulations at different temperatures. Our research group has established that noscapine is a chemotherapeutic agent for the treatment of drug resistant cancers; however, noscapine was also being used as anti-malarial, anti-stroke and cough-suppressant. This study suggests for the first time that noscapine exerts its antiviral effects by inhibiting viral protein synthesis.

Identification of Selective Inhibitors of LdDHFR Enzyme Using Pharmacoinformatic Methods

Dihydrofolate reductase (DHFR) is a well-known enzyme of the folate metabolic pathway and it is a validated drug target for leishmaniasis. However, only a few leads are reported against Leishmania donovani DHFR (LdDHFR), and thus, there is a need to identify new inhibitors. In this article, pharmacoinformatic tools such as molecular docking, virtual screening, absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling, and molecular dynamics (MD) simulations were utilized to identify potential LdDHFR inhibitors. Initially, a natural DHFR substrate (dihydrofolate), a classical DHFR inhibitor (methotrexate), and a potent LdDHFR inhibitor, that is, "5-(3-(octyloxy)benzyl)pyrimidine-2,4-diamine" (LEAD) were docked in the active site of the LdDHFR and MD simulated to understand the binding mode characteristics of the substrates/inhibitors in the LdDHFR. The shape of the LEAD molecule was used as a query for shape-based virtual screening, while the three-dimensional structure of LdDHFR was utilized for docking-based virtual screening. In silico ADMET factors were also considered during virtual screening. These two screening processes yielded 25 suitable hits, which were further validated for their selectivity toward LdDHFR using molecular docking and prime molecular mechanics/generalized born surface area analysis in the human DHFR (HsDHFR). Best six hits, which were selective and energetically favorable for the LdDHFR, were chosen for MD simulations. The MD analysis showed that four of the hits exhibited very good binding affinity for LdDHFR with respect to HsDHFR, and two hits were found to be more selective than the reported potent LdDHFR inhibitor. The present study thus identifies hits that can be further designed and modified as potent LdDHFR inhibitors.

Selective Inhibition of HDAC1 by Macrocyclic Polypeptide for the Treatment of Glioblastoma: A Binding Mechanistic Analysis Based on Molecular Dynamics

Glioblastoma (GBM) is the most common and aggressive intracranial malignant brain tumor, and the abnormal expression of HDAC1 is closely correlated to the progression, recurrence and metastasis of GBM cells, making selective inhibition of HDAC1 a promising strategy for GBM treatments. Among all available selective HDAC1 inhibitors, the macrocyclic peptides have gained great attention due to their remarkable inhibitory selectivity on HDAC1. However, the binding mechanism underlying this selectivity is still elusive, which increases the difficulty of designing and synthesizing the macrocyclic peptide-based anti-GBM drug. Herein, multiple computational approaches were employed to explore the binding behaviors of a typical macrocyclic peptide FK228 in both HDAC1 and HDAC6. Starting from the docking conformations of FK228 in the binding pockets of HDAC1&6, relatively long MD simulation (500 ns) shown that the hydrophobic interaction and hydrogen bonding of E91 and D92 in the Loop2 of HDAC1 with the Cap had a certain traction effect on FK228, and the sub-pocket formed by Loop1 and Loop2 in HDAC1 could better accommodate the Cap group, which had a positive effect on maintaining the active conformation of FK228. While the weakening of the interactions between FK228 and the residues in the Loop2 of HDAC6 during the MD simulation led to the large deflection of FK228 in the binding site, which also resulted in the decrease in the interactions between the Linker region of FK228 and the previously identified key amino acids (H134, F143, H174, and F203). Therefore, the residues located in Loop1 and Loop2 contributed in maintaining the active conformation of FK228, which would provide valuable hints for the discovery and design of novel macrocyclic polypeptide HDAC inhibitors.