3D-QSAR studies of Dipeptidyl peptidase IV inhibitors using a docking based alignment

Dipeptidyl peptidase IV (DPP-IV) deactivates the incretin hormones GLP-1 and GIP by cleaving the penultimate proline or alanine from the N-terminal (P1-position) of the peptide. Inhibition of this enzyme will prevent the degradation of the incretin hormones and maintain glucose homeostasis; this makes it an attractive target for the development of drugs for diabetes. This paper reports 3D-QSAR analysis of several DPP-IV inhibitors, which were aligned by the receptor-based technique. The conformation of the molecules in the active site was obtained through docking methods. The QSAR models were generated on two training sets composed of 74 and 25 molecules which included phenylalanine, thiazolidine, and fluorinated pyrrolidine analogs. The 3D-QSAR models are robust with statistically significant r², q², and values. The CoMFA and CoMSIA models were used to design some new inhibitors with several fold higher binding affinity. Figure The CoMFA contours around molecule D1T155 (a) steric contours - favored (green); disfavored (yellow) (b) electrostatic contours - electropositive (blue); electronegative (red)     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

Comparative QSAR studies on peptide deformylase inhibitors

Comparative quantitative structure–activity relationship (QSAR) analyses of peptide deformylase (PDF) inhibitors were performed with a series of previously published (British Biotech Pharmaceuticals, Oxford, UK) reverse hydroxamate derivatives having antibacterial activity against Escherichia coli PDF, using 2D and 3D QSAR methods, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and hologram QSAR (HQSAR). Statistically reliable models with good predictive power were generated from all three methods (CoMFA r ² = 0.957, q ² = 0.569; CoMSIA r ² = 0.924, q ² = 0.520; HQSAR r ² = 0.860, q ² = 0.578). The predictive capability of these models was validated by a set of compounds that were not included in the training set. The models based on CoMFA and CoMSIA gave satisfactory predictive r ² values of 0.687 and 0.505, respectively. The model derived from the HQSAR method showed a low predictability of 0.178 for the test set. In this study, 3D prediction models showed better predictive power than 2D models for the test set. This might be because 3D information is more important in the case of datasets containing compounds with similar skeletons. Superimposition of CoMFA contour maps in the active site of the PDF crystal structure showed a meaningful correlation between receptor–ligand binding and biological activity. The final QSAR models, along with information gathered from 3D contour and 2D contribution maps, could be useful for the design of novel active inhibitors of PDF. Figure Superimposition of comparative molecular field analysis (CoMFA) contour plot in the active site of peptide deformylase (PDF)     

3D-QSAR and docking studies of 3-arylquinazolinethione derivatives as selective estrogen receptor modulators

3D-QSAR and molecular docking analysis were performed to explore the interaction of estrogen receptors (ERα and ERβ) with a series of 3-arylquinazolinethione derivatives. Using the conformations of these compounds revealed by molecular docking, CoMFA analysis resulted in the first quantitative structure-activity relationship (QSAR) and first quantitative structure-selectivity relationship (QSSR) models predicting the inhibitory activity against ERβ and the selectivity against ERá. The q² and R² values, along with further testing, indicate that the obtained 3D-QSAR and 3D-QSSR models will be valuable in predicting both the inhibitory activity and selectivity of 3-arylquinazolinethione derivatives for these protein targets. A set of 3D contour plots drawn based on the 3D-QSAR and 3D-QSSR models reveal modifications of substituents at C2 and C5 of the quinazoline which my be useful to improve both the activity and selectivity of ERβ/ ERα. Results showed that both the steric and electrostatic factors should appropriately be taken into account in future rational design and development of more active and more selective ERβ inhibitors for the therapeutic treatment of osteoporosis. Figure Structures of ERβ binding with compounds 1aar, 1ax and 1aag obtained from molecular docking     

Molecular dynamics study of the initial stages of catalyzed single-wall carbon nanotubes growth: force field development

Effective force fields for Ni-C interactions developed by Yamaguchi and Maruyama for the formation of metallofullerenes are modified to simulate the catalyzed growth of single-wall carbon nanotubes on Ni_n clusters with n >20, and the reactive empirical bond order Brenner potential for C-C interactions is also revised to include the effect of the metal atoms on such interactions. Figure Force field parameters for carbon-metal interactions obtained from DFT calculations in small clusters.     

Homology modeling and examination of the effect of the D92E mutation on the H5N1 nonstructural protein NS1 effector domain

Virulent H5N1 strains of influenza virus often harbor a D92E point mutation in the nonstructural protein NS1. This crucial mutation has been correlated with increased virulence and/or cytokine resistance, but the structural implications of such a change are still unclear. Furthermore, NS1 protein could also be a potential target for the development of novel antiviral agents against H5N1 strains. Therefore, a reasonable 3D model of H5N1 NS1 is important for the understanding of the molecular basis of increased virulence and the design of novel antiviral agents. Based on the crystal structure of a non-H5N1 NS1 protein, a model of H5N1 NS1 was developed by homology modeling, molecular mechanics and molecular dynamics simulations. It was found that the D92E mutation could result in weakened interactions of the carboxylate side chain with other phosphorylated residues, thereby activating phosphorylation of NS1. Figure Superposition of snapshots picked from the two molecular dynamic (MD) trajectories: a H5N1 NS1 homology model and b non-H5N1 NS1 crystal structure after 0 (green ribbon), 5 (blue ribbon) and 10 ns (pink ribbon) MD simulation     

Molecular docking guided 3D-QSAR CoMFA analysis of N-4-Pyrimidinyl-1H-indazol-4-amine inhibitors of leukocyte-specific protein tyrosine kinase

Inhibition of leukocyte-specific protein tyrosine kinase (Lck) activity offers one of the approaches for the treatment of T-cell mediated inflammatory disorders including rheumatoid arthritis, transplant rejection and inflammatory bowel disease. To explore the relationship between the structures of the N-4 Pyrimidinyl-1H-indazol-4-amines and their Lck inhibition, 3D-QSAR study using CoMFA analysis have been performed on a dataset of 42 molecules. The bioactive conformation of the template molecule, selected as the most potent molecule 23 from the series was obtained by performing molecular docking at the ATP binding site of Lck, which is then used to build the rest of the molecules in the series. The constructed CoMFA model is robust with of 0.603 and conventional r² of 0.983. The predictive power of the developed model was obtained using a test set of 10 molecules, giving predictive correlation coefficient of 0.921. CoMFA contour analysis was performed to obtain useful information about the structural requirements for the Lck inhibitors which could be utilized in its future design. Figure CoMFA steric contour map. Sterically favored areas (contribution level 80%) are represented by green polyhedra. Sterically disfavored areas (contribution level 20%) are represented by yellow polyhedra. The active molecule 23 shown in capped sticks. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Docking studies suggest ligand-specific δ-opioid receptor conformations

An automated docking procedure was used to study binding of a series of δ-selective ligands to three models of the δ-opioid receptor. These models are thought to represent the three ligand-specific receptor conformations. Docking results are in agreement with point mutation studies and suggest that different ligands—agonists and antagonists—may bind to the same binding site under different receptor conformations. Docking to different receptor models (conformations) also suggests that by changing to a receptor-specific conformation, the receptor may open or close different binding sites to other ligands. Figure  Ligands 5 (green) and 6 (orange) in bindingpocket BP1 of the R1 δ-opioid receptor model     

Aromaticity and electronic properties of Heterosuperbenzene (Heterohexabenzocoronene)

A global electrophilicity parameter and the aromaticity of some heterocyclic polyaromatic hydrocarbons were evaluated on the basis of DFT calculations. The substitution of carbon atoms by nitrogen atoms dramatically changes the global electrophilicity of the molecules, with the fully substituted molecule being the most electrophilic with a reactivity very close to that of fullerene. Figure Fully substituted heterohexabenzocoronene Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Novel distance-based atom-type topological indices DAI for QSPR/QSAR studies of alcohols

In this work, we propose a distance-based atom-type topological index (DAI) for quantitative structure-property/activity relationship (QSPR/QSAR) studies. The newly constructed index, which codes the structural environment of each atom type in a molecule, can be calculated simply. These atom-type topological indices, along with our recently proposed Lu index, were used to construct QSPR/QSAR models for several representative physical properties and biological activities of several data sets of alcohols with a range of non-hydrogen atoms by using multiple linear regression (MLR) analysis. The efficiency of these indices is verified by high quality QSPR models. The results indicate that the combined use of Lu and DAI indices promises to be a useful method for QSPR/QSAR analysis of complex compounds.      

CHIH-DFT determination of the molecular structure and IR and UV spectra of solanidine

Solanidine is the steroidal aglycon of some potato glycoalkaloids and a very important precursor for the synthesis of hormones and some pharmacologically active compounds. In this work, we make use of a new chemistry model within Density Functional Theory, called CHIH-DFT, to calculate the molecular structure of solanidine, as well to predict its infrared and ultraviolet spectra. The calculated values are compared with the experimental data available for this molecule as a means of validation of our proposed chemistry model. Figure Molecular structure of solanidine calculated with the CHIH(small) model chemistry     

Molecular models of protein kinase 6 from Plasmodium falciparum

Cyclin-dependent kinases (CDKs) have been identified as potential targets for development of drugs, mainly against cancer. These studies generated a vast library of chemical inhibitors of CDKs, and some of these molecules can also inhibit kinases identified in the Plasmodium falciparum genome. Here we describe structural models for Protein Kinase 6 from P. falciparum (PfPK6) complexed with Roscovitine and Olomoucine. These models show clear structural evidence for differences observed in the inhibition, and may help designing inhibitors for PfPK6 generating new potential drugs against malaria. Figure Ribbon diagram of PfPK6 complexed with a roscovitine and b olomoucine     

Evaluation of the chemical reactivity in lignin precursors using the Fukui function

The hydroxycinnamyl alcohols: p-coumarol, coniferol and sinapol are considered the basic units and precursors of lignins models. In this work, the specific reactivity of these molecules was studied. We investigate their intrinsic chemical reactivity in terms of the Fukui function, applying the principle of hard and soft acids and bases (HSAB) in the framework of the density functional theory (DFT). Comparisons of their nucleophilic, electrophilic and free radical reactivity show their most probably sites to form linkages among them. It is found that the most reactive sites, for reactions involving free radicals, are the carbons at the β-position in the p-coumarol and sinapol molecules, whilst the regions around the carbon-oxygen bond of the phenoxyl group are the most reactive in coniferol. Figure Isocontour plots for the free radical form of the Fukui function f ⁰ (r), showing the reactive sites toward electron-rich/poor reactants     

Topochemical models for the prediction of voltage-gated sodium channel binding activity of hydantoins and related non-hydantoins

The relationship of Wiener’s topochemical index—a distance based topochemical index, molecular connectivity topochemical index—an adjacency based topochemical index and eccentric connectivity topochemical index—an adjacency-cum-distance based topochemical index with sodium channel binding activity has been investigated. A dataset comprising 50 hydantoins and related non-hydantoins was selected. The dataset was divided equally into training and test sets. The values of the three topochemical indices for all the compounds present in both the training and test sets were computed using an in-house computer program. The resulting data was analyzed and suitable models were developed after identification of the active ranges in the training set. Subsequently, a biological activity was assigned to each compound involved in the training set using these models, which was then compared with the reported sodium channel binding activity. An accuracy of prediction of the order of >99% was observed using the proposed models. Cross-validation of these models using the test set revealed an exceptionally high accuracy of ∼95%. 3,5-disubstituted-5-phenylhydantoin     

3D–QSAR studies of orvinol analogs as κ-opioid agonists

Orvinols are potent analgesics that target opioid receptors. However, their analgesic mechanism remains unclear and no significant preference for subtype opioid receptor has been achieved. In order to find new orvinols that target the κ-receptor, comparative 3D–QSAR studies were performed on 26 orvinol analogs using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best predictions for the κ-receptor were obtained with the CoMFA standard model (q ²=0.686, r ²=0.947) and CoMSIA model combined steric, electrostatic, hydrophobic, and hydrogen bond donor/acceptor fields (q ²=0.678, r ²=0.914). The models built were further validated by a test set made up of seven compounds, leading to predictive r ² values of 0.672 for CoMFA and 0.593 for CoMSIA. The study could be helpful for designing and prepare new category κ-agonists from orvinols.      

Proton transfer in some periodic molecular systems

The electronic structure of representative hydrogen bonded systems: hydrogen cyanide, imidazole and malonic acid have been studied at the non-empirical level. The role of the dimensionality on the potential barrier for the proton transfer has been examined. It was shown that it depends on the crystal structure and only in some cases like hydrogen cyanide or imidazole the relevant crystals may be considered as one-dimensional. However, for more complicated crystallographic structures, e.g. malonic acid, the evaluated barrier is strongly dependent on the dimensionality taken into account in our calculations. Figure Density of states for the imidazole dimer: MD-EQ- standard proton position (black line), MD-TR- proton position in the middle of the hydrogen bond (red line)     

Insights into catalytic activity of industrial enzyme Co-nitrile hydratase. Docking studies of nitriles and amides

Nitrile hydratase (NHase) is an enzyme containing non-corrin Co³⁺ in the non-standard active site. NHases from Pseudonocardia thermophila JCM 3095 catalyse hydration of nitriles to corresponding amides. The efficiency of the enzyme is 100 times higher for aliphatic nitriles then aromatic ones. In order to understand better this selectivity dockings of a series of aliphatic and aromatic nitriles and related amides into a model protein based on an X-ray structure were performed. Substantial differences in binding modes were observed, showing better conformational freedom of aliphatic compounds. Distinct interactions with postranslationally modified cysteines present in the active site of the enzyme were observed. Modeling shows that water molecule activated by a metal ion may easily directly attack the docked acrylonitrile to transform this molecule into acryloamide. Thus docking studies provide support for one of the reaction mechanisms discussed in the literature. Figure Crystalographic structure of Pseudonocardia thermophila JCM 3095 nitrile hydratase (a) and the non-standard active site (b)