Docking and quantitative structure-activity relationship of oxadiazole derivates as inhibitors of GSK3β

Computational Modeling to Explain Why 5,5-Diarylpentadienamides are TRPV1 Antagonists

Several years ago, the crystallographic structures of the transient receptor potential vanilloid 1 (TRPV1) in the presence of agonists and antagonists were reported, providing structural information about its chemical activation and inactivation. TRPV1’s activation increases the transport of calcium and sodium ions, leading to the excitation of sensory neurons and the perception of pain. On the other hand, its antagonistic inactivation has been explored to design analgesic drugs. The interactions between the antagonists 5,5-diarylpentadienamides (DPDAs) and TRPV1 were studied here to explain why they inactivate TRPV1. The present work identified the structural features of TRPV1–DPDA complexes, starting with a consideration of the orientations of the ligands inside the TRPV1 binding site by using molecular docking. After this, a chemometrics analysis was performed (i) to compare the orientations of the antagonists (by using LigRMSD), (ii) to describe the recurrent interactions between the protein residues and ligand groups in the complexes (by using interaction fingerprints), and (iii) to describe the relationship between topological features of the ligands and their differential antagonistic activities (by using a quantitative structure–activity relationship (QSAR) with 2D autocorrelation descriptors). The interactions between the DPDA groups and the residues Y511, S512, T550, R557, and E570 (with a recognized role in the binding of classic ligands), and the occupancy of isoquinoline or 3-hydroxy-3,4-dihydroquinolin-2(1H)-one groups of the DPDAs in the vanilloid pocket of TRPV1 were clearly described. Based on the results, the structural features that explain why DPDAs inactivate TRPV1 were clearly exposed. These features can be considered for the design of novel TRPV1 antagonists.

Theoretical Studies on the Selectivity Mechanisms of Glycogen Synthase Kinase 3β (GSK3β) with Pyrazine ATP-competitive Inhibitors by 3DQSAR, Molecular Docking, Molecular Dynamics Simulation and Free Energy Calculations

Background

Glycogen synthase kinase-3 (GSK3) is associated with various key biological processes and has been considered as an important therapeutic target for the treatment of many diseases. Great efforts have been made on the development of GSK3 inhibitors, especially ATP-competitive GSK3β inhibitor, but it is still a great challenge to develop selective GSK3β inhibitors because of the high sequence homology with other kinases.

Objective

In order to reveal the selectivity mechanisms of GSK3β inhibition at the molecular level, a series of ATP-competitive GSK3β inhibitor was analyzed by a systematic computational method, combining 3D-QSAR, molecular docking, molecular dynamic simulations and free energy calculations.

Methods

Firstly, 3D-QSAR with CoMFA was built to explore the general structure activity relationships. Secondly, CDOCKER and Flexible docking were employed to predicted the reasonable docking poses of all studied inhibitors. And then, both GSK3β and CDK2 complexes were selected to conduct molecular dynamics simulations. Finally, the free energy calculations were employed to find the key selective-residues.

Results

CoMFA model suggested the steric, hydrophobic fields play key roles in the bioactivities of inhibitors, and the binding mechanisms were well analyzed through molecular docking. The binding free energies predicted are in good agreement with the experimental bioactivities and the free energy calculations showed that the binding of GSK3β/inhibitors was mainly contributed from hydrogen bonding and hydrophobic interaction.

Conclusion

Some key residues for selective binding were highlighted, which may afford important guidance for the rational design of novel ATP-competitive GSK3β inhibitors.

Design, synthesis, structure, in vitro cytotoxic activity evaluation and docking studies on target enzyme GSK-3β of new indirubin-3'-oxime derivatives

The addition of chalcone and amine components into indirubin-3′-oxime resulted in 15 new derivatives with high yields. Structures of new derivatives were also elucidated through 1D, 2D-NMR and HR-MS(ESI) spectra and X-ray crystallography. All designed compounds were screened for cytotoxic activity against four human cancer cell lines (HepG2, LU-1, SW480 and HL-60) and one human normal kidney cell line (HEK-293). Compound 6f exhibited the most marked cytotoxicity meanwhile cytotoxicity of compounds 6e, 6h and 6l was more profound toward cancer cell lines than toward normal cell. These new derivatives were further analyzed via molecular docking studies on GSK-3β enzyme. Docking analysis shows that most of the derivatives exhibited potential inhibition activity against GSK-3β with characteristic interacting residues in the binding site. The fast pulling of ligand scheme was then employed to refine the binding affinity and mechanism between ligands and GSK-3β enzyme. The computational results are expected to contribute to predicting enzyme target of the trial inhibitors and their possible interaction, from which the design of new cytotoxic agents could be created in the future.

Structural Requirements of N-alpha-Mercaptoacetyl Dipeptide (NAMdP) Inhibitors of Pseudomonas Aeruginosa Virulence Factor LasB: 3D-QSAR, Molecular Docking, and Interaction Fingerprint Studies

The zinc metallopeptidase Pseudomonas elastase (LasB) is a virulence factor of Pseudomonas aeruginosa (P. aeruginosa), a pathogenic bacterium that can cause nosocomial infections. The present study relates the structural analysis of 118 N-alpha-mercaptoacetyl dipeptides (NAMdPs) as LasB inhibitors. Field-based 3D-QSAR and molecular docking methods were employed to describe the essential interactions between NAMdPs and LasB binding sites, and the chemical features that determine their differential activities. We report a predictive 3D-QSAR model that was developed according to the internal and external validation tests. The best model, including steric, electrostatic, hydrogen bond donor, hydrogen bond acceptor, and hydrophobic fields, was found to depict a three-dimensional map with the local positive and negative effects of these chemotypes on the LasB inhibitory activities. Furthermore, molecular docking experiments yielded bioactive conformations of NAMdPs inside the LasB binding site. The series of NAMdPs adopted a similar orientation with respect to phosphoramidon within the LasB binding site (crystallographic reference), where the backbone atoms of NAMdPs are hydrogen-bonded to the LasB residues N112, A113, and R198, similarly to phosphoramidon. Our study also included a deep description of the residues involved in the protein-ligand interaction patterns for the whole set of NAMdPs, through the use of interaction fingerprints (IFPs).

Omega hydroxylated JA-Ile is an endogenous bioactive jasmonate that signals through the canonical jasmonate signaling pathway

Jasmonates are fatty acid derivatives that control several plant processes including growth, development and defense. Despite the chemical diversity of jasmonates, only jasmonoyl-L-isoleucine (JA-Ile) has been clearly characterized as the endogenous ligand of the jasmonate co-receptors (COI1-JAZs) in higher plants. Currently, it is accepted that ω-hydroxylation of JA-Ile leads to inactivation of the molecule. This study shows that ω-hydroxylated JA-Ile (12-OH-JA-Ile) retains bioactivity and signals through the canonical JA-pathway. The results suggest that 12-OH-JA-Ile differentially activates a subset of JA-Ile co-receptors that may control and/or modulate particular jasmonate dependent responses. It is proposed that after a strong immune response mediated by JA-Ile, the ω-hydroxylated form modulates JA-Ile activated processes thereby improving plant resilience.

Docking, Interaction Fingerprint, and Three-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) of Sigma1 Receptor Ligands, Analogs of the Neuroprotective Agent RC-33

The human Sigma1 receptor (S1R), which has been identified as a target with an important role in neuropsychological disorders, was first crystallized 3 years ago. Since S1R structure has no relation with another previous crystallized structures, the presence of the new crystal is an important hallmark for the design of agonists and antagonists against this important target. Some years ago, our group identified RC-33, a potent and selective S1R agonist, endowed with neuroprotective properties. In this work, drawing on new structural information, we studied the interactions of RC-33 and its analogs with the S1R binding site by using computational methods such as docking, interaction fingerprints, and receptor-guided alignment three dimensional quantitative structure–activity relationship (3D-QSAR). We found that RC-33 and its analogs adopted similar orientations within S1R binding site, with high similitude with orientations of the crystallized ligands; such information was used for identifying the residues involved in chemical interactions with ligands. Furthermore, the structure-activity relationship of the studied ligands was adequately described considering classical QSAR tests. All relevant aspects of the interactions between the studied compounds and S1R were covered here, through descriptions of orientations, binding interactions, and features that influence differential affinities. In this sense, the present results could be useful in the future design of novel S1R modulators.

Rational Design of Novel Phosphoinositide 3-Kinase Gamma (PI3Kγ) Selective Inhibitors: A Computational Investigation Integrating 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Phosphoinositide 3-kinase gamma (PI3Kγ) draws an increasing attention due to its link with deadly cancer, chronic inflammation and allergy. But the development of PI3Kγ selective inhibitors is still a challenging endeavor because of the high sequence homology with the other PI3K isoforms. In order to acquire valuable information about the interaction mechanism between potent inhibitors and PI3Kγ, a series of PI3Kγ isoform-selective inhibitors were analyzed by a systematic computational method, combining 3D-QSAR, molecular docking, molecular dynamic (MD) simulations, free energy calculations and decomposition. The general structure-activity relationships were revealed and some key residues relating to selectivity and high activity were highlighted. It provides precious guidance for rational virtual screening, modification and design of selective PI3Kγ inhibitors. Finally, ten novel inhibitors were optimized and P10 showed satisfactory predicted bioactivity, demonstrating the feasibility to develop potent PI3Kγ inhibitors through this computational modeling and optimization.

Mycobacterium tuberculosis serine/threonine protein kinases: structural information for the design of their specific ATP-competitive inhibitors

In the last decades, human protein kinases (PKs) have been relevant as targets in the development of novel therapies against many diseases, but the study of Mycobacterium tuberculosis PKs (MTPKs) involved in tuberculosis pathogenesis began much later and has not yet reached an advanced stage of development. To increase knowledge of these enzymes, in this work we studied the structural features of MTPKs, with focus on their ATP-binding sites and their interactions with inhibitors. PknA, PknB, and PknG are the most studied MTPKs, which were previously crystallized; ATP-competitive inhibitors have been designed against them in the last decade. In the current work, reported PknA, PknB, and PknG inhibitors were extracted from literature and their orientations inside the ATP-binding site were proposed by using docking method. With this information, interaction fingerprints were elaborated, which reveal the more relevant residues for establishing chemical interactions with inhibitors. The non-crystallized MTPKs PknD, PknF, PknH, PknJ, PknK, and PknL were also studied; their three-dimensional structural models were developed by using homology modeling. The main characteristics of MTPK ATP-binding sites (the non-crystallized and crystallized MTPKs, including PknE and PknI) were accounted; schemes of the main polar and nonpolar groups inside their ATP-binding sites were constructed, which are suitable for a major understanding of these proteins as antituberculotic targets. These schemes could be used for establishing comparisons between MTPKs and human PKs in order to increase selectivity of MTPK inhibitors. As a key tool for guiding medicinal chemists interested in the design of novel MTPK inhibitors, our work provides a map of the structural elements relevant for the design of more selective ATP-competitive MTPK inhibitors.

Insights into the Structural Requirements of 2(S)-Amino-6-Boronohexanoic Acid Derivatives as Arginase I Inhibitors: 3D-QSAR, Docking, and Interaction Fingerprint Studies

Human arginase I (hARGI) is an important enzyme involved in the urea cycle; its overexpression has been associated to cardiovascular and cerebrovascular diseases. In the last years, several congeneric sets of hARGI inhibitors have been reported with possible beneficial roles for the cardiovascular system. At the same time, crystallographic data have been reported including hARGI–inhibitor complexes, which can be considered for the design of novel inhibitors. In this work, the structure–activity relationship (SAR) of Cα substituted 2(S)-amino-6-boronohexanoic acid (ABH) derivatives as hARGI inhibitors was studied by using a three-dimensional quantitative structure–activity relationships (3D-QSAR) method. The predictivity of the obtained 3D-QSAR model was demonstrated by using internal and external validation experiments. The best model revealed that the differential hARGI inhibitory activities of the ABH derivatives can be described by using steric and electrostatic fields; the local effects of these fields in the activity are presented. In addition, binding modes of the above-mentioned compounds inside the hARGI binding site were obtained by using molecular docking. It was found that ABH derivatives adopted the same orientation reported for ABH within the hARGI active site, with the substituents at Cα exposed to the solvent with interactions with residues at the entrance of the binding site. The hARGI residues involved in chemical interactions with inhibitors were identified by using an interaction fingerprints (IFPs) analysis.

Flavonoids as CDK1 Inhibitors: Insights in Their Binding Orientations and Structure-Activity Relationship

In the last years, the interactions of flavonoids with protein kinases (PKs) have been described by using crystallographic experiments. Interestingly, different orientations have been found for one flavonoid inside different PKs and different chemical substitutions lead to different orientations of the flavonoid scaffold inside one PK. Accordingly, orientation predictions of novel analogues could help to the design of flavonoids with high PK inhibitory activities. With this in mind, we studied the binding modes of 37 flavonoids (flavones and chalcones) inside the cyclin-dependent PK CDK1 using docking experiments. We found that the compounds under study adopted two different orientations into the active site of CDK1 (orientations I and II in the manuscript). In addition, quantitative structure-activity relationship (QSAR) models using CoMFA and CoMSIA methodologies were constructed to explain the trend of the CDK1 inhibitory activities for the studied flavonoids. Template-based and docking-based alignments were used. Models developed starting from docking-based alignment were applied for describing the whole dataset and compounds with orientation I. Adequate R2 and Q2 values were obtained by each method; interestingly, only hydrophobic and hydrogen bond donor fields describe the differential potency of the flavonoids as CDK1 inhibitors for both defined alignments and subsets. Our current application of docking and QSAR together reveals important elements to be drawn for the design of novel flavonoids with increased PK inhibitory activities.

Is It Reliable to Use Common Molecular Docking Methods for Comparing the Binding Affinities of Enantiomer Pairs for Their Protein Target?

Molecular docking is a computational chemistry method which has become essential for the rational drug design process. In this context, it has had great impact as a successful tool for the study of ligand–receptor interaction modes, and for the exploration of large chemical datasets through virtual screening experiments. Despite their unquestionable merits, docking methods are not reliable for predicting binding energies due to the simple scoring functions they use. However, comparisons between two or three complexes using the predicted binding energies as a criterion are commonly found in the literature. In the present work we tested how wise is it to trust the docking energies when two complexes between a target protein and enantiomer pairs are compared. For this purpose, a ligand library composed by 141 enantiomeric pairs was used, including compounds with biological activities reported against seven protein targets. Docking results using the software Glide (considering extra precision (XP), standard precision (SP), and high-throughput virtual screening (HTVS) modes) and AutoDock Vina were compared with the reported biological activities using a classification scheme. Our test failed for all modes and targets, demonstrating that an accurate prediction when binding energies of enantiomers are compared using docking may be due to chance. We also compared pairs of compounds with different molecular weights and found the same results.