Development of a Chirality-Sensitive Flexibility Descriptor for 3+3D-QSAR

Recent advances in chiral selectors immobilization and chiral mobile phase additives in liquid chromatographic enantio-separations: A review

For decades now, the separation of chiral enantiomers of drugs has been gaining the interest and attention of researchers. In 1991, the first guidelines for development of chiral drugs were firstly released by the US-FDA. Since then, the development in chromatographic enantioseparation tools has been fast and variable, aiming at creating a suitable environment where the physically and chemically identical enantiomers can be separated. Among those tools, the immobilization of chiral selectors (CS) on different stationary phases and the chiral mobile phase additives (CMPA) which have been progressed and studied extensively. This review article highlights the major advances in immobilization of CS together with their different recognition mechanisms as well as CMPA as a cheaper and successful alternative for chiral stationary phases. Moreover, the role of molecular modeling tool as a pre-step in the choice of CS for evaluating possible interactions with different ligands has been pointed up. Illustrations of reported methods and updates for immobilized CS and CMPA have been included.

Modelling the enantiorecognition of structurally diverse pharmaceuticals on O-substituted polysaccharide-based stationary phases

This study aims to develop models to predict the retention, separation and elution sequence of the enantiomers of structurally diverse pharmaceuticals. More specifically, Quantitative Structure Retention Relationships (QSRR) models are built that describe the relationship between molecular descriptors and retention. Eighteen structurally diverse chiral mixtures, each consisting of a pair of enantiomers, were analyzed on two polysaccharide chiral stationary phases, Chiralcel OD-RH (cellulose tris(3,5-dimethylphenylcarbamate)) and Lux amylose-2 (amylose tris(5-chloro-2-methylphenylcarbamate)), applying either a basic or an acidic mobile phase, and their retention factor and elution sequence were determined. Both achiral and, in-house defined, chiral descriptors were used as descriptive variables to build the models. Linear regression techniques, i.e. stepwise multiple linear regression (sMLR) and partial least squares (PLS) regression, were applied to model the retention or separation as a function of the descriptors. In a first step, models were built with only achiral descriptors to model the global retention of both enantiomers of a chiral molecule. Subsequently, models were built with only chiral descriptors to predict the enantioseparation and elution sequence, and finally, models were considered with both descriptor types to predict the retention, the separation and the elution sequence of the enantiomers. The global retention was predicted well by the sMLR models with only achiral descriptors. The models with only chiral descriptors were not found suitable to predict the enantioseparation and elution sequence. Finally, the models containing both chiral and achiral descriptors allowed predicting the retention well, but their ability to predict the elution sequence and separation of the enantiomers differed widely for the chromatographic systems considered.

Beyond Rotatable Bond Counts: Capturing 3D Conformational Flexibility in a Single Descriptor

A new molecular descriptor, nConf20, based on chemical connectivity, is presented which captures the accessible conformational space of a molecule. Currently the best available two-dimensional descriptors for quantifying the flexibility of a particular molecule are the rotatable bond count (RBC) and the Kier flexibility index. We present a descriptor which captures this information by sampling the conformational space of a molecule using the RDKit conformer generator. Flexibility has previously been identified as a key feature in determining whether a molecule is likely to crystallize or not. For this application, nConf20 significantly outperforms previously reported single-variable classifiers and also assists rule-based analysis of black-box machine learning classification algorithms.

Zwitterionic structures: from physicochemical properties toward computer-aided drug designs

Zwitterions, used widely in chemical, biological and medicinal fields, show distinct physicochemical properties relative to ordinary ampholytes, which largely decide their bioavailability and biological activities. In the present manuscript, these properties are discussed in order to facilitate our understanding of zwitterionic structures, followed by various examples of zwitterionic drugs and the critical role these properties play. We specifically focus our discussions on neuraminidase inhibitors (NAIs), which are used in the treatment and prevention of influenza, covering their computer-assisted design, transformation to zwitterionic isomers and interaction mechanisms of NAIs with proteins. The discovery and development of NAIs provide useful insights that may assist in the exploration of new zwitterionic drugs.

4D-QSAR study of HEPT derivatives by electron conformational-genetic algorithm method

In this work, the EC-GA method, a hybrid 4D-QSAR approach that combines the electron conformational (EC) and genetic algorithm optimization (GA) methods, was applied in order to explain pharmacophore (Pha) and predict anti-HIV-1 activity by studying 115 compounds in the class of 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio) thymine (HEPT) derivatives as non-nucleoside reverse transcriptase inhibitors (NNRTIs). The series of NNRTIs were partitioned into four training and test sets from which corresponding quantitative structure-activity relationship (QSAR) models were constructed. Analysis of the four QSAR models suggests that the three models generated from the training and test sets used in previous works yielded comparable results with those of previous studies. Model 4, the data set of which was partitioned randomly into two training and test sets with 11 descriptors, including electronical and geometrical parameters, showed good statistics both in the regression (r2(training) )= 0.867, r2test = 0.923) and cross-validation (q (2) = 0.811, q2(ext1) = 0.909, q2(ext2) = 0.909) for the training set of 80 compounds and the test set of 27 compounds. The prediction of the anti-HIV-1 activity of HEPT compounds by means of the EC-GA method allowed for a quantitatively consistent QSAR model. In addition, eight novel compounds never tested experimentally have been designed theoretically using model 4.

Probabilistic Modeling of Conformational Space for 3D Machine Learning Approaches

We present a new probabilistic encoding of the conformational space of a molecule that allows for the integration into common similarity calculations. The method uses distance profiles of flexible atom-pairs and computes generative models that describe the distance distribution in the conformational space. The generative models permit the use of probabilistic kernel functions and, therefore, our approach can be used to extend existing 3D molecular kernel functions, as applied in support vector machines, to build QSAR models. The resulting kernels are valid 4D kernel functions and reduce the dependency of the model quality on suitable conformations of the molecules. We showed in several experiments the robust performance of the 4D kernel function, which was extended by our approach, in comparison to the original 3D-based kernel function. The new method compares the conformational space of two molecules within one kernel evaluation. Hence, the number of kernel evaluations is significantly reduced in comparison to common kernel-based conformational space averaging techniques. Additionally, the performance gain of the extended model correlates with the flexibility of the data set and enables an a priori estimation of the model improvement.

Atomic local neighborhood flexibility incorporation into a structured similarity measure for QSAR

In this work, we introduce a new method to regard the geometry in a structural similarity measure by approximating the conformational space of a molecule. Our idea is to break down the molecular conformation into the local conformations of neighbor atoms with respect to core atoms. This local geometry can be implicitly accessed by the trajectories of the neighboring atoms, which are emerge by rotatable bonds. In our approach, the physicochemical atomic similarity, which can be used in structured similarity measures, is augmented by a local flexibility similarity, which gives a rough estimate of the similarity of the local conformational space. We incorporated this new type of encoding the flexibility into the optimal assignment molecular similarity approach, which can be used as a pseudokernel in support vector machines. The impact of the local flexibility was evaluated on several published QSAR data sets. This lead to an improvement of the model quality on 9 out of 10 data sets compared to the unmodified optimal assignment kernel.

Antibody recognition and conformational flexibility of a plaque-specific beta-amyloid epitope modulated by non-native peptide flanking regions

Here we report on the synthesis, antibody binding, and QSAR studies of a series of linear and cyclic peptides containing a beta-amyloid plaque-specific epitope (Abeta(4-10); FRHDSGY). In these constructs, two or three alpha- l-Ala, alpha- d-Ala, or beta-Ala residues were introduced at both N- and C-termini of the epitope as non-native flanking sequences. Cyclization of the linear Abeta(4-10) epitope peptide resulted in reduced antibody binding. However, the antibody binding could be fully compensated by insertion of alanine flanks into the corresponding cyclic peptides. These results indicate that the modification of a beta-amyloid plaque-specific epitope by combination of cyclization and flanking sequences could generate highly antigenic peptides compared to the native sequence. A novel 3D QSAR method, which explicitly handles conformational flexibility, was developed for the case of such molecular libraries. This method led to the prediction of the binding conformation for the common FRHDSGY sequence.

Multi-dimensional QSAR in drug discovery

Quantitative structure-activity relationships (QSAR) is an area of computational research that builds virtual models to predict quantities such as the binding affinity or the toxic potential of existing or hypothetical molecules. Although a wealth of experimental data emphasizes the active role of the target protein in the binding process, QSAR studies are frequently restricted to the properties of the small-molecule ligand. This review aims at discussing recent QSAR concepts exploring higher dimensions (simulation of induced fit, simultaneous exploration of alternative binding modes, and solvation scenarios), and their benefit for the drug-discovery process.

Novel Approach for the Numerical Characterization of Molecular Chirality

The use of chiral compounds as pharmaceuticals and agrochemicals continues to increase, warranting numerical characterization of chirality in order to develop structure-activity relationship models involving these compounds. Enantiomers are identical in all scalar properties and, hence, are not differentiated by topological indices and 3-D descriptors. Three distinct measures of chirality were developed to discriminate diastereomers and enantiomers. The novel topological indices treat chirality as a continuous measure, and hence we prefer to call it the Relative Chirality Index (RCI). Application of RCI in developing SAR is illustrated with the repellency data for the diastereomers of picaridin and AI3-37220.