Molecular docking and molecular dynamics simulation studies to identify potent AURKA inhibitors: assessing the performance of density functional theory, MM-GBSA and mass action kinetics calculations

Synthesis and biological assessment of novel 4H-chromene-3-carbonitrile derivatives as tyrosinase inhibitors

Excessive activity of the tyrosinase enzyme during melanogenesis results in hyperpigmentation in the skin. To address this issue, there is a need to develop effective tyrosinase inhibitors as a treatment for hyperpigmentation. In this study, we synthesized some novel 4H-chromene-3-carbonitrile compounds (6a-o) and assessed their inhibitory activities against tyrosinase, comparing them with kojic acid, which is known as a positive control. Compound 6f emerged as the most effective inhibitor, with an IC50 of 35.38 ± 2.12 µM. Kinetic studies of 6f exhibited competitive inhibition, with Ki = 16.15 µM. Molecular docking studies highlighted the importance of π-π stacking and hydrogen bonding interactions within the binding site. Molecular dynamics simulations showed that the R-enantiomer 6f exhibited superior binding stability compared to the S-enantiomer, with a lower standard deviation of RMSD and more persistent interactions with the key active site residues. These findings underscore the potential of the R-enantiomer of compound 6f as a potent tyrosinase inhibitor and provide insights for developing effective treatments for hyperpigmentation and related skin conditions.

Designing a novel drug-drug conjugate as a prodrug for breast cancer therapy: in silico insights

Breast cancer (BC) poses a significant global health threat, necessitating innovative therapeutic approaches. The ribosomal s6 kinase 2 (RSK2) has emerged as a promising target due to its roles in cell proliferation and survival. This study proposes a drug-drug conjugate prodrug comprising Methotrexate (hydrophobic) and Capecitabine (hydrophilic) for BC treatment. In silico approaches, including Molecular Docking, Molecular Dynamics Simulations, MM-PBSA, ADME, and DFT calculations were employed to evaluate the prodrug's potential. The designed MET-CAP ligand exhibits a robust docking score (-8.980 kcal/mol), superior binding affinity (-53.16 kcal/mol), and stable dynamic behavior (0.62 nm) compared to native ligands. The DFT results reveal intramolecular charge transfer in MET-CAP (HLG = 0.09 eV), indicating its potential as a BC inhibitor. ADME analysis suggests satisfactory pharmaceutically relevant properties. The results indicate that the conjugated MET-CAP ligand exhibits favorable binding characteristics, stability, and pharmaceutically relevant properties, making it a potential RSK2 inhibitor for BC therapy. The multifaceted approach provides insights into binding interactions, stability, and pharmacokinetic properties, laying the foundation for further experimental validation and potential clinical development.

A discovery of potent kaempferol derivatives as multi-target medicines against diabetes as well as bacterial infections: an in silico approach

Flavonoids demonstrate beneficial effects on human health because flavonoids contain important biological properties. Kaempferol is a flavonol, type of flavonoid found in eatable plants and in plants usually employed in ancient drugs (Moringa oleifera, Tilia spp., fern genus spp. and gingko etc.). Some medicinal studies have shown that the use of foods full of kaempferol decreases the risk of many (cancer, vascular) diseases. All the data of 50 kaempferol derivatives were collected from PubChem database. Through Schrödinger software, 3D-QSAR study was performed for 50 compounds by using method of field base. Conformer of kaempferol derivatives was docked against anti-diabetic, anti-microbial co-crystal structures and protein. To monitor the best anti-diabetic and antibacterial agent, particular kaempferol derivatives were downloaded from PubChem database. Virtual screening by molecular docking provided four lead compounds with four different proteins. These hit compounds were found to be potent inhibitor for diabetic enzymes alpha-amylase and DPP IV and had the potential to suppress DNA gyrase and dihydrofolate reductase synthesis. Molecular dynamic simulation of docked complexes evaluates the value of root mean square fluctuation by iMOD server. Kaempferol 3-O-alpha-L-(2, 3-di-Z-p-coumaroyl) rhamnoside (42) compound used as anti-diabetic and kaempferol 3-O-gentiobioside (3) as antibacterial with good results can be used for drug discovery.Communicated by Ramaswamy H. Sarma.

A biochemistry-oriented drug design: synthesis, anticancer activity, enzymes inhibition, molecular docking studies of novel 1,2,4-triazole derivatives

In this study, we researched the reactions of 5-(5-bromofuran-2-yl)-4-methyl-1,2,4-triazole-3-thiol and 5-thiophene-(3-ylmethyl)-4R-1,2,4-triazole-3-thiols with some halogen-containing compounds, a number of new compounds were synthesized (1.1-1.5 and 2.1-2.8). These compounds showed excellent to good inhibitory activities on acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes. For obtaining the effects of these compounds on AChE and BChE enzymes were determined spectrophotometrically according to Ellman. IC50 values of these enzymes were ranging between 1.63 and 17.68 nM for AChE and 8.71 and 84.02 nM for BChE. After, prostate cancer is the second leading cause of cancer-related mortality for men over the age of 65 in developed countries. Current treatment options remain limited in the treatment of advanced-stage prostate cancer leading to biochemical recurrence in almost 40% of the patients. Therefore, there is an urgent need for development of novel therapeutic tools for treatment of prostate cancer patients. In this study, we aimed at analyzing the potential of all compounds against prostate cancer cells. We found that, of the tested compounds, 2.1, 2.2 and 2.3 showed significant cytotoxic activities against PC3 prostate cancer cells, although their effect on the viability of normal prostate cells was limited. These findings suggest their selective targeting potential for prostate cancer cells and offer them as candidate therapeutic agents against prostate cancer. The inhibitory activities of some chemical compounds, such as (1.1-1.5 and 2.1-2.8) were assessed by performing the molecular docking study in the presence of AChE, BChE and prostate cancer protein. MM/GBSA methods are calculated binding free energy. Finally, ADME/T analysis was performed to examine the drug properties of the 13 studied molecules.Communicated by Ramaswamy H. Sarma.

Antioxidant enzyme activities, molecular docking studies, MM-GBSA, and molecular dynamic of chlorpyrifos in freshwater fish Capoeta umbla

Chlorpyrifos (CPF), which was started to be used in 1965, is a broad spectrum organophosphate insecticide that is used more and more day by day. Commonly used to control pests in farmland and homes, CPF is more toxic to fish than organochlorine compounds. CPF poses a serious threat to the health of humans and aquatic organisms. This paper studies the relationship between CPF exposure and antioxidant enzyme activities in gill, kidney and liver tissues of Capoeta umbla. Different time intervals (12, 24, 48, 72, and 96 h) and CPF doses (55 and 110 µg L-1) were used in the study. Spectrophotometrical measures were taken in all tissues for antioxidant enzyme activities and malondialdehyde (MDA) levels, as indices of the lipid peroxidation (LPO). A positive relationship between CPF and MDA levels was found in the study at a statistically significant level (p < 0.05). The study also found a negative relationship between CPF levels and catalase (CAT), superoxide dismutase (SOD), glutathione peroxidase (GPx), and glutathione reductase (GR) activity. Independent variables in the study can act as biomarkers of CPF exposure. The study recommends employing proper ecotoxicological risk evaluations in cases of CPF usage as a pesticide. The activities of the studied molecules against various proteins that are crystal structure of human peroxiredoxin 5 (PDB ID: 1HD2) has docking score value is -2.67, crystal structure of Bovine Xanthine Oxidase (PDB ID: 3NRZ) has docking score value is -3.76, and crystal structure of antibacterial FabH (PDB ID: 4Z8D) has docking score value is -3.16, were compared. Molecular dynamic (MD) calculations were made in 100 ns. MM/GBSA methods are calculated binding free energy. Afterwards, ADME/T analysis was performed to examine the some properties of the molecules.Communicated by Ramaswamy H. Sarma.

In Silico and In Vivo Evaluation of the Maqui Berry (Aristotelia chilensis (Mol.) Stuntz) on Biochemical Parameters and Oxidative Stress Markers in a Metabolic Syndrome Model

Metabolic syndrome (MetS) is a complex disease that includes metabolic and physiological alterations in various organs such as the heart, pancreas, liver, and brain. Reports indicate that blackberry consumption, such as maqui berry, has a beneficial effect on chronic diseases such as cardiovascular disease, obesity, and diabetes. In the present study, in vivo and in silico studies have been performed to evaluate the molecular mechanisms implied to improve the metabolic parameters of MetS. Fourteen-day administration of maqui berry reduces weight gain, blood fasting glucose, total blood cholesterol, triacylglycerides, insulin resistance, and blood pressure impairment in the diet-induced MetS model in male and female rats. In addition, in the serum of male and female rats, the administration of maqui berry (MB) improved the concentration of MDA, the activity of SOD, and the formation of carbonyls in the group subjected to the diet-induced MetS model. In silico studies revealed that delphinidin and its glycosylated derivatives could be ligands of some metabolic targets such as α-glucosidase, PPAR-α, and PPAR-γ, which are related to MetS parameters. The experimental results obtained in the study suggest that even at low systemic concentrations, anthocyanin glycosides and aglycones could simultaneously act on different targets related to MetS. Therefore, these molecules could be used as coadjuvants in pharmacological interventions or as templates for designing new multitarget molecules to manage patients with MetS.

Design, syntheses, theoretical calculations, MM-GBSA, potential anti-cancer and enzyme activities of novel Schiff base compounds

In this study, new Schiff base compounds (SB-F-OH, SB-Cl-OH and SB-Br-OH) were derived from chalcone-derived amine compounds containing halogen groups and 4-hydroxybenzaldehyde. Also, their phthalonitrile compounds (SB-F-CN, SB-Cl-CN and SB-Br-CN) have been synthesized. The structures of these compounds were elucidated by NMR, FT-IR and Mass spectroscopic methods. The quantum chemical parameters were calculated at B3LYP/6-31++g(d,p), HF/6-31++g(d,p) and M062X/6-31++g(d,p) levels. As the biological application of the synthesized compounds, (i) their inhibition properties of the synthesized compounds on Acetylcholinesterase (AChE) and Butyrylcholinesterase (BChE) metabolic enzymes were investigated, and their potential anticancer activities against neuroblastoma (NB; SH-SY5Y) and healthy fibroblast (NIH-3T3) cell lines were determined by in vitro assays. All compounds showed inhibition at nanomolar level with the Ki values in the range of 97.86 ± 30.51-516.82 ± 31.42 nM for AChE, 33.21 ± 4.45-78.50 ± 8.91 nM for BChE, respectively. It has been determined that all tested compounds have a remarkable cytotoxic effect against SH-SY5Y, and IC50 values were significantly lower than NIH-3T3 cells. The lowest IC50 value was observed in SB-Cl-OH (7.48 ± 0.86 µM) and SB-Cl-CN (7.31 ± 0.69 µM). The molecular docking of the molecules was also investigated using crystal structure of AChE enzyme protein (PDB ID: 4M0E), crystal structure of BChE protein (PDB ID: 6R6V) and SH-SY5Y cancer protein (PDB ID: 2F3F, 3PBL and 5WIV). The ADME properties of the compounds were investigated. MM/GBSA method is calculated binding free energy. Afterwards, ADME/T analysis was performed to examine the some properties of the molecules.Communicated by Ramaswamy H. Sarma.

Combination of bendamustine-azacitidine against Syk target of breast cancer: an in silico study

Breast cancer (BC) is the most serious and second leading cause of death in women worldwide. When breast cancer is diagnosed and treated early, the chance of long-term survival is up to 90%. On the other hand, 90% of BC patient deaths are due to metastasis and a lack of effective early diagnosis. The existing conventional chemotherapy provides negative feedback due to transportation barriers towards the action sites, multidrug resistance, poor bio-availability, non-specific delivery and systemic side effects on the healthy tissue. Syk protein Kinase has been reported in BC, as a tumor modulator, providing a pro-survival signal and also by restricting epithelial-mesenchymal transition, enhancing cell-cell interactions and inhibiting migration. In the present study, we explored the possibility of targeting BC by attenuating Syk protein Kinase. Hence, we have conjugated the hydrophobic Bendamustine (BEN) and hydrophilic Azacitidine (AZA) anticancer drugs to evaluate their efficacy against BC. The native drugs (BEN and AZA) and designed drug-drug conjugate (BEN-AZA) were docked with Syk protein. Then, the docked complex was performed for Binding Free Energy and Molecular Dynamics Simulations. Furthermore, DFT and ADME properties were carried out. The results revealed that the designed drug-drug conjugate has a better docking score, ΔGbind and admirable stability throughout the simulation when compared with native drugs. In DFT and ADME analyses, the designed drug-drug conjugate has shown good stereo electronic features and pharmaceutical relevant parameters than that of native drugs. The overall results suggested that the designed drug-drug conjugate may be a suitable candidate for BC treatment.Communicated by Ramaswamy H. Sarma.

Design, synthesis, molecular modeling, DFT, ADME and biological evaluation studies of some new 1,3,4-oxadiazole linked benzimidazoles as anticancer agents and aromatase inhibitors

Breast cancer is the most frequent female cancer and second cause of cancer-related deaths among women around the world. Two thirds of breast cancer patients have hormone-dependent tumors, which is very likely be treated with hormonal therapy. Aromatase is involved in the biosynthesis of estrogen thus a critical target for breast cancer. In this study, in order to identify new aromatase enzyme inhibitors, a series of benzimidazole-1,3,4-oxadiazole derivatives were synthesized and characterized by ¹H NMR, ¹³C NMR, and MS spectra analyses. In the in vitro anticancer assay, all the compounds tested anticancer activities using MTT-based assay against five cancer cell lines (MCF-7, A549, HeLa, C6, and HepG2). Among them, compound 5a exhibited the most potent activity with IC₅₀ values of 5.165 ± 0.211 μM and 5.995 ± 0.264 μM against MCF-7 and HepG2 cell lines. Compound 5a was included in the BrdU test to determine the DNA synthesis inhibition effects for both cell types. Furthermore, compound 5c was also found to be more effective than doxorubicin on the HeLa cell line. The selectivity of anticancer activity was evaluated in NIH3T3 cell line. In vitro, enzymatic inhibition assays of aromatase enzyme were performed for compound 5a acting on the MCF-7 cell line. For compound 5a, in silico molecular docking and dynamics simulations against aromatase enzyme was performed to determine possible protein-ligand interactions and stability. DFT study was performed to evaluate the quantum mechanical and electronic properties of compound 5a. Finally, the theoretical ADME properties of the potential aromatase inhibitor compound 5a were analyzed by calculations.Communicated by Ramaswamy H. Sarma.

Molecular docking and dynamic simulation approach to decipher steroidal sapogenins (genus Trillium) derived agonists for glucocorticoid receptor

Steroidal sapogenins (SS) are structural analogues of steroidal drugs, which are frequently used for the treatment of several diseases including reproductive, malignancies, neurological, and inflammation-related diseases. The glucocorticoid receptor (GR) is a nuclear receptor that regulates development, metabolism, and inflammation, in response to steroidal ligands. Therefore, GR is considered as a potential therapeutic target for steroidal agents to the treatment of inflammation-related diseases. We hypothesized that SS may act as an agonist for GR due to structural similarity with corticosteroids. In this study, we carried out in silico screening of various SS from the genus Trillium to check their potential as an agonist for GR. Our data suggest that out of 42 SS, only 7 molecules have interacted with GR. However, molecular mechanics with generalized Born and surface area (MM-GBSA) analysis revealed that only two SS (SS 38 and SS 39) molecules bind favorably to GR. Among these, SS 38 (docking score: -9.722 Kcal/mol and MM-GBSA ΔG_bind: -50.192 Kcal/mol) and SS 39 (docking score: -11.20 Kcal/mol and MM-GBSA ΔG_bind: -58.937 Kcal/mol) have best docking and MM-GBSA scores. Molecular dynamics (MD) simulation studies of SS 38, SS 39, and dexamethasone-GR complex revealed that both SS shows hydrogen bonding and hydrophobic interaction with GR over the 120 ns simulation with mild fluctuations. The current study suggests that SS 38 and SS 39 may be further explored as a potential agonist to treat several disease conditions mediated by GR.

The Essential Oil from Oliveria decumbens Vent. (Apiaceae) as Inhibitor of Breast Cancer Cell (MCF-7) Growth

Oliveria decumbens Vent. is an aromatic and medicinal plant traditionally used in Iran for the treatment of infections, gastrointestinal diseases, cancer, and inflammation. This research was aimed at investigating the pharmacological potential of O. decumbens essential oil (OEO) and its main compounds, focusing on OEO's cytotoxic effects on MCF-7 breast cancer cells. OEO was obtained by hydro-distillation, and the chemical constituents were identified using GC-MS. Thymol, carvacrol, γ-terpinene, and p-cymene were the main OEO constituents. When MCF-7 cells were treated with OEO, the expressions of genes related to apoptosis (BIM and Bcl-2), tumor suppression (PTEN), and cell growth inhibition (AURKA), were evaluated using real-time PCR. Moreover, molecular docking was used for studying in silico the interaction of OEO principal compounds with PTEN and AURKA. The expression of AURKA was significantly reduced since the OEO treatment enhanced the expression of PTEN. Through in silico molecular docking, it was revealed that thymol, carvacrol, p-cymene, and γ-terpinene can activate PTEN and thus inhibit AURKA. Additionally, the DNA fragmentation assay, acridine orange/ethidium bromide (AO/EB) double-staining assay, and real-time PCR highlighted the fact that the OEO treatment could activate apoptosis and inhibit cell proliferation. Therefore, OEO is a viable candidate to be employed in the pharmaceutical industry, specifically as a possible agent for cancer therapy.

Potent Bioactive Compounds From Seaweed Waste to Combat Cancer Through Bioinformatics Investigation

The seaweed industries generate considerable amounts of waste that must be appropriately managed. This biomass from marine waste is a rich source of high-value bioactive compounds. Thus, this waste can be adequately utilized by recovering the compounds for therapeutic purposes. Histone deacetylases (HDACs) are key epigenetic regulators established as one of the most promising targets for cancer chemotherapy. In the present study, our objective is to find the HDAC 2 inhibitor. We performed top-down in silico methodologies to identify potential HDAC 2 inhibitors by screening compounds from edible seaweed waste. One hundred ninety-three (n = 193) compounds from edible seaweeds were initially screened and filtered with drug-likeness properties using SwissADME. After that, the filtered compounds were followed to further evaluate their binding potential with HDAC 2 protein by using Glide high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP), and quantum polarized ligand docking (QPLD). One compound with higher negative binding energy was selected, and to validate the binding mode and stability of the complex, molecular dynamics (MD) simulations using Desmond were performed. The complex-binding free energy calculation was performed using molecular mechanics-generalized born surface area (MM-GBSA) calculation. Post-MD simulation analyses such as PCA, DCCM, and free energy landscape were also evaluated. The quantum mechanical and electronic properties of the potential bioactive compounds were assessed using the density functional theory (DFT) study. These findings support the use of marine resources like edible seaweed waste for cancer drug development by using its bioactive compounds. The obtained results encourage further in vitro and in vivo research. Our in silico findings show that the compound has a high binding affinity for the catalytic site of the HDAC 2 protein and has drug-likeness properties, and can be utilized in drug development against cancer.

Atomic level and structural understanding of natural ligands inhibiting Helicobacter pylori peptide deformylase through ligand and receptor based screening, SIFT, molecular dynamics and DFT - a structural computational approach

Helicobacter pylori is a Gram-negative microaerophilic gastric pathogen, responsible for the cause of peptic ulcer around half of the global population. Although several antibiotics and combination therapies have been employed for H. pylori-related gastric ulcer and cancer regiments, identifying potent inhibitors for specific targets of this bacterium will help assessing better treatment periodicity and methods to eradicate H. pylori. Herein, 1,000,000 natural compounds were virtually screened against Helicobacter pylori Peptide deformylase (HpPDF). Pharmacophore hypotheses were created using ligand and receptor-based pharmacophore modeling of GLIDE. Stringent HTVS and IFD docking protocol of GLIDE predicted leads with stable intermolecular bonds and scores. Molecular dynamics simulation of HpPDF was carried out for 100 ns using GROMACS. Hits ZINC00225109 and ZINC44896875 came up with a glide score of -9.967 kcal/mol and -12.114 kcal/mol whereas; reference compound actinonin produced a glide score of -9.730 kcal/mol. Binding energy values of these hits revealed the involvement of significant Van der Waals and Coulomb forces and the deduction of lipophilic forces that portray the deep hydrophobic residues in the S1pocket of H. pylori. The DFT analysis established the electron density-based features of the molecules and observed that the results correlate with intermolecular docking interactions. Analysis of the MD trajectories revealed the crucial residues involved in HpPDF - ligand binding and the conformational changes in the receptor. We have identified and deciphered the crucial features necessary for the potent ligand binding at catalytic site of HpPDF. The resulting ZINC natural compound hits from the study could be further employed for potent drug development.Communicated by Ramaswamy H. Sarma.

Molecular docking and dynamic simulations study for repurposing of multitarget coumarins against SARS-CoV-2 main protease, papain-like protease and RNA-dependent RNA polymerase

Proteases and RNA-Dependent RNA polymerase, major enzymes which are essential targets involved in the life and replication of SARS-CoV-2. This study aims at in silico examination of the potential ability of coumarins and their derivatives to inhibit the replication of SARS-Cov-2 through multiple targets, including the main protease, papain-like protease and RNA-Dependent RNA polymerase. Several coumarins as biologically active compounds were studied, including coumarin antibiotics and some naturally reported antiviral coumarins. Aminocoumarin antibiotics, especially coumermycin, showed a high potential to bind to the enzymes’ active site, causing possible inhibition and termination of viral life. They demonstrate the ability to bind to residues essential for triggering the crucial cascades within the viral cell. Molecular dynamics simulations for 50 ns supported these data pointing out the formation of rigid, stable Coumermycin/enzyme complexes. These findings strongly suggest the possible use of Coumermycin, Clorobiocin or Novobiocin in the fight against COVID-19, but biological evidence is still required to support such suggestions.

Computational identification of 2,4-disubstituted amino-pyrimidines as L858R/T790M-EGFR double mutant inhibitors using pharmacophore mapping, molecular docking, binding free energy calculation, DFT study and molecular dynamic simulation

Pharmacophore modelling studies have been performed for a series of 2,4-disubstituted-pyrimidines derivatives as EGFR L858R/T790M tyrosine kinase inhibitors. The high scoring AARR.15 hypothesis was selected as the best pharmacophore model with the highest survival score of 3.436 having two hydrogen bond acceptors and two aromatic ring features. Pharmacophore-based virtual screening followed by structure-based yielded the six molecules (ZINC17013227, ZINC17013215, ZINC9573324, ZINC9573445, ZINC24023331 and ZINC17013503) from the ZINC database with significant in silico predicted activity and strong binding affinity towords the EGFR L858R/T790M tyrosine kinase. In silico toxicity and cytochrome profiling indicates that all the 06 virtually screened compounds were substrate/inhibitors of the CYP-3A4 metabolizing enzyme and were non-carcinogenic and devoid of Ames mutagenesis. Density functional theory (DFT) and molecular dynamic (MD) simulation further validated the obtained hits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00113-x.

Molecular studies of immunological enzyme clumping factor B for the inhibition of Staphylococcus aureus with essential oils of Nigella sativa

Essential oils from black cumin seeds (Nigella sativa) have largely been used in the manufacturing of nutraceuticals and functional food products due to the presence of a wide variety of bioactive compounds. However, their applications in the pharmaceutical sector have recently attracted interest and started blooming. The present research elucidates the in silico and in vitro efficacies of active leads from essential oil of N sativa against the human pathogenic bacterium Staphylococcus aureus. Biofilm development has become an inevitable situation in the health care sector. Lowering the efficacies of antimicrobial drugs is one of the vital ramifications that resulted in the emergence of multidrug resistance. Clumping factor B (clfB) of S aureus plays a key role in the human immune functions during pathogenesis. Through STRING analysis, the interacting protein partners of clfB were found to regulate biofilm pathway. Therefore, eight ligands from essential oil are docked with the critical clfB protein, which revealed p-cymene, thymoquinone and carvacrol as the robust ligands with highest binding affinity. Therefore, antibiofilm potential of N sativa essential oil at in vitro states was evaluated against S aureus. Further, real time PCR analysis showed that the expression of clfB and intercellular adhesion gene (icaA and icaD) was significantly altered upon treatment with essential oil. Altogether, the findings confirmed the antibiofilm efficacy of N sativa essential oil against S aureus. Hence, the essential oil from N. sativa was envisaged to be promising candidate to treat S aureus biofilm mediated infection.

High-Throughput Screening and Quantum Mechanics for Identifying Potent Inhibitors Against Mac1 Domain of SARS-CoV-2 Nsp3

SARS-CoV-2 encodes the Mac1 domain within the large nonstructural protein 3 (Nsp3), which has an ADP-ribosylhydrolase activity conserved in other coronaviruses. The enzymatic activity of Mac1 makes it an essential virulence factor for the pathogenicity of coronavirus (CoV). They have a regulatory role in counteracting host-mediated antiviral ADP-ribosylation, which is unique part of host response towards viral infections. Mac1 shows highly conserved residues in the binding pocket for the mono and poly ADP-ribose. Therefore, SARS-CoV-2 Mac1 enzyme is considered as an ideal drug target and inhibitors developed against them can possess a broad antiviral activity against CoV. ADP-ribose-1 phosphate bound closed form of Mac1 domain is considered for screening with large database of ZINC. XP docking and QPLD provides strong potential lead compounds, that perfectly fits inside the binding pocket. Quantum mechanical studies expose that, substrate and leads have similar electron donor ability in the head regions, that allocates tight binding inside the substrate-binding pocket. Molecular dynamics study confirms the substrate and new lead molecules presence of electron donor and acceptor makes the interactions tight inside the binding pocket. Overall binding phenomenon shows both substrate and lead molecules are well-adopt to bind with similar binding mode inside the closed form of Mac1.

From genomes to molecular dynamics - A bottom up approach in extrication of SARS CoV-2 main protease inhibitors

The recent pandemic Coronavirus disease-19 outbreak had traumatized global countries since its origin in late December 2019. Though the virus originated in China, it has spread rapidly across the world due its firmly established community transmission. To successfully tackle the spread and further infection, there needs a clear multidimensional understanding of the molecular mechanisms. Henceforth, 942 viral genome sequences were analysed to predict the core genomes crucial in virus life cycle. Additionally, 35 small interfering RNA transcripts were predicted that can target specifically the viral core proteins and reduce pathogenesis. The crystal structure of Covid-19 main protease-6LU7 was chosen as an attractive target due to the factors that there were fewer mutations and whose structure had significant identity to the annotated protein sequence of the core genome. Drug repurposing of both recruiting and non recruiting drugs was carried out through molecular docking procedures to recognize bitolterol as a good inhibitor of Covid-19 protease. The study was extended further to screen antiviral phytocompounds through quantitative structure activity relationship and molecular docking to identify davidigenin, from licorice as the best novel lead with good interactions and binding energy. The docking of the best compounds in all three categories was validated with molecular dynamics simulations which implied stable binding of the drug and lead molecule. Though the studies need clinical evaluations, the results are suggestive of curbing the pandemic.

Microsecond MD Simulation and Multiple-Conformation Virtual Screening to Identify Potential Anti-COVID-19 Inhibitors Against SARS-CoV-2 Main Protease

The recent pandemic outbreak of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), raised global health and economic concerns. Phylogenetically, SARS-CoV-2 is closely related to SARS-CoV, and both encode the enzyme main protease (M^pro/3CL^pro), which can be a potential target inhibiting viral replication. Through this work, we have compiled the structural aspects of M^pro conformational changes, with molecular modeling and 1-μs MD simulations. Long-scale MD simulation resolves the mechanism role of crucial amino acids involved in protein stability, followed by ensemble docking which provides potential compounds from the Traditional Chinese Medicine (TCM) database. These lead compounds directly interact with active site residues (His41, Gly143, and Cys145) of M^pro, which plays a crucial role in the enzymatic activity. Through the binding mode analysis in the S1, S1′, S2, and S4 binding subsites, screened compounds may be functional for the distortion of the oxyanion hole in the reaction mechanism, and it may lead to the inhibition of M^pro in SARS-CoV-2. The hit compounds are naturally occurring compounds; they provide a sustainable and readily available option for medical treatment in humans infected by SARS-CoV-2. Henceforth, extensive analysis through molecular modeling approaches explained that the proposed molecules might be promising SARS-CoV-2 inhibitors for the inhibition of COVID-19, subjected to experimental validation.

Docking-based virtual screening of TβR1 inhibitors: evaluation of pose prediction and scoring functions

To improve the reliability of virtual screening for transforming growth factor-beta type 1 receptor (TβR1) inhibitors, 2 docking methods and 11 scoring functions in Discovery Studio software were evaluated and validated in this study. LibDock and CDOCKER protocols were performed on a test set of 24 TβR1 protein-ligand complexes. Based on the root-mean-square deviation (RMSD) values (in Å) between the docking poses and co-crystal conformations, the CDOCKER protocol can be efficiently applied to obtain more accurate dockings in medium-size virtual screening experiments of TβR1, with a successful docking rate of 95%. A dataset including 281 known active and 8677 inactive ligands was used to determine the best scoring function. The receiver operating characteristic (ROC) curves were used to compare the performance of scoring functions in attributing best scores to active than inactive ligands. The results show that Ludi 1, PMF, Ludi 2, Ludi 3, PMF04, PLP1, PLP2, LigScore2, Jain and LigScore1 are better scoring functions than the random distribution model, with AUC of 0.864, 0.856, 0.842, 0.812, 0.776, 0.774, 0.769, 0.762, 0.697 and 0.660, respectively. Based on the pairwise comparison of ROC curves, Ludi 1 and PMF were chosen as the best scoring functions for virtual screening of TβR1 inhibitors. Further enrichment factors (EF) analysis also supports PMF and Ludi 1 as the top two scoring functions.

Theoretical Insights into the Anti-SARS-CoV-2 Activity of Chloroquine and Its Analogs and In Silico Screening of Main Protease Inhibitors

Corona virus disease (COVID-19) is a dangerous disease rapidly spreading all over the world today. Currently there are no treatment options for it. Drug repurposing studies explored the potency of antimalarial drugs, chloroquine and hydroxychloroquine, against SARS-CoV-2 virus. These drugs can inhibit the viral protease, called chymotrypsin-like cysteine protease, also known as Main protease (3CLpro); hence, we studied the binding efficiencies of 4-aminoquinoline and 8-aminoquinoline analogs of chloroquine. Six compounds furnished better binding energies than chloroquine and hydroxychloroquine. The interactions with the active site residues especially with Cys145 and His41, which are involved in catalytic diad for proteolysis, make these compounds potent main protease inhibitors. A regression model correlating binding energy and the molecular descriptors for chloroquine analogs was generated with R2 = 0.9039 and Q2 = 0.8848. This model was used to screen new analogs of primaquine and molecules from the Asinex compound library. The docking and regression analysis showed these analogs to be more potent inhibitors of 3CLpro than hydroxychloroquine and primaquine. The molecular dynamic simulations of the hits were carried out to determine the binding stabilities. Finally, we propose four compounds that show drug likeness toward SARS-CoV-2 that can be further validated through in vitro and in vivo studies.

Bacterial protein azurin and derived peptides as potential anti-SARS-CoV-2 agents: insights from molecular docking and molecular dynamics simulations

The current pandemic SARS-CoV-2 has wreaked havoc in the world, and neither drugs nor vaccine is available for the treatment of this disease. Thus, there is an immediate need for novel therapeutics that can combat this deadly infection. In this study, we report the therapeutic assessment of azurin and its peptides: p18 and p28 against the viral structural S-protein and non-structural 3CL^pro and PL^pro proteins. Among the analyzed complexes, azurin docked relatively well with the S2 domain of S-protein compared to the other viral proteins. The derived peptide p18 bound to the active site domain of the PL^pro protein; however, in other complexes, lesser interactions were recorded. The second azurin derived peptide p28, fared the best among the docked proteins. p28 interacted with all the three viral proteins and the host ACE-2 receptor by forming several electrostatic and hydrogen bonds with the S-protein, 3CL^pro, and PL^pro. MD simulations indicated that p28 exhibited a strong affinity to S-protein and ACE-2 receptor, indicating a possibility of p28 as a protein-protein interaction inhibitor. Our data suggest that the p28 has potential as an anti-SARS-CoV-2 agent and can be further exploited to establish its validity in the treatment of current and future SARS-CoV crisis.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening and molecular dynamics simulation of SARS-CoV-2 Guanine-N7 methyltransferase (nsp14) for identifying antiviral inhibitors against COVID-19

The recent pandemic caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) calls the whole world into a medical emergency. For tackling Coronavirus Disease 2019 (COVID-19), researchers from around the world are swiftly working on designing and identifying inhibitors against all possible viral key protein targets. One of the attractive drug targets is guanine-N7 methyltransferase which plays the main role in capping the 5′-ends of viral genomic RNA and sub genomic RNAs, to escape the host’s innate immunity. We performed homology modeling and molecular dynamic (MD) simulation, in order to understand the molecular architecture of Guanosine-P3-Adenosine-5’,5’-Triphosphate (G3A) binding with C-terminal N7-MTase domain of nsp14 from SARS-CoV-2. The residue Asn388 is highly conserved in present both in N7-MTase from SARS-CoV and SARS-CoV-2 and displays a unique function in G3A binding. For an in-depth understanding of these substrate specificities, we tried to screen and identify inhibitors from the Traditional Chinese Medicine (TCM) database. The combination of several computational approaches, including screening, MM/GBSA, MD simulations, and PCA calculations, provides the screened compounds that readily interact with the G3A binding site of homology modeled N7-MTase domain. Compounds from this screening will have strong potency towards inhibiting the substrate-binding and efficiently hinder the viral 5’-end RNA capping mechanism. We strongly believe the final compounds can become COVID-19 therapeutics, with huge international support.

Communicated by Ramaswamy H. Sarma.

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.

α-Glucosidase inhibitors from Chinese bayberry (Morella rubra Sieb. et Zucc.) fruit: molecular docking and interaction mechanism of flavonols with different B-ring hydroxylations

Inhibition of α-glucosidase alleviates postprandial high glycemic levels in diabetic or prediabetic population.