Novel amino acids indices based on quantum topological molecular similarity and their application to QSAR study of peptides

Topological Distance-Based Electron Interaction Tensor to Apply a Convolutional Neural Network on Drug-like Compounds

Deep learning (DL) models in quantitative structure-activity relationship fed the molecular structure directly to the network without using human-designed descriptors by representing molecule as a graph or string (e.g., SMILES code). However, these two representations were oversimplification of real molecules to reflect chemical properties of molecular structures. Given that the choice of molecular representation determines the architecture of the DL model to apply, a novel way of molecular representation can open a way to apply diverse DL networks developed and used in other fields. A topological distance-based electron interaction (TDEi) tensor has been developed in this study inspired by the quantum mechanical model of the molecule, which defines a molecule with electrons and protons. In the TDEi tensor, the atomic orbital (AO) of each atom is represented by an electron configuration (EC) vector, which is a bit string based on the presence and absence of electrons in each AO according to spin indicated by positive and negative signs. Interactions between EC vectors were calculated based on the topological distance between atoms in a molecule. As a molecular structure was translated into 3D array, CNN models (modified VGGNet) were applied using a TDEi tensor to predict four physicochemical properties of drug-like compound datasets: MP (275,131), Lipop (4193), Esol (1127), and Freesolv (639). Models achieved good prediction accuracy. PCA showed that a stronger correlation was observed between the extracted features and the target endpoint as features were extracted from the deeper layer.

In Silico Rational Design and Virtual Screening of Bioactive Peptides Based on QSAR Modeling

Predicting the bioactivity of peptides is an important challenge in drug development and peptide research. In this study, numerical descriptive vectors (NDVs) for peptide sequences were calculated based on the physicochemical properties of amino acids (AAs) and principal component analysis (PCA). The resulted NDV had the same length as the peptide sequence, so that each entry of NDV corresponded to one AA in the sequence. They were then applied to quantitative structure-activity relationship (QSAR) analysis of angiotensin-converting enzyme (ACE) inhibitor dipeptides, bitter-tasting dipeptides, and nonameric binding peptides of the human leukocyte antigens (HLA-A*0201). Multiple linear regression was used to construct the QSAR models. For each peptide set, a proper subset of physicochemical properties was chosen by the ant colony optimization algorithm. The leave-one-out cross-validation (q _loo ²) values were 0.855, 0.936, and 0.642 and the root-mean-square errors (RMSEs) were 0.450, 0.149, and 0.461. Our results revealed that the new numerical descriptive vector can afford extensive characterization of peptide sequence so that it can be easily employed in peptide QSAR studies. Moreover, the proposed numerical descriptive vectors were able to determine hot spot residues in the peptides under study.

Quantitative sequence-activity modeling of ACE peptide originated from milk using ACC-QTMS amino acid indices

Up to now, numerous peptides/hydrolysates derived from casein and whey protein have shown angiotensin-I-converting enzyme (ACE) inhibitory. In this research, quantum topological molecular similarity (QTMS) indices of amino acids were utilized in quantitative sequence-activity modeling (QSAM) to predict the activity of a set of milk-driven peptides with ACE inhibition. Since the derived peptides have not the same number of residues, we overcame this issue by auto cross covariance (ACC) methodology. Then, some QSAMs were built to predict the pIC₅₀ value of ACE peptides derived from Bovine Casein and Whey. The model established an acceptable relationship between the selected variables and the pIC₅₀ of the peptides. To estimate the performance of the developed models, casein and whey proteins from human, goat, bovine and sheep were virtually broken by trypsin and chymotrypsin enzymes and the ACE activity of the resultant virtual peptides were predicted and some new ACE peptides were proposed.

Assessment of cytotoxicity of choline chloride-based natural deep eutectic solvents against human HEK-293 cells: A QSAR analysis

Deep eutectic solvents (DESs) are a new generation of solvents. To consider them as green solvents, investigation of their toxicity is essential. In this work, the cytotoxicity of a number of natural deep eutectic solvents (NADESs) against HEK-293 human embryonic kidney cells was evaluated by MTT assay. The NADESs were prepared with choline chloride (ChCl) as hydrogen-bond acceptor (HBA) and different sugar alcohols as hydrogen-bond donor (HBD) constituents. They showed IC₅₀ values in the range of 3.52-75.46 mM. These results were used to evaluate the effect of structural parameters on the cytotoxicity of the studied NADESs by using quantitative structure activity relationship (QSAR) analysis. A three-parameter linear model was obtained between - log(IC₅₀) as a dependent variable and structural descriptors as independent variables. Rotatable bond number (RBN), mean atomic van der Waals volume (Mv) and the interaction of second power carbon numbers with the molar ratio of HBA to HBD in each NADES (C² Ratio), were three major parameters. The statistical model covered about 76.4% and 69.8% variance of data in training and leave-one-out cross-validation, respectively. This work, as the first study on the QSAR analysis of DESs, can provide a good perspective for designing greener novel DESs.

An Introduction to the Basic Concepts in QSAR-Aided Drug Design

The need for the development of new drugs to combat existing and newly identified conditions is unavoidable. One of the important tools used in the advanced drug development pipeline is computer-aided drug design. Traditionally, to find a drug many ligands were synthesized and evaluated for their effectiveness using suitable bioassays and if all other drug-likeness features were met, the candidate(s) would possibly reach the market. Although this approach is still in use in advanced format, computational methods are an indispensable component of modern drug development projects. One of the methods used from very early days of rationalizing the drug design approaches is Quantitative Structure-Activity Relationship (QSAR). This chapter overviews QSAR modeling steps by introducing molecular descriptors, mathematical model development for relating biological activities to molecular structures, and model validation. At the end, several successful cases where QSAR studies were used extensively are presented.

A New Descriptor of Amino Acids-SVGER and its Applications in Peptide QSAR

In the study of peptide quantitative structure activity relationship (QSAR), a new descriptor of amino acids (SVGER) was calculated. It was applied in two peptides which are angiotensin converting enzyme inhibitors and bitter tasting threshold of di-peptide. QSAR models were built by stepwise multiple regression-multiple linear regression (SMR-MLR) and stepwise multiple regression-partial least square regression (SMR-PLS). In the SMR-MLR models for angiotensin converting enzyme inhibitors, the squared cross-validation correlation coefficient (Q_LOO² ) was 0.907, squared correlation coefficient between predicted and observed activities (R_cum² ) was 0.977 and external multiple correlation coefficient (Q_ext² ) was 0.867. The corresponding data for the bitter tasting threshold of di-peptide were 0.802, 0.966, 0.719. While in the SMR-PLS model, Q_LOO² , R_cum² and Q_ext² were 0.804, 0.915, 0.858 for angiotensin converting enzyme inhibitors and 0.782, 0.881, 0.747 for bitter tasting threshold of di-peptide. Our results showed that descriptor SVGER can afford good account of relationships between activity and structure of peptide drugs.

Variable selection in multivariate calibration based on clustering of variable concept

Recently we have proposed a new variable selection algorithm, based on clustering of variable concept (CLoVA) in classification problem. With the same idea, this new concept has been applied to a regression problem and then the obtained results have been compared with conventional variable selection strategies for PLS. The basic idea behind the clustering of variable is that, the instrument channels are clustered into different clusters via clustering algorithms. Then, the spectral data of each cluster are subjected to PLS regression. Different real data sets (Cargill corn, Biscuit dough, ACE QSAR, Soy, and Tablet) have been used to evaluate the influence of the clustering of variables on the prediction performances of PLS. Almost in the all cases, the statistical parameter especially in prediction error shows the superiority of CLoVA-PLS respect to other variable selection strategies. Finally the synergy clustering of variable (sCLoVA-PLS), which is used the combination of cluster, has been proposed as an efficient and modification of CLoVA algorithm. The obtained statistical parameter indicates that variable clustering can split useful part from redundant ones, and then based on informative cluster; stable model can be reached.

An Introduction to the Basic Concepts in QSAR-Aided Drug Design

The need for the development of new drugs to combat existing and newly identified conditions is unavoidable. One of the important tools used in the advanced drug development pipeline is computer-aided drug design. Traditionally, to find a drug many ligands were synthesized and evaluated for their effectiveness using suitable bioassays and if all other drug-likeness features were met, the candidate(s) would possibly reach the market. Although this approach is still in use in advanced format, computational methods are an indispensable component of modern drug development projects. One of the methods used from very early days of rationalizing the drug design approaches is Quantitative Structure-Activity Relationship (QSAR). This chapter overviews QSAR modeling steps by introducing molecular descriptors, mathematical model development for relating biological activities to molecular structures, and model validation. At the end, several successful cases where QSAR studies were used extensively are presented.

A pipeline for improved QSAR analysis of peptides: physiochemical property parameter selection via BMSF, near-neighbor sample selection via semivariogram, and weighted SVR regression and prediction

In this paper, we present a pipeline to perform improved QSAR analysis of peptides. The modeling involves a double selection procedure that first performs feature selection and then conducts sample selection before the final regression analysis. Five hundred and thirty-one physicochemical property parameters of amino acids were used as descriptors to characterize the structure of peptides. These high-dimensional descriptors then go through a feature selection process given by the binary matrix shuffling filter (BMSF) to obtain a set of important low-dimensional features. Each descriptor that passes the BMSF filtering also receives a weight defined through its contribution to reduce the estimation error. These selected features served as the predictors for subsequent sample selection and modeling. Based on the weighted Euclidean distances between samples, a common range was determined with high-dimensional semivariogram and then used as a threshold to select the near-neighbor samples from the training set. For each sample to be predicted, the QSAR model was established using SVR with the weighted, selected features based on the exclusive set of near-neighbor training samples. Prediction was conducted for each test sample accordingly. The performances of this pipeline are tested with the QSAR analysis of angiotensin-converting enzyme inhibitors and HLA-A*0201 data sets. Improved prediction accuracy was obtained in both applications. This pipeline can optimize the QSAR modeling from both the feature selection and sample selection perspectives. This leads to improved accuracy over single selection methods. We expect this pipeline to have extensive application prospect in the field of regression prediction.

PAAQD: Predicting immunogenicity of MHC class I binding peptides using amino acid pairwise contact potentials and quantum topological molecular similarity descriptors

Prediction of peptide immunogenicity is a promising approach for novel vaccine discovery. Conventionally, epitope prediction methods have been developed to accelerate the process of vaccine production by searching for candidate peptides from pathogenic proteins. However, recent studies revealed that peptides with high binding affinity to major histocompatibility complex molecules (MHCs) do not always result in high immunogenicity. Therefore, it is promising to predict the peptide immunogenicity rather than epitopes in order to discover new vaccines more effectively. To this end, we developed a novel T-cell reactivity predictor which we call PAAQD. Nonapeptides were encoded numerically, using combining information of amino acid pairwise contact potentials (AAPPs) and quantum topological molecular similarity (QTMS) descriptors. Encoded data were used in the construction of our classification model. Our numerical experiments suggested that the predictive performance of PAAQD is at least comparable with POPISK, one of the pioneering techniques for T-cell reactivity prediction. Also, our experiment suggested that the first and eighth positions of nonapeptides are the most important for immunogenicity and most of the anchor residues in epitope prediction were not important in T-cell reactivity prediction. The R implementation of PAAQD is available at http://pirun.ku.ac.th/~fsciiok/PAAQD.rar.