An in silico pharmacological approach toward the discovery of potent inhibitors to combat drug resistance HIV‐1 protease variants

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

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

Drug repurposing toward the inhibition of RNA-dependent RNA polymerase of various flaviviruses through computational study

Numerous pathogens affecting human is present in the flavivirus family namely west nile, dengue, yellow fever, and zika which involves in development of global burden and distressing the environment economically. Till date, no approved drugs are available for targeting these viruses. The threat which urged the identification of small molecules for the inhibition of these viruses is the spreading of serious viral diseases. The recent outbreak of zika and dengue infections postured a solemn risk to worldwide public well-being. RNA-dependent RNA polymerase (RdRp) is the supreme adaptable enzymes of all the RNA viruses which is responsible for the replication and transcription of genome among the structural and nonstructural proteins of flaviviruses. It is understood that the RdRp of the flaviviruses are similar stating that the japanese encephalitis and west nile shares 70% identity with zika whereas the dengue serotype 2 and 3 shares the identity of 76% and 81%, respectively. In this study, we investigated the binding site of four flaviviral RdRp and provided insights into various interaction of the molecules using the computational approach. Our study helps in recognizing the potent compounds that could inhibit the viral protein as a common inhibitor. Additionally, with the conformational stability analysis, we proposed the possible mechanism of inhibition of the identified common small molecule toward RdRp of flavivirus. Finally, this study could be an initiative for the identification of common inhibitors and can be explored further for understanding the mechanism of action through in vitro studies for the study on efficacy.

Molecular Dynamics Study on Selected Bioactive Phytochemicals as Potential Inhibitors of HIV-1 Subtype C Protease

Acquired immunodeficiency syndrome (AIDS), one of the deadliest global diseases, is caused by the Human Immunodeficiency Virus (HIV). To date, there are no known conventional drugs that can cure HIV/AIDS, and this has prompted continuous scientific efforts in the search for novel and potent anti-HIV therapies. In this study, molecular dynamics simulation (MDS) and computational techniques were employed to investigate the inhibitory potential of bioactive compounds from selected South African indigenous plants against HIV-1 subtype C protease (HIVpro). Of the eight compounds (CMG, MA, UA, CA, BA, UAA, OAA and OA) evaluated, only six (CMG (-9.9 kcal/mol), MA (-9.3 kcal/mol), CA (-9.0 kcal/mol), BA (-8.3 kcal/mol), UAA (-8.5 kcal/mol), and OA (-8.6 kcal/mol)) showed favourable activities against HIVpro and binding landscapes like the reference FDA-approved drugs, Lopinavir (LPV) and Darunavir (DRV), with CMG and MA having the highest binding affinities. Using the structural analysis (root-mean-square deviation (RMSD), fluctuation (RMSF), and radius of gyration (RoG) of the bound complexes with HIVpro after 350 ns, structural evidence was observed, indicating that the six compounds are potential lead candidates for inhibiting HIVpro. This finding was further corroborated by the structural analysis of the enzyme-ligand complexe systems, where structural mechanisms of stability, flexibility, and compactness of the study metabolites were established following binding with HIVpro. Furthermore, the ligand interaction plots revealed that the metabolites interacted hydrophobically with the active site amino residues, with identification of other key residues implicated in HIVpro inhibition for drug design. Overall, this is the first computational report on the anti-HIV-1 activities of CMG and MA, with efforts on their in vitro and in vivo evaluations underway. Judging by the binding affinity, the degree of stability, and compactness of the lead metabolites (CMG, MA, CA, BA, OA, and UAA), they could be concomitantly explored with conventional HIVpro inhibitors in enhancing their therapeutic activities against the HIV-1 serotype.

Atom-based 3D-QSAR and DFT analysis of 5-substituted 2-acylaminothiazole derivatives as HIV-1 latency-reversing agents

HIV-1 latency consists of viral DNA; integrated inside the host genome; it remains transcriptional silent. Combined Antiretroviral Therapy (cART) and the host immune system fail to recognize the latency cells or reservoirs, representing a major barrier to eradicating the HIV-1 infection. The Shock and Kill emerged as a promising strategy to target these cells using Latency reversal agents (LRAs); partially succeeded in producing viral mRNA but failed to reduce the size of reservoirs. In this Context, 2-acylaminothiazole class derivatives appeared as promising HIV-1 latency-reversing agents. In this study, we have developed an atom-based 3 D-QSAR model by utilizing the 49 active compounds of the 5-substituted 2-acylaminothiazoles derivatives. These compounds are further randomly divided into training (37) and test (12) datasets, yielding statistically significant R² (0.90) and Q² (0.85) results, respectively. The internal and external validation of the model shows highly robust and reliable results. Next, the model was visualized to check the favourable and unfavourable groups in terms of hydrogen bond donor, electron-withdrawing and hydrophobic group on the most active compound 96 and least active compound 30. The investigated model reveals the structural insights required for obtaining more leads that are potent. Finally, DFT calculations on the most and least active compounds were performed to support the atom-based 3 D-QSAR model. Overall, this study will aid in understanding the minimum structural requirement and functional group required for screening the novel potent leads as HIV-1 latency reversal agents.Communicated by Ramaswamy H. Sarma.

Structural insights into the RNA interaction with Yam bean Mosaic virus (coat protein) from Pachyrhizus erosus using bioinformatics approach

Plants are constantly threatened by a virus infection, i.e., Potyviruses, the second largest genus of plant viruses which results in several million-dollar losses in various essential crops globally. Yam bean (Pachyrhizus erosus) is considered to be one of the essential tuberous legume crops holding a great potential source of starch. Yam Bean Mosaic Virus (YBMV) of Potyvirus group belonging to the family potyviridae affects Yam bean and several angiosperms both in the tropical and sub-tropical regions causing large economical losses in crops. In this study, we attempted to understand the sequence-structure relationship and mode of RNA binding mechanism in YBMV CP using in silico integrative modeling and all-atoms molecular dynamics (MD) simulations. The assembly of coat protein (CP) subunits from YBMV and the plausible mode of RNA binding were compared with the experimental structure of CP from Watermelon mosaic virus potyvirus (5ODV). The transmembrane helix region is present in the YBMV CP sequence ranging from 76 to 91 amino acids. Like the close structural-homolog, 24 CPs monomeric sub-units formed YBMV a conserved fold. Our computational study showed that ARG124, ARG155, and TYR151 orient towards the inner side of the virion, while, THR122, GLN125, SER92, ASP94 reside towards the outer side of the virion. Despite sharing very low sequence similarity with CPs from other plant viruses, the strongly conserved residues Ser, Arg, and Asp within the RNA binding pocket of YBMV CP indicate the presence of a highly conserved RNA binding site in CPs from different families. Using several bioinformatics tools and comprehensive analysis from MD simulation, our study has provided novel insights into the RNA binding mechanism in YBMV CP. Thus, we anticipate that our findings from this study will be useful for the development of new therapeutic agents against the pathogen, paving the way for researchers to better control this destructive plant virus.

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.

E-pharmacophore based screening to identify potential HIV-1 gp120 and CD4 interaction blockers for wild and mutant types

HIV-1 gp120 provides a multistage viral entry process through the conserved CD4 binding site. Hunting of potential blockers can diminish the interaction of gp120 with the CD4 host receptor leading to the suppression of HIV-1 infection. Structure-based pharmacophore virtual screening followed by binding free energy calculation, molecular dynamics (MD) simulation and density functional theory (DFT) calculation is applied to discriminate the potential blockers from six small molecule databases. Five compounds from six databases exhibited vital interactions with key residues ASP368, GLU370, ASN425, MET426, TRP427 and GLY473 of gp120, involved in the binding with CD4, host receptor. Most importantly, compound NCI-254200 displayed strong communication with key residues of wild type and drug resistance single mutant gp120 (M426L and W427V) even in the dynamic condition, evidenced from MD simulation. This investigation provided a potential compound NCI-254200 which may show inhibitory activity against HIV-1 gp120 variant interactions with CD4 host cell receptors.

Identification of immucillin analogue natural compounds to inhibit Helicobacter pylori MTAN through high throughput virtual screening and molecular dynamics simulation

Abstract

One in every two humans is having Helicobacter pylori (H. pylori) in stomach causing gastric ulcer. Emergence of several drugs in eliminating H. pylori has paved way for emergence of multidrug resistance in them. This resistance is thriving and thereby necessitating the need of a potent drug. Identifying a potential target for medication is crucial. Bacterial 5'-methylthioadenosine/S-enosyl homocysteine nucleosidase (MTAN) is a multifunctional enzyme that controls seven essential metabolic pathways. It functions as a catalyst in the hydrolysis of the N-ribosidic bond of adenosine-based metabolites: S-adenosylhomocysteine (SAH), 5'-methylthioadenosine (MTA), 5'-deoxyadenosine (5'-DOA), and 6-amino-6-deoxyfutalosine. H. pylori unlike other bacteria and humans utilises an alternative pathway for menaquinone synthesis. It utilises Futosiline pathway for menaquinone synthesis which are obligatory component in electron transport pathway. Therefore, the enzymes functioning in this pathway represent them-self as a point of attack for new medications. We targeted MTAN protein of H. pylori to find out a potent natural hit to inhibit its growth. A comparative analysis was made with potent H. pylori MTAN (HpMTAN) known inhibitor, 5'-butylthio-DADMe-Immucillin-A (BuT-DADMe-ImmA) and ZINC natural subset database. Optimized ligands from the ZINC natural database were virtually screened using ligand based pharmacophore hypothesis to obtain the most efficient and potent inhibitors for HpMTAN. The screened leads were evaluated for their therapeutic likeness. Furthermore, the ligands that passed the test were subjected for MM-GBSA with MTAN to reveal the essential features that contributes selectivity. The results showed that Van der Waals contributions play a central role in determining the selectivity of MTAN. Molecular dynamics (MD) studies were carried out for 100 ns to assess the stability of ligands in the active site. MD analysis showed that binding of ZINC00490333 with MTAN is stable compared to reference inhibitor molecule BuT-DADMe-ImmA. Among the natural inhibitors screened after various docking procedures ZINC00490333 has highest binding score for HpMTAN (- 13.987). The ZINC inhibitor was successful in reproducing the BuT-DADMe-ImmA interactions with HpMTAN. Hence we suggest that ZINC00490333 compound may represent as a good lead in designing novel potent inhibitors of HpMTAN. This in silico approach indicates the potential of this molecule for advancing a further step in gastric ulcer treatment.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00081-2.

Computational analysis identifies druggable mutations in human rBAT mediated Cystinuria

Culex quinquefasciatus Cqm1 protein acts as the receptor for Lysinibacillus sphaericus mosquito-larvicidal binary (BinAB) toxin that is used worldwide for mosquito control. We found amino acid transporter protein, rBAT, as phylogenetically closest Cqm1 homolog in humans. The present study reveals large evolutionary distance between Cqm1 and rBAT, and rBAT ectodomain lacks the sequence motif which serves as binding-site for the BinAB toxin. Thus, BinAB toxin can be expected to remain safe for humans. rBAT (heavy subunit; SLC3A1) and catalytic b^0,+AT (light subunit; SLC7A9), linked by single disulfide bond, mediate renal reabsorption of cystine and dibasic amino acids in Na⁺ independent manner. Mutations in rBAT cause type I Cystinuria disease which shows global prevalence, and rBAT can be thought as an important pharmacological target. However, 3D structures of rBAT and b^0,+AT, the two components of b^0,+ heteromeric amino acid transporter systems, are not available. We constructed a reliable homology model of rBAT using Cqm1 coordinates and that of transmembrane b^0,+AT subunit using LAT1 coordinates. Mapping of pathogenic mutations onto rBAT ectodomain revealed their scattered distribution throughout the rBAT protein. Further, our computational simulations-based scoring of several known deleterious mutations of rBAT revealed that mutations those do not compromise the protein fold and stability, are localized on the same face of the molecule. These residues are expected to interact with the b^0,+AT transporter. The present study thus identifies druggable sites on rBAT that could be targeted for the treatment of type I Cystinuria.Communicated by Ramaswamy H. Sarma.

Current Computational Approaches for the Development of Anti-HIV Inhibitors: An Overview

BACKGROUND

Today, HIV-1 infection has become an extensive problem to public health and a greater challenge to all working researchers throughout the world. Since the beginning of HIV-1 virus, several antiviral therapeutic agents have been developed at various stages to combat HIV-1 infection. But, many of antiviral drugs are on the platform of drug resistance and toxicology issues, needs an urgent constructive investigation for the development of productive and protective therapeutics to make an improvement of individual life suffering with viral infection. As developing a novel agent is very costly, challenging and time taking route in the recent times.

METHODS

The review summarized about the modern approaches of computational aided drug discovery to developing a novel inhibitor within a short period of time and less cost.

RESULTS

The outcome suggests on the premise of reported information that the computational drug discovery is a powerful technology to design a defensive and fruitful therapeutic agents to combat HIV-1 infection and recover the lifespan of suffering one.

CONCLUSION

Based on survey of the reported information, we concluded that the current computational approaches is highly supportive in the progress of drug discovery and controlling the viral infection.

An Innovative Sequence-to-Structure-Based Approach to Drug Resistance Interpretation and Prediction: The Use of Molecular Interaction Fields to Detect HIV-1 Protease Binding-Site Dissimilarities

In silico methodologies have opened new avenues of research to understanding and predicting drug resistance, a pressing health issue that keeps rising at alarming pace. Sequence-based interpretation systems are routinely applied in clinical context in an attempt to predict mutation-based drug resistance and thus aid the choice of the most adequate antibiotic and antiviral therapy. An important limitation of approaches based on genotypic data exclusively is that mutations are not considered in the context of the three-dimensional (3D) structure of the target. Structure-based in silico methodologies are inherently more suitable to interpreting and predicting the impact of mutations on target-drug interactions, at the cost of higher computational and time demands when compared with sequence-based approaches. Herein, we present a fast, computationally inexpensive, sequence-to-structure-based approach to drug resistance prediction, which makes use of 3D protein structures encoded by input target sequences to draw binding-site comparisons with susceptible templates. Rather than performing atom-by-atom comparisons between input target and template structures, our workflow generates and compares Molecular Interaction Fields (MIFs) that map the areas of energetically favorable interactions between several chemical probe types and the target binding site. Quantitative, pairwise dissimilarity measurements between the target and the template binding sites are thus produced. The method is particularly suited to understanding changes to the 3D structure and the physicochemical environment introduced by mutations into the target binding site. Furthermore, the workflow relies exclusively on freeware, making it accessible to anyone. Using four datasets of known HIV-1 protease sequences as a case-study, we show that our approach is capable of correctly classifying resistant and susceptible sequences given as input. Guided by ROC curve analyses, we fined-tuned a dissimilarity threshold of classification that results in remarkable discriminatory performance (accuracy ≈ ROC AUC ≈ 0.99), illustrating the high potential of sequence-to-structure-, MIF-based approaches in the context of drug resistance prediction. We discuss the complementarity of the proposed methodology to existing prediction algorithms based on genotypic data. The present work represents a new step toward a more comprehensive and structurally-informed interpretation of the impact of genetic variability on the response to HIV-1 therapies.

Constitutive Inflammatory Cytokine Storm: A Major Threat to Human Health

Immune cells secrete small protein molecules that aim for cell-cell communications. These small molecules are called cytokines. Targeting cancer cells with administration of bispecific antibodies and natural extracts results in elevated circulating levels of inflammatory cytokines, including interferon-γ and interleukin (IL)-6, which lead to cell toxicity. Sustained release of cytokines due to immunotherapy or hormonal issues causes various diseases. Novel T cell-engaging therapies and monoclonal antibodies cause cytokine release syndrome. Efforts are being carried out to maximize the chance for therapeutic benefit from immunotherapy while minimizing the risk for life-threatening complications of sustained cytokine release. Neurodegeneration and cardiac diseases are the prominent diseases caused by inflammatory cytokines. The phenomenon is called cytokine storm. Cytokines can act antagonistically or synergistically. Constitutive expression of proinflammatory cytokines such as IL-3 and IL-6 causes organ damage and unbearable pain. In this review, we will discuss the regulators of cytokine release, its types, its implications on human health, and treatment.