Exploring the potential of 2-arylbenzimidazole scaffolds as novel α-amylase inhibitors: QSAR, molecular docking, simulation and pharmacokinetic studies.

Previous studies have shown that 2-arylbenzimidazole derivatives have a strong anti-diabetic effect. To further explore this potential, we develop new analogues of the compound using ligand-based drug design and tested their inhibitory and binding properties through QSAR analyses, molecular docking, dynamic simulations and pharmacokinetic studies. By using quantitative structure activity relationship and ligand-based modification, a highly precise predictive model and design of potent compounds was developed from the derivatives of 2-arylbenzimidazoles. Molecular docking and simulation studies were then conducted to identify the optimal binding poses and pharmacokinetic profiles of the newly generated therapeutic drugs. DFT was employed to optimize the chemical structures of 2-arylbenzimidazole derivatives using B3LYP/6-31G* as the basis set. The model with the highest R2trng set, R2adj, Q2cv, and R2test sets (0.926, 0.912, 0.903, and 0.709 respectively) was chosen to predict the inhibitory activities of the derivatives. Five analogues designed using ligand-based strategy had higher activity than the hit molecule. Additionally, the designed molecules had more favorable MolDock scores than the hit molecule and acarbose and simulation studies confirm on their stability and binding affinities towards the protein. The ADME and druglikeness properties of the analogues indicated that they are safe to consume orally and have a high potential for total clearance. The results of this study showed that the suggested analogues could act as α-amylase inhibitors, which could be used as a basis for the creation of new drugs to treat type 2 diabetes mellitus.

Identification of novel peptide inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR): molecular docking and MD simulation studies.

The presence of drug-resistant variants of Plasmodium parasites within the population has presented a substantial obstacle to the eradication of Malaria. As a result, numerous research groups have directed their efforts towards creating new medication candidates that specifically target parasites. In this study, our main objective was to identify tri-peptide inhibitors for Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) with the aim of finding a new peptide that exhibits superior binding properties compared to the current inhibitor, WR99210. In order to achieve this objective, a virtual library consisting of 8000 tripeptides was generated and subjected to computational screening against wild-type PfDHFR. The purpose of this screening was to discover the most effective binders at the active site. The four most optimal tripeptides identified (Trp-Trp-Glu, Trp-Phe-Tyr, Phe-Trp-Trp, Tyr-Trp-Trp) exhibited significant non-covalent interactions inside the active site of PfDHFR and had binding energies ranging from -9.5 to -9.0 kcal/mol and WR99210 had a binding energy of -6.2 kcal/mol. A 250 ns Molecular Dynamics (MD) simulation was performed to investigate the kinetic and thermodynamic characteristics of the protein-ligand complexes. The Root Mean Square Deviation (RMSD) values for the optimal tripeptides fell within the allowed range, indicating the stability of the ligands inside the protein complex. The Ki value for the most effective tripeptide was 0.3482 µM, whereas WR99210 had a Ki value of 1.02 µM. This article presents the initial discovery of peptide inhibitors targeting PfDHFR. In this text, we provide a comprehensive explanation of the interactions that occur between peptides and the enzyme.Communicated by Ramaswamy H. Sarma.

Ligand based-design of potential schistosomiasis inhibitors through QSAR, homology modeling, molecular dynamics, pharmacokinetics, and DFT studies.

Objectives: Schistosomiasis, a neglected tropical disease, is a leading cause of mortality in affected geographic areas. Currently, because no vaccine for schistosomiasis is available, control measures rely on widespread administration of the drug praziquantel (PZQ). The mass administration of PZQ has prompted concerns regarding the emergence of drug resistance. Therefore, new therapeutic targets and potential compounds are necessary to combat schistosomiasis. Methods: Twenty-four potent derivatives of PZQ were optimized via density functional theory (DFT) at the B3LYP/6-31G∗ level. Quantitative structureactivity relationship (QSAR) models were generated and statistically validated, and a lead candidate was selected to develop therapeutic options with improved efficacy against schistosomiasis. The biological and binding energies of the designed compounds were evaluated. In addition, molecular dynamics; drug-likeness; absorption, distribution, metabolism, excretion, and toxicity (ADMET); and DFT studies were performed on the newly designed compounds. Results: Five QSAR models were generated, among which model 1 had favorable validation parameters (R2train: 0.957, R2adj: 0.941, LOF: 0.101, Q2cv: 0.906, and R2test: 0.783) and was chosen to identify a lead candidate. Other statistical parameters for the chosen model included variance inflation factor values ranging from 1.242 to 1.678, and a Y-scrambling coefficient (cRp2) of 0.747. Five new compounds were designed with improved predicted activity (ranging from 5.081 to 7.022) surpassing those of both the lead compound and PZQ (predicted pEC50 of 5.545). Molecular dynamics simulation revealed high binding affinity of the proposed compounds toward the target receptor. ADMET and drug-likeness assessments indicated adherence to Lipinski's rule of five criteria, thereby suggesting pharmacological and oral safety. In addition, DFT analysis indicated resistance to electronic alteration during chemical reactions. Conclusion: The proposed compounds exhibited potential drug characteristics, thus indicating their suitability for further investigation to enhance schistosomiasis treatment options.

Chemical library design, QSAR modeling and molecular dynamics simulations of naturally occurring coumarins as dual inhibitors of MAO-B and AChE.

Coumarins are a highly privileged scaffold in medicinal chemistry. It is present in many natural products and is reported to display various pharmacological properties. A large plethora of compounds based on the coumarin ring system have been synthesized and were found to possess biological activities such as anticonvulsant, antiviral, anti-inflammatory, antibacterial, antioxidant as well as neuroprotective properties. Despite the wide activity spectrum of coumarins, its naturally occurring derivatives are yet to be investigated in detail. In the current study, a chemical library was created to assemble all chemical information related to naturally occurring coumarins from the literature. Additionally, a multi-stage virtual screening combining QSAR modeling, molecular docking, and ADMET prediction was conducted against monoamine oxidase B and acetylcholinesterase, two relevant targets known for their neuroprotective properties and 'disease-modifying' potential in Parkinson's and Alzheimer's disease. Our findings revealed ten coumarin derivatives that may act as dual-target drugs against MAO-B and AChE. Two coumarin candidates were selected from the molecular docking study: CDB0738 and CDB0046 displayed favorable interactions for both proteins as well as suitable ADMET profiles. The stability of the selected coumarins was assessed through 100 ns molecular dynamics simulations which revealed promising stability through key molecular interactions for CDB0738 to act as dual inhibitor of MAO-B and AChE. However, experimental studies are necessary to evaluate the bioactivity of the proposed candidate. The current results may generate an increasing interest in bioprospecting naturally occurring coumarins as potential candidates against relevant macromolecular targets by encouraging virtual screening studies against our chemical library.Communicated by Ramaswamy H. Sarma.

In-silico molecular modelling studies of some camphor imine based compounds as anti-influenza A (H1N1) pdm09 virus agents.

The advent of influenza A (H1N1) drug-resistant strains led to the search quest for more potent inhibitors of the influenza A virus, especially in this devastating COVID-19 pandemic era. Hence, the present research utilized some molecular modelling strategies to unveil new camphor imine-based compounds as anti-influenza A (H1N1) pdm09 agents. The 2D-QSAR results revealed GFA-MLR (R2train = 0.9158, Q2=0.8475) and GFA-ANN (R2train = 0.9264, Q2=0.9238) models for the anti-influenza A (H1N1) pdm09 activity prediction which have passed the QSAR model acceptability thresholds. The results from the 3D-QSAR studies also revealed CoMFA (R2train =0.977, Q2=0.509) and CoMSIA_S (R2train =0.976, Q2=0.527) models for activity predictions. Based on the notable information derived from the 2D-QSAR, 3D-QSAR, and docking analysis, ten (10) new camphor imine-based compounds (22a-22j) were designed using the most active compound 22 as the template. Furthermore, the high predicted activity and binding scores of compound 22j were further justified by the high reactive sites shown in the electrostatic potential maps and other quantum chemical calculations. The MD simulation of 22j in the active site of the influenza hemagglutinin (HA) receptor confirmed the dynamic stability of the complex. Moreover, the appraisals of drug-likeness and ADMET properties of the proposed compounds showed zero violation of Lipinski's criteria with good pharmacokinetic profiles. Hence, the outcomes in this work recommend further in-depth in vivo and in-vitro investigations to validate these theoretical findings.Communicated by Ramaswamy H. Sarma.

Immunoinformatics and reverse vaccinology approach in designing a novel highly immunogenic multivalent peptide-based vaccine against the human monkeypox virus.

Background: Monkeypox is a highly infectious zoonotic disease, often resulting in complications ranging from respiratory illnesses to vision loss. The escalating global incidence of its cases demands prompt attention, as the absence of a proven post-exposure treatment underscores the criticality of developing an effective vaccine. Methods: Interactions of the viral proteins with TLR2 and TLR4 were investigated to assess their immunogenic potentials. Highly immunogenic proteins were selected and subjected to epitope mapping for identifying B-cell and MHC class I and II epitopes. Epitopes with high antigenicity were chosen, considering global population coverage. A multi-target, multi-epitope vaccine peptide was designed, incorporating a beta-defensin 2 adjuvant, B-cell epitopes, and MHC class I and II epitopes. Results: The coordinate structure of the engineered vaccine was modeled and validated. In addition, its physicochemical properties, antigenicity, allergenicity, and virulence traits were evaluated. Molecular docking studies indicated strong interactions between the vaccine peptide and the TLR2 receptor. Furthermore, molecular dynamics simulations and immune simulation studies reflected its potent cytosolic stability and robust immune response dynamics induced by the vaccine. Conclusion: This study explored an innovative structure-guided approach in the use of immunoinformatics and reverse vaccinology in pursuit of a novel multi-epitope vaccine against the highly immunogenic monkeypox viral proteins. The simulation studies indicated the engineered vaccine candidate to be promising in providing prophylaxis to the monkeypox virus; nevertheless, further in vitro and in vivo investigations are required to prove its efficacy.

Molecular modelling studies of substituted indole derivatives as novel influenza a virus inhibitors.

The emergence of drug-resistant strains motivate researchers to find new innovative anti-IAV candidates with a different mode of action. In this work, molecular modelling strategies, such as 2D-QSAR, 3D-QSAR, molecular docking, molecular dynamics, FMOs, and ADMET were applied to some substituted indoles as IAV inhibitors. The best-developed 2D-QSAR models, MLR (Q2 = 0.7634, R2train = 0.8666) and ANN[4-3-1] (Q2 = 0.8699, R2train = 0.8705) revealed good statistical validation for the inhibitory response predictions. The 3D-QSAR models, CoMFA (Q2 = 0.504, R2train = 0.805) and CoMSIA/SEDHA (Q2 = 0.619, R2train = 0.813) are selected as the best 3D models following the global thresholds. In addition, the contour maps generated from the CoMFA and CoMSIA models illustrate the relationship between the molecular fields and the inhibitory effects of the studied molecules. The results of the studies led to the design of five new molecules (24a-e) with enhanced anti-IAV activities and binding potentials using the most active molecule (24) as the template scaffold. The conformational stability of the best-designed molecules with the NA protein showed hydrophobic and H-bonds with the key residues from the molecular dynamics simulations of 100 ns. Furthermore, the global reactivity indices from the DFT calculations portrayed the relevance of 24c in view of its smaller band gap as also justified by our QSAR and molecular simulation studies.Communicated by Ramaswamy H. Sarma.

Exploring natural products as multi-target-directed drugs for Parkinson's disease: an in-silico approach integrating QSAR, pharmacophore modeling, and molecular dynamics simulations.

Parkinson's disease is a neurodegenerative disorder characterized by the progressive loss of dopaminergic neurons in the midbrain. Current treatments provide limited symptomatic relief without halting disease progression. A multi-targeting approach has shown potential benefits in treating neurodegenerative diseases. In this study, we employed in silico approaches to explore the COCONUT natural products database and identify novel drug candidates with multi-target potential against relevant Parkinson's disease targets. QSAR models were developed to screen for potential bioactive molecules, followed by a hybrid virtual screening approach involving pharmacophore modeling and molecular docking against MAO-B, AA2AR, and NMDAR. ADME evaluation was performed to assess drug-like properties. Our findings revealed 22 candidates that exhibited the desired pharmacophoric features. Particularly, two compounds: CNP0121426 and CNP0242698 exhibited remarkable binding affinities, with energies lower than -10 kcal/mol and promising interaction profiles with the chosen targets. Furthermore, all the ligands displayed desirable pharmacokinetic properties for brain-targeted drugs. Lastly, molecular dynamics simulations were conducted on the lead candidates, belonging to the dihydrochalcone and curcuminoid class, to evaluate their stability over a 100 ns timeframe and compare their dynamics with reference complexes. Our findings revealed the curcuminoid CNP0242698 to have an overall better stability with the three targets compared to the dihydrochalcone, despite the high ligand RMSD, the curcuminoid CNP0242698 showed better protein stability, implying ligand exploration of different orientations. Similarly, AA2AR exhibited higher stability with CNP0242698 compared to the reference complex, despite the high initial ligand RMSD due to the bulkier active site. In NMDAR, CNP0242698 displayed good stability and less fluctuations implying a more restricted conformation within the smaller active site of NMDAR. These results may serve as lead compounds for the development and optimization of natural products as multi-target disease-modifying natural remedies for Parkinson's disease patients. However, experimental assays remain necessary to validate these findings.Communicated by Ramaswamy H. Sarma.

Identification of high-affinity pyridoxal kinase inhibitors targeting cancer therapy: an integrated docking and molecular dynamics simulation approach.

Pyridoxal kinase (PDXK) is a vitamin B6-dependent transferase enzyme encoded by the PDXK gene, crucial for leukemic cell proliferation. Disruption of its activity causes altered metabolism and reduced levels of nucleotides and polyamines. PDXK and pyridoxal 5'-phosphate (PLP) are overexpressed in various carcinomas, making them promising targets for drug design against cancer. Targeting PDXK may hold promise as a therapeutic approach for cancer treatment. This study focused on discovering potential inhibitors that could selectively interrupt the binding of pyridoxal phosphate (PLP) to pyridoxal kinase (PDXK). A commercially available library of 7,28,747 natural and druglike compounds was virtually screened using a molecular docking approach to target the substrate binding pocket of PDXK. Six promising inhibitors were identified, and all-atom molecular dynamics simulations were conducted on the PDXK-ligand complexes for 100 ns to assess their binding conformational stability. The simulation results indicated that the binding of ZINC095099376, ZINC01612996, ZINC049841390, ZINC095098959, ZINC01482077, and ZINC03830976 induced a slight structural change and stabilized the PDXK structure. This analysis provided valuable information about the critical residues involved in the PDXK-PLP complex formation and can be utilized in designing specific and effective PDXK inhibitors. According to this study, these compounds could be developed as anticancer agents targeting PDXK as a potential candidate for further study.Communicated by Ramaswamy H. Sarma.

Repurposing of drugs against methyltransferase as potential Zika virus therapies.

In recent years, the outbreak of infectious disease caused by Zika Virus (ZIKV) has posed a major threat to global public health, calling for the development of therapeutics to treat ZIKV disease. Several possible druggable targets involved in virus replication have been identified. In search of additional potential inhibitors, we screened 2895 FDA-approved compounds using Non-Structural Protein 5 (NS5) as a target utilizing virtual screening of in-silco methods. The top 28 compounds with the threshold of binding energy -7.2 kcal/mol value were selected and were cross-docked on the three-dimensional structure of NS5 using AutoDock Tools. Of the 2895 compounds screened, five compounds (Ceforanide, Squanavir, Amcinonide, Cefpiramide, and Olmesartan_Medoxomil) ranked highest based on filtering of having the least negative interactions with the NS5 and were selected for Molecular Dynamic Simulations (MDS) studies. Various parameters such as RMSD, RMSF, Rg, SASA, PCA and binding free energy were calculated to validate the binding of compounds to the target, ZIKV-NS5. The binding free energy was found to be -114.53, -182.01, -168.19, -91.16, -122.56, and -150.65 kJ mol-1 for NS5-SFG, NS5-Ceforanide, NS5-Squanavir, NS5-Amcinonide, NS5-Cefpiramide, and NS5-Ol_Me complexes respectively. The binding energy calculations suggested Cefpiramide and Olmesartan_Medoxomil (Ol_Me) as the most stable compounds for binding to NS5, indicating a strong rationale for their use as lead compounds for development of ZIKV inhibitors. As these drugs have been evaluated on pharmacokinetics and pharmacodynamics parameters only, in vitro and in vivo testing and their impact on Zika viral cell culture may suggest their clinical trials on ZIKV patients.

Modeling of chitosan modified PLGA atorvastatin-curcumin conjugate (AT-CU) nanoparticles, overcoming the barriers associated with PLGA: An approach for better management of atherosclerosis.

Conjugate drugs are evolving into potent techniques in the drug development process for enhancing the biopharmaceutical, physicochemical, and pharmacokinetic properties. Atorvastatin (AT) is the first line of treatment for coronary atherosclerosis; however its therapeutic efficacy is limited because of its poor solubility and fast pass metabolism. Curcumin (CU) is evidenced in several crucial signaling pathways linked to lipid regulation and inflammation. To enhance the therapeutic efficacy and physical properties of AT and CU, a new conjugate derivative (AT-CU) was synthesized and assessed by in silico, in vitro characterizations, and in vivo efficacy through mice model. Although the biocompatibility and biodegradability of Polylactic-co-Glycolic Acid (PLGA) in nanoparticles are well documented, burst release is a common issue with this polymer. Hence the current work used chitosan as a drug release modifier to the PLGA nanoparticles. The chitosan-modified PLGA AT-CU nanoparticles were prepaid by single emulsion and solvent evaporation technique. With raising the concentration of chitosan the particle size grew from 139.2 nm to 197.7 nm, the zeta potential rose from -20.57 mV to 28.32 mV, and the drug encapsulation efficiency improved from 71.81% to 90.57%. At 18 h, the burst release of AT-CU from PLGA nanoparticles was seen, hitting abruptly 70.8%. For chitosan-modified PLGA nanoparticles, the burst release pattern was significantly reduced which could be due to the adsorption of the drug on the surface of chitosan. The efficiency of the ideal formulation i.e F4 (chitosan/PLGA = 0.4) in treating atherosclerosis was further strongly evidenced by in vivo investigation.

Drugs swapping in coronavirus strains: a structural biology view.

Coronavirus belongs to the coronaviridae family, having a single-stranded RNA as genetic material of 26-42 kb in size. The first coronavirus infection emerged in 2002, caused by SARS-CoV1. Since then, genome sequences and three-dimensional structures of crucial proteins and enzymes of the virus have been studied in detail. The novel coronavirus (nCoV) outbreak has caused the COVID19 pandemic, which is responsible for the deaths of millions of people worldwide. The nCoV was later renamed as SARS-CoV2. The details of most of the COV proteins are available at the atomic and molecular levels. The entire genome is made up of 12 open reading frames that code for 27 different proteins. The spike surface glycoprotein, the envelope protein, the nucleocapsid protein, and the membrane protein are the four structural proteins which are required for virus attachment, entrance, assembly, and pathogenicity. The remaining proteins encoded are called non-structural (NSPs) and support the survival of the virus. Several non-structural proteins are also validated targets for drug development against coronavirus and are being used for drug design purposes. To perform a comparative study, sequences and three-dimensional structures of four crucial viral enzymes, Mpro, PLpro, RdRp, and EndoU from SARS-CoV1 and SARS-CoV2 variants were analyzed. The key structural elements and ligands recognizing amino acid residues were found to be similar in enzymes from both strains. The significant sequences and structural resemblance also suggest that a drug developed either for SARS-CoV1 or SARS-CoV2 using these enzymes may also have the potential to cross-react.Communicated by Ramaswamy H. Sarma.

Exploring potential inhibitor of SARS-CoV2 replicase from FDA approved drugs using insilico drug discovery methods.

Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV2) is responsible for fetal pneumonia called COVID19. SARS-CoV2 emerged in Wuhan, Hubei Province of China in December 2019. The COVID19 pandemic has now gripped the entire world with more than 70 million cases and over 1.5 million deaths so far. There no treatment option for COVID19 is in term of a drug or vaccine is currently available. Therefore drug repurposing may only provide a quick method for utilizing existing drugs for a therapeutic option. The virus genome contains several non-structural proteins (NSP) which serve as target for designing of antiviral agents. NSP9 of SARS-CoV2 encodes for a replicase enzyme which is essential for the virus replication in the host cell. In search of potent inhibitors, we have screened FDA approved drugs against NSP9 using in silico methods. Five drugs fluspirilene, troglitazone, alvesco, dihydroergotoxine and avodart were found to have highest affinities with the replicase. The molecular dynamics simulation (MDS) studies demonstrated strong drugs binding and stable NSP9-drugs complexes formation. The findings are also strongly supported by root-mean-square deviation, root-mean-square fluctuation, radius of gyration, and hydrogen bond analysis of the complexes. Principal component analysis showed the stable conformation of NSP9 upon drug binding. It could be inferred that these five drugs individually or in combinations may be used as potential inhibitors of NSP9 of SARS-CoV-2 after exploring their in vivo antiviral potential.Communicated by Ramaswamy H. Sarma.

In silico identification and validation of natural antiviral compounds as potential inhibitors of SARS-CoV-2 methyltransferase.

The novel Coronavirus disease 2019 (COVID-19) is potentially fatal and caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Due to the unavailability of any proven treatment or vaccination, the outbreak of COVID-19 is wreaking havoc worldwide. Hence, there is an urgent need for therapeutics targeting SARS-CoV-2. Since, botanicals are an important resource for several efficacious antiviral agents, natural compounds gaining significant attention for COVID-19 treatment. In the present study, methyltranferase (MTase) of the SARS-CoV-2 is targeted using computational approach. The compounds were identified using molecular docking, virtual screening and molecular dynamics simulation studies. The binding mechanism of each compound was analyzed considering the stability and energetic parameter using in silico methods. We have found four natural antiviral compounds Amentoflavone, Baicalin, Daidzin and Luteoloside as strong inhibitors of methyltranferase of SARS-CoV-2. ADMET prediction and target analysis of the selected compounds showed favorable results. MD simulation was performed for four top-scored molecules to analyze the stability, binding mechanism and energy requirements. MD simulation studies indicated energetically favorable complex formation between MTase and the selected antiviral compounds. Furthermore, the structural effects on these substitutions were analyzed using the principles of each trajectories, which validated the interaction studies. Our analysis suggests that there is a very high probability that these compounds may have a good potential to inhibit Methyltransferase (MTase) of SARS-CoV-2 and to be used in the treatment of COVID-19. Further studies on these natural compounds may offer a quick therapeutic choice to treat COVID-19.Communicated by Ramaswamy H. Sarma.

Structure-activity relationship (SAR) and molecular dynamics study of withaferin-A fragment derivatives as potential therapeutic lead against main protease (Mpro) of SARS-CoV-2.

The spread of novel coronavirus SARS-CoV-2 has directed to a state of an unprecedented global pandemic. Many synthetic compounds and FDA-approved drugs have been significantly inhibitory against the virus, but no SARS-CoV-2 solution has been identified. However, small molecule fragment-based derivatives of potent phytocompounds may serve as promising inhibitors against SARS-CoV-2. In the pursuit of exploring novel SARS-CoV-2 inhibitors, we generated small molecule fragment derivatives from potent phytocompounds using neural networking and machine learning-based tools, which can cover unexplored regions of the chemical space that still retain lead-like properties. Out of 300 derivative molecules from withaferin-A, hesperidin, and baicalin, 30 were screened out with synthetic accessibility scores > 4 having the best ADME properties. The withaferin-A derivative molecules 61 and 64 exhibited a significant binding affinity of - 7.84 kcal/mol and - 7.94 kcal/mol. The docking study reveals that withaferin-A mol 61 forms 5 polar H-bonds with the Mpro where amino acids involved are GLU166, THR190, CYS145, MET165, and GLN152 and upon QSAR analysis showed a minimal predicted IC50 value of 7762.47 nM. Furthermore, the in silico cytotoxicity predictions, pharmacophore modeling, and molecular dynamics simulation studies have resulted in predicting the highly potent small molecule derivative from withaferin-A (phytocompound from Withania somnifera) to be the potential inhibitor of SARS-CoV 2 protease (Mpro) and a promising future lead candidate against COVID-19. The rationale of choosing withaferin-A from Withania somnifera (Ashwagandha) was propelled by the innumerous applications of Ashwagandha for the treatment of various antiviral diseases, common cold, and fever since time immemorial. Graphical abstract.

Identification of potential inhibitors of SARS-COV-2 endoribonuclease (EndoU) from FDA approved drugs: a drug repurposing approach to find therapeutics for COVID-19.

SARS-CoV-2 is causative agent of COVID-19, which is responsible for severe social and economic disruption globally. Lack of vaccine or antiviral drug with clinical efficacy suggested that drug repurposing approach may provide a quick therapeutic solution to COVID-19. Nonstructural protein-15 (NSP15) encodes for an uridylate-specific endoribonuclease (EndoU) enzyme, essential for virus life cycle and an attractive target for drug development. We have performed in silico based virtual screening of FDA approved compounds targeting EndoU in search of COVID-19 drugs from commercially available approved molecules. Two drugs Glisoxepide and Idarubicin used for treatment for diabetes and leukemia, respectively, were selected as stronger binder of EndoU. Both the drugs bound to the active site of the viral endonuclease by forming attractive intermolecular interactions with catalytically essential amino acid residues, His235, His250, and Lys290. Molecular dynamics simulation studies showed stable conformation dynamics upon drugs binding to endoU. The binding free energies for Glisoxepide and Idarubicin were calculated to be -141 ± 11 and -136 ± 16 kJ/mol, respectively. The IC50 were predicted to be 9.2 µM and 30 µM for Glisoxepide and Idarubicin, respectively. Comparative structural analysis showed the stronger binding of EndoU to Glisoxepide and Idarubicin than to uridine monophosphate (UMP). Surface area calculations showed buried are of 361.8Å2 by Glisoxepide which is almost double of the area occupied by UMP suggesting stronger binding of the drug than the ribonucleotide. However, further studies on these drugs for evaluation of their clinical efficacy and dose formulations may be required, which may provide a quick therapeutic option to treat COVID-19. Communicated by Ramaswamy H. Sarma.

Identification of hot spot residues on serine-arginine protein kinase-1 by molecular dynamics simulation studies.

Serine-arginine protein kinase-1 (SRPK1) is a highly specific kinase that recognizes serine-arginine dipeptide repeats and phosphorylates SR rich splicing factor ASF/SF2 in a cell-cycle regulated manner. SRPK1 processively phosphorylates serine residues on its substrate ASF/SF2. Elevated expression pattern of both SRPK1 and ASF/SF2 and their association with various carcinomas have established SRPK1 as a potent target for drug design against cancers. In order to develop specific inhibitors the binding of ASF/SF2 to SRPK1 is desired to be selectively interrupted. We have performed molecular dynamics simulation studies on crystal structure of SRPK1 complex with ASF/SF2. The ASF/SF2 acquired a stable binding on the surface of SRPK1 with strong attractive forces. Analysis revealed that there was no major position shifting of the core ß-sheet region within the catalytic site of SRPK1 when present in the state of ASF/SF2 bound in comparison to apo form. Global motions of SRPK1 indicated that major stable structural changes occurred after the substrate binding. The interactions between SRPK1 and ASF/SF2 were examined and calculated during molecular dynamics simulation of 1 µs. Molecular dynamics study indicated Arg84, Lys85, Leu86, Lys174, Tyr227 and Leu479 residues of SRPK1 as essential hot spots involved in the stable binding with substrate. Structural analysis of the binding affinity and hot spot investigation provided significant information on ASF/SF2 binding which may also be considered for designing of the novel specific inhibitors of SRPK1 for the applications in cancer therapy.Communicated by Ramaswamy H. Sarma.

Identification of a novel and potent small molecule inhibitor of SRPK1: mechanism of dual inhibition of SRPK1 for the inhibition of cancer progression.

Protein kinases are the family of attractive enzyme targets for drug design with relevance to cancer biology. Serine arginine protein kinase 1 (SRPK1) is responsible for the phosphorylation of serine/arginine (SR)-rich proteins. Alternative Splicing Factor/Splicing Factor 2 (ASF/SF2) involved in mRNA editing. ASF/SF2 is over expressed in many cancers and plays crucial roles in the cell survival. Phosphorylation of ASF/SF2 is decisive for its functions in cancer. In search of potential anticancer therapeutic agents for attenuating phosphorylation of ASF/SF2, we have explored specific and potential inhibitors of SRPK1 from natural and drug like compounds databases using in-silico methods. Compound ZINC02154892 (C02) was found to be the most potent inhibitor for SRPK1. In-vitro molecular and cell biology studies have shown C02 as a potent and specific inhibitor of phosphorylation of ASF/SF2 and cell survival in leukemic cell line. Structural analysis of SRPK1 with compound C02 revealed a unique pattern of binding targeting ATP binding site along with inhibiting recruitment of ASF/SF2 by SRPK1. The possibilities of compound C02 to be used as a lead compound paving way for the development of potent and specific inhibitors of SRPK1 for designing of novel potential anticancer inhibitor is inferred from the current studies.

Structural studies on dihydrouridine synthase A (DusA) from Pseudomonas aeruginosa.

Dihydrouridination is one of the abundant modifications in tRNA editing. The presence of dihydrouridine is attributed to tRNA stability desired for the efficient gene translation process. The conversion of uridine to dihydrouridine is catalyzed by flavine containing enzyme called dihydrouridine synthase (Dus). We report first ever information about DusA enzyme from Pseudomonas aeruginosa in form of structural and functional studies. The gene coding for DusA from P. aeruginosa (PADusA) was cloned, expressed and purified, using recombinant DNA technology methods. Thermal and chemical stability of PADusA was determined with respect to temperature and urea-induced equilibrium unfolding experiments, with monitoring the change of ellipticity at 200-260 nm by Circular Dichroism (CD) spectroscopy. Unfolding studies revealed that PADusA has acquired a stable tertiary structure fold with a Tm value of 46.2 °C and Cm of 2.7 M for urea. The enzyme contains 43% α-helices and 16% ß-strands. The three dimensional structure of PADusA was modeled using insilico methods. In order to understand the mechanism of substrate recognition and catalysis, tRNA and puromycin were docked on PADusA structure and their binding was analyzed. The structural features suggested that PADusA may also form a novel target for structure based drug design of antimicrobial agents.