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 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.     

Preparation of l-proline based aeruginosin 298-A analogs: Optimization of the P1-moiety

Aeruginosins are a family of naturally occurring oligopeptides that share a common bicyclic amino acid core structure. Many compounds in the family are inhibitors of serine proteases, such as thrombin and trypsin. Thrombin is an important enzyme in the blood coagulation cascade, and is a promising target for anticoagulant drug development. In order to understand the structure–activity relationship (SAR) and to find selective thrombin inhibitors, we synthesized a series of aeruginosin 298-A analogs, in which the P₂ bicyclic amino acid was replaced by a l-proline residue. The structure optimization was focused on modification of the P₁ position. In choosing the P₁ group, an effort was made to avoid using the highly basic guanidine groups present in nearly all naturally occurring aeruginosins. The synthesis and enzyme assays of these aeruginosin analogs against thrombin and trypsin are reported. We found that several compounds with neutral P₁ groups exhibit excellent selectivity over trypsin and good potency against thrombin. The SAR data of the P₁ groups obtained here can be used in preparing other thrombin inhibitors with better selectivity against trypsin.     

Discovery of novel purine nucleoside derivatives as phosphodiesterase 2 (PDE2) inhibitors: Structure-based virtual screening,optimization and biological evaluation

Phosphodiesterase 2 (PDE2) has received much attention for the potential treatment of the central nervous system (CNS) disorders and pulmonary hypertension. Herein, we identified that clofarabine (4), an FDA-approved drug, displayed potential PDE2 inhibitory activity (IC₅₀?=?3.12?±?0.67?μM) by structure-based virtual screening and bioassay. Considering the potential therapeutic benefit of PDE2, a series of purine nucleoside derivatives based on the structure and binding mode of 4 were designed, synthesized and evaluated, which led to the discovery of the best compound 14e with a significant improvement of inhibitory potency (IC₅₀?=?0.32?±?0.04?μM). Further molecular docking and molecular dynamic (MD) simulations studies revealed that 5′-benzyl group of 14e could interact with the unique hydrophobic pocket of PDE2 by forming extra van der Waals interactions with hydrophobic residues such as Leu770, Thr768, Thr805 and Leu809, which might contribute to its enhancement of PDE2 inhibition. These potential compounds reported in this article and the valuable structure-activity relationships (SARs) might bring significant instruction for further development of potent PDE2 inhibitors.     

Molecular modeling studies,synthesis and biological evaluation of dabigatran analogues as thrombin inhibitors

In this work, 48 thrombin inhibitors based on the structural scaffold of dabigatran were analyzed using a combination of molecular modeling techniques. We generated three-dimensional quantitative structure–activity relationship (3D-QSAR) models based on three alignments for both comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) to highlight the structural requirements for thrombin protein inhibition. In addition to the 3D-QSAR study, Topomer CoMFA model also was established with a higher leave-one-out cross-validation q² and a non-cross-validation r², which suggest that the three models have good predictive ability. The results indicated that the steric, hydrophobic and electrostatic fields play key roles in QSAR model. Furthermore, we employed molecular docking and re-docking simulation explored the binding relationship of the ligand and the receptor protein in detail. Molecular docking simulations identified several key interactions that were also indicated through 3D-QSAR analysis. On the basis of the obtained results, two compounds were designed and predicted by three models, the biological evaluation in vitro (IC₅₀) demonstrated that these molecular models were effective for the development of novel potent thrombin 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.     

Stereoselectivity of phosphotriesterase with paraoxon derivatives: a computational study

The bacterial enzyme phosphotriesterase (PTE) exhibits stereoselectivity toward hydrolysis of chiral substrates with a preference for the S_p enantiomer. In this work, docking analysis and two explicit-solvent molecular dynamics (MD) simulations were performed to characterize and differentiate the structural dynamics of PTE bound to the S_p and R_p paraoxon derivative enantiomers (R_p-1 and S_p-1) hydrolyzed with distinct catalytic efficiencies. Comparative analysis of the molecular trajectories for PTE bound to R_p-1 and S_p-1 suggested that substrate binding induced conformational changes in the loops near the active site. After 100 ns of MD simulation, the Zn β²⁺ metal ion formed hexacoordinated- and tetracoordinated geometries in the S_p-1-PTE and R_p-1-PTE ensembles, respectively. Simulation results further showed that the hydrogen bond between Asp301 and His254 occurred with a higher probability after S_p-1 binding to PTE (47.5%) than that after R_p-1 binding (22.2%). These results provide a qualitative and molecular-level explanation for the 10 orders of magnitude increase in the catalytic efficiency of PTE toward the S_p enantiomer of paraoxon.     

Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

New substituted aminopyrimidine derivatives as BACE1 inhibitors: in silico design,synthesis and biological assays

We report in this work new substituted aminopyrimidine derivatives acting as inhibitors of the catalytic site of BACE1. These compounds were obtained from a molecular modeling study. The theoretical and experimental study reported here was carried out in several steps: docking analysis, Molecular Dynamics (MD) simulations, Quantum Theory Atom in Molecules (QTAIM) calculations, synthesis and bioassays and has allowed us to propose some compounds of this series as new inhibitors of the catalytic site of BACE1. The QTAIM study has allowed us to obtain an excellent correlation between the electronic densities and the experimental data of IC₅₀. Also, using combined techniques (MD simulations and QTAIM calculations) enabled us to describe in detail the molecular interactions that stabilize the different L-R complexes. In addition, our results allowed us to determine what portion of these compounds should be changed in order to increase their affinity with the BACE1. Another interesting result is that a sort of synergism was observed when the effects of these new catalytic site inhibitors were combined with Ac-Tyr5-Pro6-Tyr7-Asp8-Ile9-Pro10-Leu11-NH₂, which we have recently reported as a modulator of BACE1 acting on its exosite.     

Theoretical studies on the structures,intra- and inter-molecular hydrogen bonding interactions in HNF and HNF–H2O clusters in the gaseous,aqueous and solid phases

Hydrazimium nitroformate ([N₂H₅]⁺[C(NO₂)₃]^? , HNF) is an ionic oxidiser used in solid propellants. Its properties are easily affected by H₂O because of its hygroscopicity. In this article, density functional theory (DFT) and molecular dynamics (MD) were employed to study the isolated HNF molecule and the HNF–H₂O cluster in gas phase and in the aqueous solution. Three stable conformations were obtained for HNF in the gas phase and in the aqueous solution, respectively, and each conformation can form several different HNF–H₂O clusters. Irrespective of whether it is in gas phase or in solution, intramolecular hydrogen bond interactions and other interactions (e.g. the binding energy, the dispersion energy, the second-order perturbation energy and the energy gap between frontier orbitals) of HNF are weaker in the clusters than in the isolated state. The initial decomposition energy of the cluster is lower than that of the isolated HNF molecule in both gaseous and aqueous phases, while the dissociation processes are the same. Molecular dynamic simulations showed that the clustered H₂O elongates and weakens the C–NO₂ bond in the solid HNF–H₂O cluster compared with that in the solid HNF. H₂O reduces and weakens intramolecular N–HΛO bonds too, and O–HΛN is the dominant intermolecular hydrogen bond between HNF and H₂O.     

Inhibitory effects of epigallocatechin gallate and its glucoside on the human intestinal maltase inhibition

Human intestinal maltase (HMA) is an ??-glucosidase responsible for the hydrolysis of ??-1,4-linkages from the non-reducing end of malto-oligosaccharides. HMA has become an important target in the treatment of type-2 diabetes. In this study, epigallocatechin gallate (EGCG) and EGCG glucoside (EGCG-G1) were identified as inhibitors of HMA by an in vitro assay with IC₅₀ of 20 ± 1.0 and 31.5 ± 1.0 ??M, respectively. A Lineweaver-Burk plot confirmed that EGCG and EGCG-G1 were competitive inhibitors of maltose substrate against HMA and inhibition kinetic constants (K _i ) calculated from a Dixon plot were 5.93 ± 0.26 and 7.88 ± 0.57 ??M, respectively. Both EGCG and EGCG-G1 bound to the active site of HMA with numerous hydrophobic and hydrogen bond interactions.     

Long time molecular dynamic simulation on the agonist binding and activation of the β2-adrenergic receptor

β₂-Adrenergic receptor (β₂AR) plays a critical role in mediating the effects of catecholamine hormones. Due to the flexibility of the structure of its active state, study of agonist–β2AR is usually performed by molecular dynamic (MD) simulation. In this study we show the representative characteristics of agonist binding and activation on β₂AR by MD simulation. The binding site and the conformational changes in the specific regions of β2AR are reasonable which confirmed the conclusion that agonist–β2AR reached its active-like state. We have studied the disruption of non-covalent intramolecular interactions, including the conserved DRY motif, the rotamer toggle switch and the ionic lock, the cytoplasmic ends of transmembranes 5 and 6, and some water-mediated hydrogen bond network regions. We conclude that agonist induced β₂AR to its active conformation, or at least the active-like conformation.     

Structural insight into the inhibition of acetylcholinesterase by 2,3,4, 5-tetrahydro-1, 5-benzothiazepines

Benzothiazepines 1–3 inhibited acetylcholinesterase (AChE; EC 3.1.1.7) enzyme in a concentration-dependent fashion with IC₅₀ values of 1.0 ± 0.002, 1.2 ± 0.005 and 1.3 ± 0.001 μM, respectively. By using linear-regression equations, Lineweaver-Burk, Dixon plots and their secondary replots were constructed which indicated that compounds 1–3 are non-competitive inhibitors of AChE with K_i values of 0.8 ± 0.04, 1.1 ± 0.002, and 1.5 ± 0.001 μM, respectively. Molecular docking studies revealed that all the compounds are completely buried inside the aromatic gorge of AChE, extending deep into the gorge of AChE. A comparison of the docking results of compounds 1–3 displayed that these compounds generally adopt the same binding mode in the active site of AChE. The superposition of the docked structures demonstrated that the non-flexible benzothiazepine always penetrate into the aromatic gorge through the six-membered ring A, which allowed the ligands to interact simultaneously with more than one subsites of the active center of AChE. The higher AChE inhibitory potential of compounds 1–3 was found to be the cumulative effect of hydrophobic contacts and π-π interactions between the ligands and AChE. The relatively high affinity of benzothiazepine 1 with AChE was found to be due to additional hydrogen bond in benzothiazepine 1-AChE complex. The results indicated that substitution of halogen and methyl groups by hydrogen at aromatic ring of the benzothiazepine decreased the affinity of these molecules towards enzyme that may be due to the polar non-polar repulsions of these moieties with the amino acid residues in the active site of AChE. The observed binding modes of benzothiazepines 1–3 in the active site of AChE explain the affinities of benzothiazepines and provide a rational basis for the structure-based drug design of benzothiazepines with improved pharmacological properties.     

Biflavones from Ginkgo biloba as inhibitors of human thrombin

Ginkgo Biloba leaf extract has been widely used for the prevention and treatment of thrombosis and cardiovascular disease in both eastern and western countries, but the bioactive constituents and the underlying mechanism of anti-thrombosis have not been fully characterized. The purpose of this study was to investigate the inhibitory effects of major constituents in Ginkgo biloba on human thrombin, a key serine protease regulating the blood coagulation cascade and the processes of thrombosis. To this end, a fluorescence-based biochemical assay was used to assay the inhibitory effects of sixteen major constituents from Ginkgo biloba on human thrombin. Among all tested natural compounds, four biflavones (ginkgetin, isoginkgetin, bilobetin and amentoflavone), and five flavonoids (luteolin, apigenin, quercetin, kaempferol and isorhamnetin) were found with thrombin inhibition activity, with the IC₅₀ values ranging from 8.05 μM to 82.08 μM. Inhibition kinetic analyses demonstrated that four biflavones were mixed inhibitors against thrombin-mediated Z-GGRAMC acetate hydrolysis, with the K_i values ranging from 4.12 μM to 11.01 μM. Molecular docking method showed that the four biflavones could occupy the active cavity with strong interactions of salt bridges and hydrogen bonds. In addition, mass spectrometry-based lysine labeling reactivity assay suggested that the biflavones could bind on human thrombin at exosite I rather than exosite II. All these findings suggested that the biflavones in Ginkgo biloba were naturally occurring inhibitors of human thrombin, and these compounds could be used as lead compounds for the development of novel thrombin inhibitors with improved efficacy and high safety profiles.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

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.     

Synthesis of CC,CN coupled novel substituted dibutyl benzothiazepinone derivatives and evaluation of their thrombin inhibitory activity

The formation of a thrombus is a key event in thromboembolic disorders, that contribute to high mortality and morbidity in affected patients. In the present study, we synthesized a library of novel substituted 3,3-dibutyl-8-methoxy-2,3-dihydrobenzo [b] [1,4] thiazepin-4(5H)-one derivatives which were tested for their platelet aggregation and thrombin inhibitory activity. Among the tested compounds, 3,3-dibutyl-7-(2-chlorophenyl)-8-methoxy-2,3-dihydrobenzo[b] [1,4]thiazepin-4(5H)-one (DCT) displayed the maximum thrombin inhibition with an IC₅₀ value of 3.85 μM and thus DCT was chosen for further studies. Next, the effect of DCT on primary hemostasis was evaluated using agonist-induced platelet aggregation model. The lead compound inhibited the collagen- or ADP- or thrombin-induced platelet aggregation in a dose-dependent manner. Furthermore, DCT prolonged the process of clot formation when evaluating plasma re-calcification time (320 ± 11 sec at 5 µg DCT), activated partial thromboplastin time (58.0 ± 0.01 sec at 2 µg), and prothrombin time (14.7 ± 0.01 sec at 5 µg). Molecular docking studies suggested that the benzothiazepinones evaluated here consistently display hydrogen bonding with Ser214 of thrombin, which is similar to that of the co-crystallized ligand (1-(2R)-2-amino-3-phenyl-propanoyl-N-(2,5dichlorophenyl)methylpyrrolidine-2-carboxamide). DCT displayed additional hydrogen bonding to Ser195 and π-π interactions between its methoxyphenyl groups and Trp60, thereby providing a structural rationale for the observed biological effect.     

Computational approach for the assessment of inhibitory potency against beta-amyloid aggregation

Beta-amyloid (Aβ) plaques are one of the hallmarks of Alzheimer’s disease. Their presence in the brain leads to neurodegeneration and memory decline. Therefore, search for new drugs able to decrease formation of such deposits is of great interest. Our previously developed multifunctional compounds inhibited transformation of monomers into fibrils. Herein, we describe the computational approach for the assessment of inhibitory activity against Aβ aggregation. The influence of novel inhibitors on amyloid Aβ_17-42 was studied by employing of molecular docking and all-atom molecular dynamics simulations. We found that the number of intermolecular backbone hydrogen bonds at the end of 100 ns MD simulation was correlated with the level of anti-aggregation potency of studied compounds. Such data may be successfully applied to in silico design of novel inhibitors of Aβ aggregation.     

Crystal packing driven by metal-ligand interactions, hydrogen bonds, π-π stacking and anion-π stacking

Two polymeric copper(II) compounds have been synthesized with the ligand bis(pyrimidin-2yl)amine (dipm), by means of coordination bonds, Watson-Crick type hydrogen bond interactions and π-π stacking. Both supramolecules hold the same cationic building block, namely [Cu(dipm)₂(H₂O)₂]²⁺. However, their crystal structure significantly differs and this variation apparently arises from the nature of the anion which induces dissimilar crystal packings. The crystal growth is driven by several synergistic intermolecular interactions, i.e., coordination and hydrogen bonds, π-π and anion-π stacks.