Dynamics and flexibility of human aromatase probed by FTIR and time resolved fluorescence spectroscopy

The thin line between monooxygenases and peroxygenases. P450s, UPOs, MMOs, and LPMOs: A brick to bridge fields of expertise

Many scientific fields, although driven by similar purposes and dealing with similar technologies, often appear so isolated and far from each other that even the vocabularies to describe the very same phenomenon might differ. Concerning the vast field of biocatalysis, a special role is played by those redox enzymes that employ oxygen-based chemistry to unlock transformations otherwise possible only with metal-based catalysts. As such, greener chemical synthesis methods and environmentally-driven biotechnological approaches were enabled over the last decades by the use of several enzymes and ultimately resulted in the first industrial applications. Among what can be called today the environmental biorefinery sector, biomass transformation, greenhouse gas reduction, bio-gas/fuels production, bioremediation, as well as bulk or fine chemicals and even pharmaceuticals manufacturing are all examples of fields in which successful prototypes have been demonstrated employing redox enzymes. In this review we decided to focus on the most prominent enzymes (MMOs, LPMO, P450 and UPO) capable of overcoming the ∼100 kcal mol-1 barrier of inactivated CH bonds for the oxyfunctionalization of organic compounds. Harnessing the enormous potential that lies within these enzymes is of extreme value to develop sustainable industrial schemes and it is still deeply coveted by many within the aforementioned fields of application. Hence, the ambitious scope of this account is to bridge the current cutting-edge knowledge gathered upon each enzyme. By creating a broad comparison, scientists belonging to the different fields may find inspiration and might overcome obstacles already solved by the others. This work is organised in three major parts: a first section will be serving as an introduction to each one of the enzymes regarding their structural and activity diversity, whereas a second one will be encompassing the mechanistic aspects of their catalysis. In this regard, the machineries that lead to analogous catalytic outcomes are depicted, highlighting the major differences and similarities. Finally, a third section will be focusing on the elements that allow the oxyfunctionalization chemistry to occur by delivering redox equivalents to the enzyme by the action of diverse redox partners. Redox partners are often overlooked in comparison to the catalytic counterparts, yet they represent fundamental elements to better understand and further develop practical applications based on mono- and peroxygenases.

Novel D-modified heterocyclic androstane derivatives as potential anticancer agents: Synthesis, characterization, in vitro and in silico studies

Cancer remains a major health concern worldwide. The most frequently diagnosed types of cancer are caused by abnormal production or action of steroid hormones. In the present study, the synthesis and structural characterization of new heterocyclic androstane derivatives with D-homo lactone, 17α-(pyridine-2''-ylmethyl) or 17(E)-(pyridine-2''-ylmethylidene) moiety are presented. All compounds were evaluated for their anti-proliferative activity against HeLa cervical cancer cell line and non-cancerous kidney MDCK cells, where A-homo lactam compound 9A showed the greatest selectivity. Based on in vitro binding assays, N-formyl lactam compound 18 appeared to be the strong and isoform-selective ligand for ERα, while compound 9A displayed binding affinity for the GR-LBD, but also inhibited aldo-keto reductase 1C4 enzyme. Out of four selected compounds, methylpyrazolo derivative 13 showed potential for aromatase binding, while in silico studies provided insight into experimentally confirmed protein-ligand interactions.

Design, synthesis, biological activity evaluation and structure-activity relationships of new steroidal aromatase inhibitors. The case of C-ring and 7β substituted steroids

In this work, new steroidal aromatase inhibitors (AIs) were designed, synthesized, and tested. In one approach, C-ring substituted steroids namely those functionalized at C-11 position with an α or β hydroxyl group or with a carbonyl group as well as C-9/C-11 steroidal olefins and epoxides were studied. It was found that the carbonyl group at C-11 is more beneficial for aromatase inhibition than the hydroxyl group, and that the C-ring epoxides were more potent than the C-ring olefins, leading to the discovery of a very strong AI, compound 7, with an IC₅₀ of 0.011 μM, better than Exemestane, the steroidal AI in clinical use, which presents an IC₅₀ of 0.050 μM. In another approach, we explored the biological activity of A-ring C-1/C-2 steroidal olefins and epoxides in relation to aromatase inhibition and compared it with the biological activity of C-ring C-9/C-11 steroidal olefins and epoxides. On the contrary to what was observed for the C-ring olefins and epoxides, the A-ring epoxides were less potent than A-ring olefins. Finally, the effect of 7β-methyl substitution on aromatase inhibition was compared with 7α-methyl substitution, showing that 7β-methyl is better than 7α-methyl substitution. Molecular modelling studies showed that the 7β-methyl on C-7 seems to protrude into the opening to the access channel of aromatase in comparison to the 7α-methyl. This comparison led to find the best steroidal AI (12a) of this work with IC₅₀ of 0.0058 μM. Compound 12a showed higher aromatase inhibition capacity than two of the three AIs currently in clinical use.

Machine Learning-Assisted Prediction of the Biological Activity of Aromatase Inhibitors and Data Mining to Explore Similar Compounds

Designing molecules for drugs has been a hot topic for many decades. However, it is hard and expensive to find a new molecule. Thus, the cost of the final drug is also increased. Machine learning can provide the fastest way to predict the biological activity of druglike molecules. In the present work, machine learning models are trained for the prediction of the biological activity of aromatase inhibitors. Data was collected from the literature. Molecular descriptors are calculated to be used as independent features for model training. The results showed that the R ² values for linear regression, random forest regression, gradient boosting regression, and bagging regression are 0.58, 0.84, 0.77, and 0.80, respectively. Using these models, it is possible to predict the activity of new molecules in a short period of time and at a reasonable cost. Furthermore, Tanimoto similarity is used for similarity analysis, as well as a chemical database is mined to search for similar molecules. Nonetheless, this study provides a framework for repurposing other effective drug molecules to prevent cancer.

Depicting the proton relay network in human aromatase: New insights into the role of the alcohol‐acid pair

Human aromatase is the cytochrome P450 catalyzing the conversion of androgens into estrogens in a three steps reaction essential to maintain steroid hormones balance. Here we report the capture and spectroscopic characterization of its compound I (Cpd I), the main reactive species in cytochromes P450. The typical spectroscopic transitions indicating the formation of Cpd I are detected within 0.8 s when mixing aromatase with meta‐chloroperoxybenzoic acid. The estrogen product is obtained from the same reaction mixture, demonstrating the involvement of Cpd I in aromatization reaction. Site‐directed mutagenesis is applied to the acid‐alcohol pair D309 and T310 and to R192, predicted to be part of the proton relay network. Mutants D309N and R192Q do not lead to Cpd I with an associated loss of activity, confirming that these residues are involved in proton delivery for Cpd I generation. Cpd I is captured for T310A mutant and shows 2.9‐ and 4.4‐fold faster rates of formation and decay, respectively, compared to wild‐type (WT). However, its activity is lower than the WT and a larger amount of H₂O₂ is produced during catalysis, indicating that T310 has an essential role in proton gating for generation of Cpd 0 and Cpd I and for their stabilization. The data provide new evidences on the role of threonine belonging to the conserved “acid‐alcohol” pair and known to be crucial for oxygen activation in cytochromes P450.

Molecular Lego of Human Cytochrome P450: The Key Role of Heme Domain Flexibility for the Activity of the Chimeric Proteins

The cytochrome P450 superfamily are heme-thiolate enzymes able to carry out monooxygenase reactions. Several studies have demonstrated the feasibility of using a soluble bacterial reductase from Bacillus megaterium, BMR, as an artificial electron transfer partner fused to the human P450 domain in a single polypeptide chain in an approach known as ‘molecular Lego’. The 3A4-BMR chimera has been deeply characterized biochemically for its activity, coupling efficiency, and flexibility by many different biophysical techniques leading to the conclusion that an extension of five glycines in the loop that connects the two domains improves all the catalytic parameters due to improved flexibility of the system. In this work, we extend the characterization of 3A4-BMR chimeras using differential scanning calorimetry to evaluate stabilizing role of BMR. We apply the ‘molecular Lego’ approach also to CYP19A1 (aromatase) and the data show that the activity of the chimeras is very low (<0.003 min−1) for all the constructs tested with a different linker loop length: ARO-BMR, ARO-BMR-3GLY, and ARO-BMR-5GLY. Nevertheless, the fusion to BMR shows a remarkable effect on thermal stability studied by differential scanning calorimetry as indicated by the increase in Tonset by 10 °C and the presence of a cooperative unfolding process driven by the BMR protein domain. Previously characterized 3A4-BMR constructs show the same behavior of ARO-BMR constructs in terms of thermal stabilization but a higher activity as a function of the loop length. A comparison of the ARO-BMR system to 3A4-BMR indicates that the design of each P450-BMR chimera should be carefully evaluated not only in terms of electron transfer, but also for the biophysical constraints that cannot always be overcome by chimerization.

Assessment of Five Pesticides as Endocrine-Disrupting Chemicals: Effects on Estrogen Receptors and Aromatase

Pesticides are widely applied all over the world, and pesticide exposure can induce different biological effects posing a possible threat to human health. Due to their effects on the endocrine system, some pesticides are classified as endocrine disruptors. The aim of the study is to assess the interference of five pesticides on estrogen biosynthesis and estrogen signaling. Three neonicotinoid insecticides (Acetamiprid, Clothianidin, and Thiamethoxam), a carbamate insecticide (Methiocarb) and a herbicide (Oxadiazon) were tested. The effect of pesticides on estrogen biosynthesis was studied through an ELISA assay using a recombinant form of human aromatase, the enzyme that catalyzes the transformation of androgens to estrogens. Moreover, the effect of pesticides on estrogen signaling was assessed using a gene reporter assay on MELN cells, which measures estrogen receptor-mediated estrogenic activity. The results of the ELISA assay showed that the pesticides did not alter aromatase activity (no interference with estrogen biosynthesis), while the results of the gene reporter assay showed that only Methiocarb was able to alter estrogen signaling at high doses. The estrogenic activity of Methiocarb, expressed as 17β-estradiol equivalency factor (EEF), was equal to 8.0 × 10⁻⁸. In conclusion, this study suggested that Methiocarb should be considered a potential endocrine disruptor.

Mechanisms of allosteric and mixed mode aromatase inhibitors

Aromatase (CYP19) catalyzes the last biosynthetic step of estrogens in mammals and is a primary drug target for hormone-related breast cancer. However, treatment with aromatase inhibitors is often associated with adverse effects and drug resistance. In this study, we used virtual screening targeting a predicted cytochrome P450 reductase binding site on aromatase to discover four novel non-steroidal aromatase inhibitors. The inhibitors have potencies comparable to the noncompetitive tamoxifen metabolite, endoxifen. Our two most potent inhibitors, AR11 and AR13, exhibit both mixed-type and competitive-type inhibition. The cytochrome P450 reductase-CYP19 coupling interface likely acts as a transient binding site. Our modeling shows that our inhibitors bind better at different sites near the catalytic site. Our results predict the location of multiple ligand binding sites on aromatase. The combination of modeling and experimental results supports the important role of the reductase binding interface as a low affinity, promiscuous ligand binding site. Our new inhibitors may be useful as alternative chemical scaffolds that may show different adverse effects profiles than current clinically used aromatase inhibitors.

Polymorphism on human aromatase affects protein dynamics and substrate binding: spectroscopic evidence

Human aromatase is a member of the cytochrome P450 superfamily, involved in steroid hormones biosynthesis. In particular, it converts androgen into estrogens being therefore responsible for the correct sex steroids balance. Due to its capacity in producing estrogens it has also been considered as a promising target for breast cancer therapy. Two single-nucleotide polymorphisms (R264C and R264H) have been shown to alter aromatase activity and they have been associated to an increased or decreased risk for estrogen-dependent pathologies. Here, the effect of these mutations on the protein dynamics is investigated by UV/FTIR and time resolved fluorescence spectroscopy. H/D exchange rates were measured by FTIR for the three proteins in the ligand-free, substrate- and inhibitor-bound forms and the data indicate that the wild-type enzyme undergoes a conformational change leading to a more compact tertiary structure upon substrate or inhibitor binding. Indeed, the H/D exchange rates are decreased when a ligand is present. In the variants, the exchange rates in the ligand-free and -bound forms are similar, indicating that a structural change is lacking, despite the single amino acid substitution is located in the peripheral shell of the protein molecule. Moreover, the fluorescence lifetimes data show that the quenching effect on tryptophan-224 observed upon ligand binding in the wild-type, is absent in both variants. Since this residue is located in the catalytic pocket, these findings suggest that substrate entrance and/or retention in the active site is partially compromised in both mutants. A contact network analysis demonstrates that the protein structure is organized in two main clusters, whose connectivity is altered by ligand binding, especially in correspondence of helix-G, where the amino acid substitutions occur. Our findings demonstrate that SNPs resulting in mutations on aromatase surface modify the protein flexibility that is required for substrate binding and catalysis. The cluster analysis provides a rationale for such effect, suggesting helix G as a possible target for aromatase inhibition.

Molecular and Structural Evolution of Cytochrome P450 Aromatase

Aromatase is the cytochrome P450 enzyme converting androgens into estrogen in the last phase of steroidogenesis. As estrogens are crucial in reproductive biology, aromatase is found in vertebrates and the invertebrates of the genus Branchiostoma, where it carries out the aromatization reaction of the A-ring of androgens that produces estrogens. Here, we investigate the molecular evolution of this unique and highly substrate-selective enzyme by means of structural, sequence alignment, and homology modeling, shedding light on its key role in species conservation. The alignments led to the identification of a core structure that, together with key and unique amino acids located in the active site and the substrate recognition sites, has been well conserved during evolution. Structural analysis shows what their roles are and the reason why they have been preserved. Moreover, the residues involved in the interaction with the redox partner and some phosphorylation sites appeared late during evolution. These data reveal how highly substrate-selective cytochrome P450 has evolved, indicating that the driving forces for evolution have been the optimization of the interaction with the redox partner and the introduction of phosphorylation sites that give the possibility of modulating its activity in a rapid way.

Molecular Basis for Endocrine Disruption by Pesticides Targeting Aromatase and Estrogen Receptor

The intensive use of pesticides has led to their increasing presence in water, soil, and agricultural products. Mounting evidence indicates that some pesticides may be endocrine disrupting chemicals (EDCs), being therefore harmful for the human health and the environment. In this study, three pesticides, glyphosate, thiacloprid, and imidacloprid, were tested for their ability to interfere with estrogen biosynthesis and/or signaling, to evaluate their potential action as EDCs. Among the tested compounds, only glyphosate inhibited aromatase activity (up to 30%) via a non-competitive inhibition or a mixed inhibition mechanism depending on the concentration applied. Then, the ability of the three pesticides to induce an estrogenic activity was tested in MELN cells. When compared to 17β-estradiol, thiacloprid and imidacloprid induced an estrogenic activity at the highest concentrations tested with a relative potency of 5.4 × 10⁻¹⁰ and 3.7 × 10⁻⁹, respectively. Molecular dynamics and docking simulations predicted the potential binding sites and the binding mode of the three pesticides on the structure of the two key targets, providing a rational for their mechanism as EDCs. The results demonstrate that the three pesticides are potential EDCs as glyphosate acts as an aromatase inhibitor, whereas imidacloprid and thiacloprid can interfere with estrogen induced signaling.

Effector role of cytochrome P450 reductase for androstenedione binding to human aromatase

Cytochromes P450 constitute a large superfamily of monooxygenases involved in many metabolic pathways. Most of them are not self-sufficient and need a reductase protein to provide the electrons necessary for catalysis. It was shown that the redox partner plays a role in the modulation of the structure and function of some bacterial P450 enzymes. Here, the effect of NADPH-cytochrome reductase (CPR) on human aromatase (Aro) is studied for what concerns its role in substrate binding. Pre-steady-state kinetic experiments indicate that both the substrate binding rates and the percentage of spin shift detected for aromatase are increased when CPR is present. Moreover, aromatase binds the substrate through a conformational selection mechanism, suggesting a possible effector role of CPR. The thermodynamic parameters for the formation of the CPR-Aro complex were studied by isothermal titration calorimetry. The dissociation constant of the complex formation is 4.5 folds lower for substrate-free compared to the substrate-bound enzyme. The enthalpy change observed when the CPR-Aro complex forms in the absence of the substrate are higher than in its presence, indicating that more interactions are formed/broken in the former case. Taken together, our data confirm that CPR has a role in promoting aromatase conformation optimal for substrate binding.

Identification of endocrine disrupting chemicals acting on human aromatase

Human aromatase is the cytochrome P450 catalysing the conversion of androgens into estrogens playing a key role in the endocrine system. Due to this role, it is likely to be a target of the so-called endocrine disrupting chemicals, a series of compounds able to interfere with the hormone system with toxic effects. If on one side the toxicity of some compounds such as bisphenol A is well known, on the other side the toxic concentrations of such compounds as well as the effect of the many other molecules that are in contact with us in everyday life still need a deep investigation. The availability of biological assays able to detect the interaction of chemicals with key molecular targets of the endocrine system represents a possible solution to identify potential endocrine disrupting chemicals. Here the so-called alkali assay previously developed in our laboratory is applied to test the effect of different compounds on the activity of human aromatase. The assay is based on the detection of the alkali product that forms upon strong alkali treatment of the NADP⁺ released upon enzyme turnover. Here it is applied on human aromatase and validated using anastrozole and sildenafil as known aromatase inhibitors. Out of the small library of compounds tested, resveratrol and ketoconazole resulted to inhibit aromatase activity, while bisphenol A and nicotine were found to exert an inhibitory effect at relatively high concentrations (100μM), and other molecules such as lindane and four plasticizers did not show any significant effect. These data are confirmed by quantification of the product estrone in the same reaction mixtures through ELISA. Overall, the results show that the alkali assay is suitable to screen for molecules that interfere with aromatase activity. As a consequence it can also be applied to other molecular targets of EDCs that use NAD(P)H for catalysis in a high throughput format for the fast screening of many different compounds as endocrine disrupting chemicals. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Nanodiscs in Membrane Biochemistry and Biophysics

Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.

Impact of R264C and R264H polymorphisms in human aromatase function

The cytochrome P450 aromatase is involved in the last step of sex hormones biosynthesis by converting androgens into estrogens. The human enzyme is highly polymorphic and literature data correlate aromatase single nucleotide polymorphisms to the onset of pathologies such as breast cancer and neurodegenerative diseases. The aims of this study were i) to study the influence of the mutations R264C and R264H on the structure-function of the enzyme also upon phosphorylation by selected kinases and ii) to compare the activity of the variants to that of aromatase wild type in two different cell lines. Far-UV circular dichroism spectroscopy, thermal denaturation experiments and CO-binding assay showed that the two polymorphic variants are correctly folded. Steady-state kinetics experiments showed that rArom R264C and R264H exhibit a 1.5 and 3.4 folds lower catalytic efficiency, respectively, when compared to the wild type protein. Since R264 is part of the consensus motif of PKA and PKG1, phosphorylation experiments were performed to study the effect on aromatase function. Phosphorylation by PKA caused a decrease in activity by 36.2%, 49.3% and 27.9% in the wild type, R264C and R264H proteins respectively. Phosphorylation by PKG1 was also found to decrease the activity by 30.3%, 30.5% and 15.4% in the wild type, R264C and R264H proteins respectively. Experiments performed on the three full-length proteins expressed in human MCF-7 breast cancer cells and rat ST14A neuronal cells showed that, depending on the cell line used, the activity of the proteins is different, implicating different cellular mechanisms modulating aromatase activity. This work demonstrate that R264 polymorphism causes an intrinsic alteration of aromatase activity together with a different consensus for phosphorylation by different kinases, indicating that estrogen production can be different when such mutations are present. These findings are significant in understanding the onset and treatment of pathologies in which aromatase has been shown to be involved.

Effect of sildenafil on human aromatase activity: From in vitro structural analysis to catalysis and inhibition in cells

Aromatase catalyses the conversion of androgens into estrogens and is a well-known target for breast cancer therapy. As it has been suggested that its activity is affected by inhibitors of phosphodiesterase-5, this work investigates the potential interaction of sildenafil with aromatase. This is carried out both at molecular level through structural and kinetics assays applied to the purified enzyme, and at cellular level using neuronal and breast cancer cell lines. Sildenafil is found to bind to aromatase with a K_D of 0.58±0.05μM acting as a partial and mixed inhibitor with a maximal inhibition of 35±2%. Hyperfine sublevel correlation spectroscopy and docking studies show that sildenafil binds to the heme iron via its 6th axial water ligand. These results also provide information on the starting molecular scaffold for the development of new generations of drugs designed to inhibit aromatase as well as phosphodiesterase-5, a new emerging target for breast cancer therapy.

FTIR spectral signature of anticancer drugs. Can drug mode of action be identified?

Infrared spectroscopy has brought invaluable information about proteins and about the mechanism of action of enzymes. These achievements are difficult to transpose to living organisms as all biological molecules absorb in the mid infrared, with usually a high degree of overlap. Deciphering the contribution of each enzyme is therefore almost impossible. On the other hand, small changes in the infrared spectra of cells induced by environmental conditions or drugs may provide an accurate signature of the metabolic shift experienced by the cell as a response to a change in the growth medium. The present paper aims at reviewing the contribution of infrared spectroscopy to the description of small chemical changes that occur in cells when they are exposed to a drug. In particular, this review will focus on cancer cells and anti-cancer drugs. Results accumulated so far tend to demonstrate that infrared spectroscopy could be a very accurate descriptor of the mode of action of anticancer drugs. If confirmed, such a segmentation of potential drugs according to their "mode of action" will be invaluable for the discovery of new therapeutic molecules. This article is part of a Special Issue entitled: Physiological Enzymology and Protein Functions.

Evidence for an elevated aspartate pK(a) in the active site of human aromatase

Aromatase (CYP19A1), the enzyme that converts androgens to estrogens, is of significant mechanistic and therapeutic interest. Crystal structures and computational studies of this enzyme shed light on the critical role of Asp(309) in substrate binding and catalysis. These studies predicted an elevated pK(a) for Asp(309) and proposed that protonation of this residue was required for function. In this study, UV-visible absorption, circular dichroism, resonance Raman spectroscopy, and enzyme kinetics were used to study the impact of pH on aromatase structure and androstenedione binding. Spectroscopic studies demonstrate that androstenedione binding is pH-dependent, whereas, in contrast, the D309N mutant retains its ability to bind to androstenedione across the entire pH range studied. Neither pH nor mutation perturbed the secondary structure or heme environment. The origin of the observed pH dependence was further narrowed to the protonation equilibria of Asp(309) with a parallel set of spectroscopic studies using exemestane and anastrozole. Because exemestane interacts with Asp(309) based on its co-crystal structure with the enzyme, its binding is pH-dependent. Aromatase binding to anastrozole is pH-independent, consistent with the hypothesis that this ligand exploits a distinct set of interactions in the active site. In summary, we assign the apparent pK(a) of 8.2 observed for androstenedione binding to the side chain of Asp(309). To our knowledge, this work represents the first experimental assignment of a pK(a) value to a residue in a cytochrome P450. This value is in agreement with theoretical calculations (7.7-8.1) despite the reliance of the computational methods on the conformational snapshots provided by crystal structures.

Computational approaches elucidate the allosteric mechanism of human aromatase inhibition: a novel possible route to Small-molecule regulation of CYP450s activities?

Human aromatase (HA) is a P450 cytochrome (CYP) with an essential role in estrogen biosynthesis. Since more than 70% of breast cancers are positive for estrogenic receptor (ER), the reduction of estrogen physiological concentrations through HA inhibition is one of most important therapeutic strategies against this cancer type. Recently, experimental evidence showed that selected taxmoxifen metabolites, which are typically used as estrogen receptor modulators (SERMs), inhibit HA through an allosteric mechanism. In this work, we present a computational protocol to (i) characterize the structural framework and (ii) define the atomistic details of the determinants for the noncompetitive inhibition mechanism. Our calculations identify two putative binding sites able to efficiently bind all tamoxifen metabolites. Analysis of long-scale molecular dynamics simulations reveal that endoxifen, the most effective noncompetitive inhibitor, induces significant enzyme rigidity by binding in one of the possible peripheral sites. The consequence of this binding event is the suppression of one of the functional enzymatic collective motions associated with breathing of the substrate access channel. Moreover, an internal dynamics-based alignment of HA with six other human cytochromes shows that this collective motion is common to other members of the CYP450 protein family. On this basis, our findings may thus be of help for the development of new (pan)inhibitors for the therapeutic treatment of cancer, targeting and modulating the activity of HA and of estrogen receptor, and may also stimulate the development of new drug design strategies for chemoprevention and chemoprotection via allosteric inhibition of CYP450 proteins.