Active site acidic residues and structural analysis of modelled human aromatase: A potential drug target for breast cancer

SummaryThis study sheds new light on the role of acidic residues present in the active site cavity of human aromatase. Eight acidic residues (E129, D222, E245, E302, D309, E379, D380 and D476) lining the cavity are identified and studied using comparative modeling, docking, molecular dynamics as well as statistical techniques. The structural environment of these acidic residues is studied to assess the stability of the corresponding carboxylate anions. Results indicate that the environment of the residues E245, E302 and D222 is most suitable for carboxylate ion formation in the uncomplexed form. However, the stability of D309, D222 and D476 anions is seen to increase on complexation to steroidal substrates. In particular, the interaction between D309 and T310, which assists proton transfer, is found to be formed following androgen/nor-androgen complexation. The residue D309 is found to be clamped in the presence of substrate which is not observed in the case of the other residues although they exhibit changes in properties following substrate binding. Information entropic analysis indicates that the residues D309, D222 and D476 have more conformational flexibility compared to E302 and E245 prior to substrate binding. Interaction similar to that between D476 and D309, which is expected to assist androgen aromatization, is proposed between E302 and E245. The inhibition of aromatase activity by 4-hydroxy androstenedione (formestane) is attributed to a critical hydrogen bond formation between the hydroxy moiety and T310/D309 as well as the large distance from D476. The results corroborate well with earlier site directed mutagenesis studies.

[1]  S. Neidle,et al.  A detailed molecular model for human aromatase , 1993, The Journal of Steroid Biochemistry and Molecular Biology.

[2]  B Honig,et al.  Extracting hydrophobic free energies from experimental data: relationship to protein folding and theoretical models. , 1991, Biochemistry.

[3]  E. Simpson,et al.  Functional domains of human aromatase cytochrome P450 characterized by linear alignment and site-directed mutagenesis. , 1993, Molecular endocrinology.

[4]  S. Rault,et al.  Evidence for new non-steroidal human aromatase inhibitors and comparison with equine aromatase inhibition for an understanding of the mammalian active site , 1998 .

[5]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

[6]  D. Zhou,et al.  Functional domains of aromatase cytochrome P450 inferred from comparative analyses of amino acid sequences and substantiated by site-directed mutagenesis experiments. , 1992, The Journal of biological chemistry.

[7]  Andrew C. R. Martin,et al.  Assessment of comparative modeling in CASP2 , 1997, Proteins.

[8]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[9]  E. Simpson,et al.  A three‐dimensional model of aromatase cytochrome P450 , 1995, Protein science : a publication of the Protein Society.

[10]  Sason Shaik,et al.  Mechanism of oxidation reactions catalyzed by cytochrome p450 enzymes. , 2004, Chemical reviews.

[11]  D T Jones,et al.  Protein secondary structure prediction based on position-specific scoring matrices. , 1999, Journal of molecular biology.

[12]  F. Gurd,et al.  pH-dependent processes in proteins. , 1985, CRC critical reviews in biochemistry.

[13]  Shiuan Chen,et al.  Structure-function studies of human aromatase , 1993, The Journal of Steroid Biochemistry and Molecular Biology.

[14]  J. Antosiewicz,et al.  Empirical relationships between protein structure and carboxyl pKa values in proteins , 2002, Proteins.

[15]  Chris de Graaf,et al.  Cytochrome p450 in silico: an integrative modeling approach. , 2005, Journal of medicinal chemistry.

[16]  Shridhar R. Gadre.  Information entropy and Thomas-Fermi theory , 1984 .

[17]  Anna L. Blobaum,et al.  P450 active site architecture and reversibility: inactivation of cytochromes P450 2B4 and 2B4 T302A by tert-butyl acetylenes. , 2005, Biochemistry.

[18]  Matthias Keil,et al.  Identifification of Substrate Channels and Protein Cavities , 1998 .

[19]  D E McRee,et al.  Microsomal cytochrome P450 2C5: comparison to microbial P450s and unique features. , 2000, Journal of inorganic biochemistry.

[20]  C David Stout,et al.  Structure of mammalian cytochrome P450 2C5 complexed with diclofenac at 2.1 A resolution: evidence for an induced fit model of substrate binding. , 2003, Biochemistry.

[21]  T. Poulos,et al.  A site-directed mutagenesis study of human placental aromatase. , 1992, The Journal of biological chemistry.

[22]  E. Cavalieri,et al.  Genotoxic metabolites of estradiol in breast: potential mechanism of estradiol induced carcinogenesis , 2003, The Journal of Steroid Biochemistry and Molecular Biology.

[23]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[24]  D. Eisenberg,et al.  Assessment of protein models with three-dimensional profiles , 1992, Nature.

[25]  G. N. Sastry,et al.  Aromatase inhibitors: a new paradigm in breast cancer treatment. , 2004, Current medicinal chemistry. Anti-cancer agents.

[26]  Shiuan Chen,et al.  Kinetic properties of aromatase mutants Pro308Phe, Asp309Asn, and Asp309Ala and their interactions with aromatase inhibitors , 1992, The Journal of Steroid Biochemistry and Molecular Biology.

[27]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[28]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[29]  G. Séralini,et al.  Study of substrate specificity of human aromatase by site directed mutagenesis. , 2002, European journal of biochemistry.

[30]  G Narahari Sastry,et al.  Structural and conformational changes concomitant with the E1-E2 transition in H(+)K(+)-ATPase: a comparative protein modeling study. , 2004, Biochemical and biophysical research communications.

[31]  M. Sternberg,et al.  Enhanced genome annotation using structural profiles in the program 3D-PSSM. , 2000, Journal of molecular biology.

[32]  A. Cavalli,et al.  Linking CoMFA and protein homology models of enzyme-inhibitor interactions: an application to non-steroidal aromatase inhibitors. , 2000, Bioorganic & medicinal chemistry.

[33]  J. Thijssen.  Local biosynthesis and metabolism of oestrogens in the human breast. , 2004, Maturitas.

[34]  M. Akhtar,et al.  The impact of aromatase mechanism on other P450s , 1997, The Journal of Steroid Biochemistry and Molecular Biology.

[35]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[36]  K. Korzekwa,et al.  Evaluation of the mechanism of aromatase cytochrome P450. A site-directed mutagenesis study. , 2001, European journal of biochemistry.

[37]  R. Fine,et al.  A CFF 91-based Continuum Solvation Model: Solvation Free Energies of Small Organic Molecules and Conformations of the Alanine Dipeptide in Solution , 1994 .

[38]  S. Pochapsky,et al.  Detection of a high-barrier conformational change in the active site of cytochrome P450cam upon binding of putidaredoxin. , 2005, Journal of the American Chemical Society.

[39]  Pieter F. W. Stouten,et al.  A molecular mechanics/grid method for evaluation of ligand–receptor interactions , 1995, J. Comput. Chem..

[40]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[41]  C. Laughton,et al.  Binding characteristics of seven inhibitors of human aromatase: a site-directed mutagenesis study. , 1996, Cancer research.