Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Insights into the inhibition of xanthine oxidase by curcumin

As a natural pigment, curcumin exhibits multiple biological activities. Previous studies have investigated the inhibition of xanthine oxidase (XO) by curcumin. In the present work, based on the molecular docking simulations, it is interesting to find that parent curcumin binds weakly to XO, while its degradation products, for example, trans-6-(4′-hydroxy-3′-methoxyphenyl)-2,4-dioxo-5-hexenal, exhibit effective inhibitory activities against XO. The findings shed new light on the underlying mechanisms of curcumin in inhibiting XO and also have potential implication that both parent curcumin and its degradation products should be taken into account when exploring the mechanisms of curcumin’s biological activities.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Pharmacophore feature prediction and molecular docking approach to identify novel anti‐HCV protease inhibitors

Discovering a potential drug for HCV treatment is a challenging task in the field of drug research. This study initiates with computational screening and modeling of promising ligand molecules. The foremost modeling method involves the identification of novel compound and its molecular interaction based on pharmacophore features. A total of 197 HCV compounds for NS3/4A protein target were screened for our study. The pharmacophore models were generated using PHASE module implemented in Schrodinger suite. The pharmacophore features include one hydrogen bond acceptor, one hydrogen bond donor, and three hydrophobic sites. As a result, based on mentioned hypothesis the model ADHHH.159 corresponds to the CID 59533233. Furthermore, docking was performed using maestro for all the 197 compounds. Among these, the CID 59533313 and 59533233 possess the best binding energy of ?11.75 and ?10.40 kcal/mol, respectively. The interactions studies indicated that the CID complexed with the NS3/4A protein possess better binding affinity with the other compounds. Further the compounds were subjected to calculate the ADME properties. Therefore, it can be concluded that these two compounds could be a potential alternative drug for the development of HCV.     

Homology modeling, molecular dynamics, e-pharmacophore mapping and docking study of Chikungunya virus nsP2 protease

To date, no suitable vaccine or specific antiviral drug is available to treat Chikungunya viral (CHIKV) fever. Hence, it is essential to identify drug candidates that could potentially impede CHIKV infection. Here, we present the development of a homology model of nsP2 protein based on the crystal structure of the nsP2 protein of Venezuelan equine encephalitis virus (VEEV). The protein modeled was optimized using molecular dynamics simulation; the junction peptides of a nonstructural protein complex were then docked in order to investigate the possible protein–protein interactions between nsP2 and the proteins cleaved by nsP2. The modeling studies conducted shed light on the binding modes, and the critical interactions with the peptides provide insight into the chemical features needed to inhibit the CHIK virus infection. Energy-optimized pharmacophore mapping was performed using the junction peptides. Based on the results, we propose the pharmacophore features that must be present in an inhibitor of nsP2 protease. The resulting pharmacophore model contained an aromatic ring, a hydrophobic and three hydrogen-bond donor sites. Using these pharmacophore features, we screened a large public library of compounds (Asinex, Maybridge, TOSLab, Binding Database) to find a potential ligand that could inhibit the nsP2 protein. The compounds that yielded a fitness score of more than 1.0 were further subjected to Glide HTVS and Glide XP. Here, we report the best four compounds based on their docking scores; these compounds have IDs of 27943, 21362, ASN 01107557 and ASN 01541696. We propose that these compounds could bind to the active site of nsP2 protease and inhibit this enzyme. Furthermore, the backbone structural scaffolds of these four lead compounds could serve as building blocks when designing drug-like molecules for the treatment of Chikungunya viral fever.     

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     

Pharmacophore and structure-based drug design,molecular dynamics and admet/tox studies to design novel potential pad4 inhibitors

We have used docking (GLIDE), pharmacophore modeling (Discovery Studio), long trajectory molecular dynamics (Discovery Studio) and ADMET/Tox (QikProp and DEREK) to investigate PAD4 in order to determine potential novel inhibitors and hits. We have carried out virtual screening in the ZINC natural compounds database. Pharmacokinetics and Toxicity of the best hits were assessed using databases implemented in softwares that create models based on chemical structures taking into account consideration about the toxicophoric groups. A wide variety of pharmaceutical relevant properties are determined in order to make decisions about molecular suitability. After screening and analysis, the 6 most promising PAD4 inhibitors are suggested, with strong interactions (pi-stacking, hydrogen bonds, hydrophobic contacts) and suitable pharmacotherapeutic profile as well.     

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.     

Hierarchical virtual screening of the dual MMP-2/HDAC-6 inhibitors from natural products based on pharmacophore models and molecular docking

The dual-target inhibitors tend to improve the response rate in treating tumors, comparing with the single-target inhibitors. Matrix metalloproteinase-2 (MMP-2) and histone deacetylase-6 (HDAC-6) are attractive targets for cancer therapy. In this study, the hierarchical virtual screening of dual MMP-2/HDAC-6 inhibitors from natural products is investigated. The pharmacophore model of MMP-2 inhibitors is built based on ligands, but the pharmacophore model of HDAC-6 inhibitors is built based on the experimental crystal structures of multiple receptor–ligand complexes. The reliability of these two pharmacophore models is validated subsequently. The hierarchical virtual screening, combining these two different pharmacophore models of MMP-2 and HDAC-6 inhibitors with molecular docking, is carried out to identify the dual MMP-2/HDAC-6 inhibitors from a database of natural products. The four potential dual MMP-2/HDAC-6 inhibitors of natural products, STOCK1 N-46177, STOCK1 N-52245, STOCK1 N-55477, and STOCK1 N-69706, are found. The studies of binding modes show that the screened four natural products can simultaneously well bind with the MMP-2 and HDAC-6 active sites by different kinds of interactions, to inhibit the MMP-2 and HDAC-6 activities. In addition, the ADMET properties of screened four natural products are assessed. These found dual MMP-2/HDAC-6 inhibitors of natural products could serve as the lead compounds for designing the new dual MMP-2/HDAC-6 inhibitors having higher biological activities by carrying out structural modifications and optimizations in the future studies.     

Virtual screening to identify <Emphasis Type="Italic">Leishmania braziliensis N</Emphasis>-myristoyltransferase inhibitors: pharmacophore models,docking, and molecular dynamics

Leishmaniasis is caused by several protozoa species belonging to genus Leishmania that are hosted by humans and other mammals. Millions of new cases are recorded every year and the drugs available on the market do not show satisfactory efficacy and safety. A hierarchical virtual screening approach based on the pharmacophore model, molecular docking, and molecular dynamics was conducted to identify possible Leishmania braziliensis N-misristoyltransferase (LbNMT) inhibitors. The adopted pharmacophore model had three main features: four hydrophobic centers, four hydrogen-bond acceptor atoms, and one positive nitrogen center. The molecules (n=15,000) were submitted to alignment with the pharmacophore model and only 27 molecules aligned to model. Six molecules were submitted to molecular docking, using receptor PDB ID 5A27. After docking, the ZINC35426134 was a top-ranked molecule (? 64.61 kcal/mol). The molecule ZINC35426134 shows hydrophobic interactions with Phe82, Tyr209, Val370, and Leu391 and hydrogen bonds with Asn159, Tyr318, and Val370. Molecular dynamics simulations were performed with the protein in its APO and HOLO forms for 37 ns in order to assess the stability of the protein–ligand complex. Results showed that the HOLO form was more stable than the APO one, and it suggests that the ZINC35426134 binding stabilizes the enzyme. Therefore, the selected molecule has the potential to meet the herein proposed target.     

Combined pharmacophore,virtual screening and molecular dynamics studies to identify Bruton’s tyrosine kinase inhibitors

Bruton’s tyrosine Kinase (BTK) is a cytoplasmic, non-receptor tyrosine kinase expressed in hematopoietic cells. BTK plays a critical role in many cellular signalling pathways making it a potential target to treat autoimmune diseases and cancer. BTK signalling is important for the production of arthritis-associated antibodies, and inhibiting BTK will help the system to block the production of disease-associated antibodies. In this study, we have implemented ligand-based pharmacophore modelling and virtual screening against natural compounds followed by molecular docking, density functional theory and molecular dynamics studies for 50 ns. Four compounds with high affinity towards BTK were identified, and it could be used as a potent lead molecule for designing BTK inhibitor.     

Pharmacoinformatics analysis to identify inhibitors of Mtb-ASADH

Aspartate-semialdehyde dehydrogenase (ASADH; EC 1.2.1.11) is a key enzyme in the biosynthesis of essential amino acids in prokaryotes and fungi, inhibition of ASADH leads to the development of novel antitubercular agents. In the present work, a combined structure and ligand-based pharmacophore modeling, molecular docking, and molecular dynamics (MD) approaches were employed to identify potent inhibitors of mycobacterium tuberculosis (Mtb)-ASADH. The structure-based pharmacophore hypothesis consists of three hydrogen bond acceptor (HBA), two negatively ionizable, and one positively ionizable center, while ligand-based pharmacophore consists of additional one HBA and one hydrogen bond donor features. The validated pharmacophore models were used to screen the chemical databases (ZINC and NCI). The screened hits were subjected to ADME and toxicity filters, and subsequently to the molecular docking analysis. Best-docked 25 compounds carry the characteristics of highly electronegative functional groups (–COOH and –NO₂) on both sides and exhibited the H-bonding interactions with highly conserved residues Arg99, Arg249, and His256. For further validation of docking results, MD simulation studies were carried out on two representative compounds NSC51108 and ZINC04203124. Both the compounds remain bound to the key active residues of Mtb-ASADH during the MD simulations. These identified hits can be further used for lead optimization and in the design more potent inhibitors against Mtb-ASADH.     

Exploration and validation of diphosphate‐based Plasmodium LytB inhibitors using computational approaches

LytB or IspH is an indispensable enzyme and a suitable drug target of Plasmodium falciparum that participate in isoprenoid biosynthesis of nonmevalonate pathway (MEP). Recently, we have investigated the structural dynamics of Plasmodium LytB and proposed some novel diphosphate‐based inhibitors using molecular modeling and docking studies. Here, we have tried to characterize those previously screened molecules by quantitative structure activity relationships and pharmacophore‐based analyses, as well as validated the dynamics of their interactions with LytB protein. Five total compounds having PubChem CID 516 , 125696 , 165275 , 448012 , and 9921431 were predicted with significant inhibitory activity by quantitative structure activity relationships and pharmacophore models. Again, the molecular dynamics simulation results showed that these five compounds are able to form stable complexes with the receptor through many direct and water mediated interactions. The binding free energies calculated by Poisson‐Boltzmann surface area method resulted within the range between ?99.77 and ? 43.74 kcal/mol, which favoured their profound inhibitory affinity. Residues of LytB like His41, His74, Ser222, Ser223, and Asn224 in LytB were the main protagonists in contributing the majority of interaction energies to the ligands. Finally, the ADMET, toxicity, and drug‐likeness scores also affirm these compounds to be considered for further development of new antimalarial inhibitor in the future.     

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Evaluation of the combination mode and features of p38 MAPK inhibitors: construction of different pharmacophore models and molecular docking

Abnormal expression of tumour necrosis factor-α (TNF-α) can lead to various pathological reactions, such as arthritis, psoriasis, krone disease, etc. p38 mitogen-activated protein kinase (p38 MAPK) is an important signal transduction enzyme that plays important roles in influencing the release of intracellular TNF-α factor. It is very meaningful to study the targeting kinase with specific inhibitors in the treatment of related diseases. In order to achieve a deeper insight, it is necessary to analyse the structural characteristics and the action mode of the p38 MAPK inhibitors in the active site. In the study, a ligand-based common feature pharmacophore model and the receptor structure-based pharmacophore model were constructed, respectively. Their common chemical features consisted of the hydrophobic groups (H) and the hydrogen bond acceptors (A), and kept the consistency of spatial structure distribution. Then, the molecular docking and molecular dynamics simulation were performed with the eight training set compounds. The binding characteristics of molecules binding were described in the topological region of the active site. Finally, the structure–activity relationship (SAR) was obtained by analysing docking results with the different pharmacophore models. This research leads to the proposal of an interaction model in the p38 MAPK active site and provides guidance for the screening and design of more potent and selective p38 MAPK inhibitors.     

Inhibitory effects of flavonoids on aldehyde oxidase activity

Flavonoids are an important group of natural compounds that can interfere with the activity of some enzymes. In this study, effects of various flavonoids on aldehyde oxidase (AO) activity were evaluated in vitro. AO was partially purified from guinea pig liver. The effects of 12 flavonoids from three subclasses of flavon-3-ol, flavan-3-ol and flavanone on the oxidation of vanillin and phenanthridine as substrates of AO and xanthine as a substrate of xanthine oxidase (XO) were investigated spectrophotometrically. Among the 12 flavonoids, myricetin and quercetin were the most potent inhibitors of both AO and XO. In general, the oxidation of vanillin was more inhibited by flavonoids than that of phenanthridine. Almost all of the flavonoids inhibited AO activity more potently than XO, which was more evident with non-planner flavanols. A planner structure seems to be essential for a potent inhibitory effect and any substitution by sugar moieties reduces the inhibitory effects. This study could provide a new insight into AO natural inhibitors with potential to lead to some food-drug interactions.     

Inhibitory effects of flavonoids on aldehyde oxidase activity

Flavonoids are an important group of natural compounds that can interfere with the activity of some enzymes. In this study, effects of various flavonoids on aldehyde oxidase (AO) activity were evaluated in vitro. AO was partially purified from guinea pig liver. The effects of 12 flavonoids from three subclasses of flavon-3-ol, flavan-3-ol and flavanone on the oxidation of vanillin and phenanthridine as substrates of AO and xanthine as a substrate of xanthine oxidase (XO) were investigated spectrophotometrically. Among the 12 flavonoids, myricetin and quercetin were the most potent inhibitors of both AO and XO. In general, the oxidation of vanillin was more inhibited by flavonoids than that of phenanthridine. Almost all of the flavonoids inhibited AO activity more potently than XO, which was more evident with non-planner flavanols. A planner structure seems to be essential for a potent inhibitory effect and any substitution by sugar moieties reduces the inhibitory effects. This study could provide a new insight into AO natural inhibitors with potential to lead to some food-drug interactions.