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.

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.