Theoretical studies of the mechanism of N-hydroxylation of primary aromatic amines by cytochrome P450 1A2: radicaloid or anionic?

Primary aromatic and heteroaromatic amines are notoriously known as potential mutagens and carcinogens. The major event of the mechanism of their mutagenicity is N-hydroxylation by P450 enzymes, primarily P450 1A2 (CYP1A2), which leads to the formation of nitrenium ions that covalently modify nucleobases of DNA. Energy profiles of the NH bond activation steps of two possible mechanisms of N-hydroxylation of a number of aromatic amines by CYP1A2, radicaloid and anionic, are studied by dispersion-corrected DFT calculations. The classical radicaloid mechanism is mediated by H-atom transfer to the electrophilic ferryl-oxo intermediate of the P450 catalytic cycle (called Compound I or Cpd I), whereas the alternative anionic mechanism involves proton transfer to the preceding nucleophilic ferrous-peroxo species. The key structural features of the catalytic site of human CYP1A2 revealed by X-ray crystallography are maintained in calculations. The obtained DFT reaction profiles and additional calculations that account for nondynamical electron correlation suggest that Cpd I has higher thermodynamic drive to activate aromatic amines than the ferrous-peroxo species. Nevertheless, the anionic mechanism is demonstrated to be consistent with a variety of experimental observations. Thus, energy of the proton transfer from aromatic amines to the ferrous-peroxo dianion splits aromatic amines into two classes with different mutagenicity mechanisms. Favorable or slightly unfavorable barrier-free proton transfer is inherent in compounds that undergo nitrenium ion mediated mutagenicity. Monocyclic electron-rich aromatic amines that do not follow this mutagenicity mechanism show significantly unfavorable proton transfer. Feasibility of the entire anionic mechanism is demonstrated by favorable Gibbs energy profiles of both chemical steps, NH bond activation, and NO bond formation. Taken together, results suggest that the N-hydroxylation of aromatic amines in CYP1A2 undergoes the anionic mechanism. Possible reasons for the apparent inability of Cpd I to activate aromatic amines in CYP1A2 are discussed.

Experimental testing of quantum mechanical predictions of mutagenicity: aminopyrazoles.

A computational method for predicting the likelihood of aromatic amines being active in the Ames test for mutagenicity was trialed on a set of aminopyrazoles. A virtual array of compounds was generated from the available sets of hydrazines and α-cyanoaldehydes (or ketones) and quantum mechanical calculations used to compute a probability of being active in the Ames test. The compounds selected for synthesis and testing were not based on the predictions and so spanned the range of predicted probabilities. The subsequently generated results of the Ames test were in good correspondence with the predictions and confirm this approach as a useful means of predicting likely mutagenic risk.

Theoretical studies of chemical reactivity of metabolically activated forms of aromatic amines toward DNA.

The metabolism of aromatic and heteroaromatic amines (ArNH2) results in nitrenium ions (ArNH⁺) that modify nucleobases of DNA, primarily deoxyguanosine (dG), by forming dG-C8 adducts. The activated amine nitrogen in ArNH⁺ reacts with the C8 of dG, which gives rise to mutations in DNA. For the most mutagenic ArNH2, including the majority of known genotoxic carcinogens, the stability of ArNH⁺ is of intermediate magnitude. To understand the origin of this observation as well as the specificity of reactions of ArNH⁺ with guanines in DNA, we investigated the chemical reactivity of the metabolically activated forms of ArNH2, that is, ArNHOH and ArNHOAc, toward 9-methylguanine by DFT calculations. The chemical reactivity of these forms is determined by the rate constants of two consecutive reactions leading to cationic guanine intermediates. The formation of ArNH⁺ accelerates with resonance stabilization of ArNH⁺, whereas the formed ArNH⁺ reacts with guanine derivatives with the constant diffusion-limited rate until the reaction slows down when ArNH⁺ is about 20 kcal/mol more stable than PhNH⁺. At this point, ArNHOH and ArNHOAc show maximum reactivity. The lowest activation energy of the reaction of ArNH⁺ with 9-methylguanine corresponds to the charge-transfer π-stacked transition state (π-TS) that leads to the direct formation of the C8 intermediate. The predicted activation barriers of this reaction match the observed absolute rate constants for a number of ArNH⁺. We demonstrate that the mutagenic potency of ArNH2 correlates with the rate of formation and the chemical reactivity of the metabolically activated forms toward the C8 atom of dG. On the basis of geometric consideration of the π-TS complex made of genotoxic compounds with long aromatic systems, we propose that precovalent intercalation in DNA is not an essential step in the genotoxicity pathway of ArNH2. The mechanism-based reasoning suggests rational design strategies to avoid genotoxicity of ArNH2 primarily by preventing N-hydroxylation of ArNH2.

Explanation for main features of structure-genotoxicity relationships of aromatic amines by theoretical studies of their activation pathways in CYP1A2.

Aromatic and heteroaromatic amines (ArNH(2)) represent a class of potential mutagens that after being metabolically activated covalently modify DNA. Activation of ArNH(2) in many cases starts with N-hydroxylation by P450 enzymes, primarily CYP1A2. Poor understanding of structure-mutagenicity relationships of ArNH(2) limits their use in drug discovery programs. Key factors that facilitate activation of ArNH(2) are revealed by exploring their reaction intermediates in CYP1A2 using DFT calculations. On the basis of these calculations and extensive analysis of structure-mutagenicity data, we suggest that mutagenic metabolites are generated by ferric peroxo intermediate, (CYP1A2)Fe(III)-OO(-), in a three-step heterolytic mechanism. First, the distal oxygen of the oxidant abstracts proton from H-bonded ArNH(2). The subsequent proximal protonation of the resulting (CYP1A2)Fe(III)-OOH weakens both the O-O and the O-H bonds of the oxidant. Heterolytic cleavage of the O-O bond leads to N-hydroxylation of ArNH(-) via S(N)2 mechanism, whereas cleavage of the O-H bond results in release of hydroperoxy radical. Thus, our proposed reaction offers a mechanistic explanation for previous observations that metabolism of aromatic amines could cause oxidative stress. The primary drivers for mutagenic potency of ArNH(2) are (i) binding affinity of ArNH(2) in the productive binding mode within the CYP1A2 substrate cavity, (ii) resonance stabilization of the anionic forms of ArNH(2), and (iii) exothermicity of proton-assisted heterolytic cleavage of N-O bonds of hydroxylamines and their bioconjugates. This leads to a strategy for designing mutagenicity free ArNH(2): Structural alterations in ArNH(2), which disrupt geometric compatibility with CYP1A2, hinder proton abstraction, or strongly destabilize the nitrenium ion, in this order of priority, prevent genotoxicity.

Boronic acids-a novel class of bacterial mutagen.

Boronic acids and their esters are important building blocks in organic syntheses including those for drug substances and for which, as far as it can be determined, there are no published reports of testing for genotoxicity. A number of boronic acids have now been tested in this laboratory using Salmonella typhimurium strains TA1535, TA1537, TA98 and TA100 and Escherichia coli strain WP2uvrA(pKM101). Twelve of the 13 structures presented here were found to be mutagenic. All the compounds except one were active only in TA100 and/or WP2uvrA(pKM101), did not require S9 activation and produced relatively weak responses, i.e. no more than seven times the concurrent solvent-control values at >1000µg/plate. The single exception was also weakly mutagenic for TA1537 in the presence of S9. Results with two compounds mutagenic for both TA100 and WP2uvrA(pKM101) showed no evidence of DNA-adduct formation detectable by (32)P-postlabelling. It appears that boronic acids represent a novel class of bacterial mutagen that may not act by direct covalent binding to DNA. However, their mechanism of action remains to be elucidated and it cannot yet be determined whether or not they present a real genotoxic hazard.

Discovery and optimization of a series of carbazole ureas as NPY5 antagonists for the treatment of obesity.

The hypothesis that antagonists of the neuropeptide Y5 receptor would provide safe and effective appetite suppressants for the treatment of obesity has prompted vigorous research to identify suitable compounds. We discovered a series of acylated aminocarbazole derivatives (e.g., 3a) that are potent and selective Y5 antagonists, representing interesting starting points but suffering from poor bioavailability and concerns about potential toxicity as a consequence of the embedded aminocarbazole fragment. It proved relatively easy to improve the drug metabolism and pharmacokinetic (DMPK) properties by variation of the side chain (as in 4a) but difficult to eliminate the aminocarbazole fragment. For compounds in this series to have the potential to be drugs, we believed that both the compound itself and the component aniline must be free of mutagenic activity. Parallel structure-activity relationship studies looking at the effects of ring substitution have proved that it is possible by incorporation of a 4-methyl substituent to produce carbazole ureas with potent Y5 activity, comprised of carbazole anilines that in themselves are devoid of mutagenic activity in the Ames test. Compound 4o (also known as NPY5RA-972) is highly selective with respect to Y1, Y2, and Y4 receptors (and also to a diverse range of unrelated receptors and enzymes), with an excellent DMPK profile including central nervous system penetration. NPY5RA-972 (4o) is a highly potent Y5 antagonist in vivo but does not block neuropeptide Y-induced feeding nor does it reduce feeding in rats, suggesting that the Y5 receptor alone has no significant role in feeding in these models.

Rationally designing safer anilines: the challenging case of 4-aminobiphenyls.

We describe how we have been able to design 4-aminobiphenyls that are nonmutagenic (inactive in the Ames test). No such 4-aminobiphenyls were known to us, but insights provided by quantum mechanical calculations have permitted us to design and synthesize some examples. Importantly, the quantum mechanical calculations could be combined with predictions of other properties of the compounds that contained the 4-aminobiphenyls so that these remained druglike. Having found compounds that are not active, the calculations can provide insight into which factors (electronic and conformational in this case) are important. The calculations provided SAR-like information that was able guide the design of further examples of 4-aminobiphenyls that are not active in the Ames test.