Lena Ripa

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.

Increasing Selectivity of CC Chemokine Receptor 8 Antagonists by Engineering Nondesolvation Related Interactions with the Intended and Off-Target Binding Sites

The metabolic stability and selectivity of a series of CCR8 antagonists against binding to the hERG ion channel and cytochrome Cyp2D6 are studied by principal component analysis. It is demonstrated that an efficient way of increasing metabolic stability and selectivity of this series is to decrease compound lipophilicity by engineering nondesolvation related attractive interactions with CCR8, as rationalized by three-dimensional receptor models. Although such polar interactions led to increased compound selectivity, such a strategy could also jeopardize the DMPK profile of compounds. However, once increased potency is found, the lipophilicity can be readjusted by engineering hydrophobic substituents that fit to CCR8 but do not fit to hERG. Several such lipophilic fragments are identified by two-dimensional fragment-based QSAR analysis. Electrophysiological measurements and site-directed mutagenesis studies indicated that the repulsive interactions of these fragments with hERG are caused by steric hindrances with residue F656.

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

The metabolism of aromatic and heteroaromatic amines (ArNH₂) 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 ArNH₂, 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 ArNH₂, 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 ArNH₂ 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 ArNH₂. The mechanism-based reasoning suggests rational design strategies to avoid genotoxicity of ArNH₂ primarily by preventing N-hydroxylation of ArNH₂.

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.

Theoretical studies of the second step of the nitric oxide synthase reaction: Electron tunneling prevents uncoupling

Nitric oxide (NO·) is a messenger molecule with diverse physiological roles including host defense, neurotransmission and vascular function. The synthesis of NO· from l-arginine is catalyzed by NO-synthases and occurs in two steps through the intermediary Nω-hydroxy-l-arginine (NHA). In both steps the P450-like reaction cycle is coupled with the redox cycle of the cofactor tetrahydrobiopterin (H4B). The mechanism of the second step is studied by Density Functional Theory calculations to ascertain the canonical sequence of proton and electron transfer (PT and ET) events. The proposed mechanism is controlled by the interplay of two electron donors, H4B and NHA. Consistent with experimental data, the catalytic cycle proceeds through the ferric-hydroperoxide complex (Cpd 0) and the following aqua-ferriheme resting state, and involves interim partial oxidation of H4B. The mechanism starts with formation of Cpd 0 from the ferrous-dioxy reactant complex by PT from the C-ring heme propionate coupled with hole transfer to H4B through the highest occupied π-orbital of NHA as a bridge. This enables PT from NHA+· to the proximal oxygen leading to the shallow ferriheme-H2O2 oxidant. Subsequent Fenton-like peroxide bond cleavage triggered by ET from the NHA-derived iminoxy-radical leads to the protonated Cpd II diradicaloid singlet stabilized by spin delocalization in H4B, and the closed-shell coordination complex of HO- with iminoxy-cation. The complex is converted to the transient C-adduct, which releases intended products upon PT to the ferriheme-HO- complex coupled with ET to the H4B+·. Deferred ET from the substrate or undue ET from/to the cofactor leads to side products.