Molecular docking and MD simulation studies of 4-thiazol-N-(pyridin-2-yl)pyrimidin-2-amine derivatives as novel inhibitors targeted to CDK2/4/6.

PURPOSE: Nowadays, cyclin-dependent kinase 4/6 (CDK4/6) inhibitors have been approved for treating metastatic breast cancer and have achieved inspiring curative effects. But some discoveries have indicated that CDK 4/6 are not the requisite factors in some cell types because CDK2 partly compensates for the inhibition of CDK4/6. Thus, it is urgent to design CDK2/4/6 inhibitors for significantly enhancing their potency. This study aims to explore the mechanism of the binding of CDK2/4/6 kinases and their inhibitors to design novel CDK2/4/6 inhibitors for significantly enhancing their potency in different kinds of cancers. MATERIALS AND METHODS: A series of 72 disparately functionalized 4-substituted N-phenylpyrimidin-2-amine derivatives exhibiting potent inhibitor activities against CDK2, CDK4 and CDK6 were collected to apply to this research. The total set of these derivatives was divided into a training set (54 compounds) and a test set (18 compounds). The derivatives were constructed through the sketch molecule module in SYBYL 6.9 software. A Powell gradient algorithm and Tripos force field were used to calculate the minimal structural energy and the minimized structure was used as the initial conformation for molecular docking. By the means of 3D-QSAR models, partial least squares (PLS) analysis, molecular dynamics (MD) simulations and binding free energy calculations, we can find the relationship between structure and biological activity. RESULTS: In this study, we used molecular docking, 3D-QSAR and molecular dynamics simulation methods to comprehensively analyze the interaction and structure-activity relationships of 72 new CDK2/4/6 inhibitors. We used detailed statistical data to reasonably verify the constructed 3D-QSAR models for three receptors (q2 of CDK2 = 0.714, R2pred = 0.764, q2 = 0.815; R2pred of CDK4 = 0.681, q2 = 0.757; R2pred of CDK6 = 0.674). MD simulations and decomposition energy analysis validated the reasonability of the docking results and identified polar interactions as crucial factors that influence the different bioactivities of the studied inhibitors of CDK2/4/6 receptors, especially the electrostatic interactions of Lys33/35/43 and Asp145/158/163. The nonpolar interaction with Ile10/12/19 was also critical for the differing potencies of the CDK2/4/6 inhibitors. We concluded that the following probably enhanced the bioactivity against CDK2/4/6 kinases: (1) electronegative groups at the N1-position and electropositive and moderate-sized groups at ring E; (2) electrogroups featured at R2; (3) carbon atoms at the X-position or ring C replaced by a benzene ring; and (4) an electrogroup as R4. CONCLUSION: Previous studies, to our knowledge, only utilized a single approach of 3D-QSAR and did not integrate this method with other sophisticated techniques such as molecular dynamics simulations to discover new potential inhibitors of CDK2, CDK4, or CDK6. So we applied the intergenerational technology, such as 3D-QSAR technology, molecular docking simulation techniques, molecular dynamics simulations and MMPBSA19/MMGBSA20-binding free energy calculations to statistically explore the correlations between the structure with biological activities. The constructed 3D-QSAR models of the three receptors were reasonable and confirmed by the excellent statistical data. We hope the results obtained from this work will provide some useful references for the development of novel CDK2/4/6 inhibitors.

The molecular basis of drug selectivity for α5 subunit-containing GABAA receptors.

α5 subunit-containing γ-aminobutyric acid type A (GABAA) receptors represent a promising drug target for neurological and neuropsychiatric disorders. Altered expression and function contributes to neurodevelopmental disorders such as Dup15q and Angelman syndromes, developmental epilepsy and autism. Effective drug action without side effects is dependent on both α5-subtype selectivity and the strength of the positive or negative allosteric modulation (PAM or NAM). Here we solve structures of drugs bound to the α5 subunit. These define the molecular basis of binding and α5 selectivity of the ß-carboline, methyl 6,7-dimethoxy-4-ethyl-ß-carboline-3-carboxylate (DMCM), type II benzodiazepine NAMs, and a series of isoxazole NAMs and PAMs. For the isoxazole series, each molecule appears as an 'upper' and 'lower' moiety in the pocket. Structural data and radioligand binding data reveal a positional displacement of the upper moiety containing the isoxazole between the NAMs and PAMs. Using a hybrid molecule we directly measure the functional contribution of the upper moiety to NAM versus PAM activity. Overall, these structures provide a framework by which to understand distinct modulator binding modes and their basis of α5-subtype selectivity, appreciate structure-activity relationships, and empower future structure-based drug design campaigns.

A general chemical crosslinking strategy for structural analyses of weakly interacting proteins applied to preTCR-pMHC complexes.

T lymphocytes discriminate between healthy and infected or cancerous cells via T-cell receptor-mediated recognition of peptides bound and presented by cell-surface-expressed major histocompatibility complex molecules (MHCs). Pre-T-cell receptors (preTCRs) on thymocytes foster development of αßT lymphocytes through their ß chain interaction with MHC displaying self-peptides on thymic epithelia. The specific binding of a preTCR with a peptide-MHC complex (pMHC) has been identified previously as forming a weak affinity complex with a distinct interface from that of mature αßTCR. However, a lack of appropriate tools has limited prior efforts to investigate this unique interface. Here we designed a small-scale linkage screening protocol using bismaleimide linkers for determining residue-specific distance constraints between transiently interacting protein pairs in solution. Employing linkage distance restraint-guided molecular modeling, we report the oriented solution docking geometry of a preTCRß-pMHC interaction. The linkage model of preTCRß-pMHC complex was independently verified with paramagnetic pseudocontact chemical shift (PCS) NMR of the unlinked protein mixtures. Using linkage screens, we show that the preTCR binds with differing affinities to peptides presented by MHC in solution. Moreover, the C-terminal peptide segment is a key determinant in preTCR-pMHC recognition. We also describe the process for future large-scale production and purification of the linked constructs for NMR, X-ray crystallography, and single-molecule electron microscopy studies.

Nearest-neighbor NMR spectroscopy: categorizing spectral peaks by their adjacent nuclei.

Methyl-NMR enables atomic-resolution studies of structure and dynamics of large proteins in solution. However, resonance assignment remains challenging. The problem is to combine existing structural informational with sparse distance restraints and search for the most compatible assignment among the permutations. Prior classification of peaks as either from isoleucine, leucine, or valine reduces the search space by many orders of magnitude. However, this is hindered by overlapped leucine and valine frequencies. In contrast, the nearest-neighbor nuclei, coupled to the methyl carbons, resonate in distinct frequency bands. Here, we develop a framework to imprint additional information about passively coupled resonances onto the observed peaks. This depends on simultaneously orchestrating closely spaced bands of resonances along different magnetization trajectories, using principles from control theory. For methyl-NMR, the method is implemented as a modification to the standard fingerprint spectrum (the 2D-HMQC). The amino acid type is immediately apparent in the fingerprint spectrum. There is no additional relaxation loss or an increase in experimental time. The method is validated on biologically relevant proteins. The idea of generating new spectral information using passive, adjacent resonances is applicable to other contexts in NMR spectroscopy.

Cytotoxicity and Anti-inflammatory Properties of Apigenin-Derived Isolaxifolin.

The rare flavonoid isolaxifolin, a potent insecticide, has been touted as a potential grain-protecting agent. In order to assess any impact of this natural product on human health and to explore its various other biological properties, we have established a semisynthesis from the simpler but structurally related and more abundant natural product apigenin. The five-step reaction sequence has provided, for the first time, sufficient material for an in-depth evaluation of the cytotoxic properties of the title natural product. The impact of isolaxifolin on certain pro-inflammatory cytokines in murine macrophage RAW 264.7 cells has also been examined. Such studies have revealed that isolaxifolin displays no toxic effects toward normal cells while displaying greater cytotoxicities against certain cancer cell lines than its synthetic precursor apigenin. Furthermore, unlike apigenin, isolaxifolin only reduced NO, TNF-α, and IL-6 secretions in LPS-induced RAW 264.7 cells in a rather modest and dose-independent manner.

Using tert-Butyl Groups in a Ligand To Identify Its Binding Site on a Protein.

Few methods allow determining the binding site of tightly binding ligands. We show that ligands containing a tert-butyl (e.g., Boc) group produce easily observable nuclear Overhauser effects (NOE) with the target protein even when the tert-butyl group is not highly solvent exposed. NOEs with methyl groups of the target protein are readily assigned by selectively isotope labeling, presenting a practical and quick way to pinpoint the location of the ligand without any prior specific nuclear magnetic resonance assignments of the protein. The approach works for nonexchanging ligands as well as for weakly binding ligands.

Probing the stereoselectivity of OleD-catalyzed glycosylation of cardiotonic steroids.

The glycosyltransferase OleD variant as a catalyst for the glycosylation of four pairs of epimers of cardiotonic steroids (CTS) are assessed. The results of this study demonstrated that the OleD-catalyze glycosylation of CTS is significantly influenced by the configuration at C-3 and the A/B fusion mode. 3ß-OH and A/B ring cis fusion are favoured by OleD (ASP). An epoxide ring at C-14 and C-15 further increases the bioconversion rate; while an acetyl group at C-16 and lactone ring type at C-17 did not influence the biotransformation. A high conversion rate corresponded to a low K m value. A molecular docking simulation showed that filling of hydrophobic pocket II and interaction with residue Tyr115 may play an important role in the glycosylation reactions catalyzed by OleD glycosyltransferases. Furthermore, the glycosylation products showed a stronger inhibitory activity for Na+, K+-ATPase than the corresponding aglycones. This study provides the first stereoselective properties for OleD (ASP) catalyzed glycosylation.

Sensitive NMR Approach for Determining the Binding Mode of Tightly Binding Ligand Molecules to Protein Targets.

Structure-guided drug design relies on detailed structural knowledge of protein-ligand complexes, but crystallization of cocomplexes is not always possible. Here we present a sensitive nuclear magnetic resonance (NMR) approach to determine the binding mode of tightly binding lead compounds in complex with difficult target proteins. In contrast to established NMR methods, it does not depend on rapid exchange between bound and free ligand or on stable isotope labeling, relying instead on a tert-butyl group as a chemical label. tert-Butyl groups are found in numerous protein ligands and deliver an exceptionally narrow and tall (1)H NMR signal. We show that a tert-butyl group also produces outstandingly intense intra- and intermolecular NOESY cross-peaks. These enable measurements of pseudocontact shifts generated by lanthanide tags attached to the protein, which in turn allows positioning of the ligand on the protein. Once the ligand has been located, assignments of intermolecular NOEs become possible even without prior resonance assignments of protein side chains. The approach is demonstrated with the dengue virus NS2B-NS3 protease in complex with a high-affinity ligand containing a tert-butyl group.

O-tert-Butyltyrosine, an NMR tag for high-molecular-weight systems and measurements of submicromolar ligand binding affinities.

O-tert-Butyltyrosine (Tby) is an unnatural amino acid that can be site-specifically incorporated into proteins using established orthogonal aminoacyl-tRNA synthetase/tRNA systems. Here we show that the tert-butyl group presents an outstanding NMR tag that can readily be observed in one-dimensional (1)H NMR spectra without any isotope labeling. Owing to rapid bond rotations and the chemical equivalence of the protons of a solvent-exposed tert-butyl group from Tby, the singlet resonance from the tert-butyl group generates an easily detectable narrow signal in a spectral region with limited overlap with other methyl resonances. The potential of the tert-butyl (1)H NMR signal in protein research is illustrated by the observation and assignment of two resonances in the Bacillus stearothermophilus DnaB hexamer (320 kDa), demonstrating that this protein preferentially assumes a 3-fold rather than 6-fold symmetry in solution, and by the quantitative measurement of the submicromolar dissociation constant Kd (0.2 µM) of the complex between glutamate and the Escherichia coli aspartate/glutamate binding protein (DEBP, 32 kDa). The outstanding signal height of the (1)H NMR signal of the Tby tert-butyl group allows Kd measurements using less concentrated protein solutions than usual, providing access to Kd values 1 order of magnitude lower than established NMR methods that employ direct protein detection for Kd measurements.

The dengue virus NS2B-NS3 protease retains the closed conformation in the complex with BPTI.

The C-terminal ß-hairpin of NS2B (NS2Bc) in the dengue virus NS2B-NS3 protease is required for full enzymatic activity. In crystal structures without inhibitor and in the complex with bovine pancreatic trypsin inhibitor (BPTI), NS2Bc is displaced from the active site. In contrast, nuclear magnetic resonance (NMR) studies in solution only ever showed NS2Bc in the enzymatically active closed conformation. Here we demonstrate by pseudocontact shifts from a lanthanide tag that NS2Bc remains in the closed conformation also in the complex with BPTI. Therefore, the closed conformation is the best template for drug discovery.

Binding mode of the activity-modulating C-terminal segment of NS2B to NS3 in the dengue virus NS2B-NS3 protease.

The two-component dengue virus NS2B-NS3 protease (NS2B-NS3pro) is an established drug target but inhibitor design is hampered by uncertainties about its 3D structure in solution. Crystal structures reported very different conformations for the functionally important C-terminal segment of the NS2B cofactor (NS2Bc), indicating open and closed conformations in the absence and presence of inhibitors, respectively. An earlier NMR study in solution indicated that a closed state is the preferred conformation in the absence of an artificial linker engineered between NS2B and NS3pro. To obtain direct structural information on the fold of unlinked NS2B-NS3pro in solution, we tagged NS3pro with paramagnetic tags and measured pseudocontact shifts by NMR to position NS2Bc relative to NS3pro. NS2Bc was found to bind to NS3pro in the same way as reported in a previously published model and crystal structure of the closed state. The structure is destabilized, however, by high ionic strength and basic pH, showing the importance of electrostatic forces to tie NS2Bc to NS3pro. Narrow NMR signals previously thought to represent the open state are associated with protein degradation. In conclusion, the closed conformation of the NS2B-NS3 protease is the best model for structure-guided drug design.

Synthesis, molecular docking, and biofilm formation inhibitory activity of 5-substituted 3,4-dihalo-5h-furan-2-one derivatives on Pseudomonas aeruginosa.

Pseudomonas aeruginosa (P. aeruginosa) colonize on most wounds and live as biofilm, which causes antibiotic resistance and wounds unhealed. To investigate the effects of 5-substituted 3,4-dihalo-5H-furan-2-one compounds on biofilm formation of P. aeruginosa, a set of 5-(aryl-1'-hydroxy-methyl)- or 5-(aryl-2-methylene)-3,4-dihalo-5H-furan-2-one compounds were designed and synthesized. Their inhibitory activities on biofilm formation of P. aeruginosa were studied by MIC assay, quantitative analysis of biofilm inhibition, and observation of biofilm formation with SEM. It was found that compounds 2i, 3f, 3i showed remarkable effects of biofilm formation inhibition on P. aeruginosa. Furthermore, molecular docking was performed to identify the key structural features of these compounds with the binding site of LasR receptor.

Molecular modeling studies on phosphonic acid-containing thiazole derivatives: design for fructose-1,6-bisphosphatase inhibitors.

Presently, an in silico modeling was carried out on a series of 63 phosphonic acid-containing thiazole derivatives as fructose-1,6-bisphosphatase (FBPase) inhibitors using CoMFA/CoMSIA and molecular docking methods. The CoMFA and CoMSIA models using 51 molecules in the training set gave r (cv) (2) values of 0.675 and 0.619, r ( 2 ) values of 0.985 and 0.979, respectively. The systemic external validation indicated that our CoMFA and CoMSIA models possessed high predictive powers with r (0) (2) values of 0.995 and 0.994, r (m(test)) (2) values of 0.887 and 0.860, respectively. The 3D contour maps of the CoMFA and CoMSIA provided smooth and interpretable explanation of the structure-activity relationship for the inhibitors. Molecular docking studies revealed that a phosphonic group was essential for binding to the AMP binding site, and some key features were also identified. The analyses of the 3D contour plots and molecular docking results permitted interesting conclusions about the effects of different substituent groups at different positions of the common scaffold, which might guide the design of novel FBPase inhibitors with higher activity and bioavailability. A set of 60 new analogues were designed by utilizing the results revealed in the present study, and were predicted with significantly improved potencies in the developed models. The findings can be quite useful to aid the designing of new fructose-1,6-biphophatase inhibitors with improved biological response.

3D-QSAR and molecular docking studies of azaindole derivatives as Aurora B kinase inhibitors.

The Aurora kinases have been regarded as attractive targets for the development of new anticancer agents. Recently a series of azaindole derivatives with Aurora B inhibitory activities were reported. To explore the relationship between the structures of substituted azaindole derivatives and their inhibition of Aurora B, 3D-QSAR and molecular docking studies were performed on a dataset of 41 compounds. 3D-QSAR, including CoMFA and CoMSIA, were applied to identify the key structures impacting their inhibitory potencies. The CoMSIA model showed better results than CoMFA, with r(2)(cv) value of 0.575 and r(2) value of 0.987. 3D contour maps generated from CoMFA and CoMSIA along with the docking binding structures provided enough information about the structural requirements for better activity. Based on the structure-activity relationship revealed by the present study, we have designed a set of novel Aurora B inhibitors that showed excellent potencies in the developed models. Thus, our results allowed us to design new derivatives with desired activities.

Molecular modelling studies on d-annulated benzazepinones as VEGF-R2 kinase inhibitors using docking and 3D-QSAR.

Chemotypes comprising the d-annulated 1,3-dihydro-2H-1-benzazepin-2-one scaffold derived from the paullone structure were found to be potent vascular endothelial growth factor receptor 2 (VEGF-R2) kinase inhibitors. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies were performed on a series of d-annulated benzazepinones with VEGF-R2 kinase inhibition activities. The comparative molecular field analysis and comparative molecular similarity indices analysis models using 32 molecules in the training set gave r(2)(cv) values of 0.811 and 0.769, r(2) values of 0.962 and 0.953, respectively. 3D contour maps generated from the two models revealed that the electron-withdrawing groups at R(1) and the bulky, electron-withdrawing as well as hydrogen bond donor groups at R(2) position are favourable; the bulky, hydrogen bond acceptor substituent at R(3) and the minor groups at R(4) position may benefit the potency. We have designed a series of novel VEGF-R2 inhibitors by utilizing the SAR results revealed in the present study, which were predicted with excellent potencies in the developed models. The results may aid in designing of potential VEGF-R2 inhibitors with better activities.

Understanding the structure-activity relationship of betulinic acid derivatives as anti-HIV-1 agents by using 3D-QSAR and docking.

Novel anti-HIV-1 agents derived from betulinic acid have been greatly concerned. 3D-QSAR and molecular docking studies were applied to rationalize the structural requirements responsible for the anti-HIV activity of these compounds. The CoMFA and CoMSIA models resulted from 28 molecules gave r(cv)² values of 0.599 and 0.630, r² values of 0.994 and 0.958, respectively. To estimate the predictive ability of the 3D-QSAR model, an external validation was employed. Based on the contour maps generated from both CoMFA and CoMSIA, we have identified some key features in the betulinic acid derivatives that are responsible for the anti-HIV activity. Molecular docking was used to explore the binding mode between these derivatives and HIV gp120. We have therefore designed a series of novel betulinic acid derivatives by utilizing the SAR results revealed in the present study, which were predicted with excellent potencies in the developed models. The results provide a valuable method to design new betulinic acid derivatives as anti-HIV-1 agents.

Molecular modeling studies on imidazo[4,5-b]pyridine derivatives as Aurora A kinase inhibitors using 3D-QSAR and docking approaches.

3D-QSAR and docking studies were performed on sixty imidazo[4,5-b]pyridine derivatives as Aurora A kinase inhibitors. The CoMFA and CoMSIA models using forthy-eight molecules in the training set, gave r(cv)(2) values of 0.774 and 0.800, r(2) values of 0.975 and 0.977, respectively. The external validation indicated that both CoMFA and CoMSIA models possessed high predictive powers with r(pred)(2) values of 0.933 and 0.959, r(m)(2) values of 0.883 and 0.915, respectively. 3D contour maps generated from the two models along with docking binding structures have identified several key structural requirements responsible for the activity. A set of thirty new analogues were proposed by utilizing the results revealed in the present study, and were predicted with significantly improved potencies in the developed models.

3D-QSAR studies and molecular docking on [5-(4-amino-1H-benzoimidazol-2-yl)-furan-2-yl]-phosphonic acid derivatives as fructose-1,6-biphophatase inhibitors.

Fructose-1,6-biphophatase has been regarded as a novel therapeutic target for the treatment of type 2 diabetes mellitus (T2DM). 3D-QSAR and docking studies were performed on a series of [5-(4-amino-1H-benzoimidazol-2-yl)-furan-2-yl]-phosphonic acid derivatives as fructose-1,6-biphophatase inhibitors. The CoMFA and CoMSIA models using thirty-seven molecules in the training set gave r (cv) (2) values of 0.614 and 0.598, r (2) values of 0.950 and 0.928, respectively. The external validation indicated that our CoMFA and CoMSIA models possessed high predictive powers with r (0) (2) values of 0.994 and 0.994, r (m) (2) values of 0.751 and 0.690, respectively. Molecular docking studies revealed that a phosphonic group was essential for binding to the receptor, and some key features were also identified. A set of forty new analogues were designed by utilizing the results revealed in the present study, and were predicted with significantly improved potencies in the developed models. The findings can be quite useful to aid the designing of new fructose-1,6-biphophatase inhibitors with improved biological response.

3D-QSAR and docking studies on pyrazolo[4,3-h]qinazoline-3-carboxamides as cyclin-dependent kinase 2 (CDK2) inhibitors.

3D-QSAR and docking studies were performed on a series of pyrazolo[4,3-h]quinazoline-3-carboxamides as CDK2/CyA inhibitors. The CoMFA and CoMSIA models using 54 molecules in the training set, gave r(cv)(2) values of 0.644 and 0.507, r(2) values of 0.959 and 0.951, respectively. 3D contour maps generated from the two models were applied to identify features important for the activity and better understand the interaction between the inhibitors and the receptor. Molecular docking was employed to explore the binding mode between these compounds and the receptor, as well as help understanding the structure-activity relationship revealed by CoMFA and CoMSIA. The results provide a useful guideline for the rational design of novel CDKs inhibitors.