3D QSAR studies,pharmacophore modeling,and virtual screening of diarylpyrazole–benzenesulfonamide derivatives as a template to obtain new inhibitors,using human carbonic anhydrase II as a model protein

A 3D-QSAR modeling was performed on a series of diarylpyrazole-benzenesulfonamide derivatives acting as inhibitors of the metalloenzyme carbonic anhydrase (CA, EC 4.2.1.1). The compounds were collected from two datasets with the same scaffold, and utilized as a template for a new pharmacophore model to screen the ZINC database of commercially available derivatives. The datasets were divided into training, test, and validation sets. As the first step, comparative molecular field analysis (CoMFA), CoMFA region focusing and comparative molecular similarity indices analysis (CoMSIA) in parallel with docking studies were applied to a set of 41 human (h) CA II inhibitors. The validity and the prediction capacity of the resulting models were evaluated by leave-one-out (LOO) cross-validation approach. The reliability of the model for the prediction of possibly new CA inhibitors was also tested.     

3D-QSAR and molecular docking studies of 2-pyrimidinecarbonitrile derivatives as inhibitors against falcipain-3

The three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on a series of falcipain-3 inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. A training set containing 42 molecules served to establish the QSAR models. The optimum CoMFA and CoMSIA models obtained for the training set were statistically significant with cross-validated correlation coefficients r(cv)(2) (q(2)) of 0.549 and 0.608, and conventional correlation coefficients (r(2)) of 0.976 and 0.932, respectively. An independent test set of 12 molecules validated the external predictive power of both models with predicted correlation coefficients (r(pred)(2)) for CoMFA and CoMSIA as 0.697 and 0.509, respectively. The docking of inhibitors into falcipain-3 active site using GOLD software revealed the vital interactions and binding conformation of the inhibitors. The CoMFA and CoMSIA field contour maps agree well with the structural characteristics of the binding pocket of falcipain-3 active site, which suggests that the information rendered by 3D-QSAR models and the docking interactions can provide guidelines for the development of improved falcipain-3 inhibitors.     

A combined approach of docking and 3D QSAR study of beta-ketoacyl-acyl carrier protein synthase III (FabH) inhibitors

The enzyme FabH catalyzes the initial step of fatty acid biosynthesis via a type II fatty acid synthase. The pivotal role of this essential enzyme combined with its unique structural features and ubiquitous occurrence in bacteria has made it an attractive new target for the development of antibacterial and antiparasitic compounds. Three-dimensional quantitative structure-activity relationship (3D QSAR) studies such as comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) and docking simulations were conducted on a series of potent benzoylaminobenzoic acids. Docking studies were employed to position the inhibitors into the FabH active site to determine the probable binding conformation. A reasonable correlation between the predicated binding free energy and the inhibitory activity was found. CoMFA and CoMSIA were performed based on the docking conformations, giving q(2) of 0.637 and 0.697 for CoMFA and CoMSIA models, respectively. The predictive ability of the models was validated using a set of compounds that were not included in the training set and progressive scrambling test. Mapping the 3D QSAR models to the active site of FabH related that some important amino acid residues are responsible for protein-inhibitor interaction. These results should be applicable to the prediction of the activities of new FabH inhibitors, as well as providing structural understanding.     

Molecular simulation of a series of benzothiazole PI3Kα inhibitors: probing the relationship between structural features, anti-tumor potency and selectivity

The phosphatidylinositol 3-kinase α (PI3Kα) was genetically validated as a promising therapeutic target for developing novel anticancer drugs. In order to explore the structure-activity correlation of benzothiazole series as inhibitors of PI3Kα, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) were performed on 61 promising molecules to build 3D-QSAR models based on both the ligand- and receptor-based methods. The best CoMFA and CoMSIA models had a cross-validated coefficient r(cv)(2) of 0.618 and 0.621, predicted correlation coefficient r(pred) (2) of 0.812 and 0.83, respectively, proving their high correlative and predictive abilities on both the training and test sets. In addition, docking analysis and molecular dynamics simulation (MD) were also applied to elucidate the probable binding modes of these inhibitors at the ATP binding pocket. Based on the contour maps and MD results, some key structural factors responsible for the activity of this series of compounds were revealed as follows: (1) Ring-A has a strong preference for bulky hydrophobic or aromatic groups; (2) Electron-withdrawing groups at the para position of ring-B and hydrophilic substituents in ring-B region may benefit the potency; (3) A polar substituent like -NHSO(2)- between ring-A and ring-B can enhance the activity of the drug by providing hydrogen bonding interaction with the protein target. The satisfactory results obtained from this work strongly suggest that the developed 3D-QSAR models and the obtained PI3Kα inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and be helpful in future PI3Kα inhibitor design.     

A novel 3D-QSAR comparative molecular field analysis (CoMFA) model of imidazole and quinazolinone functionalized p38 MAP kinase inhibitors

In this study we describe a new comparative molecular field analysis (CoMFA) model of dihydroquinazolinone and tetrasubstituted imidazole compounds with p38 MAPK inhibitory activity. A series of 51 (a training set of 40 and a test set of 11) dihydroquinazolinone [Bioorg. Med. Chem. Lett. 2003, 13, 277.] and tetrasubstituted imidazole [J. Med. Chem. 1999, 42, 2180.] derivatives known as p38 mitogen-activated protein kinase (p38 MAPK) selective inhibitors was studied by quantitative structure-activity relationship (3D-QSAR) analysis using comparative molecular field analysis. The 3D-QSAR models were generated and evaluated by a scheme that combines a genetic algorithm (GA) optimization with partial least squares (PLS) regression and by crossvalidation using the leave-one-out technique. The model was able to efficiently predict the activities of the compounds of the test set, suggesting that it can be used for the planning of new p38 MAPK inhibitor candidates useful to treat chronic inflammatory states.     

Structural and dynamics analysis of matrix metalloproteinases MMP-2 complexed with chemically modified tetracyclines (CMTs)

Matrix metalloproteinases (MMPs) play a critical role in physiological processes and pathological conditions such tumor invasion and metastasis. In recent years, a number of MMP inhibitors have been proposed, including the chemically modified tetracyclines (CMTs), which have been evaluated in preclinical cancer models showing promising results. This work provides insights into the structure and dynamics of the MMP-2 catalytic domain complexed with seven CMT (CMT-n), based on the analysis of molecular dynamics trajectories in solution. The comparative analysis of various relevant molecular aspects of the different complexes of MMP-2 and CMT-n derivatives was performed aiming to elucidate the effect of ligands on the enzyme structure. These include the radial distribution function of the water molecules around the catalytic zinc, the solvent accessible surface area for the inhibitors and the root-mean-square fluctuation for all amino acid residues. The results help to understand the differences in the binding modes of related compounds and, therefore, add to further design of novel tetracycline-based inhibitors for MMP enzymes.     

Comparative molecular field analysis (CoMFA) models of phenylethanolamine N-methyltransferase (PNMT) and the alpha2-adrenoceptor: the development of new, highly selective inhibitors of PNMT

As a guide to the development of new and more selective inhibitors of phenylethanolamine N-methyltransferase (PNMT) vs the alpha2-adrenoceptor, we have performed a comparative molecular field analysis (CoMFA) on a series of 80 benzylamine analogues. Using the models obtained, we have proposed a series of 3-trifluoromethyl-1,2,3,4-tetrahydroisoquinolines and predicted the activity of other analogues.     

Two- and three-dimensional quantitative structure-activity relationships for a series of purine nucleoside phosphorylase inhibitors

Comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis, and hologram quantitative structure-activity relationship (HQSAR) studies were conducted on a series of 52 training set inhibitors of calf spleen purine nucleoside phosphorylase (PNP). Significant cross-validated correlation coefficients (CoMFA, q(2)=0.68; CoMSIA, q(2)=0.66; and HQSAR, q(2)=0.70) were obtained, indicating the potential of the models for untested compounds. The models were then used to predict the inhibitory potency of 16 test set compounds that were not included in the training set, and the predicted values were in good agreement with the experimental results. The final QSAR models along with the information gathered from 3D contour and 2D contribution maps should be useful for the design of novel inhibitors of PNP having improved potency.     

Using multidimensional patterns of amino acid attributes for QSAR analysis of peptides

On the basis of exploratory factor analysis, six multidimensional patterns of 516 amino acid attributes, namely, factor analysis scales of generalized amino acid information (FASGAI) involving hydrophobicity, alpha and turn propensities, bulky properties, compositional characteristics, local flexibility and electronic properties, are proposed to represent structures of 48 bitter-tasting dipeptides and 58 angiotensin-converting enzyme inhibitors. Characteristic parameters related to bioactivities of the peptides studied are selected by genetic algorithm, and quantitative structure–activity relationship (QSAR) models are constructed by partial least square (PLS). Our results by a leave-one-out cross validation are compared with the previously known structure representation method and are shown to give slightly superior or comparative performance. Further, two data sets are divided into training sets and test sets to validate the characterization repertoire of FASGAI. Performance of the PLS models developed by training samples by a leave-one-out cross validation and external validation for test samples are satisfying. These results demonstrate that FASGAI is an effective representation technique of peptide structures, and that FASGAI vectors have many preponderant characteristics such as straightforward physicochemical information, high characterization competence and easy manipulation. They can be further applied to investigate the relationship between structures and functions of various peptides, even proteins.     

Three-dimensional quantitative structure activity relationships (3-D-QSAR) of antihyperglycemic agents.

A three-dimensional quantitative structure activity relationship study (3-D-QSAR) was performed on a set of thiazolidinedione antihyperglycemic agents using the comparative molecular field analysis (CoMFA) method. The CoMFA models were derived from a training set of 53 compounds. Fifteen compounds, which were not used in model generation were used to validate the CoMFA models. All the compounds were superimposed to the template structure by atom-based and shape-based strategies. The SYBYL QSAR rigid body field fit was also used for aligning the ligands. A total of twelve different alignments were generated. The resulting models exhibited a good cross-validated r2cv values (0.624-0.764) and the conventional r2 values (0.689-0.921). A more robust cross-validation test using cross-validation by 2 groups (leave half out method) was performed 100 times to ascertain the predictiveness of the CoMFA models. The mean of r2cv values from 100 runs ranged from 0.611-0.690. Few models exhibited good external predictivity. These models were then used to define a hypothetical receptor model for antihyperglycemic agents.     

3D-QSAR studies on thieno[3,2-d]pyrimidines as phosphodiesterase IV inhibitors

Cyclic nucleotide phosphodiesterase IV (PDE IV) inhibitors find utility in asthma and Chronic Obstructive Pulmonary Disease (COPD) therapy. A series of 29 thieno[3,2-d]pyrimidines with affinity for PDE IV was subjected to three dimensional quantitative structure activity relationship (3D-QSAR) studies using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Both CoMFA and CoMSIA provided statistically valid models with good correlative and predictive power. The incorporation of hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields showed insignificant improvement in CoMSIA model. The 3D-QSAR models provide information for predicting the affinity of related compounds and designing more potent inhibitors.     

Pharmacophore-based 3D-QSAR as a predictive method for the QSAR analysis on a series of potent and selective inhibitors for three kinases of RTK family

For targets belonging to the same family of receptors, inhibitors often act at more than one biological target and produce a synergistic effect. Separate pharmacophore-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were developed from our data-set for the kinase insert domain-containing receptor (KDR), cKit and Flt1 inhibitors. These models showed excellent internal predictability and consistency; validation using test-set compounds yielded a good predictive power for the pIC₅₀-value. The field contour maps (CoMFA and CoMSIA) corresponding to the KDR, cKit and Flt1 kinase subtypes reflected the characteristic similarities and differences between these types. These maps provided valuable information to facilitate structural modifications of the inhibitor to increase selectivity of the KDR over cKit and Flt1.     

Matrix metalloproteinases and their endogenous regulators in squamous cervical carcinoma (A review of own results)

Own results of long-term studies of expression of matrix metalloproteinases (MMPs) and their endogenous regulators examined in fibroblasts transformed by oncogene E7 HPV16 (TF), immortalized fibroblasts (IF), cell lines associated with HPV16 and HPV18, and tumor tissue samples from patients with squamous cervical carcinoma (SCC) associated with HPV16 have been summarized. Transfection of fibroblasts with the E7 HPV16 oncogen was accompanied by induction of collagenase (MMP-1, MMP-14) and gelatinase (MMP-9) gene expression and the increase in catalytic activity of these MMP, while gelatinase MMP-2 expression remained unchanged. MMP expression correlated with the tumorigenic of transformed clones. Expression of MMP-9 was found only in TF. In TF expression mRNA TIMP-1 decreased, while expression of the genatinase inhibitor, TIMP-2, increased. Collagenase activity and expression of the MMP-14 (collagenase) mRNA increased, while gelatinase activity remained unchanged. The destructive potential of TF is associated with induction of collagenases, gelatinase MMP-9 and decreased levels of MMP inhibitors. MMP-9 may serve as a TF marker. Invasive potential of cell lines associated with HPV18 (HeLa and S4-1) was more pronounced than that of cell lines associated with HPV16 (SiHa and Caski). In most cell lines mRNA levels of collagenases MMP-1 and MMP-14 and the activator (uPA) increased, while gelatinase MMP-2 mRNA and tissue inhibitors mRNAs changed insignificantly. MMP-2 activity significantly increased in Caski and HeLa cell lines, while MMP-9 expression in these cell lines was not detected. The comparative study of expression MMP of and their endogenous regulators performed using SCC tumor samples associated with HPV16 has shown that the invasive and metastatic potentials of tumor tissue in SCC is obviously associated with increased expression of collagenases MMP-1, MMP-14 and gelatinase MMP-9, as well as decreased expression of inhibitors (TIMP-1 and TIMP-2), and to a lesser extent with increased expression of MMP-2. MMP-1 and MMP-9 can serve as markers of invasive and metastatic potential of the SCC tumor. The morphologically normal tissue adjacent to the tumor tissue is characterized by significant expression of MMP-1, MMP-2, and MMP-9. This also contributes to the increased destructive potential of the tumor.     

Molecular docking and 3D-QSAR studies on inhibitors of DNA damage signaling enzyme human PARP-1

Poly (ADP-ribose) polymerase-1 (PARP-1) operates in a DNA damage signaling network. Molecular docking and three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on human PARP-1 inhibitors. Docked conformation obtained for each molecule was used as such for 3D-QSAR analysis. Molecules were divided into a training set and a test set randomly in four different ways, partial least square analysis was performed to obtain QSAR models using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Derived models showed good statistical reliability that is evident from their r2, q2(loo) and r2(pred) values. To obtain a consensus for predictive ability from all the models, average regression coefficient r2(avg) was calculated. CoMFA and CoMSIA models showed a value of 0.930 and 0.936, respectively. Information obtained from the best 3D-QSAR model was applied for optimization of lead molecule and design of novel potential inhibitors.     

Structure-based and ligand-based drug design for microsomal prostaglandin E synthase-1 inhibitors

Microsomal prostaglandin E synthase-1 (mPGES-1) has been regarded as an attractive drug for inflammation-related diseases. In search of new mPGES-1 inhibitors, we performed virtual screening using our traditional Chinese medicine and natural products database (http://tcm.cmu.edu.tw/) and constructed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) using a training set of 30 experimentally tested mPGES-1 inhibitors. The CoMFA and CoMSIA models derived were statistically significant with cross-validated coefficient values of 0.808 for CoMFA and 0.829 for CoMSIA and non-cross-validated coefficient values of 0.829 for CoMFA and 0.980 for CoMSIA. Docking and de novo evolution design gave three top derivatives, 2-O-caffeoyl tartaric acid-Evo_2, glucogallin-Evo_1 and 3-O-feruloylquinic acid-Evo_7 that have higher binding affinities than the control, glutathione. These three derivatives have interactions with Arg70, Arg73, Arg110, Arg126 and Arg38, which all are mPGES-1 key active site residues. In addition, these derivatives fit well into the CoMFA and CoMSIA models, with hydrophobic, hydrophilic and electropositive substructures mapped onto corresponding contour plots. Hence, we suggest that these three de novo compounds could be a starting basis for new mPGES-1 inhibitors.     

The isolation, characterization, and molecular cloning of a 75-kDa gelatinase B-like enzyme, a member of the matrix metalloproteinase (MMP) family. An avian enzyme that is MMP-9-like in its cell expression pattern but diverges from mammalian gelatinase B in sequence and biochemical properties

We have isolated a novel 75-kDa gelatinase from a chicken macrophage cell line, HD11. Biochemical and immunological characterization of the purified enzyme demonstrated that it is distinct from the chicken 72-kDa gelatinase A (MMP-2). The enzyme is capable of specific gelatin binding and rapid gelatin cleavage. Incubation with an organomercurial compound (p-aminophenylmercuric acetate) induces proteolytic processing and activation of this enzyme, and the resultant gelatinolytic activity is sensitive to both zinc chelators and tissue inhibitors of metalloproteinases. A full-length cDNA for the enzyme has been cloned, and sequence analysis demonstrated that the enzyme possesses the characteristic multidomain structure of an MMP gelatinase including a cysteine switch prodomain, three fibronectin type II repeats, a catalytic zinc binding region, and a hemopexin-like domain. The 75-kDa gelatinase is produced by phorbol ester-treated chicken bone marrow cells, monocytes, and polymorphonuclear leukocytes, cell types that charac- teristically produce the 92-kDa mammalian gelatinase B (MMP-9). The absence of a 90-110-kDa gelatinase in these cell types indicates that the 75-kDa gelatinase is likely the avian counterpart of gelatinase B. However, the protein is only 59% identical to human gelatinase B, whereas all previously cloned chicken MMP homologues are 75-90% identical to their human counterparts. In addition, the new 75-kDa chicken gelatinase lacks the type V collagen domain that is found in all mammalian gelatinase Bs. Furthermore, the secreted enzyme appears structurally distinct from known gelatinase Bs and the activated enzyme can cleave fibronectin, which is not a substrate for mammalian gelatinase B. Thus the results of this study indicate that a second MMP gelatinase exists in chickens, and although it is MMP-9/gelatinase B-like in its overall domain structure and expression pattern, it appears to be biochemically divergent from mammalian gelatinase B.     

Molecular design of new aggrecanases-2 inhibitors

Aggrecanases-2 is a very important potential drug target for the treatment of osteoarthritis. In this study, a series of known aggrecanases-2 inhibitors was analyzed by the technologies of three-dimensional quantitative structure–activity relationships (3D-QSAR) and molecular docking. Two 3D-QSAR models, which based on comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods, were established. Molecular docking was employed to explore the details of the interaction between inhibitors and aggrecanases-2 protein. According to the analyses for these models, several new potential inhibitors with higher activity predicted were designed, and were supported by the simulation of molecular docking. This work propose the fast and effective approach to design and prediction for new potential inhibitors, and the study of the interaction mechanism provide a better understanding for the inhibitors binding into the target protein, which will be useful for the structure-based drug design and modifications.