Physicochemical Properties and Elimination of the Activity of Anti-Nutritional Serine Protease Inhibitors from Mulberry Leaves

Mulberry leaf is an excellent protein resource that can be used as feed additive for livestock and poultry. Nevertheless, the use of mulberry leaves in animal diets is limited by its protease inhibitors, tannic acid and other anti-nutritional factors. This study systematically analyzed the type and activity of serine protease inhibitors (SPIs) from the leaves of 34 mulberry varieties, aiming to reveal the physicochemical properties and inactivation mechanism of SPIs. The types and activities of trypsin inhibitors (TIs) and chymotrypsin inhibitors (CIs) exhibited polymorphisms among different mulberry varieties. The highest number of types of inhibitors was detected in Jinshi, with six TIs (TI-1~TI-6) and six CIs (CI-1~CI-6). TIs and CIs exhibited strong thermal and acid–base stability. High-temperature and high-pressure treatment could reduce the activities of TIs and CIs to a certain extent. β-mercaptoethanol treatment could completely abolish TIs and CIs, suggesting that the disulfide bridges were critical for their inhibitory activities. The Maillard reaction could effectively eliminate the inhibitory activities of TI-1~TI-4 and CI-1~CI-4. This study reveals the physicochemical properties and inactivation mechanisms of the anti-nutritional SPIs from mulberry leaves, which is helpful to exploit mulberry-leaf food with low-activity SPIs, promote the development and utilization of mulberry-leaf resources in animal feed and provide reference for mulberry breeding with different functions.     

Unified QSAR & network‐based computational chemistry approach to antimicrobials. II. Multiple distance and triadic census analysis of antiparasitic drugs complex networks

In the previous work, we reported a multitarget Quantitative Structure‐Activity Relationship (mt‐QSAR) model to predict drug activity against different fungal species. This mt‐QSAR allowed us to construct a drug–drug multispecies Complex Network (msCN) to investigate drug–drug similarity (González‐Díaz and Prado‐Prado, J Comput Chem 2008, 29, 656). However, important methodological points remained unclear, such as follows: (1) the accuracy of the methods when applied to other problems; (2) the effect of the distance type used to construct the msCN; (3) how to perform the inverse procedure to study species–species similarity with multidrug resistance CNs (mdrCN); and (4) the implications and necessary steps to perform a substructural Triadic Census Analysis (TCA) of the msCN. To continue the present series with other important problem, we developed here a mt‐QSAR model for more than 700 drugs tested in the literature against different parasites (predicting antiparasitic drugs). The data were processed by Linear Discriminate Analysis (LDA) and the model classifies correctly 93.62% (1160 out of 1239 cases) in training. The model validation was carried out by means of external predicting series; the model classified 573 out of 607, that is, 94.4% of cases. Next, we carried out the first comparative study of the topology of six different drug–drug msCNs based on six different distances such as Euclidean, Chebychev, Manhattan, etc. Furthermore, we compared the selected drug–drug msCN and species–species mdsCN with random networks. We also introduced here the inverse methodology to construct species–species msCN based on a mt‐QSAR model. Last, we reported the first substructural analysis of drug–drug msCN using Triadic Census Analysis (TCA) algorithm. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Quantitative structure-activity relationship studies of boron-containing dipeptide proteasome inhibitors using calculated mathematical descriptors

Topological indices (TIs) and atom pairs (APs) were used to develop quantitative structure-activity relationship (QSAR) models of a set of 58 dipeptide boronic acids which are potent inhibitors of proteasome and have found applications in the treatment of various types of cancers. Of the three linear regression methods used for QSAR development, viz., principal components regression (PCR), partial least square (PLS), and ridge regression (RR), the last method gave the most satisfactory models whereas the remaining two methods yielded poor models. RR results obtained in this paper using TIs and APs are comparable to the CoMFA and CoMSIA results reported in the literature with the same set of compounds.     

MMM-QSAR recognition of ribonucleases without alignment: comparison with an HMM model and isolation from Schizosaccharomyces pombe, prediction, and experimental assay of a new sequence

The study of type III RNases constitutes an important area in molecular biology. It is known that the pac1+ gene encodes a particular RNase III that shares low amino acid similarity with other genes despite having a double-stranded ribonuclease activity. Bioinformatics methods based on sequence alignment may fail when there is a low amino acidic identity percentage between a query sequence and others with similar functions (remote homologues) or a similar sequence is not recorded in the database. Quantitative structure-activity relationships (QSAR) applied to protein sequences may allow an alignment-independent prediction of protein function. These sequences of QSAR-like methods often use 1D sequence numerical parameters as the input to seek sequence-function relationships. However, previous 2D representation of sequences may uncover useful higher-order information. In the work described here we calculated for the first time the spectral moments of a Markov matrix (MMM) associated with a 2D-HP-map of a protein sequence. We used MMMs values to characterize numerically 81 sequences of type III RNases and 133 proteins of a control group. We subsequently developed one MMM-QSAR and one classic hidden Markov model (HMM) based on the same data. The MMM-QSAR showed a discrimination power of RNAses from other proteins of 97.35% without using alignment, which is a result as good as for the known HMM techniques. We also report for the first time the isolation of a new Pac1 protein (DQ647826) from Schizosaccharomyces pombe strain 428-4-1. The MMM-QSAR model predicts the new RNase III with the same accuracy as other classical alignment methods. Experimental assay of this protein confirms the predicted activity. The present results suggest that MMM-QSAR models may be used for protein function annotation avoiding sequence alignment with the same accuracy of classic HMM models.     

具有除草活性的大环内酯类衍生物的定量构效关系

研究了一系列结构新颖的具有除草活性的大环内酯衍生物的定量构效关系(QSAR). 构建的比较分子力场分析(CoMFA)、比较分子近似指数分析(CoMSIA)和全息定量构效关系(HQSAR)分子模型的交叉验证系数r2cv均大于0.5, 非交叉验证系数r2都超过0.8, 表明获取的QSAR模型具有可信的预测能力. 对CoMFA、CoMSIA模型的三维(3D)等势图分析, 发现除了立体场和静电场外, 疏水场和氢键受体场也是影响大环内酯类化合物除草活性的重要因素. 构建的HQSAR模型的原子贡献图提示的结构改造信息与三维QSAR的结果基本一致. 利用CoMFA、CoMSIA模型提供的信息,对目前已合成的活性最高化合物B1-3进行分子结构改造, 预测结果发现部分化合物可能具有更好的除草活性.     

A combined molecular modeling study on gelatinases and their potent inhibitors

Zinc‐dependent matrix metalloproteinase (MMP) family is considered to be an attractive target because of its important role in many physiological and pathological processes. In the present work, a molecular modeling study combining protein‐, ligand‐ and complex‐based computational methods was performed to analyze a new series of β‐N‐biaryl ether sulfonamide hydroxamates as potent inhibitors of gelatinase A (MMP‐2) and gelatinase B (MMP‐9). Firstly, the similarities and differences between the binding sites of MMP‐2 and MMP‐9 were analyzed through sequence alignment and structural superimposition. Secondly, in order to extract structural features influencing the activities of these inhibitors, quantitative structure‐activity relationship (QSAR) models using genetic algorithm‐multiple linear regression (GA‐MLR), comparative molecular field (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were developed. The proposed QSAR models could give good predictive ability for the studied inhibitors. Thirdly, docking study was employed to further explore the binding mode between the ligand and protein. The results from all the above analyses could provide the information about the similarities and differences of the binding mode between the MMP‐2, MMP‐9 and their potent inhibitors. The obtained results can provide very useful information for the design of new potential inhibitors. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

GTI-space: the space of generalized topological indices

A new extension of the generalized topological indices (GTI) approach is carried out to represent “simple” and “composite” topological indices (TIs) in an unified way. This approach defines a GTI-space from which both simple and composite TIs represent particular subspaces. Accordingly, simple TIs such as Wiener, Balaban, Zagreb, Harary and Randić connectivity indices are expressed by means of the same GTI representation introduced for composite TIs such as hyper-Wiener, molecular topological index (MTI), Gutman index and reverse MTI. Using GTI-space approach we easily identify mathematical relations between some composite and simple indices, such as the relationship between hyper-Wiener and Wiener index and the relation between MTI and first Zagreb index. The relation of the GTI-space with the sub-structural cluster expansion of property/activity is also analysed and some routes for the applications of this approach to QSPR/QSAR are also given.     

Homology model of the rainbow trout estrogen receptor (rtERalpha) and docking of endocrine disrupting chemicals (EDCs)

A model for rainbow trout (Oncorhynchus mykiss) estrogen receptor (rtERa) was built by homology with the human estrogen receptor (hERa). A high level of sequence conservation between the two receptors was found with 64% and 80% of identity and similarity, respectively. Selected endocrine disrupting chemicals were docked into the ligand binding domain (LBD) of rtERa and the corresponding free binding energies Delta(DeltaG(bind)) values were calculated. A Quantitative Structure-Activity Relationship (QSAR) model between the relative binding affinity data and the Delta(DeltaG(bind)) values was derived in order to predict which further organic pollutants are likely to bind to rtERa.     

芳基磺酰胺类丙酮酸激酶M2激动剂的基于多复合物的药效团模型和定量构效关系

丙酮酸激酶M2(PKM2)是肿瘤治疗中最具发展潜力的靶点之一. 本文以一系列丙酮酸激酶M2-激动剂复合物的晶体结构为基础, 采用基于多复合物的药效团(MCBP)方法产生了PKM2 的药效团模型. 并使用该药效团模型产生了62个芳基磺酰胺类PKM2激动剂的活性构象和分子叠合, 通过三维定量构效关系(3D-QSAR)方法研究了该类PKM2激动剂与PKM2蛋白的相互作用, 并建立了相关预测模型. 比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)模型的交互验证相关系数q²分别为0.545 和0.653, 非交互验证相关系数r²分别为0.966和0.987. 本研究为进一步结构优化、设计和合成新型的PKM2激动剂提供了理论依据.     

Optimal neighbor selection in molecular similarity: comparison of arbitrary versus tailored prediction spaces

Three classes of arbitrary quantitative molecular similarity analysis (QMSA) methods have been computed using atom pairs (APs), topological indices (TIs), and principal components (PCs) derived from topological indices. Tailored QMSA models have been developed from TIs selected through ridge regression. K-nearest neighbor (kNN) based estimation has been applied to all of the methods to estimate normal vapor pressure (p(vap)) and water solubility (sol) for a set of 194 chemicals. Results show that the tailored QMSA methods are superior to arbitrary similarity methods in estimating both of these properties for the given set of chemicals.     

Design,synthesis, pharmacological evaluation and in silico ADMET prediction of novel substituted benzimidazole derivatives as angiotensin II–AT1 receptor antagonists based on predictive 3D QSAR models

In this study we designed novel substituted benzimidazole derivatives and predicted their absorption, distribution, metabolism, excretion and toxicity (ADMET) properties, based on a predictive 3D QSAR study on 132 substituted benzimidazoles as AngII–AT₁ receptor antagonists. The two best predicted compounds were synthesized and evaluated for AngII–AT₁ receptor antagonism. Three different alignment tools for comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used. The best 3D QSAR models were obtained using the rigid body (Distill) alignment method. CoMFA and CoMSIA models were found to be statistically significant with leave-one-out correlation coefficients (q²) of 0.630 and 0.623, respectively, cross-validated coefficients (r²_cv) of 0.651 and 0.630, respectively, and conventional coefficients of determination (r²) of 0.848 and 0.843, respectively. 3D QSAR models were validated using a test set of 24 compounds, giving satisfactory predicted results (r²_pred) of 0.727 and 0.689 for the CoMFA and CoMSIA models, respectively. We have identified some key features in substituted benzimidazole derivatives, such as lipophilicity and H-bonding at the 2- and 5-positions of the benzimidazole nucleus, respectively, for AT₁ receptor antagonistic activity. We designed 20 novel substituted benzimidazole derivatives and predicted their activity. In silico ADMET properties were also predicted for these designed molecules. Finally, the compounds with best predicted activity were synthesized and evaluated for in vitro angiotensin II–AT₁ receptor antagonism.     

Design of topological indices: computer-oriented approach

A novel method is suggested for constructing topological indices (TIs) of molecular graphs which models human logic. This method is described in terms of a block scheme, consisting of the mutually connected elementary blocks. In each block the simple transformations of a molecular graph are fulfilled. A variant of the transformation is selected from the list of possible variants. Every TI is obtained as a result of the sequential execution of a number of operations, corresponding to some ‘walk’ on the block scheme. This walk can be selected both randomly and by the investigator. The suggested method can serve as a basis for the development of the respective computer program which may be used for the automatic construction of any number of TIs of differing nature. By this process one can also obtain the TIs that are unlikely to be constructed manually, due to their complexity. The set of obtained TIs may be used for building the structure–property models. In the case of an unsatisfactory result the obtained set of TIs may be changed using the described generator of TIs. A number of examples of application of the suggested approach for the building QSAR/QSPR models is given.     

QSAR models of cytochrome P450 enzyme 1A2 inhibitors using CoMFA,CoMSIA and HQSAR

Quantitative structure–activity relationship (QSAR) studies were conducted on an in-house database of cytochrome P450 enzyme 1A2 inhibitors using the comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA) and hologram QSAR (HQSAR) approaches. The database consisted of 36 active molecules featuring varied core structures. The model based on the naphthalene substructure alignment incorporating 19 molecules yielded the best model with a CoMFA cross validation value q² of 0.667 and a Pearson correlation coefficient r² of 0.976; a CoMSIA q² value of 0.616 and r² value of 0.985; and a HQSAR q² value of 0.652 and r² value of 0.917. A second model incorporating 34 molecules aligned using the benzene substructure yielded an acceptable CoMFA model with q² value of 0.5 and r² value of 0.991. Depending on the core structure of the molecule under consideration, new CYP1A2 inhibitors will be designed based on the results from these models.