Catalytic Mechanism of Cytochrome P450 for 5′-Hydroxylation of Nicotine: Fundamental Reaction Pathways and Stereoselectivity

Abstract

A series of computational methods were used to study how cytochrome P450 2A6 (CYP2A6) interacts with (S)-(−)-nicotine, demonstrating that the dominant molecular species of (S)-(−)-nicotine in CYP2A6 active site exists in the free base state (with two conformations, SRt and SRc), despite of the fact that the protonated state is dominant for the free ligand in solution. The computational results reveal that the dominant pathway of nicotine metabolism in CYP2A6 is through nicotine free base oxidation. Further, first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations were carried out to uncover the detailed reaction pathways for the CYP2A6-catalyzed nicotine 5′-hydroxlylation reaction. In the determined CYP2A6-(S)-(−)-nicotine binding structures, the oxygen of Compound I (Cpd I) can abstract a hydrogen from either the trans-5′- or the cis-5′-position of (S)-(−)-nicotine. CYP2A6-catalyzed (S)-(−)-nicotine 5′-hydroxylation consists of two reaction steps, i.e. the hydrogen transfer from the 5′-position of (S)-(−)-nicotine to the oxygen of Cpd I (the H-transfer step), followed by the recombination of the (S)-(−)-nicotine moiety with the iron-bound hydroxyl group to generate the 5′-hydroxynicotine product (the O-rebound step). The H-transfer step is rate-determining. The 5′-hydroxylation proceeds mainly with the stereoselective loss of the trans-5′-hydrogen, i.e. the 5′-hydrogen trans to the pyridine ring. The calculated overall stereoselectivity of ~97% favoring the trans-5′-hydroxylation is close to the observed stereoselectivity of 89~94%. This is the first time to demonstrate that a CYP substrate exists dominantly in one protonation state (cationic species) in solution, but uses its less-favorable protonation state (neutral free base) to perform the enzymatic reaction.

Recent density functional theory model calculations of drug metabolism by cytochrome P450

“…ab initio, DFT, or semiempirical method) treatment of atoms involved in the active site with a MM description of the native protein environment [135][136][137]. QM/MM approach is becoming increasingly popular for the study of biochemical reactions [138][139][140].…”

Mechanistic insights into the catalytic reaction of plant allene oxide synthase (pAOS) via QM and QM/MM calculations

“…The mechanism of P450 enzymes, with oxidation mechanisms and inhibition has received particular attention [30][31][32][33][34][35][36]. Plant AOS and related atypical P450s in the prostaglandin H 2 to prostacyclin pathway have also been studied extensively.…”

Tautomeric transformations and reactivity of isoindole and sila-indole: A computational study

“…[28][29][30][31] The study of solvent e®ect on the stability of isomers indicates, (I)-isomer is still the most stable one among all the isomers. On the other hand, ÁE 1 values decrease and ÁE 2 values increase, as polarity of solvent increase Among the transition states, TSb state changes considerably with polarity of solvent.…”