Conformational coverage by a genetic algorithm: saturation of conformational space.

The molecular modeling is traditionally based on analysis of minimum energy conformers. Such simplifying assumptions could doom to failure the modeling studies given the significant variation of the geometric and electronic characteristics across the multitude of energetically reasonable conformers representing the molecules. Moreover, it has been found that the lowest energy conformers of chemicals are not necessarily the active ones with respect to various endpoints. Hence, the selection of active conformers appears to be as important as the selection of molecular descriptors in the modeling process. In this respect, we have developed effective tools for conformational analysis based on a genetic algorithm (GA), published in J. Chem. Inf. Comput. Sci. (1994, 34, 234; 1999, 39 (6), 997) and J. Chem. Inf. Model. (2005, 45 (2), 283). This paper presents a further improvement of the evolutionary algorithm for conformer generation minimizing the sensitivity of conformer distributions from the effect of smoothing parameter and improving the reproducibility of conformer distributions given the nondeterministic character of the genetic algorithm (GA). The ultimate goal of the saturation is to represent the conformational space of chemicals with an optimal number of conformers providing a stable conformational distribution which cannot be further perturbed by the addition of new conformers. The generation of stable conformational distributions of chemicals by a limited number of conformers will improve the adequacy of the subsequent molecular modeling analysis. The impact of the saturation procedure on conformer distributions in a specific structural space is illustrated by selected examples. The effect of the procedure on similarity assessment between chemicals is discussed.

Identifying the structural requirements for chromosomal aberration by incorporating molecular flexibility and metabolic activation of chemicals.

Modeling the potential of chemicals to induce chromosomal damage has been hampered by the diversity of mechanisms which condition this biological effect. The direct binding of a chemical to DNA is one of the underlying mechanisms that is also responsible for bacterial mutagenicity. Disturbance of DNA synthesis due to inhibition of topoisomerases and interaction of chemicals with nuclear proteins associated with DNA (e.g., histone proteins) were identified as additional mechanisms leading to chromosomal aberrations (CA). A comparative analysis of in vitro genotoxic data for a large number of chemicals revealed that more than 80% of chemicals that elicit bacterial mutagenicity (as indicated by the Ames test) also induce CA; alternatively, only 60% of chemicals that induce CA have been found to be active in the Ames test. In agreement with this relationship, a battery of models is developed for modeling CA. It combines the Ames model for bacterial mutagenicity, which has already been derived and integrated into the Optimized Approach Based on Structural Indices Set (OASIS) tissue metabolic simulator (TIMES) platform, and a newly derived model accounting for additional mechanisms leading to CA. Both models are based on the classical concept of reactive alerts. Some of the specified alerts interact directly with DNA or nuclear proteins, whereas others are applied in a combination of two- or three-dimensional quantitative structure-activity relationship models assessing the degree of activation of the alerts from the rest of the molecules. The use of each of the alerts has been justified by a mechanistic interpretation of the interaction. In combination with a rat liver S9 metabolism simulator, the model explained the CA induced by metabolically activated chemicals that do not elicit activity in the parent form. The model can be applied in two ways: with and without metabolic activation of chemicals.

Identification of the structural requirements for mutagenicity by incorporating molecular flexibility and metabolic activation of chemicals I: TA100 model.

Traditional attempts to model genotoxicity data have been limited to congeneric data sets, primarily because the mechanism of action was ignored, and frequently, the chemicals required metabolism to the active species. In this exercise, the COmmon REactivity PAtterns (COREPA) approach was used to delineate the structural requirements for eliciting mutagenicity in terms of ranges of descriptors associated with three-dimensional molecular structures. The database used to build the mutagenicity model includes 1196 structurally diverse chemicals tested in the Ames assay by the National Toxicology Program. This manuscript describes the development of the TA100 model that predicts the results of mutagenicity testing using only the Ames TA100 strain. The TA100 model was developed using 148 chemicals that tested positive in TA100 strain without rat liver enzymes (S-9) and 188 chemicals that tested positive in TA100 strain with rat liver enzymes. A decision tree was developed by first comparing the reactivity profile of chemicals that were positive in TA100 without rat liver enzymes to the reactivity profile of the remaining 1048 chemicals. This approach correctly identified 82% of the primary acting mutagens and 94% of the nonmutagens in the training set. The 188 chemicals in the training set that are positive only in the presence of metabolic activation would pass through the decision tree as negative. The next step was to identify the chemicals that are positive only in the presence of metabolic activation. To accomplish this, a series of hierarchically ordered metabolic transformations were used to develop an S-9 metabolism simulator that was applied to each of the 1048 chemicals. The potential metabolites were then screened through the decision tree to identify reactive mutagens. This model correctly identified 77% of the metabolically activated chemicals in a training set. A computer system that applies the COREPA models and predicts mutagenicity of chemicals, including their metabolic activation, was developed. Each prediction is accompanied by a probabilistic estimate of the chemical being in the structural domain covered by the training set.