Genomic coordinates of predicted binding pocket I would like to know the genomic coordinates of a proteins ligand-binding pocket.  For example, on this PDB page you can see that there is xray crystallography data for how CYP2C9 interacts with two ligands. I would like to know which amino acids form the space on the protein that binds the ligands, then translate (de-translate?) that back to the nucleotide sequence. Is anyone familiar with either a database that has this information or how to extract the ligand interaction data from the PDB files?

Thank you
It turns out there is a tool that does this beautifully.  I used the proteinToGenome function from the ensembldb package.

The code looks like this:

# Define the Ensemble database object for the 
edbX &lt;- filter(EnsDb.Hsapiens.v75, filter = ~ seq_name == "X") 
# Set the first and last amino acid you want to convert and the ensembl protein id
GAGE10_prt &lt;- IRanges(start = 5, end = 9, names = "ENSP00000385415")
# Perform the conversion.  GAGE10_prt is an object that has the coordinates
GAGE10_prt &lt;- proteinToGenome(GAGE10_prt, edbX)