Toward prediction of class II mouse major histocompatibility complex peptide binding affinity: in silico bioinformatic evaluation using partial least squares, a robust multivariate statistical technique

The accurate identification of T-cell epitopes remains a principal goal of bioinformatics within immunology. As the immunogenicity of peptide epitopes is dependent on their binding to major histocompatibility complex (MHC) molecules, the prediction of binding affinity is a prerequisite to the reliable prediction of epitopes. The iterative self-consistent (ISC) partial-least-squares (PLS)-based additive method is a recently developed bioinformatic approach for predicting class II peptide-MHC binding affinity. The ISC-PLS method overcomes many of the conceptual difficulties inherent in the prediction of class II peptide-MHC affinity, such as the binding of a mixed population of peptide lengths due to the open-ended class II binding site. The method has applications in both the accurate prediction of class II epitopes and the manipulation of affinity for heteroclitic and competitor peptides. The method is applied here to six class II mouse alleles (I-Ab, I-Ad, I-Ak, I-As, I-Ed, and I-Ek) and included peptides up to 25 amino acids in length. A series of regression equations highlighting the quantitative contributions of individual amino acids at each peptide position was established. The initial model for each allele exhibited only moderate predictivity. Once the set of selected peptide subsequences had converged, the final models exhibited a satisfactory predictive power. Convergence was reached between the 4th and 17th iterations, and the leave-one-out cross-validation statistical terms--q2, SEP, and NC--ranged between 0.732 and 0.925, 0.418 and 0.816, and 1 and 6, respectively. The non-cross-validated statistical terms r2 and SEE ranged between 0.98 and 0.995 and 0.089 and 0.180, respectively. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made freely available online (http://www.jenner.ac.uk/MHCPred).     

HLA-A*0201限制性CTL表位肽的三维定量构效关系的研究

运用比较分子力场(CoMFA)和比较分子相似性指数分析(CoMFA)方法研究了50个HLA-A^*0201限制性CTL表位九肽结构与亲和性间的关系,另外15个表位九肽作为预测集用于检验模型的预测能力.结果表明采用CoMSIA得到的构效关系模型(q^2=0.628,r^2=0.997,F=840.419)要明显优于采用CoMFA得到的构效关系模型.在CoMSIA计算中,当引入疏水场时,三维构效关系模型得到明显改善,通过该三维构效关系模型,可较精确地估算预测集中15个CTL表位肽与HLA-A^*0201间的亲和力(r^2pred=0.743).通过分析分子场等值面图在空间的分布,可以观察到表位肽分子周围的立体及疏水特征对表位肽与HLA-A^*0201间结合亲和力的影响,从而为进一步对CTL表位肽进行结构改造并基于此进行治疗性疫苗分子设计提供理论基础.     

Performance of (consensus) kNN QSAR for predicting estrogenic activity in a large diverse set of organic compounds

A novel method (in the context of quantitative structure-activity relationship (QSAR)) based on the k nearest neighbour (kNN) principle, has recently been introduced for the derivation of predictive structure-activity relationships. Its performance has been tested for estimating the estrogen binding affinity of a diverse set of 142 organic molecules. Highly predictive models have been obtained. Moreover, it has been demonstrated that consensus-type kNN QSAR models, derived from the arithmetic mean of individual QSAR models were statistically robust and provided more accurate predictions than the great majority of the individual QSAR models. Finally, the consensus QSAR method was tested with 3D QSAR and log P data from a widely used steroid benchmark data set.     

A 3D QSAR study on a set of dopamine D4 receptor antagonists

The molecular alignments obtained from a previously reported pharmacophore model have been employed in a three-dimensional quantitative structure-activity relationship (3D QSAR) study, to obtain a more detailed insight into the structure-activity relationships for D(2) and D(4) receptor antagonists. The frequently applied CoMFA method and the related CoMSIA method were used. Statistically significant models have been derived with these two methods, based on a set of 32 structurally diverse D(2) and D(4) receptor antagonists. The CoMSIA and the CoMFA methods produced equally good models expressed in terms of q(2) values. The predictive power of the derived models were demonstrated to be high. Graphical interpretation of the results, provided by the CoMSIA method, brings to light important structural features of the compounds related to either low- or high-affinity D(2) or D(4) antagonism. The results of the 3D QSAR studies indicate that bulky N-substituents decrease D(2) binding, whereas D(4) binding is enhanced. Electrostatically favorable and unfavorable regions exclusive to D(2) receptor binding were identified. Likewise, certain hydrogen-bond acceptors can be used to lower D(2) affinity. These observations may be exploited for the design of novel dopamine D(4) selective antagonists.     

Performance of (consensus) kNN QSAR for predicting estrogenic activity in a large diverse set of organic compounds

A novel method (in the context of quantitative structure–activity relationship (QSAR)) based on the k nearest neighbour (kNN) principle, has recently been introduced for the derivation of predictive structure–activity relationships. Its performance has been tested for estimating the estrogen binding affinity of a diverse set of 142 organic molecules. Highly predictive models have been obtained. Moreover, it has been demonstrated that consensus-type kNN QSAR models, derived from the arithmetic mean of individual QSAR models were statistically robust and provided more accurate predictions than the great majority of the individual QSAR models. Finally, the consensus QSAR method was tested with 3D QSAR and log?P data from a widely used steroid benchmark data set.     

Combinatorial QSAR modeling of chemical toxicants tested against Tetrahymena pyriformis

Selecting most rigorous quantitative structure-activity relationship (QSAR) approaches is of great importance in the development of robust and predictive models of chemical toxicity. To address this issue in a systematic way, we have formed an international virtual collaboratory consisting of six independent groups with shared interests in computational chemical toxicology. We have compiled an aqueous toxicity data set containing 983 unique compounds tested in the same laboratory over a decade against Tetrahymena pyriformis. A modeling set including 644 compounds was selected randomly from the original set and distributed to all groups that used their own QSAR tools for model development. The remaining 339 compounds in the original set (external set I) as well as 110 additional compounds (external set II) published recently by the same laboratory (after this computational study was already in progress) were used as two independent validation sets to assess the external predictive power of individual models. In total, our virtual collaboratory has developed 15 different types of QSAR models of aquatic toxicity for the training set. The internal prediction accuracy for the modeling set ranged from 0.76 to 0.93 as measured by the leave-one-out cross-validation correlation coefficient ( Q abs2). The prediction accuracy for the external validation sets I and II ranged from 0.71 to 0.85 (linear regression coefficient R absI2) and from 0.38 to 0.83 (linear regression coefficient R absII2), respectively. The use of an applicability domain threshold implemented in most models generally improved the external prediction accuracy but at the same time led to a decrease in chemical space coverage. Finally, several consensus models were developed by averaging the predicted aquatic toxicity for every compound using all 15 models, with or without taking into account their respective applicability domains. We find that consensus models afford higher prediction accuracy for the external validation data sets with the highest space coverage as compared to individual constituent models. Our studies prove the power of a collaborative and consensual approach to QSAR model development. The best validated models of aquatic toxicity developed by our collaboratory (both individual and consensus) can be used as reliable computational predictors of aquatic toxicity and are available from any of the participating laboratories.     

Lessons learned in induced fit docking and metadynamics in the Drug Design Data Resource Grand Challenge 2

Two of the major ongoing challenges in computational drug discovery are predicting the binding pose and affinity of a compound to a protein. The Drug Design Data Resource Grand Challenge 2 was developed to address these problems and to drive development of new methods. The challenge provided the 2D structures of compounds for which the organizers help blinded data in the form of 35 X-ray crystal structures and 102 binding affinity measurements and challenged participants to predict the binding pose and affinity of the compounds. We tested a number of pose prediction methods as part of the challenge; we found that docking methods that incorporate protein flexibility (Induced Fit Docking) outperformed methods that treated the protein as rigid. We also found that using binding pose metadynamics, a molecular dynamics based method, to score docked poses provided the best predictions of our methods with an average RMSD of 2.01 Å. We tested both structure-based (e.g. docking) and ligand-based methods (e.g. QSAR) in the affinity prediction portion of the competition. We found that our structure-based methods based on docking with Smina (Spearman ρ?=?0.614), performed slightly better than our ligand-based methods (ρ?=?0.543), and had equivalent performance with the other top methods in the competition. Despite the overall good performance of our methods in comparison to other participants in the challenge, there exists significant room for improvement especially in cases such as these where protein flexibility plays such a large role.     

CoMFA, CoMSIA, and molecular hologram QSAR studies of novel neuronal nAChRs ligands-open ring analogues of 3-pyridyl ether

3-Pyridyl ethers are excellent nAChRs ligands, which show high subtype selectivity and binding affinity to alpha4beta2 nAChR. Although the quantitative structure-activity relationship (QSAR) of nAChRs ligands has been widely investigated using various classes of compounds, the open ring analogues of 3-pyridyl ethers have been less involved in these studies due to the greater flexibility of this kind of molecule. In this study, two three-dimensional QSAR techniques and one two-dimensional QSAR technique were used to correlate the molecular structure with the biological activity of 64 analogues of 3-pyridyl ethers. Three different QSAR models were established. Their performances in the QSAR studies of open ring analogues of 3-pyridyl ethers were evaluated by the statistical values in the corresponding models. All models exhibited satisfactory predictive power. Of these models, the HQSAR behaved optimally in terms of the statistical values with q2=0.845, r2=0.897. Finally, graphic interpretation of three different models provided coincident information about the interaction of the ligand-receptor complex and supplied useful guidelines for the synthesis of novel, potent ligands.     

Triplex selective 2-(2-naphthyl)quinoline compounds: origins of affinity and new design principles

A novel competition dialysis assay was used to investigate the structural selectivity of a series of substituted 2-(2-naphthyl)quinoline compounds designed to target triplex DNA. The interaction of 14 compounds with 13 different nucleic acid sequences and structures was studied. A striking selectivity for the triplex structure poly dA:[poly dT](2) was found for the majority of compounds studied. Quantitative analysis of the competition dialysis binding data using newly developed metrics revealed that these compounds are among the most selective triplex-binding agents synthesized to date. A quantitative structure-affinity relationship (QSAR) was derived using triplex binding data for all 14 compounds used in these studies. The QSAR revealed that the primary favorable determinant of triplex binding free energy is the solvent accessible surface area. Triplex binding affinity is negatively correlated with compound electron affinity and the number of hydrogen bond donors. The QSAR provides guidelines for the design of improved triplex-binding agents.