Recent Trends in Drug Design and Discovery.

INTRODUCTION: Structure-based drug design is a wide area of identification of selective inhibitors of a target of interest. From the time of the availability of three dimensional structure of the drug targets, mostly the proteins, many computational methods had emerged to address the challenges associated with drug design process. Particularly, drug-likeness, druggability of the target protein, specificity, off-target binding, etc., are the important factors to determine the efficacy of new chemical inhibitors. OBJECTIVE: The aim of the present research was to improve the drug design strategies in field of design of novel inhibitors with respect to specific target protein in disease pathology. Recent statistical machine learning methods applied for structural and chemical data analysis had been elaborated in current drug design field. METHODS: As the size of the biological data shows a continuous growth, new computational algorithms and analytical methods are being developed with different objectives. It covers a wide area, from protein structure prediction to drug toxicity prediction. Moreover, the computational methods are available to analyze the structural data of varying types and sizes of which, most of the semi-empirical force field and quantum mechanics based molecular modeling methods showed a proven accuracy towards analysing small structural data sets while statistics based methods such as machine learning, QSAR and other specific data analytics methods are robust for large scale data analysis. RESULTS: In this present study, the background has been reviewed for new drug lead development with respect specific drug targets of interest. Overall approach of both the extreme methods were also used to demonstrate with the plausible outcome. CONCLUSION: In this chapter, we focus on the recent developments in the structure-based drug design using advanced molecular modeling techniques in conjunction with machine learning and other data analytics methods. Natural products based drug discovery is also discussed.

Identification of Potential Inhibitors for Targets Involved in Dengue Fever.

Lethality due to dengue infection is a global threat. Nearly 400 million people are affected every year, which approximately costs 500 million dollars for surveillance and vector control itself. Many investigations on the structure-function relationship of proteins expressed by the dengue virus are being made for more than a decade and had come up with many reports on small molecule drug discovery. In this review, we present a detailed note on viral proteins and their functions as well as the inhibitors discovered/designed so far using experimental and computational methods. Further, the phytoconstituents from medicinal plants, specifically the extract of the papaya leaves, neem and bael, which combat dengue infection via dengue protease, helicase, methyl transferase and polymerase are summarized.

Role of disordered regions in transferring tyrosine to its cognate tRNA.

Aminoacyl tRNA synthetase (AARS) plays an important role in transferring each amino acid to its cognate tRNA. Specifically, tyrosyl tRNA synthetase (TyrRS) is involved in various functions including protection from DNA damage due to oxidative stress, protein synthesis and cell signaling and can be an attractive target for controlling the pathogens by early inhibition of translation. TyrRS has two disordered regions, which lack a stable 3D structure in solution, and are involved in tRNA synthetase catalysis and stability. One of the disordered regions undergoes disorder-to-order transition (DOT) upon complex formation with tRNA whereas the other remains disordered (DR). In this work, we have explored the importance of these disordered regions using molecular dynamics simulations of both free and RNA-complexed states. We observed that the DOT and DR regions of the first subunit acts as a flap and interact with the acceptor arm of the tRNA. The DOT-DR flap closes when tyrosine (TyrRSTyr) is present at the active site of the complex and opens in the presence of tyrosine monophosphate (TyrRSYMP). The DOT and DR regions of the second subunit interact with the anticodon stem as well as D-loop of the tRNA, which might be involved in stabilizing the complex. The anticodon loop of the tRNA binds to the structured region present in the C-terminal of the protein, which is observed to be flexible during simulations. Detailed energy calculations also show that TyrRSTyr complex has stronger binding energy between tRNA and protein compared to TyrRSYMP; on the contrary, the anticodon is strongly bound in TyrRSYMP. The results obtained in the present study provide additional insights for understanding catalysis and the involvement of disordered regions in Tyr transfer to cognate tRNA.

Identification of novel NAD(P)H dehydrogenase [quinone] 1 antagonist using computational approaches.

NAD(P)H: quinone oxidoreductase 1 (NQO1) inhibitors are proved as promising therapeutic agents against cancer. This study is to determine potent NAD(P)H-dependent NQO1 inhibitors with new scaffold. Pharmacophore-based three-dimensional (3D) QSAR model has been built based on 45 NQO1 inhibitors reported in the literature. The structure-function correlation coefficient graph represents the relationship between phase activity and phase predicted activity for training and test sets. A QSAR model statistics shows the excellent correlation of the generated model. Pharmacophore hypothesis (AARR) yielded a statistically significant 3D QSASR model with a correlation coefficient of r2 = 0.99 as well as an excellent predictive power. From the analysis of pharmacophore-based virtual screening using by SPEC database, 4093 hits were obtained and were further filtered using virtual screening filters (HTVS, SP, XP) through structure based molecular docking. Based on glide energy and docking score, seven lead compounds show better binding affinity compared to the co-crystal inhibitor. The results of induced fit docking and prime/MM-GBSA suggest that leads AN-153/J117103 and AT-138/KB09997 binding with the catalytic site. Further, to understanding the stability of identified lead compounds MD simulations were done. The lead AN-153/J117103 showed the strong binding stable of the protein-ligand complex. Also the computed drug likeness reveals potential of this compound to treat cancer. AbbreviationsNQO1NAD(P)H-quinine oxidoreductase 1CPHcommon pharmacophore hypothesisPLSpartial least squireHBDhydrogen bond donorSDstandard deviationXPextra precisionIFDinduced fit dockingMM-GBSAmolecular mechanics generalized born surface areaMDSmolecular dynamics simulationRMSDroot mean square deviationRMSFroot mean square fluctuationRMSEroot mean square errorADMEabsorption distribution metabolism excretionsCommunicated by Ramaswamy H. Sarma.

A prospective compound screening contest identified broader inhibitors for Sirtuin 1.

Potential inhibitors of a target biomolecule, NAD-dependent deacetylase Sirtuin 1, were identified by a contest-based approach, in which participants were asked to propose a prioritized list of 400 compounds from a designated compound library containing 2.5 million compounds using in silico methods and scoring. Our aim was to identify target enzyme inhibitors and to benchmark computer-aided drug discovery methods under the same experimental conditions. Collecting compound lists derived from various methods is advantageous for aggregating compounds with structurally diversified properties compared with the use of a single method. The inhibitory action on Sirtuin 1 of approximately half of the proposed compounds was experimentally accessed. Ultimately, seven structurally diverse compounds were identified.

Plant latex thrombin-like cysteine proteases alleviates bleeding by bypassing factor VIII in murine model.

Hemostasis is a tightly regulated process which maintains a fluid state of blood within the vasculature and provides thrombotic response upon tissue injury. Various scientific studies have implicated the role of plant latex proteases in hemostasis using in vitro experiments. However, in vivo models substantiate their role in hemostasis. Therefore, in the present study, the effect of plant latex thrombin-like proteases (PTLPs) on hemostasis was investigated systematically using mice tail bleeding as a preclinical model. In this direction, latex protease fractions (LPFs), which showed potent thrombin-like activity, were selected as they act directly on fibrinogen to form clot and quickly stop bleeding. Thrombin-like activity was exhibited mainly by cysteine proteases. Calotropis gigantea, Carica papaya, Jatropha curcas, Oxystelma esculentum, Tabernaemontana divaricata, and Vallaris solanacea LPFs and papain from C. papaya latex significantly reduced bleeding on a topical application in normal and aspirin administered mice. In addition, PTLPs accelerated the clotting of factor VIII deficient plasma, while, papain brought back the clotting time to normal levels acting like a bypassing agent. Further, papain failed to show activity in the presence of specific cysteine protease inhibitor iodoacetic acid; confirming protease role in all the activities exhibited. At the tested dose, PTLPs except C. gigantea did not show toxicity. Further, structural and sequence comparison between PTLPs and human thrombin revealed structural and sequence dissimilarity indicating their unique nature. The findings of the present study may open up a new avenue for considering PTLPs including papain in the treatment of bleeding wounds.

The non-enzymatic RAS effector RASSF7 inhibits oncogenic c-Myc function.

c-Myc is a proto-oncogene controlling expression of multiple genes involved in cell growth and differentiation. Although the functional role of c-Myc as a transcriptional regulator has been intensively studied, targeting this protein in cancer remains a challenge. Here, we report a trimodal regulation of c-Myc function by the Ras effector, Ras-association domain family member 7 (RASSF7), a nonenzymatic protein modulating protein-protein interactions to regulate cell proliferation. Using HEK293T and HeLa cell lines, we provide evidence that RASSF7 destabilizes the c-Myc protein by promoting Cullin4B-mediated polyubiquitination and degradation. Furthermore, RASSF7 competed with MYC-associated factor X (MAX) in the formation of a heterodimeric complex with c-Myc and attenuated its occupancy on target gene promoters to regulate transcription. Consequently, RASSF7 inhibited c-Myc-mediated oncogenic transformation, and an inverse correlation between the expression levels of the RASSF7 and c-Myc genes was evident in human cancers. Furthermore, we found that RASSF7 interacts with c-Myc via its RA and leucine zipper (LZ) domains and LZ domain peptide is sufficient to inhibit c-Myc function, suggesting that this peptide might be used to target oncogenic c-Myc. These results unveil that RASSF7 and c-Myc are functionally linked in the control of tumorigenesis and open up potential therapeutic avenues for targeting the "undruggable" c-Myc protein in a subset of human cancers.

Identification of type I and type II inhibitors of c-Yes kinase using in silico and experimental techniques.

c-Yes kinase is considered as one of the attractive targets for anti-cancer drug design. The DFG (Asp-Phe-Gly) motif present in most of the kinases will adopt active and inactive conformations, known as DFG-in and DFG-out and their inhibitors are classified into type I and type II, respectively. In the present study, two screening protocols were followed for identification of c-Yes kinase inhibitors. (i) Structure-based virtual screening (SBVS) and (ii) Structure-based (SB) and Pharmacophore-based (PB) tandem screening. In SBVS, the c-Yes kinase structure was obtained from homology modeling and seven ensembles with different active site scaffolds through molecular dynamics (MD) simulations. For SB-PB tandem screening, we modeled ligand bound active and inactive conformations. Physicochemical properties of inhibitors of Src kinase family and c-Yes kinase were used to prepare target focused libraries for screenings. Our screening procedure along with docking showed 520 probable hits in SBVS and tandem screening (120 and 400, respectively). Out of 5000 compounds identified from different computational methods, 2410 were examined using kinase inhibition assays. It includes 266 compounds (5.32%) identified from our method. We observed that 14 compounds (12%) are identified by the present method out of 168 that showed > 30% inhibition. Among them, three compounds are novel, unique, and showed good inhibition. Further, we have studied the binding of these compounds at the DFG-in and DFG-out conformations and reported the probable class (type I or type II). Hence, we suggest that these compounds could be novel drug leads for regulation of colorectal cancer.

An iterative compound screening contest method for identifying target protein inhibitors using the tyrosine-protein kinase Yes.

We propose a new iterative screening contest method to identify target protein inhibitors. After conducting a compound screening contest in 2014, we report results acquired from a contest held in 2015 in this study. Our aims were to identify target enzyme inhibitors and to benchmark a variety of computer-aided drug discovery methods under identical experimental conditions. In both contests, we employed the tyrosine-protein kinase Yes as an example target protein. Participating groups virtually screened possible inhibitors from a library containing 2.4 million compounds. Compounds were ranked based on functional scores obtained using their respective methods, and the top 181 compounds from each group were selected. Our results from the 2015 contest show an improved hit rate when compared to results from the 2014 contest. In addition, we have successfully identified a statistically-warranted method for identifying target inhibitors. Quantitative analysis of the most successful method gave additional insights into important characteristics of the method used.

Structure-function relationship of Chikungunya nsP2 protease: A comparative study with papain.

Chikungunya virus is a growing human pathogen transmitted by mosquito bite. It causes fever, chills, nausea, vomiting, joint pain, headache, and swelling in the joints. Its replication and propagation depend on the protease activity of the Chikungunya virus-nsP2 protein, which cleaves the nsP1234 polyprotein replication complex into individual functional units. The N-terminal segment of papain is structurally identical with the Chikungunya virus-nsP2 protease. Hence, molecular dynamics simulations were performed to compare molecular mechanism of these proteases. The Chikungunya virus-snP2 protease shows more conformational changes and adopts an alternate conformation. However, N-terminal segment of these two proteases has identical active site scaffold with the conserved catalytic diad. Hence, some of the non-peptide inhibitors of papain were used for induced fit docking at the active site of the nsP2 to assess the binding mode. In addition, the peptides that connect different domains/protein in Chikungunya virus poly-protein were also subjected for docking. The overall results suggest that the active site scaffold is the same in both the proteases and a possibility exists to experimentally assess the efficacy of some of the papain inhibitors to inhibit the Chikungunya virus-nsP2.

Celastrol modulates inflammation through inhibition of the catalytic activity of mediators of arachidonic acid pathway: Secretory phospholipase A2 group IIA, 5-lipoxygenase and cyclooxygenase-2.

Elevated production of arachidonic acid (AA)-derived pro-inflammatory eicosanoids due to the concerted action of secretory phospholipase A2 group IIA (sPLA2IIA), 5-lipoxygenase (5-LOX) and cyclooxygenase-2 (COX-2) is a common feature of many inflammatory disorders. Hence, modulation of the bioactivity of these 3 enzymes is an important strategy to control inflammation. However, the failure of drugs specific for an individual enzyme (sPLA2IIA-, 5-LOX- or COX-2) and the success of 5-LOX/COX-2 dual inhibitors in effectively controlling inflammation in clinical trials prompted us to evaluate a common inhibitor for sPLA2IIA, 5-LOX and COX-2 enzymes. Celastrol, a quinone methide triterpene, was selected in this regard through molecular docking studies. We provide the first evidence for celastrol's ability to inhibit the catalytic activity of sPLA2IIA, 5-LOX and COX-2 enzymes. Celastrol significantly inhibited the catalytic activity of sPLA2IIA (IC50=6µM) in vitro, which is independent of substrate and calcium concentration. In addition, celastrol inhibited the catalytic activities of 5-LOX (IC50=5µM) and COX-2 (IC50=20µM) in vitro; sPLA2IIA-induced edema and carrageenan-induced edema in mice; and lipopolysaccharide-stimulated production of PGE2 in human neutrophils. Thus, celastrol modulates inflammatory responses by targeting multiple enzymes of AA pathway.

Dimethyl ester of bilirubin exhibits anti-inflammatory activity through inhibition of secretory phospholipase A2, lipoxygenase and cyclooxygenase.

Overproduction of arachidonic acid (AA) mediated by secretory phospholipase A2 group IIA (sPLA2IIA) is a hallmark of many inflammatory disorders. AA is subsequently converted into pro-inflammatory eicosanoids through 5-lipoxygenase (5-LOX) and cyclooxygenase-1/2 (COX-1/2) activities. Hence, inhibition of sPLA2IIA, 5-LOX and COX-1/2 activities is critical in regulating inflammation. We have previously reported unconjugated bilirubin (UCB), an endogenous antioxidant, as sPLA2IIA inhibitor. However, lipophilic UCB gets conjugated in liver with glucuronic acid into hydrophilic conjugated bilirubin (CB). Since hydrophobicity is pre-requisite for sPLA2IIA inhibition, conjugation reduces the efficacy of UCB. In this regard, UCB was chemically modified and derivatives were evaluated for sPLA2IIA, 5-LOX and COX-1/2 inhibition. Among the derivatives, BD1 (dimethyl ester of bilirubin) exhibited âˆ¼ 3 fold greater inhibitory potency towards sPLA2IIA compared to UCB. Both UCB and BD1 inhibited human 5-LOX and COX-2 activities; however only BD1 inhibited AA induced platelet aggregation. Molecular docking studies demonstrated BD1 as better inhibitor of aforesaid enzymes than UCB and other endogenous antioxidants. These data suggest that BD1 exhibits strong anti-inflammatory activity through inhibition of AA cascade enzymes which is of great therapeutic importance.

Identification of potential inhibitors based on compound proposal contest: Tyrosine-protein kinase Yes as a target.

A search of broader range of chemical space is important for drug discovery. Different methods of computer-aided drug discovery (CADD) are known to propose compounds in different chemical spaces as hit molecules for the same target protein. This study aimed at using multiple CADD methods through open innovation to achieve a level of hit molecule diversity that is not achievable with any particular single method. We held a compound proposal contest, in which multiple research groups participated and predicted inhibitors of tyrosine-protein kinase Yes. This showed whether collective knowledge based on individual approaches helped to obtain hit compounds from a broad range of chemical space and whether the contest-based approach was effective.

Progressive Hemorrhage and Myotoxicity Induced by Echis carinatus Venom in Murine Model: Neutralization by Inhibitor Cocktail of N,N,N',N'-Tetrakis (2-Pyridylmethyl) Ethane-1,2-Diamine and Silymarin.

Viperbite is often associated with severe local toxicity, including progressive hemorrhage and myotoxicity, persistent even after the administration of anti-snake venom (ASV). In the recent past, investigations have revealed the orchestrated actions of Zn2+ metalloproteases (Zn2+MPs), phospholipase A2s (PLA2s) and hyaluronidases (HYs) in the onset and progression of local toxicity from the bitten site. As a consequence, venom researchers and medical practitioners are in deliberate quest of potent molecules alongside ASV to tackle the brutal local manifestations induced by aforesaid venom toxins. Based on these facts, we have demonstrated the protective efficacy of inhibitor cocktail containing equal ratios of N,N,N',N'-tetrakis (2-pyridylmethyl) ethane-1,2-diamine (TPEN) and silymarin (SLN) against progressive local toxicity induced by Echis carinatus venom (ECV). In our previous study we have shown the inhibitory potentials of TPEN towards Zn2+MPs of ECV (IC50: 6.7 µM). In this study we have evaluated in vitro inhibitory potentials of SLN towards PLA2s (IC50: 12.5 µM) and HYs (IC50: 8 µM) of ECV in addition to docking studies. Further, we have demonstrated the protection of ECV induced local toxicity with 10 mM inhibitor cocktail following 15, 30 min (for hemorrhage and myotoxicity); 60 min (for hemorrhage alone) of ECV injection in murine model. The histological examination of skin and thigh muscle sections taken out from the site of ECV injection substantiated the overall protection offered by inhibitor cocktail. In conclusion, the protective efficacy of inhibitor cocktail is of high interest and can be administered locally alongside ASV to treat severe local toxicity.

Structure based approach for understanding organism specific recognition of protein-RNA complexes.

BACKGROUND: Protein-RNA interactions perform diverse functions within the cell. Understanding the recognition mechanism of protein-RNA complexes has been a challenging task in molecular and computational biology. In earlier works, the recognition mechanisms have been studied for a specific complex or using a set of non-redundant complexes. In this work, we have constructed 18 sets of same protein-RNA complexes belonging to different organisms from Protein Data Bank (PDB). The similarities and differences in each set of complexes have been revealed in terms of various sequence and structure based features such as root mean square deviation, sequence homology, propensity of binding site residues, variance, conservation at binding sites, binding segments, binding motifs of amino acid residues and nucleotides, preferred amino acid-nucleotide pairs and influence of neighboring residues for binding. RESULTS: We found that the proteins of mesophilic organisms have more number of binding sites than thermophiles and the binding propensities of amino acid residues are distinct in E. coli, H. sapiens, S. cerevisiae, thermophiles and archaea. Proteins prefer to bind with RNA using a single residue segment in all the organisms while RNA prefers to use a stretch of up to six nucleotides for binding with proteins. We have developed amino acid residue-nucleotide pair potentials for different organisms, which could be used for predicting the binding specificity. Further, molecular dynamics simulation studies on aspartyl tRNA synthetase complexed with aspartyl tRNA showed specific modes of recognition in E. coli, T. thermophilus and S. cerevisiae. CONCLUSION: Based on structural analysis and molecular dynamics simulations we suggest that the mode of recognition depends on the type of the organism in a protein-RNA complex.

Inhibitory potential of three zinc chelating agents against the proteolytic, hemorrhagic, and myotoxic activities of Echis carinatus venom.

Viperbites undeniably cause local manifestations such as hemorrhage and myotoxicity involving substantial degradation of extracellular matrix (ECM) at the site of envenomation and lead to progressive tissue damage and necrosis. The principle toxin responsible is attributed to snake venom metalloproteases (SVMPs). Treatment of such progressive tissue damage induced by SVMPs has become a challenging task for researchers and medical practitioners who are in quest of SVMPs inhibitors. In this study, we have evaluated the inhibitory potential of three specific zinc (Zn(2+)) chelating agents; N,N,N',N'-tetrakis (2-pyridylmethyl) ethane-1,2-diamine (TPEN), diethylene triamine pentaacetic acid (DTPA), tetraethyl thiuram disulfide (TTD) on Echis carinatus venom (ECV) induced hemorrhage and myotoxicity. Amongst them, TPEN has high affinity for Zn(2+) and revealed potent inhibition of ECV metalloproteases (ECVMPs) in vitro (IC50: 6.7 µM) compared to DTPA and TTD. The specificity of TPEN towards Zn(2+) was confirmed by spectral and docking studies. Further, TPEN, DTPA, and TTD completely blocked the hemorrhagic and myotoxic activities of ECV in a dose dependent manner upon co-injection; whereas, only TPEN successfully neutralized hemorrhage and myotoxicity following independent injection. Histological examinations revealed that TPEN effectively prevents degradation of dermis and basement membrane surrounding the blood vessels in mouse skin sections. TPEN also prevents muscle necrosis and accumulation of inflammatory cells at the site of ECV injections. In conclusion, a high degree of structural and functional homology between mammalian MMPs and SVMPs suggests that specific Zn(2+) chelators currently in clinical practice could be potent first aid therapeutic agents in snakebite management, particularly for local tissue damage.

Identification of novel inhibitors of Daboia russelli phospholipase A2 using the combined pharmacophore modeling approach.

Crystal structures available for Daboia russelli venom PLA(2) confirm that it undergoes dimerization with asymmetry and hence difference in the conformation of active site of the two subunits. The active site of subunit A is open and that of subunit B is closed. Pharmacophore models were generated based on the interaction of different types of inhibitors with their preferred subsites in the active site of subunit A. Particularly, the features responsible for recognizing subsites 1-3 and those of subsites 4-6 were combined as these two are involving in inflammation and anticoagulation processes, respectively. Pharmacophore model was edited to make the geometry suitable for the active site of both the subunits A and B. Final model is validated and subjected for screening a library of druglike compounds. Eight compounds were shortlisted and subjected for molecular docking and dynamics simulation to assess their binding mode with both the subunits. Based on the hydrophobic interactions and binding free energy, four compounds were selected for further biochemical assay. The overall results suggest that two compounds can bind both the subunits of PLA(2) of Daboia russelli venom in spite of its aggregated form and other two inhibit structurally very similar Naja naja PLA(2).