Cooperative ordering in homeodomain-DNA recognition: solution structure and dynamics of the MATa1 homeodomain

Structural characterization of PHOX2B and its DNA interaction shed light on the molecular basis of the +7Ala variant pathogenicity in CCHS

An expansion of poly-alanine up to +13 residues in the C-terminus of the transcription factor PHOX2B underlies the onset of congenital central hypoventilation syndrome (CCHS). Recent studies demonstrated that the alanine tract expansion influences PHOX2B folding and activity. Therefore, structural information on PHOX2B is an important target for obtaining clues to elucidate the insurgence of the alanine expansion-related syndrome and also for defining a viable therapy. Here we report by NMR spectroscopy the structural characterization of the homeodomain (HD) of PHOX2B and HD + C-terminus PHOX2B protein, free and in the presence of the target DNA. The obtained structural data are then exploited to obtain a structural model of the PHOX2B-DNA interaction. In addition, the variant +7Ala, responsible for one of the most frequent forms of the syndrome, was analysed, showing different conformational proprieties in solution and a strong propensity to aggregation. Our data suggest that the elongated poly-alanine tract would be related to disease onset through a loss-of-function mechanism. Overall, this study paves the way for the future rational design of therapeutic drugs, suggesting as a possible therapeutic route the use of specific anti-aggregating molecules capable of preventing variant aggregation and possibly restoring the DNA-binding activity of PHOX2B.

Chimeric Sex-Determining Chromosomal Regions and Dysregulation of Cell-Type Identity in a Sterile Zygosaccharomyces Allodiploid Yeast

Allodiploidization is a fundamental yet evolutionarily poorly characterized event, which impacts genome evolution and heredity, controlling organismal development and polyploid cell-types. In this study, we investigated the sex determination system in the allodiploid and sterile ATCC 42981 yeast, a member of the Zygosaccharomyces rouxii species complex, and used it to study how a chimeric mating-type gene repertoire contributes to hybrid reproductive isolation. We found that ATCC 42981 has 7 MAT-like (MTL) loci, 3 of which encode α-idiomorph and 4 encode a-idiomorph. Two phylogenetically divergent MAT expression loci were identified on different chromosomes, accounting for a hybrid a/α genotype. Furthermore, extra a-idimorph-encoding loci (termed MTLa copies 1 to 3) were recognized, which shared the same MATa1 ORFs but diverged for MATa2 genes. Each MAT expression locus was linked to a HML silent cassette, while the corresponding HMR loci were located on another chromosome. Two putative parental sex chromosome pairs contributed to this unusual genomic architecture: one came from an as-yet-undescribed taxon, which has the NCYC 3042 strain as a unique representative, while the other did not match any MAT-HML and HMR organizations previously described in Z. rouxii species. This chimeric rearrangement produces two copies of the HO gene, which encode for putatively functional endonucleases essential for mating-type switching. Although both a and α coding sequences, which are required to obtain a functional cell-type a1-α2 regulator, were present in the allodiploid ATCC 42981 genome, the transcriptional circuit, which regulates entry into meiosis in response to meiosis-inducing salt stress, appeared to be turned off. Furthermore, haploid and α-specific genes, such as MATα1 and HO, were observed to be actively transcribed and up-regulated under hypersaline stress. Overall, these evidences demonstrate that ATCC 42981 is unable to repress haploid α-specific genes and to activate meiosis in response to stress. We argue that sequence divergence within the chimeric a1-α2 heterodimer could be involved in the generation of negative epistasis, contributing to the allodiploid sterility and the dysregulation of cell identity.

Sex-determination system in the diploid yeast Zygosaccharomyces sapae

Sexual reproduction and breeding systems are driving forces for genetic diversity. The mating-type (MAT) locus represents a mutation and chromosome rearrangement hotspot in yeasts. Zygosaccharomyces rouxii complex yeasts are naturally faced with hostile low water activity (aw) environments and are characterized by gene copy number variation, genome instability, and aneuploidy/allodiploidy. Here, we investigated sex-determination system in Zygosaccharomyces sapae diploid strain ABT301(T), a member of the Z. rouxii complex. We cloned three divergent mating type-like (MTL) α-idiomorph sequences and designated them as ZsMTLα copies 1, 2, and 3. They encode homologs of Z. rouxii CBS 732(T) MATα2 (amino acid sequence identities spanning from 67.0 to 99.5%) and MATα1 (identity range 81.5-99.5%). ABT301(T) possesses two divergent HO genes encoding distinct endonucleases 100% and 92.3% identical to Z. rouxii HO. Cloning of MATA: -idiomorph resulted in a single ZsMTLA: locus encoding two Z. rouxii-like proteins MATA: 1 and MATA: 2. To assign the cloned ZsMTLα and ZsMTLA: idiomorphs as MAT, HML, and HMR cassettes, we analyzed their flanking regions. Three ZsMTLα loci exhibited the DIC1-MAT-SLA2 gene order canonical for MAT expression loci. Furthermore, four putative HML cassettes were identified, two containing the ZsMTLα copy 1 and the remaining harboring ZsMTLα copies 2 and 3. Finally, the ZsMTLA: locus was 3'-flanked by SLA2, suggesting the status of MAT expression locus. In conclusion, Z. sapae ABT301(T) displays an aααα genotype missing of the HMR silent cassette. Our results demonstrated that mating-type switching is a hypermutagenic process in Z. rouxii complex that generates genetic diversity de novo. This error-prone mechanism could be suitable to generate progenies more rapidly adaptable to hostile environments.

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.

Quality control of MATa1 splicing and exon skipping by nuclear RNA degradation

The MATa1 gene encodes a transcriptional repressor that is an important modulator of sex-specific gene expression in Saccharomyces cerevisiae. MATa1 contains two small introns, both of which need to be accurately excised for proper expression of a functional MATa1 product and to avoid production of aberrant forms of the repressor. Here, we show that unspliced and partially spliced forms of the MATa1 mRNA are degraded by the nuclear exonuclease Rat1p, the nuclear exosome and by the nuclear RNase III endonuclease Rnt1p to prevent undesired expression of non-functional a1 proteins. In addition, we show that mis-spliced forms of MATa1 in which the splicing machinery has skipped exon2 and generated exon1–exon3 products are degraded by the nuclear 5′–3′ exonuclease Rat1p and by the nuclear exosome. This function for Rat1p and the nuclear exosome in the degradation of exon-skipped products is also observed for three other genes that contain two introns (DYN2, SUS1, YOS1), identifying a novel nuclear quality control pathway for aberrantly spliced RNAs that have skipped exons.

The solution structure of the native K50 Bicoid homeodomain bound to the consensus TAATCC DNA-binding site

The solution structure of the homeodomain of the Drosophila morphogenic protein Bicoid (Bcd) complexed with a TAATCC DNA site is described. Bicoid is the only known protein that uses a homeodomain to regulate translation, as well as transcription, by binding to both RNA and DNA during early Drosophila development; in addition, the Bcd homeodomain can recognize an array of different DNA sites. The dual functionality and broad recognition capabilities signify that the Bcd homeodomain may possess unique structural/dynamic properties. Bicoid is the founding member of the K50 class of homeodomain proteins, containing a lysine residue at the critical 50th position (K50) of the homeodomain sequence, a residue required for DNA and RNA recognition; Bcd also has an arginine residue at the 54th position (R54), which is essential for RNA recognition. Bcd is the only known homeodomain with the K50/R54 combination of residues. The Bcd structure indicates that this homeodomain conforms to the conserved topology of the homeodomain motif, but exhibits a significant variation from other homeodomain structures at the end of helix 1. A key result is the observation that the side-chains of the DNA-contacting residues K50, N51 and R54 all show strong signs of flexibility in the protein-DNA interface. This finding is supportive of the adaptive-recognition theory of protein-DNA interactions.

Arm-domain interactions can provide high binding cooperativity

Peptidyl arms extending from one protein domain to another protein domain mediate many important interactions in biology. A well-studied example of this type of protein-protein interaction occurs between the yeast homeodomain proteins, MAT alpha2 and MAT a1, which form a high-affinity heterodimer on DNA. The carboxyl-terminal arm extending from MAT alpha2 to MAT a1 has been proposed to produce an allosteric conformational change in the a1 protein that generates a very large increase in the DNA binding affinity of a1. Although early studies lent some support to this model, a more recent crystal structure determination of the free a1 protein argues against any allosteric change. This note presents a thermodynamic argument that accounts for the proteins' binding behavior, so that allosteric conformational changes are not required to explain the large affinity increase. The analysis presented here should be useful in analyzing binding behavior in other systems involving arm interactions.

Crystal structures of engrailed homeodomain mutants: implications for stability and dynamics

We report the crystal structures and biophysical characterization of two stabilized mutants of the Drosophila Engrailed homeodomain that have been engineered to minimize electrostatic repulsion. Four independent copies of each mutant occupy the crystal lattice, and comparison of these structures illustrates variation that can be partly ascribed to networks of correlated conformational adjustments. Central to one network is leucine 26 (Leu26), which occupies alternatively two side chain rotameric conformations (-gauche and trans) and different positions within the hydrophobic core. Similar sets of conformational substates are observed in other Engrailed structures and in another homeodomain. The pattern of structural adjustments can account for NMR relaxation data and sequence co-variation networks in the wider homeodomain family. It may also explain the dysfunction associated with a P26L mutation in the human ARX homeodomain protein. Finally, we observe a novel dipolar interaction between a conserved tryptophan and a water molecule positioned along the normal to the indole ring. This interaction may explain the distinctive fluorescent properties of the homeodomain family.

Insights into binding cooperativity of MATa1/MATalpha2 from the crystal structure of a MATa1 homeodomain-maltose binding protein chimera

The Yeast MATa1 and MATalpha2 are homeodomain proteins that bind DNA cooperatively to repress transcription of cell type specific genes. The DNA affinity and specificity of MATa1 in the absence of MATalpha2, however, is very low. MATa1 is converted to a higher affinity DNA-binding protein by its interaction with the C-terminal tail of MATalpha2. To understand why MATa1 binds DNA weakly by itself, and how the MATalpha2 tail affects the affinity of MATa1 for DNA, we determined the crystal structure of a maltose-binding protein (MBP)-a1 chimera whose DNA binding behavior is similar to MATa1. The overall MATa1 conformation in the MBP-a1 structure, which was determined in the absence of alpha2 and DNA, is similar to that in the a1/alpha2/DNA structure. The sole difference is in the C-terminal portion of the DNA recognition helix of MATa1, which is flexible in the present structure. However, these residues are not in a location likely to be affected by binding of the MATalpha2 tail. The results argue against conformational changes in a1 induced by the tail of MATalpha2, suggesting instead that the MATalpha2 tail energetically couples the DNA binding of MATalpha2 and MATa1.

Engineered improvements in DNA-binding function of the MATa1 homeodomain reveal structural changes involved in combinatorial control

We have engineered enhanced DNA-binding function into the a1 homeodomain by making changes in a loop distant from the DNA-binding surface. Comparison of the free and bound a1 structures suggested a mechanism linking van der Waals stacking changes in this loop to the ordering of a final turn in the DNA-binding helix of a1. Inspection of the protein sequence revealed striking differences in amino acid identity at positions 24 and 25 compared to related homeodomain proteins. These positions lie in the loop connecting helix-1 and helix-2, which is involved in heterodimerization with the alpha 2 protein. A series of single and double amino acid substitutions (a1-Q24R, a1-S25Y, a1-S25F and a1-Q24R/S25Y) were engineered, expressed and purified for biochemical and biophysical study. Calorimetric measurements and HSQC NMR spectra confirm that the engineered variants are folded and are equally or more stable than the wild-type a1 homeodomain. NMR analysis of a1-Q24R/S25Y demonstrates that the DNA recognition helix (helix-3) is extended by at least one turn as a result of the changes in the loop connecting helix-1 and helix-2. As shown by EMSA, the engineered variants bind DNA with enhanced affinity (16-fold) in the absence of the alpha 2 cofactor and the variant alpha 2/a1 heterodimers bind cognate DNA with specificity and affinity reflective of the enhanced a1 binding affinity. Importantly, in vivo assays demonstrate that the a1-Q24R/S25Y protein binds with fivefold greater affinity than wild-type a1 and is able to partially suppress defects in repression by alpha 2 mutants. As a result of these studies, we show how subtle differences in residues at a surface distant from the functional site code for a conformational switch that allows the a1 homeodomain to become active in DNA binding in association with its cofactor alpha 2.

The DNA-binding domain in the Bacillus subtilis transition-state regulator AbrB employs significant motion for promiscuous DNA recognition

AbrB is a Bacillus subtilis protein responsible for regulating a diverse array of unrelated genes during periods of sub-optimal growth conditions. DNA binding by AbrB is unique in that sequence recognition is specific, yet no obvious consensus sequence of bound promoter regions is apparent. The N-terminal domain is a recently characterized representative of a novel class of DNA-binding proteins that possess a looped-hinge helix DNA-binding topology. Although the structural characterization of this DNA-binding topology contributed to an understanding of the architectural basis for recognition of DNA target sequences, specific mechanisms responsible for promiscuity in DNA sequence recognition still were not apparent. Analysis of (15)N backbone relaxation parameters shows that dynamic motion of regions directly linked to DNA binding show concerted motion on the microsecond-millisecond timescale. Furthermore, dynamic motion of the hinge region suggests that the DNA-binding region is capable of conformational orientations that allow it to accommodate DNA sequence variability in the cognate binding sites.