In-Silico molecular docking and simulation studies on novel chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage as vital inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase

The structural motifs of chalcones, flavones, and triazoles with varied substitutions have been studied for the antimalarial activity. In this study, 25 novel derivatives of chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage are docked with Plasmodium falciparum dihydroorotate dehydrogenase to establish their inhibitory activity against Plasmodium falciparum. The best binding conformation of the ligands at the catalytic site of dihydroorotate dehydrogenase are selected to characterize the best bound ligand using the best consensus score and the number of hydrogen bond interactions. The ligand namely (2E)-3-(4-{[1-(3-chloro-4-fluorophenyl)-1H-1, 2, 3-triazol-4-yl]methoxy}-3-methoxyphenyl-1-(2-hydroxy-4,6-dimethoxyphenyl)prop-2-en-1-one, is one the among the five best docked ligands, which interacts with the protein through nine hydrogen bonds and with a consensus score of five. To refine and confirm the docking study results, the stability of complexes is verified using Molecular Dynamics Simulations, Molecular Mechanics /Poisson–Boltzmann Surface Area free binding energy analysis, and per residue contribution for the binding energy. The study implies that the best docked Plasmodium falciparum dihydroorotate dehydrogenase–ligand complex is having high negative binding energy, most stable, compact, and rigid with nine hydrogen bonds. The study provides insight for the optimization of chalcone and flavone hybrids with 1, 2, 3-triazole linkage as potent inhibitors.     

Ligand‐based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore‐based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide‐quantum‐polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques. Copyright © 2014 John Wiley & Sons, Ltd.     

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment,hybrid QM/MM calculation,and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1′ pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work.

Communicated by Ramaswamy H. Sarma     

Assessment of Plasmodium falciparum PfMDR1 transport rates using Fluo‐4

Mutations in the multidrug resistance transporter of Plasmodium falciparum PfMDR1 have been implicated to play a significant role in the emergence of worldwide drug resistance, yet the molecular and biochemical mechanisms of this transporter are not well understood. Although it is generally accepted that drug resistance in P. falciparum is partly associated with PfMDR1 transport activity situated in the membrane of the digestive vacuole, direct estimates of the pump rate of this transport process in the natural environment of the intact host–parasite system have never been analysed. The fluorochrome Fluo‐4 is a well‐documented surrogate substrate of PfMDR1 and has been found to accumulate by actively being transported into the digestive vacuole of several parasitic strains. In the present study, we designed an approach to use Fluo‐4 fluorescence uptake as a measure of compartmental Fluo‐4 concentration accumulation in the different compartments of the host–parasite system. We performed a ‘reverse Fluo‐4 imaging' approach to relate fluorescence intensity to changes in dye concentration rather than Ca²⁺ fluctuations and were able to calculate the overall rate of transport for PfMDR1 in Dd2 parasites. With this assay, we provide a powerful method to selectively measure the effect of PfMDR1 mutations on substrate transport kinetics. This will be of high significance for future compound screening to test for new drugs in resistant P. falciparum strains.     

Predicting binding modes of reversible peptide‐based inhibitors of falcipain‐2 consistent with structure–activity relationships

Falcipain‐2 (FP‐2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide‐based inhibitors of FP‐2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP‐2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near‐native ligand conformations when applied to a validation set of protein‐ligand structures. Overall, we proposed substrate‐like binding modes of the studied compounds fulfilling the structural requirements for FP‐2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666–1683. © 2017 Wiley Periodicals, Inc.     

Molecular modeling reveals binding interface of γ‐tubulin with GCP4 and interactions with noscapinoids

The initiation of microtubule assembly within cells is guided by a cone shaped multi‐protein complex, γ‐tubulin ring complex (γTuRC) containing γ‐tubulin and atleast five other γ‐tubulin‐complex proteins (GCPs), i.e., GCP2, GCP3, GCP4, GCP5, and GCP6. The rim of γTuRC is a ring of γ‐tubulin molecules that interacts, via one of its longitudinal interfaces, with GCP2, GCP3, or GCP4 and, via other interface, with α/β?tubulin dimers recruited for the microtubule lattice formation. These interactions however, are not well understood in the absence of crystal structure of functional reconstitution of γTuRC subunits. In this study, we elucidate the atomic interactions between γ‐tubulin and GCP4 through computational techniques. We simulated two complexes of γ‐tubulin‐GCP4 complex (we called dimer1 and dimer2) for 25 ns to obtain a stable complex and calculated the ensemble average of binding free energies of ?158.82 and ?170.19 kcal/mol for dimer1 and ?79.53 and ?101.50 kcal/mol for dimer2 using MM‐PBSA and MM‐GBSA methods, respectively. These highly favourable binding free energy values points to very robust interactions between GCP4 and γ‐tubulin. From the results of the free‐energy decomposition and the computational alanine scanning calculation, we identified the amino acids crucial for the interaction of γ‐tubulin with GCP4, called hotspots. Furthermore, in the endeavour to identify chemical leads that might interact at the interface of γ‐tubulin‐GCP4 complex; we found a class of compounds based on the plant alkaloid, noscapine that binds with high affinity in a cavity close to γ‐tubulin‐GCP4 interface compared with previously reported compounds. All noscapinoids displayed stable interaction throughout the simulation, however, most robust interaction was observed for bromo‐noscapine followed by noscapine and amino‐noscapine. This offers a novel chemical scaffold for γ‐tubulin binding drugs near γ‐tubulin‐GCP4 interface. Proteins 2015; 83:827–843. © 2015 Wiley Periodicals, Inc.     

Exploration of new drug-like inhibitors for serine/threonine protein phosphatase 5 of Plasmodium falciparum: a docking and simulation study

Protein phosphorylation is an important mechanism that implicates in physiology of any organism including parasitic protozoa. Metallic protein Ser/Thr protein phosphatase 5 (PP5) controls various cellular signaling pathways of Plasmodium falciparum. The structure and inhibitory mechanism of PP5 in P. falciparum is not known. In fact, no experimental structural data are available for P. falciparum Ser/Thr protein phosphatase 5 (PfPP5) till date. Hence, we have proposed computer-generated model of catalytic subunit of PfPP5 and its inhibitory mechanism was analyzed. A set of 42 known natural inhibitors of protein phosphate family were docked against metal-binding catalytic site of PfPP5 and we found that cantharidin and its derivatives shows better binding energy among them. Similarity search was performed by taking these compounds as lead compounds against PubChem and ChemBank. The search result provides 3703 similar compounds; out of which 2245 qualified the Lipinski rule of five. Further, virtual screening of these compounds was performed and selected top 25 were selected on the basis of binding energy. In continuation, rigid and flexible docking of these screened compounds was performed to get the insight of interactions. Finally, top 5 compounds were verified for ADMET properties, and then, all are subjected to MD simulations for 25 ns in order to validate their stability. Compounds CBI: 3554182, CID: 23561913, and CID: 21168680 showed most stable binding, although some of hydrogen bonds pairing varied throughout simulation. These finding would be helpful to the medicinal chemists for the development of antimalarial drugs to combat this deadly disease.     

A compound‐based computational approach for the accurate determination of hot spots

A plethora of both experimental and computational methods have been proposed in the past 20 years for the identification of hot spots at a protein–protein interface. The experimental determination of a protein–protein complex followed by alanine scanning mutagenesis, though able to determine hot spots with much precision, is expensive and has no guarantee of success while the accuracy of the current computational methods for hot‐spot identification remains low. Here, we present a novel structure‐based computational approach that accurately determines hot spots through docking into a set of proteins homologous to only one of the two interacting partners of a compound capable of disrupting the protein–protein interaction (PPI). This approach has been applied to identify the hot spots of human activin receptor type II (ActRII) critical for its binding toward Cripto‐I. The subsequent experimental confirmation of the computationally identified hot spots portends a potentially accurate method for hot‐spot determination in silico given a compound capable of disrupting the PPI in question. The hot spots of human ActRII first reported here may well become the focal points for the design of small molecule drugs that target the PPI. The determination of their interface may have significant biological implications in that it suggests that Cripto‐I plays an important role in both activin and nodal signal pathways.     

Identification of Plasmodium falciparum apicoplast-targeted tRNA-guanine transglycosylase and its potential inhibitors using comparative genomics,molecular modelling,docking and simulation studies

tRNA modifications play an important role in the proper folding of tRNA and thereby determine its functionality as an adaptor molecule. Notwithstanding the centrality of this basic process in translation, a major gap in the genomics of Plasmodium falciparum is unambiguous identification of enzymes catalysing the various tRNA modifications. In this study, tRNA-modifying enzymes of P. falciparum were annotated using homology-based approach. Based on the presence of these identified enzymes, the modifications were compared with those of prokaryotic and eukaryotic organisms. Through sequence comparison and phylogenetic analysis, we have identified P. falciparum apicoplast tRNA-guanine 34 transglycosylase (TGT, EC: 2.4.2.29), which shows evidence of its prokaryotic origin. The docking analysis of the modelled TGT structures revealed that binding of quinazolinone derivatives is more favourable with P. falciparum apicoplast TGT as compared to human TGT. Molecular dynamic simulation and molecular mechanics/generalized Born surface area analysis of the complex confirmed the greater binding affinity of the ligand in the binding pocket of P. falciparum TGT protein. Further, evolutionary patterning analysis identified the amino acids of P. falciparum apicoplast TGT that are under purifying selection pressure and hence can be good inhibitor-targeting sites. Based on these computational studies, we suggest that P. falciparum apicoplast tRNA-guanine 34 transglycosylase can be a promising drug target.     

A microtitre chemiluminescence sensor for detection of pyrethroids based on dual‐dummy‐template molecularly imprinted polymer and computational simulation

The residues of pyrethroids in foods of animal origin are dangerous to the consumers, so this study presented a chemiluminescence sensor for determination of pyrethroids in chicken samples. A dual‐dummy‐template molecularly imprinted polymer capable of recognizing 10 pyrethroids was synthesized. The results of computation simulation showed that the specific 3D conformations of the templates had important influences on the polymer' recognition ability. The polymer was used to prepare a sensor on conventional 96‐well microplates, and the sample solution was added into the wells for direct absorption. The absorbed analytes were initiated with the bis(2,4,6‐trichlorophenyl)oxalate–H₂O₂–imidazole system, and the chemiluminescence intensity was used for analyte quantification. Results showed that one assay was finished within 12 min, and this sensor could be reused four times. The limits of detection for the 10 analytes were in the range o0.3–6.0 pg/ml, and the recoveries from the standards of fortified blank chicken samples were in the range 70.5–99.7%.     

Design of character‐based DNA barcode motif for species identification: A computational approach and its validation in fishes

The DNA barcodes are generally interpreted using distance‐based and character‐based methods. The former uses clustering of comparable groups, based on the relative genetic distance, while the latter is based on the presence or absence of discrete nucleotide substitutions. The distance‐based approach has a limitation in defining a universal species boundary across the taxa as the rate of mtDNA evolution is not constant throughout the taxa. However, character‐based approach more accurately defines this using a unique set of nucleotide characters. The character‐based analysis of full‐length barcode has some inherent limitations, like sequencing of the full‐length barcode, use of a sparse‐data matrix and lack of a uniform diagnostic position for each group. A short continuous stretch of a fragment can be used to resolve the limitations. Here, we observe that a 154‐bp fragment, from the transversion‐rich domain of 1367 COI barcode sequences can successfully delimit species in the three most diverse orders of freshwater fishes. This fragment is used to design species‐specific barcode motifs for 109 species by the character‐based method, which successfully identifies the correct species using a pattern‐matching program. The motifs also correctly identify geographically isolated population of the Cypriniformes species. Further, this region is validated as a species‐specific mini‐barcode for freshwater fishes by successful PCR amplification and sequencing of the motif (154 bp) using the designed primers. We anticipate that use of such motifs will enhance the diagnostic power of DNA barcode, and the mini‐barcode approach will greatly benefit the field‐based system of rapid species identification.     

Identification of acetylcholinesterase inhibitors using homogenous cell‐based assays in quantitative high‐throughput screening platforms

Acetylcholinesterase (AChE) is an enzyme responsible for metabolism of acetylcholine, a neurotransmitter associated with muscle movement, cognition, and other neurobiological processes. Inhibition of AChE activity can serve as a therapeutic mechanism, but also cause adverse health effects and neurotoxicity. In order to efficiently identify AChE inhibitors from large compound libraries, homogenous cell‐based assays in high‐throughput screening platforms are needed. In this study, a fluorescent method using Amplex Red (10‐acetyl‐3,7‐dihydroxyphenoxazine) and the Ellman absorbance method were both developed in a homogenous format using a human neuroblastoma cell line (SH‐SY5Y). An enzyme‐based assay using Amplex Red was also optimized and used to confirm the potential inhibitors. These three assays were used to screen 1368 compounds, which included a library of pharmacologically active compounds (LOPAC) and 88 additional compounds from the Tox21 program, at multiple concentrations in a quantitative high‐throughput screening (qHTS) format. All three assays exhibited exceptional performance characteristics including assay signal quality, precision, and reproducibility. A group of inhibitors were identified from this study, including known (e.g. physostigmine and neostigmine bromide) and potential novel AChE inhibitors (e.g. chelerythrine chloride and cilostazol). These results demonstrate that this platform is a promising means to profile large numbers of chemicals that inhibit AChE activity.     

A comparative study on the heparin‐binding proteomes of Toxoplasma gondii and Plasmodium falciparum

Toxoplasma gondii is an obligatory intracellular apicomplexan parasite which exploits host cell surface components in cell invasion and intracellular parasitization. Sulfated glycans such as heparin and heparan sulfate have been reported to inhibit cell invasion by T. gondii and other apicomplexan parasites such as Plasmodium falciparum. The aim of this study was to investigate the heparin‐binding proteome of T. gondii. The parasite‐derived components were affinity‐purified on the heparin moiety followed by MS fingerprinting of the proteins. The heparin‐binding proteins of T. gondii and P. falciparum were compared based on functionality and affinity to heparin. Among the proteins identified, the invasion‐related parasite ligands derived from tachyzoite/merozoite surface and the secretory organelles were prominent. However, the profiles of the proteins were different in terms of affinity to heparin. In T. gondii, the proteins with highest affinity to heparin were the intracellular components with functions of parasite development contrasted to that of P. falciparum, of which the rhoptry‐derived proteins were prominently identified. The profiling of the heparin‐binding proteins of the two apicomplexan parasites not only explained the mechanism of heparin‐mediated host cell invasion inhibition, but also, to a certain extent, revealed that the action of heparin on the parasite extended after endocytosis.     

Binding of B‐cell maturation antigen to B‐cell activating factor induces survival of multiple myeloma cells by activating Akt and JNK signaling pathways

B‐cell maturation antigen (BCMA) is expressed on normal and malignant plasma cells and represents a potential target for therapeutic intervention. In this study, we characterized the mechanism underlying the protein kinase B (Akt) and c‐Jun N‐terminal kinase (JNK) pathways and BCMA interactions in regulating multiple myeloma (MM) cell survival. It was found that the expression levels of B cell‐activating factor (BAFF) and BCMA were increased in MM cells as compared with those in normal controls. The proliferation of U266 cells was induced by recombinant human BAFF (rhBAFF) and could also be decreased by BCMA siRNA. The expression of Bcl‐2 protein was up‐regulated, and Bax protein was down‐regulated after rhBAFF treatment, which could be reversed by BCMA siRNA. Similarly, the protein p‐JNK and p‐Akt were activated by rhBAFF and could be changed by BCMA siRNA. In addition, the BCMA mRNA and protein expression levels were decreased after treatment with Akt and JNK pathway inhibitors. These results suggest that Akt and JNK pathways are involved in the regulation of BCMA. A novel BAFF/BCMA signalling pathway in MM may be a new therapeutic target for MM. Copyright © 2016 John Wiley & Sons, Ltd.     

Fullerene‐based inhibitors of HIV‐1 protease

A series of Fmoc‐Phe(4‐aza‐C₆₀)‐OH of fullerene amino acid derived peptides have been prepared by solid phase peptide synthesis, in which the terminal amino acid, Phe(4‐aza‐C₆₀)‐OH, is derived from the dipolar addition to C₆₀ of the Fmoc‐Nα‐protected azido amino acids derived from phenylalanine: Fmoc‐Phe(4‐aza‐C₆₀)‐Lys₃‐OH ( 1 ), Fmoc‐Phe(4‐aza‐C₆₀)‐Pro‐Hyp‐Lys‐OH ( 2 ), and Fmoc‐Phe(4‐aza‐C₆₀)‐Hyp‐Hyp‐Lys‐OH ( 3 ). The inhibition constant of our fullerene aspartic protease PRIs utilized FRET‐based assay to evaluate the enzyme kinetics of HIV‐1 PR at various concentrations of inhibitors. Simulation of the docking of the peptide Fmoc‐Phe‐Pro‐Hyp‐Lys‐OH overestimated the inhibition, while the amino acid PRIs were well estimated. The experimental results show that C₆₀‐based amino acids are a good base structure in the design of protease inhibitors and that their inhibition can be improved upon by the addition of designer peptide sequences. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

Self‐assembling study of sarcolipin and its mutants in multiple molecular dynamic simulations

The Sarcolipin (SLN) is a single trans‐membrane protein that can self‐assembly to dimer and oligomer for playing importantphysiological function. In this work, we addressed the dimerization of wild type SLN (wSLN) and its mutants (mSLNs) – I17A and I20A, using both coarse‐grained (CG) and atomistic (AT) molecular dynamics (MD) simulations. Our results demonstrated that wSLN homodimer assembled as a left‐handed helical complex, while mSLNs heterodimers assembled as right‐handed complexes. Analysis of residue‐residue contacts map indicated that isoleucine (Ile)‐leucione (Leu) zipper domain played an important role in dimerization. The potential of mean force (PMF) demonstrated that wSLN homodimer was more stable than mSLNs heterodimers. Meanwhile, the mSLNs heterodimers preferred right‐handed rather than left‐handed helix. AT‐MD simulations for wSLN and mSLNs were also in line with CG‐MD simulations. These results provided the insights for understanding the mechanisms of SLNs self‐assembling. Proteins 2017; 85:1065–1077. © 2017 Wiley Periodicals, Inc.