El (speTri1), dolphin (turTru1), sheep (Ovis aries), panda (ailMel1), cat (CAT), armadillo (ARMA), tenrec (TENREC), African elephant (BROAD E1), opossum (monDom5), alpaca (vicPac1), sloth (choHof1), Tasmanian devil (DEVIL7.0) and platypus (Oana-5.0) were deduced from the genome sequences at the Ensembl database [22]. To test for positive selection of codons, we aligned the available PR LBD DNA sequences of all 35 mammalian species and performed an evolutionary codon analysis using the Selecton Server [23,24]. Analysis was performed based on the phylogenetic tree of mammalian evolution and choosing the mechanistic empirical combination model (MEC) [25] with standard settings. Significance was determined by comparing the MEC test to the M8a null model. Positive selection was regarded to be significant, when the AIC score of MEC was lower than the AIC score of M8a.Competitive Binding AssayBacterial lysates were diluted 1:25 with assay buffer. 500 ml of diluted lysate was incubated with 1 nM [3H]-progesterone (Amersham Biosciences) and increasing concentrations of progesterone or DHP (both Sigma) for 16 h at 4uC. Unbound steroids were removed by adding 200 ml of 2 NoritA and 0.2 dextran (both Serva) in assay buffer to all tubes except the total-activity control, which received 200 ml assay buffer. Following 10 min incubation on ice, all tubes were centrifuged for 10 min at 3,000 rpm at 4uC. 100 ml duplicates of each supernatant were transferred into scintillation vials, mixed with 3 ml RiaLuma scintillation fluid (J.T. Baker) and MedChemExpress TBHQ measured in a scintillation counter. Specific binding was calculated by subtracting non-specific binding ([3H]-progesterone completely displaced by adding 1 mM progesterone) from the amount of [3H]-progesterone bound to the receptor. SigmaPlot 8.0 (Systat software) was used to 18325633 plot binding data, with IC50 values calculated using the one site competition setting.Molecular ModelingAll modeling procedures were performed on a dual-processor 3.06 GHz LINUX workstation or a LINUX cluster (4863.2 GHz processors). The X-ray structure of human PR LBD with progesterone (PDB code 1a28) was used as a start structure [16]. Elephant PR LBD was constructed via homology modeling using the SYBYL 7.2 software and the B chain as a starting structure. Both receptor models (ligand extracted) were energy-minimized and used for docking of the energy-minimized BTZ-043 web ligands DHP and progesterone (automated docking and scoring with Surflex-Dock [17] SYBYL 7.2). The highest-scoring docking poses were minimized and used as input structures for subsequent molecular-dynamics (MD) simulations performed using SANDER (AMBER7). The simulations were performed under periodic-boundary conditions in a TIP3 water box (WATBOX216, minimal thickness ?20 A) with 2 fs time steps using weak-coupling temperature scaling and SHAKE. Ions were added for charge neutralization. The solvent was equilibrated under NTP conditions for 50 ps (isotropic pressure scaling) to reach constant box density. The production trajectories (constant volume) were run without any 11967625 restraints for 5220 ns. The development of the potential energy and of relative center-of-mass rms deviation of the Ca atoms from the start structure was monitored. Only the parts of the trajectories in which both values reached a steady state were subjected to further evaluation.Results The G722A Substitution Changes the Ligand Specificity of the PRIn order to identify the molecular background of the altered bi.El (speTri1), dolphin (turTru1), sheep (Ovis aries), panda (ailMel1), cat (CAT), armadillo (ARMA), tenrec (TENREC), African elephant (BROAD E1), opossum (monDom5), alpaca (vicPac1), sloth (choHof1), Tasmanian devil (DEVIL7.0) and platypus (Oana-5.0) were deduced from the genome sequences at the Ensembl database [22]. To test for positive selection of codons, we aligned the available PR LBD DNA sequences of all 35 mammalian species and performed an evolutionary codon analysis using the Selecton Server [23,24]. Analysis was performed based on the phylogenetic tree of mammalian evolution and choosing the mechanistic empirical combination model (MEC) [25] with standard settings. Significance was determined by comparing the MEC test to the M8a null model. Positive selection was regarded to be significant, when the AIC score of MEC was lower than the AIC score of M8a.Competitive Binding AssayBacterial lysates were diluted 1:25 with assay buffer. 500 ml of diluted lysate was incubated with 1 nM [3H]-progesterone (Amersham Biosciences) and increasing concentrations of progesterone or DHP (both Sigma) for 16 h at 4uC. Unbound steroids were removed by adding 200 ml of 2 NoritA and 0.2 dextran (both Serva) in assay buffer to all tubes except the total-activity control, which received 200 ml assay buffer. Following 10 min incubation on ice, all tubes were centrifuged for 10 min at 3,000 rpm at 4uC. 100 ml duplicates of each supernatant were transferred into scintillation vials, mixed with 3 ml RiaLuma scintillation fluid (J.T. Baker) and measured in a scintillation counter. Specific binding was calculated by subtracting non-specific binding ([3H]-progesterone completely displaced by adding 1 mM progesterone) from the amount of [3H]-progesterone bound to the receptor. SigmaPlot 8.0 (Systat software) was used to 18325633 plot binding data, with IC50 values calculated using the one site competition setting.Molecular ModelingAll modeling procedures were performed on a dual-processor 3.06 GHz LINUX workstation or a LINUX cluster (4863.2 GHz processors). The X-ray structure of human PR LBD with progesterone (PDB code 1a28) was used as a start structure [16]. Elephant PR LBD was constructed via homology modeling using the SYBYL 7.2 software and the B chain as a starting structure. Both receptor models (ligand extracted) were energy-minimized and used for docking of the energy-minimized ligands DHP and progesterone (automated docking and scoring with Surflex-Dock [17] SYBYL 7.2). The highest-scoring docking poses were minimized and used as input structures for subsequent molecular-dynamics (MD) simulations performed using SANDER (AMBER7). The simulations were performed under periodic-boundary conditions in a TIP3 water box (WATBOX216, minimal thickness ?20 A) with 2 fs time steps using weak-coupling temperature scaling and SHAKE. Ions were added for charge neutralization. The solvent was equilibrated under NTP conditions for 50 ps (isotropic pressure scaling) to reach constant box density. The production trajectories (constant volume) were run without any 11967625 restraints for 5220 ns. The development of the potential energy and of relative center-of-mass rms deviation of the Ca atoms from the start structure was monitored. Only the parts of the trajectories in which both values reached a steady state were subjected to further evaluation.Results The G722A Substitution Changes the Ligand Specificity of the PRIn order to identify the molecular background of the altered bi.