Backbone dynamics for the wild type and a double H52R/T56W mutant of the vnd/NK-2 homeodomain from Drosophila melanogaster

Thermodynamic and structural determinants of differential Pdx1 binding to elements from the insulin and IAPP promoters

In adult mammals, the production of insulin and other peptide hormones, such as the islet amyloid polypeptide (IAPP), is limited to β-cells due to tissue-specific expression of a set of transcription factors, the best known of which is pancreatic duodenal homeobox protein 1 (Pdx1). Like many homeodomain transcription factors, Pdx1 binds to a core DNA recognition sequence containing the tetranucleotide 5'-TAAT-3'; its consensus recognition element is 5'-CTCTAAT(T/G)AG-3'. Currently, a complete thermodynamic profile of Pdx1 binding to near-consensus and native promoter sequences has not been established, obscuring the mechanism of target site selection by this critical transcription factor. Strikingly, while Pdx1 responsive elements in the human insulin promoter conform to the pentanucleotide 5'-CTAAT-3' sequence, the Pdx1 responsive elements in the human iapp promoter all contain a substitution to 5'-TTAAT-3'. The crystal structure of Pdx1 bound to the consensus nucleotide sequence does not explain how Pdx1 identifies this natural variation, if it does at all. Here we report a combination of isothermal calorimetric titrations, NMR spectroscopy, and extensive multi-microsecond molecular dynamics calculations of Pdx1 that define its interactions with a panel of natural promoter elements and consensus-derived sequences. Our results show a small preference of Pdx1 for a C base 5' relative to the core TAAT promoter element. Molecular mechanics calculations, corroborated by experimental NMR data, lead to a rational explanation for sequence discrimination at this position. Taken together, our results suggest a molecular mechanism for differential Pdx1 affinity to elements from the insulin and iapp promoter sequences.

Local folding and misfolding in the PBX homeodomain from a three-state analysis of CPMG relaxation dispersion NMR data

NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments represent a powerful approach for characterizing protein internal motions and for gaining insight into fundamental biological processes such as protein folding, catalysis, and allostery. In most cases, CPMG data are analyzed assuming that the protein exchanges between two different conformational states. Systems exchanging among more than two states are far more challenging to characterize by CPMG NMR. For example, in the case of three-state exchange in the fast time scale regime, it is difficult to uniquely connect the parameters extracted from CPMG analyses with the physical parameters of most interest, intercoversion rates, populations, and chemical shift differences for exchanging states. We have developed a grid search selection procedure that allows these physical parameters to be uniquely determined from CPMG data, based on additional information, which in this study comprises ligand-induced chemical shift perturbations. We applied this approach to the PBX homeodomain (PBX-HD), a three-helix protein with a C-terminal extension that folds into a fourth helix upon binding to DNA. We recently showed that the C-terminal extension transiently folds, even in the absence DNA, in a process that is likely tied to the cooperative binding of PBX-HD to DNA and other homeodomains. Using the grid search selection procedure, we found that PBX-HD undergoes exchange between three different conformational states, a major form in which the C-terminal extension is unfolded, the previously identified state in which the C-terminal extension forms a fourth helix, and an additional state in which the C-terminal extension is misfolded.

Structural and biophysical insights into the ligand-free Pitx2 homeodomain and a ring dermoid of the cornea inducing homeodomain mutant

The homeodomain-containing transcription factor Pitx2 (pituitary homeobox protein 2) is present in many developing embryonic tissues, including the heart. Its homeodomain is responsible for the recognition and binding to target DNA sequences and thus constitutes a major functional unit in the Pitx2 protein. Nuclear magnetic resonance techniques were employed to determine the solution structure of the native Pitx2 homeodomain and a R24H mutant that causes autosomal dominantly inherited ring dermoid of the cornea syndrome. The structures reveal that both isoforms possess the canonical homeodomain fold. However, the R24H mutation results in a 2-fold increase in DNA binding affinity and a 5 °C decrease in thermal stability, while changing the dynamic environment of the homeodomain only locally. When introduced into full-length Pitx2c, the mutation results in an only 25% loss of transactivation activity. Our data correlate well with clinical observations suggesting a milder deficiency for the R24H mutation compared to those of other Pitx2 homeodomain mutations.

Model-free analysis for large proteins at high magnetic field strengths

Protein backbone dynamics is often characterized using model-free analysis of three sets of (15)N relaxation data: longitudinal relaxation rate (R1), transverse relaxation rate (R2), and (15)N-{H} NOE values. Since the experimental data is limited, a simplified model-free spectral density function is often used that contains one Lorentzian describing overall rotational correlation but not one describing internal motion. The simplified spectral density function may be also used in estimating the overall rotational correlation time, by making the R2/R1 largely insensitive to internal motions, as well as used as one of the choices in the model selection protocol. However, such approximation may not be valid for analysis of relaxation data of large proteins recorded at high magnetic field strengths since the contribution to longitudinal relaxation from the Lorentzian describing the overall rotational diffusion of the molecule is comparably small relative to that describing internal motion. Here, we quantitatively estimate the errors introduced by the use of the simplified spectral density in model-free analysis for large proteins at high magnetic field strength.