Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

Molecular docking and molecular dynamics simulation approaches for identifying new lead compounds as potential AChE inhibitors

To provide hints for the design of novel acetylcholinesterase (AChE) inhibitors with higher potency and specificity, the binding modes of the (RS, S)-17b and (RS, R)-17b enantiomers on AChE were chosen to investigate by molecular docking and molecular dynamics simulation. The results show that the binding modes of (RS, S)-17b and (RS, R)-17b are clearly different from each other. In particular, the (RS, S)-17b and (RS, R)-17b enantiomers tend to be planar and bend conformations to interact with AChE, respectively. Furthermore, based on the binding mode on AChE and structure modification of (RS, S)-17b, two novel inhibitors (1 and 2) with higher inhibitory activity were designed. Our design strategy suggests that the number of N and O atoms should be increased, the 5, 6-dimethoxy should be transformed into ring and the indanone moiety should be ring-opening, which would result in generating potent and selective AChE inhibitors.     

Novel urushiol derivatives as HDAC8 inhibitors: rational design,virtual screening,molecular docking and molecular dynamics studies

Three series of novel urushiol derivatives were designed by introducing a hydroxamic acid moiety into the tail of an alkyl side chain and substituents with differing electronic properties or steric bulk onto the benzene ring and alkyl side chain. The compounds’ binding affinity toward HDAC8 was screened by Glide docking. The highest-scoring compounds were processed further with molecular docking, MD simulations, and binding free energy studies to analyze the binding modes and mechanisms. Ten compounds had Glide scores of ?8.2 to ?10.2, which revealed that introducing hydroxy, carbonyl, amino, or methyl ether groups into the alkyl side chain or addition of –F, –Cl, sulfonamide, benzamido, amino, or hydroxy substituents on the benzene ring could significantly increase binding affinity. Molecular docking studies revealed that zinc ion coordination, hydrogen bonding, and hydrophobic interactions contributed to the high calculated binding affinities of these compounds toward HDAC8. MD simulations and binding free energy studies showed that all complexes possessed good stability, as characterized by low RMSDs, low RMSFs of residues, moderate hydrogen bonding and zinc ion coordination and low values of binding free energies. Hie147, Tyr121, Phe175, Hip110, Phe119, Tyr273, Lys21, Gly118, Gln230, Leu122, Gly269, and Gly107 contributed favorably to the binding; and Van der Waals and electrostatic interactions provided major contributions to the stability of these complexes. These results show the potential of urushiol derivatives as HDAC8 binding lead compounds, which have great therapeutic potential in the treatment of various malignancies, neurological disorders, and human parasitic diseases.     

Febrifugine analogues as Leishmania donovani trypanothione reductase inhibitors: binding energy analysis assisted by molecular docking,ADMET and molecular dynamics simulation

Visceral leishmaniasis affects people from 70 countries worldwide, mostly from Indian, African and south American continent. The increasing resistance to antimonial, miltefosine and frequent toxicity of amphotericin B drives an urgent need to develop an antileishmanial drug with excellent efficacy and safety profile. In this study we have docked series of febrifugine analogues (n = 8813) against trypanothione reductase in three sequential docking modes. Extra precision docking resulted into 108 ligands showing better docking score as compared to two reference ligand. Furthermore, 108 febrifugine analogues and reference inhibitor clomipramine were subjected to ADMET, QikProp and molecular mechanics, the generalized born model and solvent accessibility study to ensure the toxicity caused by compounds and binding-free energy, respectively. Two best ligands (FFG7 and FFG2) qualifying above screening parameters were further subjected to molecular dynamics simulation. Conducting these studies, here we confirmed that 6-chloro-3-[3-(3-hydroxy-2-piperidyl)-2-oxo-propyl]-7-(4-pyridyl) quinazolin-4-one can be potential drug candidate to fight against Leishmania donovani parasites.     

Atomic insight into designed carbamate-based derivatives as acetylcholine esterase (AChE) inhibitors: a computational study by multiple molecular docking and molecular dynamics simulation

Over 100 variants have been designed and studied, using multiple docking methods such as Autodock Vina, ArgusLab, Molegro Virtual Docker, and Hex-Cuda, to study the effect of alteration in the structure of carbamate-based acetylcholyne esterase (AChE) inhibitors. Sixteen selected systems were then subjected to 14 ns molecular dynamics (MD) simulations. Results from all the docking methods are in agreement. Variants that involved biphenyl substituents possess the most negative binding energies in the ?37.64 to ?39.31 kJ mol^?1 range due to their π–π interactions with AChE aromatic residues. The root mean square deviation values showed that all of these components achieved equilibration after 6 ns. Gyration radius (R_g) and solvent accessibility surface area were calculated to further investigate the AChE conformational changes in the presence of these components. MD simulation results suggested that these components might interact with AChE, possibly with no major changes in AChE secondary and tertiary structures.     

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.     

New potential inhibitors of mTOR: a computational investigation integrating molecular docking,virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.     

Analysis of Coxiela burnetti dihydrofolate reductase via in silico docking with inhibitors and molecular dynamics simulation

Coxiella burnetii is a gram-negative bacterium able to infect several eukaryotic cells, mainly monocytes and macrophages. It is found widely in nature with ticks, birds, and mammals as major hosts. C. burnetii is also the biological warfare agent that causes Q fever, a disease that has no vaccine or proven chemotherapy available. Considering the current geopolitical context, this fact reinforces the need for discovering new treatments and molecular targets for drug design against C. burnetii. Among the main molecular targets against bacterial diseases reported, the enzyme dihydrofolate reductase (DHFR) has been investigated for several infectious diseases. In the present work, we applied molecular modeling techniques to evaluate the interactions of known DHFR inhibitors in the active sites of human and C. burnetii DHFR (HssDHFR and CbDHFR) in order to investigate their potential as selective inhibitors of CbDHFR. Results showed that most of the ligands studied compete for the binding site of the substrate more effectively than the reference drug trimethoprim. Also the most promising compounds were proposed as leads for the drug design of potential CbDHFR inhibitors.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Pharmacophore modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

Synthesis,molecular docking,binding free energy calculation and molecular dynamics simulation studies of benzothiazol-2-ylcarbamodithioates as Staphylococcus aureus MurD inhibitors

Abstract

A new series of benzothiazol-2-ylcarbamodithioate functional compounds 5a-f has been designed, synthesized and characterized by spectral data. These compounds were screened for their in vitro antibacterial activity against strains of Staphylococcus aureus (NCIM 5021, NCIM 5022 and methicillin-resistant isolate 43300), Bacillus subtilis (NCIM 2545), Escherichia coli (NCIM 2567), Klebsiella pneumoniae (NCIM 2706) and Psudomonas aeruginosa (NCIM 2036). Compounds 5a and 5d exhibited significant activity against all the tested bacterial strains. Specifically, compounds 5a and 5d showed potent activity against K. pneumoniae (NCIM 2706), while compound 5a also displayed potent activity against S. aureus (NCIM 5021). Compound 5d showed minimum IC₅₀ value of 13.37?μM against S. aureus MurD enzyme. Further, the binding interactions of compounds 5a-f in the catalytic pocket have been investigated using the extra-precision molecular docking and binding free energy calculation by MM-GBSA approach. A 30?ns molecular dynamics simulation of 5d/modeled S. aureus MurD enzyme was performed to determine the stability of the predicted binding conformation.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.     

Design,docking studies and molecular dynamics of new potential selective inhibitors of Plasmodium falciparum serine hydroxymethyltransferase

In this paper, former studies on the interactions of the natural substrate and potential inhibitors of Plasmodium falciparum serine hydroxymethyltransferase (PfSHMT) were used to design five new potential selective inhibitors to this enzyme. Results of the docking energies calculations of these structures inside the active sites of PfSHMT and human SHMT were used to select a more suitable structure as a potential selective inhibitor to PfSHMT. Further molecular dynamics studies of this molecule and 5-formyl-6-hydrofolic acid (natural substrate) docked inside these enzymes' active sites revealed important features for additional refinements of this structure and also additional residues in the PfSHMT active site to be considered further for designing selective inhibitors.     

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

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.