How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.     

Tunneling and classical paths for proton transfer in an enzyme reaction dominated by tunneling: oxidation of tryptamine by aromatic amine dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.     

Proton tunneling in aromatic amine dehydrogenase is driven by a short-range sub-picosecond promoting vibration: consistency of simulation and theory with experiment

Hydrogen transfer, an essential component of most biological reactions, is a quantum problem. However, the proposed role of compressive motion in promoting enzymatic H-transfer is contentious. Using molecular dynamics simulations and density functional theory (DFT) calculations, we show that, during proton tunneling in the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase (AADH), a sub-picosecond promoting vibration is inherent to the iminoquinone intermediate. We show by numerical modeling that this short-range vibration, with a frequency of approximately 165 cm-1, is consistent with "gating" motion in the hydrogen tunneling model of Kuznetsov and Ulstrup (Kuznetsov, A. M.; Ulstrup, J. Can. J. Chem. 1999, 77, 1085) in an enzymatic reaction with an observed protium/deuterium kinetic isotope effect that is not measurably temperature-dependent.     

Pressure effects on enzyme-catalyzed quantum tunneling events arise from protein-specific structural and dynamic changes

The rate and kinetic isotope effect (KIE) on proton transfer during the aromatic amine dehydrogenase-catalyzed reaction with phenylethylamine shows complex pressure and temperature dependences. We are able to rationalize these effects within an environmentally coupled tunneling model based on constant pressure molecular dynamics (MD) simulations. As pressure appears to act anisotropically on the enzyme, perturbation of the reaction coordinate (donor-acceptor compression) is, in this case, marginal. Therefore, while we have previously demonstrated that pressure and temperature dependences can be used to infer H-tunneling and the involvement of promoting vibrations, these effects should not be used in the absence of atomistic insight, as they can vary greatly for different enzymes. We show that a pressure-dependent KIE is not a definitive hallmark of quantum mechanical H-tunneling during an enzyme-catalyzed reaction and that pressure-independent KIEs cannot be used to exclude tunneling contributions or a role for promoting vibrations in the enzyme-catalyzed reaction. We conclude that coupling of MD calculations with experimental rate and KIE studies is required to provide atomistic understanding of pressure effects in enzyme-catalyzed reactions.     

Memetic reproduction and protolanguage evolution

An investigation into the evolution of protolanguages is described, exploring the tension between faithful reproduction of memes to ensure consistency and hence comprehensibility in a language and the need to generate new memes to encourage the diversity observed in evolution of new languages and dialects. Agents were given different genetic strategies for interpreting, memorising and imitating memes. The process of memetic reproduction over 50 generations was simulated with competing fitness criteria to optimise consistency of reproduction on the one hand, and diversity in the meme set on the other, while simulating variations in the complexity of seed meme sets. The results showed that selection among interpreting and memorising strategies is more critical for faithful reproduction of memes; in contrast, there was little difference between imitation strategies that selected memes either randomly, or based on frequency and similarity of the memes received in previous exchanges. More complex meme sets produced more variation favouring diversity, while simpler meme sets favoured consistency. Shorter memes were more stable than longer memes over time, but increasing the number of seed memes had less effect. The implications for the evolution of memetic reproduction and protolanguages are discussed.     

Hyperpolarization of Pyridyl Fentalogues by Signal Amplification By Reversible Exchange (SABRE)

Fentanyl, also known as ‘jackpot’, is a synthetic opiate that is 50–100 times more potent than morphine. Clandestine laboratories produce analogues of fentanyl, known as fentalogues to circumvent legislation regarding its production. Three pyridyl fentalogues were synthesized and then hyperpolarized by signal amplification by reversible exchange (SABRE) to appraise the forensic potential of the technique. A maximum enhancement of -168-fold at 1.4 T was recorded for the ortho pyridyl ¹H nuclei. Studies of the activation parameters for the three fentalogues revealed that the ratio of ligand loss trans to hydride and hydride loss in the complex [Ir(IMes)(L)₃(H)₂]⁺ (IMes=1,3-bis(2,4,6-trimethylphenyl)imidazole-2-ylidene) ranged from 0.52 to 1.83. The fentalogue possessing the ratio closest to unity produced the largest enhancement subsequent to performing SABRE at earth's magnetic field. It was possible to hyperpolarize a pyridyl fentalogue selectively from a matrix that consisted largely of heroin (97 : 3 heroin:fentalogue) to validate the use of SABRE as a forensic tool.     

Rational maps, monopoles and skyrmions

We discuss the similarities between BPS monopoles and skyrmions, and point to an underlying connection in terms of rational maps between Riemann spheres. This involves the introduction of a new ansatz for Skyrme fields. We use this to construct good approximations to several known skyrmions, including all the minimal energy configurations up to baryon number nine, and some new solutions such as a baryon number seventeen Skyrme field with the truncated icosahedron structure of a buckyball.

The new approach is also used to understand the low-lying vibrational modes of skyrmions, which are required for quantization. Along the way we discover an interesting Morse function on the space of rational maps which may be of use in understanding the Sen forms on the monopole moduli spaces.     

An NMR,ESR, ENDOR and TRIPLE resonance investigation of a cation radical formed from 1,2,4,5-tetramethoxybenzene

Several methods have been established for preparing cation radicals from 1,2,4,5-tetramethoxybenzene that allow highly resolved ESR spectra to be recorded. Precise values of the hyperfine coupling constants for the aromatic and methoxy protons have been obtained; the values are 0.2268±0.0004 and 0.0863±0.0002 mT, respectively, with dichloromethane as solvent. No temperature dependence is evident. TRIPLE resonance experiments showed that both coupling constants have the same sign. NMR experiments provided contact shift and line broadening measurements; these proved that both the above constants are positive and led to a value of 3.1 (±0.3)×10⁸M^?1 s^?1 at 23°C for the rate constant for electron exchange between the cation radical and the parent compound.     

Analysis of classical and quantum paths for deprotonation of methylamine by methylamine dehydrogenase.

The hydrogen-transfer reaction catalysed by methylamine dehydrogenase (MADH) with methylamine (MA) as substrate is a good model system for studies of proton tunnelling in enzyme reactions--an area of great current interest--for which atomistic simulations will be vital. Here, we present a detailed analysis of the key deprotonation step of the MADH/MA reaction and compare the results with experimental observations. Moreover, we compare this reaction with the related aromatic amine dehydrogenase (AADH) reaction with tryptamine, recently studied by us, and identify possible causes for the differences observed in the measured kinetic isotope effects (KIEs) of the two systems. We have used combined quantum mechanics/molecular mechanics (QM/MM) techniques in molecular dynamics simulations and variational transition state theory with multidimensional tunnelling calculations averaged over an ensemble of paths. The results reveal important mechanistic complexity. We calculate activation barriers and KIEs for the two possible proton transfers identified-to either of the carboxylate oxygen atoms of the catalytic base (Asp428beta)-and analyse the contributions of quantum effects. The activation barriers and tunnelling contributions for the two possible proton transfers are similar and lead to a phenomenological activation free energy of 16.5+/-0.9 kcal mol(-1) for transfer to either oxygen (PM3-CHARMM calculations applying PM3-SRP specific reaction parameters), in good agreement with the experimental value of 14.4 kcal mol(-1). In contrast, for the AADH system, transfer to the equivalent OD1 was found to be preferred. The structures of the enzyme complexes during reaction are analysed in detail. The hydrogen bond of Thr474beta(MADH)/Thr172beta(AADH) to the catalytic carboxylate group and the nonconserved active site residue Tyr471beta(MADH)/Phe169beta(AADH) are identified as important factors in determining the preferred oxygen acceptor. The protein environment has a significant effect on the reaction energetics and hence on tunnelling contributions and KIEs. These environmental effects, and the related clearly different preferences for the two carboxylate oxygen atoms (with different KIEs) in MADH/MA and AADH/tryptamine, are possible causes of the differences observed in the KIEs between these two important enzyme reactions.