Ensemble Modeling of Substrate Binding to Cytochromes P450:  Analysis of Catalytic Differences between CYP1A Orthologs,

Jahnavi C. Prasad,§ Jared V. Goldstone, Carlos J. Camacho, Sandor Vajda,*§ and John J. Stegeman*
Department of Biomedical Engineering, Boston University, Boston, Massachusetts 02215, Biology Department, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts 02543, and Department of Computational Biology, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
Biochemistry, 2007, 46 (10), pp 2640–2654
DOI: 10.1021/bi062320m
Publication Date (Web): February 15, 2007
Copyright © 2007 American Chemical Society

Abstract

Abstract Image

A novel application of modeling and docking approaches involving ensembles of homology models is used to understand structural bases underlying subtle catalytic differences between related cytochromes P450 (CYPs). Mammalian CYP1A1s and fish CYP1As are orthologous enzymes with similar substrate preferences. With some substrates (3,3‘,4,4‘-tetrachlorobiphenyl, TCB) oxidation rates differ by orders of magnitude, while others (e.g., benzo[a]pyrene; B[a]P) are oxidized at similar rates but with somewhat differing regiospecificity. These two environmental chemical substrates (TCB and B[a]P) as well as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) were docked to multiple models of rat, human, scup, and/or killifish CYP1As, based on multiple templates, retaining multiple poses from each model, giving ensembles of docked poses for each species. With TCB, more poses were observed closer to the heme in ensembles of rat or human CYP1A1 than of killifish CYP1A. Analysis of interacting residues suggested that differences in TCB pose distributions are due primarily to Leu387 and Val230 in killifish CYP1A. In silico mutations L387V and V230G enabled TCB to dock closer to the heme in killifish CYP1A. Mutating additional interacting residues (Ala127, Thr233, Asn317, and Tyr386) of killifish CYP1A to the corresponding residues of human CYP1A1 resulted in TCB pose distributions nearly identical with those of human CYP1A1. Docking of TCDD to sets of consensus models of killifish, rat, and human CYP1As showed species differences similar to those with TCB, but with further structural constraints possibly contributing to slower oxidation of TCDD. Docking B[a]P to sets of consensus models of the human and fish CYP1As yielded frequencies of substrate orientations correlating with known regiospecificities for metabolism of B[a]P by these enzymes. The results demonstrate the utility of this ensemble modeling method, which can account for uncertainty inherent in homology modeling and docking by producing statistical distributions of ligand positions.

Citing Articles

View all 14 citing articles

Citation data is made available by participants in CrossRef's Cited-by Linking service. For a more comprehensive list of citations to this article, users are encouraged to perform a search in SciFinder.

This article has been cited by 5 ACS Journal articles (5 most recent appear below).

  • Cover Image

    Functional Screening of Cytochrome P450 Activity and Uncoupling by Capillary Electrophoresis

    James Harskamp, Philip Britz-McKibbin, and Joanna Y. Wilson
    Analytical Chemistry2012 84 (2), 862-866
    • Functional Screening of Cytochrome P450 Activity and Uncoupling by Capillary Electrophoresis

      James Harskamp, Philip Britz-McKibbin, and Joanna Y. Wilson
      Analytical Chemistry2012 84 (2), 862-866

      Cytochrome P450s (CYPs) are functionally diverse monooxygenases responsible for oxidation of endogenous and xenobiotic compounds. The function of nonmammalian CYPs are largely unknown and tools for characterization limited. CYPs critical for xenobiotic ...

  • Cover Image

    Backbone Flexibility Controls the Activity and Specificity of a Protein−Protein Interface: Specificity in Snake Venom Metalloproteases

    Hannes G. Wallnoefer, Torsten Lingott, José María Gutiérrez, Irmgard Merfort and Klaus R. Liedl
    Journal of the American Chemical Society2010 132 (30), 10330-10337
    • Backbone Flexibility Controls the Activity and Specificity of a Protein−Protein Interface: Specificity in Snake Venom Metalloproteases

      Hannes G. Wallnoefer, Torsten Lingott, José María Gutiérrez, Irmgard Merfort and Klaus R. Liedl
      Journal of the American Chemical Society2010 132 (30), 10330-10337

      Protein−protein interfaces have crucial functions in many biological processes. The large interaction areas of such interfaces show complex interaction motifs. Even more challenging is the understanding of (multi)specificity in protein−protein binding. ...

  • Cover Image

    Discovery of New Inhibitors of d-Alanine:d-Alanine Ligase by Structure-Based Virtual Screening

    Andreja Kovač, Janez Konc, Blaž Vehar, Julieanne M. Bostock, Ian Chopra, Dušanka Janežič and Stanislav Gobec
    Journal of Medicinal Chemistry2008 51 (23), 7442-7448
    • Discovery of New Inhibitors of d-Alanine:d-Alanine Ligase by Structure-Based Virtual Screening

      Andreja Kovač, Janez Konc, Blaž Vehar, Julieanne M. Bostock, Ian Chopra, Dušanka Janežič and Stanislav Gobec
      Journal of Medicinal Chemistry2008 51 (23), 7442-7448

      The terminal dipeptide, d-Ala-d-Ala, of the peptidoglycan precursor UDPMurNAc-pentapetide is a crucial building block involved in peptidoglycan cross-linking. It is synthesized in the bacterial cytoplasm by the enzyme d-alanine:d-alanine ligase (Ddl). ...

  • Cover Image

    Prediction of Metabolism by Cytochrome P450 2C9: Alignment and Docking Studies of a Validated Database of Substrates

    Matthew J. Sykes, Ross A. McKinnon and John O. Miners
    Journal of Medicinal Chemistry2008 51 (4), 780-791
    • Prediction of Metabolism by Cytochrome P450 2C9: Alignment and Docking Studies of a Validated Database of Substrates

      Matthew J. Sykes, Ross A. McKinnon and John O. Miners
      Journal of Medicinal Chemistry2008 51 (4), 780-791

      A validated database of 70 molecules known to undergo biotransformation by CYP2C9 was collated. The molecular alignment program ROCS was used with the query molecule flurbiprofen as a basis for predicting the correct active site orientation of the CYP2C9 ...

  • Cover Image

    The Search for Endogenous Activators of the Aryl Hydrocarbon Receptor

    Linh P. Nguyen and Christopher A. Bradfield
    Chemical Research in Toxicology2008 21 (1), 102-116
    • The Search for Endogenous Activators of the Aryl Hydrocarbon Receptor

      Linh P. Nguyen and Christopher A. Bradfield
      Chemical Research in Toxicology2008 21 (1), 102-116

      The primary design of this perspective is to describe the major ligand classes of the aryl hydrocarbon receptor (AHR). A grander objective is to provide models that may help define the physiological activator or “endogenous ligand” of the AHR. We present ...

Tools

SciFinder Links

SciFinder subscribers:  Click to sign in | Not a SciFinder subscriber? Learn more at www.cas.org

Explore by:


Accession Codes

History

  • Published In Issue March 13, 2007
  • Received November 8, 2006
    Revised Manuscript Received December 28, 2006

Recommend & Share

  • Share on ACS NetworkACS Network
  • Add to FacebookFacebook
  • Tweet ThisTweet This
  • Add to CiteULikeCiteULike
  • Add to NewsvineNewsvine
  • Digg ThisDigg This
  • Add to DeliciousDelicious

Related Content

Other ACS content by these authors: