Singular value decomposition analysis of the torsional angles of dopamine reuptake inhibitor GBR 12909 analogs: effect of force field and charges.

Three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis of large, flexible molecules, such as the dopamine reuptake inhibitor GBR 12909 (1), is complicated by the fact that they can take on a wide range of closely-related conformations. The first step in the analysis is to classify the conformers into groups. Over 600 conformers each of a piperazine (2) and piperidine (3) analog of 1 were generated by random search conformational analysis using the Merck Molecular Force Field (MMFF94). Singular value decomposition (SVD) was used to group the conformers of 2 and 3 by the similarity of their non-ring torsional angles. SVD uncovered subtle differences in their conformer populations due to that fact that the conformers separate along different principal components, and ultimately to the fact that different torsional angles are the chief contributors to these components. The results were compared to our previous SVD analysis (Fiorentino, et al., Journal of Computational Chemistry, 2006, 27, 609-620) of conformer populations of 2 and 3 generated by the Tripos force field and Gasteiger-Hückel charges. Except for the dominant contribution of angle B3 to principal component 8 seen with both force fields, the angles which are chiefly responsible for the grouping of the conformers of 2 and 3 are different with both force fields. This illustrates that SVD is useful in identifying unique groupings of conformers in large data sets of flexible molecules-a first step in selecting representative conformers for 3D-QSAR modeling studies.

Conformational analysis of piperazine and piperidine analogs of GBR 12909: stochastic approach to evaluating the effects of force fields and solvent.

Analogs of the flexible dopamine reuptake inhibitor, GBR 12909 (1), may have potential utility in the treatment of cocaine abuse. As a first step in the 3D-QSAR modeling of the dopamine transporter (DAT)/serotonin transporter (SERT) selectivity of these compounds, we carried out conformational analyses of two analogs of 1: a piperazine (2) and a related piperidine (3). Ensembles of conformers consisting of local minima on the potential energy surface of the molecule were generated in the vacuum phase and in implicit solvent by random search conformational analysis using the Tripos and MMFF94 force fields. Some differences were noted in the conformer populations due to differences in the treatment of the tertiary amine nitrogen and ether oxygen atom types by the force fields. The force fields also differed in their descriptions of internal rotation around the C(sp³)-O(sp³) bond proximal to the bisphenyl moiety. Molecular orbital calculations at the HF/6-31G(d) and B3LYP/6-31G(d) levels of C-O internal rotation in model compound (5), designed to model the effect of the proximity of the bisphenyl group on C-O internal rotation, showed a broad region of low energy between -60° to 60° with minima at both -60° and 30° and a low rotational barrier at 0°, in closer agreement with the MMFF94 results than the Tripos results. Molecular mechanics calculations on model compound (6) showed that the MMFF94 force field was much more sensitive than the Tripos force field to the effects of the bisphenyl moiety on C-O internal rotation.

Quantitative structure-activity relationship studies of threo-methylphenidate analogs.

Complementary two-dimensional (2D) and three-dimensional (3D) Quantitative Structure-Activity Relationship (QSAR) techniques were used to derive a preliminary model for the dopamine transporter (DAT) binding affinity of 80 racemic threo-methylphenidate (MP) analogs. A novel approach based on using the atom-level E-state indices of the 14 common scaffold atoms in a sphere exclusion protocol was used to identify a test set for 2D- and 3D-QSAR model validation. Comparative Molecular Field Analysis (CoMFA) contour maps based on the structure-activity data of the training set indicate that the 2' position of the phenyl ring cannot tolerate much steric bulk and that addition of electron-withdrawing groups to the 3' or 4' positions of the phenyl ring leads to improved DAT binding affinity. In particular, the optimal substituents were found to be those whose bulk is mainly in the plane of the phenyl ring. Substituents with significant bulk above or below the plane of the ring led to decreased binding affinity. Suggested alterations to be explored in the design of new compounds are the placement at the 3' and 4' position of the phenyl ring of electron-withdrawing groups that lie chiefly in the plane of the ring, for example, halogen substituents on the 3',4'-benzo analog, 79. A complementary 2D-QSAR approach-partial least squares analysis using a reduced set of Molconn-Z descriptors-supports the CoMFA structure-activity interpretation that phenyl ring substitution is a major determinant of DAT binding affinity. The potential usefulness of the CoMFA models was demonstrated by the prediction of the binding affinity of methyl 2-(naphthalen-1-yl)-2-(piperidin-2-yl)acetate, an analog not in the original data set, to be in good agreement with the experimental value.

Feature extraction using molecular planes for fuzzy relational clustering of a flexible dopamine reuptake inhibitor.

Six rigid-body parameters (Shift, Slide, Rise, Tilt, Roll, Twist) are commonly used to describe the relative displacement and orientation of successive base pairs in a nucleic acid structure. The present work adapts this approach to describe the relative displacement and orientation of any two planes in an arbitrary molecule-specifically, planes which contain important pharmacophore elements. Relevant code from the 3DNA software package (Nucleic Acids Res. 2003, 31, 5108-5121) was generalized to treat molecular fragments other than DNA bases as input for the calculation of the corresponding rigid-body (or "planes") parameters. These parameters were used to construct feature vectors for a fuzzy relational clustering study of over 700 conformations of a flexible analogue of the dopamine reuptake inhibitor, GBR 12909. Several cluster validity measures were used to determine the optimal number of clusters. Translational (Shift, Slide, Rise) rather than rotational (Tilt, Roll, Twist) features dominate clustering based on planes that are relatively far apart, whereas both types of features are important to clustering when the pair of planes are close by. This approach was able to classify the data set of molecular conformations into groups and to identify representative conformers for use as template conformers in future Comparative Molecular Field Analysis studies of GBR 12909 analogues. The advantage of using the planes parameters, rather than the combination of atomic coordinates and angles between molecular planes used in our previous fuzzy relational clustering of the same data set (J. Chem. Inf. Model. 2005, 45, 610-623), is that the present clustering results are independent of molecular superposition and the technique is able to identify clusters in the molecule considered as a whole. This approach is easily generalizable to any two planes in any molecule.

DAT/SERT selectivity of flexible GBR 12909 analogs modeled using 3D-QSAR methods.

The dopamine reuptake inhibitor GBR 12909 (1-{2-[bis(4-fluorophenyl)methoxy]ethyl}-4-(3-phenylpropyl)piperazine, 1) and its analogs have been developed as tools to test the hypothesis that selective dopamine transporter (DAT) inhibitors will be useful therapeutics for cocaine addiction. This 3D-QSAR study focuses on the effect of substitutions in the phenylpropyl region of 1. CoMFA and CoMSIA techniques were used to determine a predictive and stable model for the DAT/serotonin transporter (SERT) selectivity (represented by pK(i) (DAT/SERT)) of a set of flexible analogs of 1, most of which have eight rotatable bonds. In the absence of a rigid analog to use as a 3D-QSAR template, six conformational families of analogs were constructed from six pairs of piperazine and piperidine template conformers identified by hierarchical clustering as representative molecular conformations. Three models stable to y-value scrambling were identified after a comprehensive CoMFA and CoMSIA survey with Region Focusing. Test set correlation validation led to an acceptable model, with q(2)=0.508, standard error of prediction=0.601, two components, r(2)=0.685, standard error of estimate=0.481, F value=39, percent steric contribution=65, and percent electrostatic contribution=35. A CoMFA contour map identified areas of the molecule that affect pK(i) (DAT/SERT). This work outlines a protocol for deriving a stable and predictive model of the biological activity of a set of very flexible molecules.

Hierarchical clustering analysis of flexible GBR 12909 dialkyl piperazine and piperidine analogs.

Pharmacophore modeling of large, drug-like molecules, such as the dopamine reuptake inhibitor GBR 12909, is complicated by their flexibility. A comprehensive hierarchical clustering study of two GBR 12909 analogs was performed to identify representative conformers for input to three-dimensional quantitative structure-activity relationship studies of closely-related analogs. Two data sets of more than 700 conformers each produced by random search conformational analysis of a piperazine and a piperidine GBR 12909 analog were studied. Several clustering studies were carried out based on different feature sets that include the important pharmacophore elements. The distance maps, the plot of the effective number of clusters versus actual number of clusters, and the novel derived clustering statistic, percentage change in the effective number of clusters, were shown to be useful in determining the appropriate clustering level. Six clusters were chosen for each analog, each representing a different region of the torsional angle space that determines the relative orientation of the pharmacophore elements. Conformers of each cluster that are representative of these regions were identified and compared for each analog. This study illustrates the utility of using hierarchical clustering for the classification of conformers of highly flexible molecules in terms of the three-dimensional spatial orientation of key pharmacophore elements.

Singular value decomposition of torsional angles of analogs of the dopamine reuptake inhibitor GBR 12909.

Analysis of large, flexible molecules, such as the dopamine reuptake inhibitor GBR 12909 (1), is complicated by the fact that they can take on a wide range of closely related conformations. The first step in the analysis is to classify the conformers into groups. Here, Singular Value Decomposition (SVD) was used to group conformations of GBR 12909 analogs by the similarity of their nonring torsional angles. The significance of the present work, the first application of SVD to the analysis of very flexible molecules, lies in the development of a novel scaling technique for circular data and in the grouping of molecular conformations using a technique that is independent of molecular alignment. Over 700 conformers each of a piperazine (2) and piperidine (3) analog of 1 were studied. Analysis of the score and loading plots showed that the conformers of 2 separate into three large groups due to torsional angles on the naphthalene side of the molecule, whereas those of 3 separate into nine groups due to torsional angles on the bisphenyl side of the molecule. These differences are due to nitrogen inversion at the unprotonated piperazinyl nitrogen of 2, which results in a different ensemble of conformers than those of 3, where no inversion is possible at the corresponding piperidinyl carbon.

Novel feature extraction technique for fuzzy relational clustering of a flexible dopamine reuptake inhibitor.

This paper describes a novel clustering methodology for classifying over 700 conformations of a flexible analogue of GBR 12909, a dopamine reuptake inhibitor that has completed phase I clinical trials as a treatment for cocaine abuse. The major aspect of the clustering methodology includes an efficient data-conditioning scheme where a systematic feature extraction procedure based on the structural properties of the molecule was used to reduce the associated feature space. This allowed region-specific clustering that focused on individual pharmacophore elements of the molecule. For clustering of the reduced feature set, the fuzzy clustering partitional method was utilized. Due to the relational nature of the feature data, fuzzy relational clustering was employed, and it successfully detected natural groups defined by rotational minima around N(sp(3))-C(sp(3)), O(sp(3))-C(sp(3)), and C(sp(3))-C(sp(2)) bonds. The proposed clustering methodology also employed several cluster validity measures, which corroborated the partitions produced by the clustering technique and agreed with the results of hierarchical clustering using the XCluster program. Representative structures which exhibited a reasonable spread of energies and showed good spatial coverage of the conformational space were identified for use as putative bioactive conformations in a future Comparative Molecular Field Analysis of GBR 12909 analogues. The clustering methodology developed here is capable of handling other computational chemistry problems, and the feature extraction technique can be easily generalized to other molecules.

Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods.

Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte > MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the phenyl ring rotational barrier for neutral MP and gave results very similar to those of the HF/6-31G* method.