Common variants in LSP1, 2q35 and 8q24 and breast cancer risk for BRCA1 and BRCA2 mutation carriers.

Genome-wide association studies of breast cancer have identified multiple single nucleotide polymorphisms (SNPs) that are associated with increased breast cancer risks in the general population. In a previous study, we demonstrated that the minor alleles at three of these SNPs, in FGFR2, TNRC9 and MAP3K1, also confer increased risks of breast cancer for BRCA1 or BRCA2 mutation carriers. Three additional SNPs rs3817198 at LSP1, rs13387042 at 2q35 and rs13281615 at 8q24 have since been reported to be associated with breast cancer in the general population, and in this study we evaluated their association with breast cancer risk in 9442 BRCA1 and 5665 BRCA2 mutation carriers from 33 study centres. The minor allele of rs3817198 was associated with increased breast cancer risk only for BRCA2 mutation carriers [hazard ratio (HR) = 1.16, 95% CI: 1.07-1.25, P-trend = 2.8 x 10(-4)]. The best fit for the association of SNP rs13387042 at 2q35 with breast cancer risk was a dominant model for both BRCA1 and BRCA2 mutation carriers (BRCA1: HR = 1.14, 95% CI: 1.04-1.25, P = 0.0047; BRCA2: HR = 1.18 95% CI: 1.04-1.33, P = 0.0079). SNP rs13281615 at 8q24 was not associated with breast cancer for either BRCA1 or BRCA2 mutation carriers, but the estimated association for BRCA2 mutation carriers (per-allele HR = 1.06, 95% CI: 0.98-1.14) was consistent with odds ratio estimates derived from population-based case-control studies. The LSP1 and 2q35 SNPs appear to interact multiplicatively on breast cancer risk for BRCA2 mutation carriers. There was no evidence that the associations vary by mutation type depending on whether the mutated protein is predicted to be stable or not.

Compound I in heme thiolate enzymes: a comparative QM/MM study.

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.

SULT1E1 and ID2 genes as candidates for inherited predisposition to breast and ovarian cancer in Jewish women.

Mutations in currently known genes account for only a subset of breast/ovarian cancer risk families. Three loci (2p, 4q, 22q) seemingly harbor breast cancer susceptibility genes. To explore their putative role in Jewish women, 46 affected women representing 22 high risk families were genotyped with D2S2211, D4S392, D22S278 and D22S283 and two flanking markers for each locus, and mutational analysis of ID2 (Chromosome 2) and SULT1E1 (Chromosome 4) genes was carried out in seemingly linked families. No ID2 gene mutations were detected in 8 women from the 4 families seemingly linked to D2S2211, whereas a missense mutation (His224Gln) in one affected woman from a single family was detected among 9 women from the 4 families linked to D4S392. This mutation was not found among 153 high risk, 98 sporadic breast/ovarian cancer patients, or 97 healthy controls. The SULT1E1 gene may need to be further explored as candidate breast cancer gene.

In silico design of a mutant of cytochrome P450 containing selenocysteine.

A mutant of P450(cam), in which the cysteine ligand was replaced by selenocysteine, was designed theoretically using hybrid QM/MM (quantum mechanical/molecular mechanical) calculations. The calculations of the active species, Se-CpdI (selenocysteine-Compound I), of the mutant enzyme indicate that Se-Cpd I will be formed faster than the wild-type species and be consumed more slowly in C-H hydroxylation. As such, our calculations suggest that Se-Cpd I can be observed unlike the elusive species of the wild-type enzyme (Denisov, I. G.; Makris, T. M.; Sligar, S. G.; Schlichting, I. Chem. Rev. 2005, 105, 2253-2277). Spectral features of Se-Cpd I were calculated and may assist such eventual characterization. The observation of Se-Cpd I will resolve the major puzzle in the catalytic cycle of a key enzyme in nature.

Does substrate oxidation determine the regioselectivity of cyclohexene and propene oxidation by cytochrome p450?

DFT and QM/MM computations of allylic C-H hydroxylation versus C=C epoxidation in cyclohexene and propene by Compound I of P450cam demonstrate that the relative barriers for the oxidative processes themselves are not good predictors of the observed selectivity. However, a kinetic expression previously developed (Kozuch, S.; Shaik, S. J. Am. Chem. Soc. 2006, 128, 3355) for catalytic cycles under steady-state conditions, predicts, in accord with experiment, that propene will undergo exclusive C=C epoxidation, while cyclohexene will undergo both reactions with a small preference for epoxidation. The model expression for the effective barrier of the cycle forms a general basis for understanding and predicting the selectivity of P450 isozymes.

Principal active species of horseradish peroxidase, compound I: a hybrid quantum mechanical/molecular mechanical study.

The active species, Compound I, of horseradish peroxidase (HRP) has been investigated by quantum mechanical/molecular mechanical (QM/MM) calculations using 10 different QM regions. In accord with experimental data, the lowest doublet and quartet states are found to be virtually degenerate, with two unpaired electrons on the FeO moiety and one localized on the porphyrin in an a(2u)-dominant orbital with a minor, but nonnegligible, a(1u) component. The proximal ligand appears to be imidazole rather than imidazolate. The hydrogen-bonding network around the FeO moiety (i.e., Arg38 and His42) has significant influence on the axial bonds and the spin density distribution in the FeO moiety. Including this network in the QM region was found to be essential for reproducing the experimental Mössbauer parameters. The protein environment shapes most of the subtle features of Compound I of HRP.

Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism.

The stereospecific cytochrome P450-catalyzed hydroxylation of the C(5)-H((5-exo)) bond in camphor has been studied theoretically by a combined quantum mechanical/molecular mechanical (QM/MM) approach. Density functional theory is employed to treat the electronic structure of the active site (40-100 atoms), while the protein and solvent environment (ca. 24,000 atoms) is described by the CHARMM force field. The calculated energy profile of the hydrogen-abstraction oxygen-rebound mechanism indicates that the reaction takes place in two spin states (doublet and quartet), as has been suggested earlier on the basis of calculations on simpler models ("two-state reactivity"). While the reaction on the doublet potential energy surface is nonsynchronous, yet effectively concerted, the quartet pathway is truly stepwise, including formation of a distinct intermediate substrate radical and a hydroxo-iron complex. Comparative calculations in the gas phase demonstrate the effect of the protein environment on the geometry and relative stability of intermediates (in terms of spin states and redox electromers) through steric constraints and electronic polarization.

Radical clock substrates, their C-H hydroxylation mechanism by cytochrome P450, and other reactivity patterns: what does theory reveal about the clocks' behavior?

There is an ongoing and tantalizing controversy regarding the mechanism of a key process in nature, C-H hydroxylation, by the enzyme cytochrome P450 (Auclaire, K.; Hu, Z.; Little, D. M.; Ortiz de Montellano, P. R.; Groves, J. T. J. Am. Chem. Soc. 2002, 124, 6020-6027. Newcomb, M.; Aebisher, D.; Shen, R.; Esala, R.; Chandrasena, P.; Hollenberg, P. F.; Coon, M. J. J. Am. Chem. Soc. 2003, 125, 6064-6065). To definitely resolve this controversy, theory must first address the actual systems that have been used by experiment, and that generated the controversy. This is done in the present paper, which constitutes the first extensive theoretical study of such two experimental systems, trans-2-phenylmethyl-cyclopropane (1) and trans-2-phenyl-iso-propylcyclopropane (4). The theoretical study of these substrates reveals that the only low energy pathway for C-H hydroxylation is the two-state rebound mechanism described originally for methane hydroxylation (Ogliaro, F.; Harris, N.; Cohen, S.; Filatov, M.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 8977-8989). The paper shows that the scenario of a two-state rebound mechanism accommodates much of the experimental data. The computational results provide a good match to experimental results concerning the very different extents of rearrangement for 1 (20-30%) vs 4 (virtually none), lead to product isotope effect for the reaction of 1, in the direction of the experimental result, and predict as well the observed metabolic switching from methyl to phenyl hydroxylation, which occurs upon deuteration of the methyl group. Furthermore, the study reveals that an intimate ion pair species involving an alkyl carbocation derived from 4 gives no rearranged products, again in accord with experiment. This coherent match between theory and experiment cannot be merely accidental; it comes close to being aproof that the actual mechanism of C-H hydroxylation involves the two-state reactivity revealed by theory. Analysis of the rearrangement modes of the carbocations derived from 1 and 4 excludes the participation of free carbocations during the hydroxylation of these substrates. Finally, the mechanistic significance of product isotope effect (different isotope effects for the rearranged and unrearranged alcohol products) is analyzed. It is shown to be a sensitive probe of two-state reactivity; the size of the intrinsic product isotope effect and its direction reveal the structural differences of the hydrogen abstraction transition states in the low-spin vs high-spin reaction manifolds.

A predictive pattern of computed barriers for C-h hydroxylation by compound I of cytochrome p450.

The communication presents DFT calculations of 10 different C-H hydroxylation barriers by the active species of the enzyme cytochrome P450. The work demonstrates the existence of an excellent barrier-bond energy correlation. The so-obtained equation of the straight line is demonstrated to be useful for predicting barriers of related C-H activation processes, as well as for assessing barrier heights within the protein environment. This facility is demonstrated be estimating the barrier of camphor hydroxylation by P450cam.

Searching for the second oxidant in the catalytic cycle of cytochrome P450: a theoretical investigation of the iron(III)-hydroperoxo species and its epoxidation pathways.

Iron(III)-hydroperoxo, [Por(CysS)Fe(III)-OOH](-), a key species in the catalytic cycle of cytochrome P450, was recently identified by EPR/ENDOR spectroscopies (Davydov, R.; Makris, T. M.; Kofman, V.; Werst, D. E.; Sligar, S. G.; Hoffman, B. M. J. Am. Chem. Soc. 2001, 123, 1403-1415). It constitutes the last station of the preparative steps of the enzyme before oxidation of an organic compound and is implicated as the second oxidant capable of olefin epoxidation (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555-3560), in addition to the penultimate active species, Compound I (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; pp 3-48). In response, we present a density functional study of a model species and its ethylene epoxidation pathways. The study characterizes a variety of properties of iron(III)-hydroperoxo, such as the O-O bonding, the Fe-S bonding, Fe-O and Fe-S stretching frequencies, its electron attachment, and ionization energies. Wherever possible these properties are compared with those of Compound I. The proton affinities for protonation on the proximal and distal oxygen atoms of iron(III)-hydroperoxo, and the effect of the thiolate ligand thereof, are determined. In accordance with previous results (Harris, D. L.; Loew, G. H. J. Am. Chem. Soc. 1998, 120, 8941-8948), iron(III)-hydroperoxo is a strong base (as compared with water), and its distal protonation leads to a barrier-free formation of Compound I. The origins of this barrier-free process are discussed using a valence bond approach. It is shown that the presence of the thiolate is essential for this process, in line with the "push effect" deduced by experimentalists (Sono, M.; Roach, M. P.; Coulter, E. D.; Dawson, J. H. Chem. Rev. 1996, 96, 2841-2887). Finally, four epoxidation pathways of iron(III)-hydroxperoxo are located, in which the species transfers oxygen to ethylene either from the proximal or from the distal sites, in both concerted and stepwise manners. The barriers for the four mechanisms are 37-53 kcal/mol, in comparison with 14 kcal/mol for epoxidation by Compound I. It is therefore concluded that iron(III)-hydroperoxo, as such, cannot be a second oxidant, in line with its significant basicity and poor electron-accepting capability. Possible versions of a second oxidant are discussed.

The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations.

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.