Structural basis for SRY-dependent 46-X,Y sex reversal: modulation of DNA bending by a naturally occurring point mutation

Skipping events impose repeated binding attempts: profound kinetic implications of protein-DNA conformational changes

The kinetics of protein-DNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. Protein-DNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRY-DNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific protein-DNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance protein-DNA recognition despite its intrinsically slow kinetics.

HMG-boxes, ribosomopathies and neurodegenerative disease

The UBTF E210K neuroregression syndrome is a predominantly neurological disorder caused by recurrent de novo dominant variants in Upstream Binding Factor, that is, essential for transcription of the ribosomal RNA genes. This unusual form of ribosomopathy is characterized by a slow decline in cognition, behavior, and sensorimotor functioning during the critical period of development. UBTF (or UBF) is a multi-HMGB-box protein that acts both as an epigenetic factor to establish "open" chromatin on the ribosomal genes and as a basal transcription factor in their RNA Polymerase I transcription. Here we review the possible mechanistic connections between the UBTF variants, ribosomal RNA gene transcription and the neuroregression syndrome, and suggest that DNA topology may play an important role.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

Evolutionary Landscape of SOX Genes to Inform Genotype-to-Phenotype Relationships

The SOX transcription factor family is pivotal in controlling aspects of development. To identify genotype-phenotype relationships of SOX proteins, we performed a non-biased study of SOX using 1890 open-reading frame and 6667 amino acid sequences in combination with structural dynamics to interpret 3999 gnomAD, 485 ClinVar, 1174 Geno2MP, and 4313 COSMIC human variants. We identified, within the HMG (High Mobility Group)- box, twenty-seven amino acids with changes in multiple SOX proteins annotated to clinical pathologies. These sites were screened through Geno2MP medical phenotypes, revealing novel SOX15 R104G associated with musculature abnormality and SOX8 R159G with intellectual disability. Within gnomAD, SOX18 E137K (rs201931544), found within the HMG box of ~0.8% of Latinx individuals, is associated with seizures and neurological complications, potentially through blood-brain barrier alterations. A total of 56 highly conserved variants were found at sites outside the HMG-box, including several within the SOX2 HMG-box-flanking region with neurological associations, several in the SOX9 dimerization region associated with Campomelic Dysplasia, SOX14 K88R (rs199932938) flanking the HMG box associated with cardiovascular complications within European populations, and SOX7 A379V (rs143587868) within an SOXF conserved far C-terminal domain heterozygous in 0.716% of African individuals with associated eye phenotypes. This SOX data compilation builds a robust genotype-to-phenotype association for a gene family through more robust ortholog data integration.

A Simulation Study of the Effect of Naturally Occurring Point Mutations on the SRY-DNA Complex

Molecular dynamics (MD) simulations were conducted in order to investigate the effect of the naturally occurring point mutations of the transcription factor (TF) sex-determining region Y (SRY) on the structure and dynamics of the SRY-DNA complex. The normal SRY, along with the two mutants I13T and G40R, comprising point mutations on the SRY chain, which have been clinically identified in patients with sex developmental disorders, were modeled as DNA complexes. Our modeling work aims at elucidating atomic-level structural determinants of the aberrant SRY-DNA complexation by means of μs-long MD. The results suggest that the observed disorders brought about by the G40R-DNA and I13T-DNA may arise predominantly from the destabilization of the complex being in accord with in vitro assays found elsewhere and from modifications of the DNA bending as revealed in this study. Comparative potential of mean force computations, over a sequence of short separation distances for the three complexes, verified a higher stability of the normal SRY-DNA. Examining the way the SRY mutations modulate the SRY-DNA complex dynamics at the microscopic level is important also toward elucidating molecular determinants of function for proteins capable of binding to DNA.

The DNA-Binding High-Mobility Group Box Domain of Sox Family Proteins Directly Interacts with RNA In Vitro

There is a growing body of evidence that a substantial number of protein domains identified as DNA-binding also interact with RNA to regulate biological processes. Several recent studies have revealed that the Sox2 transcription factor binds RNA through its high-mobility group box (HMGB) domain in vitro and in vivo. A high degree of conservation of this domain among members of the Sox family of transcription factors suggests that RNA-binding activity may be a general feature of these proteins. To address this hypothesis, we examined a subset of HMGB domains from human Sox family of proteins for their ability to bind both DNA and RNA in vitro. We observed selective, high-affinity interactions between Sox family HMGB domains and various model RNA elements, including a four-way junction RNA, a hairpin RNA with an internal bulge, G-quadruplex RNA, and a fragment of long noncoding RNA ES2, which is known to directly interact with Sox2. Importantly, the HMGB domains bind these RNA ligands significantly tighter than nonconsensus dsDNA and in some cases with affinities rivaling those of their consensus dsDNA sequences. These data suggest that RNA binding is a conserved feature of the Sox family of transcription factors with the potential to modulate unappreciated biological functions.

Structural insights into cooperative DNA recognition by the CCAAT-binding complex and its bZIP transcription factor HapX

The heterotrimeric CCAAT-binding complex (CBC) is a fundamental eukaryotic transcription factor recognizing the CCAAT box. In certain fungi, like Aspergilli, the CBC cooperates with the basic leucine zipper HapX to control iron metabolism. HapX functionally depends on the CBC, and the stable interaction of both requires DNA. To study this cooperative effect, X-ray structures of the CBC-HapX-DNA complex were determined. Downstream of the CCAAT box, occupied by the CBC, a HapX dimer binds to the major groove. The leash-like N terminus of the distal HapX subunit contacts the CBC, and via a flexible polyproline type II helix mediates minor groove interactions that stimulate sequence promiscuity. In vitro and in vivo mutagenesis suggest that the structural and functional plasticity of HapX results from local asymmetry and its ability to target major and minor grooves simultaneously. The latter feature may also apply to related transcription factors such as yeast Hap4 and distinct Yap family members.

Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations

Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.

Tenuous Transcriptional Threshold of Human Sex Determination. I. SRY and Swyer Syndrome at the Edge of Ambiguity

Male sex determination in mammals is initiated by SRY, a Y-encoded transcription factor. The protein contains a high-mobility-group (HMG) box mediating sequence-specific DNA bending. Mutations causing XY gonadal dysgenesis (Swyer syndrome) cluster in the box and ordinarily arise de novo. Rare inherited variants lead to male development in one genetic background (the father) but not another (his sterile XY daughter). De novo and inherited mutations occur at an invariant Tyr adjoining the motif's basic tail (box position 72; Y127 in SRY). In SRY-responsive cell lines CH34 and LNCaP, de novo mutations Y127H and Y127C reduced SRY activity (as assessed by transcriptional activation of principal target gene Sox9) by 5- and 8-fold, respectively. Whereas Y127H impaired testis-specific enhancer assembly, Y127C caused accelerated proteasomal proteolysis; activity was in part rescued by proteasome inhibition. Inherited variant Y127F was better tolerated: its expression was unperturbed, and activity was reduced by only twofold, a threshold similar to other inherited variants. Biochemical studies of wild-type (WT) and variant HMG boxes demonstrated similar specific DNA affinities (within a twofold range), with only subtle differences in sharp DNA bending as probed by permutation gel electrophoresis and fluorescence resonance-energy transfer (FRET); thermodynamic stabilities of the free boxes were essentially identical. Such modest perturbations are within the range of species variation. Whereas our cell-based findings rationalize the de novo genotype-phenotype relationships, a molecular understanding of inherited mutation Y127F remains elusive. Our companion study uncovers cryptic biophysical perturbations suggesting that the para-OH group of Y127 anchors a novel water-mediated DNA clamp.

Tenuous transcriptional threshold of human sex determination. II. SRY exploits water-mediated clamp at the edge of ambiguity

Y-encoded transcription factor SRY initiates male differentiation in therian mammals. This factor contains a high-mobility-group (HMG) box, which mediates sequence-specific DNA binding with sharp DNA bending. A companion article in this issue described sex-reversal mutations at box position 72 (residue 127 in human SRY), invariant as Tyr among mammalian orthologs. Although not contacting DNA, the aromatic ring seals the domain's minor wing at a solvent-exposed junction with a basic tail. A seeming paradox was posed by the native-like biochemical properties of inherited Swyer variant Y72F: its near-native gene-regulatory activity is consistent with the father's male development, but at odds with the daughter's XY female somatic phenotype. Surprisingly, aromatic rings (Y72, F72 or W72) confer higher transcriptional activity than do basic or polar side chains generally observed at solvated DNA interfaces (Arg, Lys, His or Gln). Whereas biophysical studies (time-resolved fluorescence resonance energy transfer and heteronuclear NMR spectroscopy) uncovered only subtle perturbations, dissociation of the Y72F complex was markedly accelerated relative to wild-type. Studies of protein-DNA solvation by molecular-dynamics (MD) simulations of an homologous high-resolution crystal structure (SOX18) suggest that Y72 para-OH anchors a network of water molecules at the tail-DNA interface, perturbed in the variant in association with nonlocal conformational fluctuations. Loss of the Y72 anchor among SRY variants presumably "unclamps" its basic tail, leading to (a) rapid DNA dissociation despite native affinity and (b) attenuated transcriptional activity at the edge of sexual ambiguity. Conservation of Y72 suggests that this water-mediated clamp operates generally among SRY and metazoan SOX domains.

Generation and mutational analysis of a transgenic mouse model of human SRY

SRY is the Y-chromosomal gene that determines male sex development in humans and most other mammals. After three decades of study, we still lack a detailed understanding of which domains of the SRY protein are required to engage the pathway of gene activity leading to testis development. Some insight has been gained from the study of genetic variations underlying differences/disorders of sex determination (DSD), but the lack of a system of experimentally generating SRY mutations and studying their consequences in vivo has limited progress in the field. To address this issue, we generated a mouse model carrying a human SRY transgene able to drive testis determination in XX mice. Using CRISPR-Cas9 gene editing, we generated novel genetic modifications in each of SRY's three domains (N-terminal, HMG box, and C-terminal) and performed a detailed analysis of their molecular and cellular effects on embryonic testis development. Our results provide new functional insights unique to human SRY and present a versatile and powerful system in which to functionally analyze variations of SRY including known and novel pathogenic variants found in DSD.

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

Novel SOX10 Mutations in Waardenburg Syndrome: Functional Characterization and Genotype-Phenotype Analysis

Waardenburg syndrome (WS) is a prevalent hearing loss syndrome, concomitant with focal skin pigmentation abnormalities, blue iris, and other abnormalities of neural crest-derived cells, including Hirschsprung’s disease. WS is clinically and genetically heterogeneous and it is classified into four major types WS type I, II, III, and IV (WS1, WS2, WS3, and WS4). WS1 and WS3 have the presence of dystopia canthorum, while WS3 also has upper limb anomalies. WS2 and WS4 do not have the dystopia canthorum, but the presence of Hirschsprung’s disease indicates WS4. There is a more severe subtype of WS4 with peripheral nerve and/or central nervous system involvement, namely peripheral demyelinating neuropathy, central dysmyelinating leukodystrophy, WS, and Hirschsprung’s disease or PCW/PCWH. We characterized the genetic defects underlying WS2, WS4, and the WS4-PCW/PCWH) using Sanger and whole-exome sequencing and cytogenomic microarray in seven patients from six unrelated families, including two with WS2 and five with WS4. We also performed multiple functional studies and analyzed genotype–phenotype correlations. The cohort included a relatively high frequency (80%) of individuals with neurological variants of WS4. Six novel SOX10 mutations were identified, including c.89C > A (p.Ser30^∗), c.207_8 delCG (p.Cys71Hisfs^∗62), c.479T > C (p.Leu160Pro), c.1379 delA (p.Tyr460Leufs^∗42), c.425G > C (p.Trp142Ser), and a 20-nucleotide insertion, c.1155_1174dupGCCCCACTATGGCTCAGCCT (p.Phe392Cysfs^∗117). All pathogenic variants were de novo. The results of reporter assays, western blotting, immunofluorescence, and molecular modeling supported the deleterious effects of the identified mutations and their correlations with phenotypic severity. The prediction of genotype–phenotype correlation and functional pathology, and dominant negative effect vs. haploinsufficiency in SOX10-related WS were influenced not only by site (first two vs. last coding exons) and type of mutation (missense vs. truncation/frameshift), but also by the protein expression level, molecular weight, and amino acid content of the altered protein. This in vitro analysis of SOX10 mutations thus provides a deeper understanding of the mechanisms resulting in specific WS subtypes and allows better prediction of the phenotypic manifestations, though it may not be always applicable to in vivo findings without further investigations.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

An Intergenic rs9275596 Polymorphism on Chr. 6p21 Is Associated with Multiple Sclerosis in Latvians

Background and objectives: Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system, leading to demyelination of neurons and potentially debilitating physical and mental symptoms. The disease is more prevalent in women than in men. The major histocompatibility complex (MHC) region has been identified as a major genetic determinant for autoimmune diseases, and its role in some neurological disorders including MS was evaluated. An intergenic single-nucleotide polymorphism (SNP), rs9275596, located between the HLA-DQB1 and HLA-DQA2 genes, is in significant association with various autoimmune diseases according to genome-wide association studies (GWASs). A cumulative effect of this SNP with other polymorphisms from this region was revealed. The aim of the study was to verify the data on rs9275596 association in multiple sclerosis in a case/control study of the Latvian population and to evaluate eventual functional significance of allele substitutions. Materials and Methods: rs9275596 (chr6:32713854; GRCh38.p12) was genotyped in 273 MS patients and 208 controls on main and sex-specific associations. Eventual functional significance of allele substitutions was evaluated in silico using publicly available tools. Results: The rs9275596 rare alleles were identified as a disease susceptibility factor in association with the MS main group and in affected females (p < 0.001 and p < 0.01, respectively). Risk factor genotypes with rare alleles included were associated with the MS common cohort (p < 0.002) and female cohort (odds ratio, OR = 2.24) and were identified as disease susceptible in males (OR = 2.41). It was shown that structural changes of rs9275596 affect the secondary structure of DNA. Functional significance of allele substitutions was evaluated on the eventual sequence affinity to transcription factors (TFs) and splicing signals similarity. A possible impact of the particular polymorphisms on the transcription and splicing efficiency is discussed. Conclusions: Our results suggest susceptibility of rs9275596 to multiple sclerosis in Latvians.

SOXopathies: Growing Family of Developmental Disorders Due to SOX Mutations

The SRY-related (SOX) transcription factor family pivotally contributes to determining cell fate and identity in many lineages. Since the original discovery that SRY deletions cause sex reversal, mutations in half of the 20 human SOX genes have been associated with rare congenital disorders, henceforward called SOXopathies. Mutations are generally de novo, heterozygous, and inactivating, revealing gene haploinsufficiency, but other types, including duplications, have been reported too. Missense variants primarily target the HMG domain, the SOX hallmark that mediates DNA binding and bending, nuclear trafficking, and protein-protein interactions. We here review key clinical and molecular features of SOXopathies and discuss the prospect that the disease family likely involves more SOX genes and larger clinical and genetic spectrums than currently appreciated.

DNAproDB: an interactive tool for structural analysis of DNA-protein complexes

Many biological processes are mediated by complex interactions between DNA and proteins. Transcription factors, various polymerases, nucleases and histones recognize and bind DNA with different levels of binding specificity. To understand the physical mechanisms that allow proteins to recognize DNA and achieve their biological functions, it is important to analyze structures of DNA–protein complexes in detail. DNAproDB is a web-based interactive tool designed to help researchers study these complexes. DNAproDB provides an automated structure-processing pipeline that extracts structural features from DNA–protein complexes. The extracted features are organized in structured data files, which are easily parsed with any programming language or viewed in a browser. We processed a large number of DNA–protein complexes retrieved from the Protein Data Bank and created the DNAproDB database to store this data. Users can search the database by combining features of the DNA, protein or DNA–protein interactions at the interface. Additionally, users can upload their own structures for processing privately and securely. DNAproDB provides several interactive and customizable tools for creating visualizations of the DNA–protein interface at different levels of abstraction that can be exported as high quality figures. All functionality is documented and freely accessible at http://dnaprodb.usc.edu.

Molecular and functional characterization of single-box high-mobility group B (HMGB) chromosomal protein from Aedes aegypti

High-mobility group B (HMGB) proteins have highly conserved, unique DNA-binding domains, HMG boxes, that can bind non-B-type DNA structures, such as bent, kinked and unwound structures, with high affinity. HMGB proteins also promote DNA bending, looping and unwinding. In this study, we determined the role of the Aedes aegypti single HMG-box domain protein AaHMGB; characterized its structure, spatiotemporal expression levels, subcellular localization, and nucleic acid binding activities; and compared these properties with those of its double-HMG-box counterpart protein, AaHMGB1. Via qRT-PCR, we showed that AaHMGB is expressed at much higher levels than AaHMGB1 throughout mosquito development. In situ hybridization results suggested a role for AaHMGB and AaHMGB1 during embryogenesis. Immunolocalization in the midgut revealed that AaHMGB is exclusively nuclear. Circular dichroism and fluorescence spectroscopy analyses showed that AaHMGB exhibits common features of α-helical structures and is more stably folded than AaHMGB1, likely due to the presence of one or two HMG boxes. Using several DNA substrates or single-stranded RNAs as probes, we observed significant differences between AaHMGB and AaHMGB1 in terms of their binding patterns, activity and/or specificity. Importantly, we showed that the phosphorylation of AaHMGB plays a critical role in its DNA-binding activity. Our study provides additional insight into the roles of single- versus double-HMG-box-containing proteins in nucleic acid interactions for better understanding of mosquito development, physiology and homeostasis.

Epigenetic regulation of male fate commitment from an initially bipotential system

A fundamental goal in biology is to understand how distinct cell types containing the same genetic information arise from a single stem cell throughout development. Sex determination is a key developmental process that requires a unidirectional commitment of an initially bipotential gonad towards either the male or female fate. This makes sex determination a unique model to study cell fate commitment and differentiation in vivo. We have focused this review on the accumulating evidence that epigenetic mechanisms contribute to the bipotential state of the fetal gonad and to the regulation of chromatin accessibility during and immediately downstream of the primary sex-determining switch that establishes the male fate.

The evolutionary process of mammalian sex determination genes focusing on marsupial SRYs

BACKGROUND

Maleness in mammals is genetically determined by the Y chromosome. On the Y chromosome SRY is known as the mammalian male-determining gene. Both placental mammals (Eutheria) and marsupial mammals (Metatheria) have SRY genes. However, only eutherian SRY genes have been empirically examined by functional analyses, and the involvement of marsupial SRY in male gonad development remains speculative.

RESULTS

In order to demonstrate that the marsupial SRY gene is similar to the eutherian SRY gene in function, we first examined the sequence differences between marsupial and eutherian SRY genes. Then, using a parsimony method, we identify 7 marsupial-specific ancestral substitutions, 13 eutherian-specific ancestral substitutions, and 4 substitutions that occurred at the stem lineage of therian SRY genes. A literature search and molecular dynamics computational simulations support that the lineage-specific ancestral substitutions might be involved with the functional differentiation between marsupial and eutherian SRY genes. To address the function of the marsupial SRY gene in male determination, we performed luciferase assays on the testis enhancer of Sox9 core (TESCO) using the marsupial SRY. The functional assay shows that marsupial SRY gene can weakly up-regulate the luciferase expression via TESCO.

CONCLUSIONS

Despite the sequence differences between the marsupial and eutherian SRY genes, our functional assay indicates that the marsupial SRY gene regulates SOX9 as a transcription factor in a similar way to the eutherian SRY gene. Our results suggest that SRY genes obtained the function of male determination in the common ancestor of Theria (placental mammals and marsupials). This suggests that the marsupial SRY gene has a function in male determination, but additional experiments are needed to be conclusive.

Decomposing protein-DNA binding and recognition using simplified protein models

We analyze the role of different physicochemical factors in protein/DNA binding and recognition by comparing the results from all-atom molecular dynamics simulations with simulations using simplified protein models. These models enable us to separate the role of specific amino acid side chains, formal amino acid charges and hydrogen bonding from the effects of the low-dielectric volume occupied by the protein. Comparisons are made on the basis of the conformation of DNA after protein binding, the ionic distribution around the complex and the sequence specificity. The results for four transcription factors, binding in either the minor or major grooves of DNA, show that the protein volume and formal charges, with one exception, play a predominant role in binding. Adding hydrogen bonding and a very small number of key amino acid side chains at the all-atom level yields results in DNA conformations and sequence recognition close to those seen in the reference all-atom simulations.

Impact of genetic variation on three dimensional structure and function of proteins

The Protein Data Bank (PDB; http://wwpdb.org) was established in 1971 as the first open access digital data resource in biology with seven protein structures as its initial holdings. The global PDB archive now contains more than 126,000 experimentally determined atomic level three-dimensional (3D) structures of biological macromolecules (proteins, DNA, RNA), all of which are freely accessible via the Internet. Knowledge of the 3D structure of the gene product can help in understanding its function and role in disease. Of particular interest in the PDB archive are proteins for which 3D structures of genetic variant proteins have been determined, thus revealing atomic-level structural differences caused by the variation at the DNA level. Herein, we present a systematic and qualitative analysis of such cases. We observe a wide range of structural and functional changes caused by single amino acid differences, including changes in enzyme activity, aggregation propensity, structural stability, binding, and dissociation, some in the context of large assemblies. Structural comparison of wild type and mutated proteins, when both are available, provide insights into atomic-level structural differences caused by the genetic variation.

Dual function of Ixr1 in transcriptional regulation and recognition of cisplatin-DNA adducts is caused by differential binding through its two HMG-boxes

Ixr1 is a transcriptional factor involved in the response to hypoxia, which is also related to DNA repair. It binds to DNA through its two in-tandem high mobility group box (HMG-box) domains. Each function depends on recognition of different DNA structures, B-form DNA at specific consensus sequences for transcriptional regulation, or distorted DNA, like cisplatin-DNA adducts, for DNA repair. However, the contribution of the HMG-box domains in the Ixr1 protein to the formation of different protein-DNA complexes is poorly understood. We have biophysically and biochemically characterized these interactions with specific DNA sequences from the promoters regulated by Ixr1, or with cisplatin-DNA adducts. Both HMG-boxes are necessary for transcriptional regulation, and they are not functionally interchangeable. The in-tandem arrangement of their HMG-boxes is necessary for functional folding and causes sequential cooperative binding to specific DNA sequences, with HMG-box A showing a higher contribution to DNA binding and bending than the HMG-box B. Binding of Ixr1 HMG boxes to specific DNA sequences is entropy driven, whereas binding to platinated DNA is enthalpy driven for HMG-box A and entropy driven for HMG-box B. This is the first proof that HMG-box binding to different DNA structures is associated with predictable thermodynamic differences. Based on our study, we present a model to explain the dual function of Ixr1 in the regulation of gene expression and recognition of distorted DNA structures caused by cisplatin treatment.

Protein-DNA interfaces: a molecular dynamics analysis of time-dependent recognition processes for three transcription factors

We have studied the dynamics of three transcription factor-DNA complexes using all-atom, microsecond-scale MD simulations. In each case, the salt bridges and hydrogen bond interactions formed at the protein-DNA interface are found to be dynamic, with lifetimes typically in the range of tens to hundreds of picoseconds, although some interactions, notably those involving specific binding to DNA bases, can be a hundred times longer lived. Depending on the complex studied, this dynamics may or may not lead to the existence of distinct conformational substates. Using a sequence threading technique, it has been possible to determine whether DNA sequence recognition is sensitive or not to such conformational changes, and, in one case, to show that recognition appears to be locally dependent on protein-mediated cation distributions.

Human Sex Determination at the Edge of Ambiguity: INHERITED XY SEX REVERSAL DUE TO ENHANCED UBIQUITINATION AND PROTEASOMAL DEGRADATION OF A MASTER TRANSCRIPTION FACTOR

A general problem is posed by analysis of transcriptional thresholds governing cell fate decisions in metazoan development. A model is provided by testis determination in therian mammals. Its key step, Sertoli cell differentiation in the embryonic gonadal ridge, is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in human SRY cause gonadal dysgenesis leading to XY female development (Swyer syndrome). Here, we have characterized an inherited mutation compatible with either male or female somatic phenotypes as observed in an XY father and XY daughter, respectively. The mutation (a crevice-forming substitution at a conserved back surface of the SRY high mobility group box) markedly destabilizes the domain but preserves specific DNA affinity and induced DNA bend angle. On transient transfection of diverse human and rodent cell lines, the variant SRY exhibited accelerated proteasomal degradation (relative to wild type) associated with increased ubiquitination; in vitro susceptibility to ubiquitin-independent ("default") cleavage by the 20S core proteasome was unchanged. The variant's gene regulatory activity (as assessed in a cellular model of the rat embryonic XY gonadal ridge) was reduced by 2-fold relative to wild-type SRY at similar levels of mRNA expression. Chemical proteasome inhibition restored native-like SRY expression and transcriptional activity in association with restored occupancy of a sex-specific enhancer element in principal downstream gene Sox9, demonstrating that the variant SRY exhibits essentially native activity on a per molecule basis. Our findings define a novel mechanism of impaired organogenesis, accelerated ubiquitin-directed proteasomal degradation of a master transcription factor leading to a developmental decision poised at the edge of ambiguity.

Stable isotope labeling methods for DNA

NMR is a powerful method for studying proteins and nucleic acids in solution. The study of nucleic acids by NMR is far more challenging than for proteins, which is mainly due to the limited number of building blocks and unfavorable spectral properties. For NMR studies of DNA molecules, (site specific) isotope enrichment is required to facilitate specific NMR experiments and applications. Here, we provide a comprehensive review of isotope-labeling strategies for obtaining stable isotope labeled DNA as well as specifically stable isotope labeled building blocks required for enzymatic DNA synthesis.

Structure and decoy-mediated inhibition of the SOX18/Prox1-DNA interaction

The transcription factor (TF) SOX18 drives lymphatic vessel development in both embryogenesis and tumour-induced neo-lymphangiogenesis. Genetic disruption of Sox18 in a mouse model protects from tumour metastasis and established the SOX18 protein as a molecular target. Here, we report the crystal structure of the SOX18 DNA binding high-mobility group (HMG) box bound to a DNA element regulating Prox1 transcription. The crystals diffracted to 1.75Å presenting the highest resolution structure of a SOX/DNA complex presently available revealing water structure, structural adjustments at the DNA contact interface and non-canonical conformations of the DNA backbone. To explore alternatives to challenging small molecule approaches for targeting the DNA-binding activity of SOX18, we designed a set of five decoys based on modified Prox1-DNA. Four decoys potently inhibited DNA binding of SOX18 in vitro and did not interact with non-SOX TFs. Serum stability, nuclease resistance and thermal denaturation assays demonstrated that a decoy circularized with a hexaethylene glycol linker and terminal phosphorothioate modifications is most stable. This SOX decoy also interfered with the expression of a luciferase reporter under control of a SOX18-dependent VCAM1 promoter in COS7 cells. Collectively, we propose SOX decoys as potential strategy for inhibiting SOX18 activity to disrupt tumour-induced neo-lymphangiogenesis.

Dynamics and recognition within a protein-DNA complex: a molecular dynamics study of the SKN-1/DNA interaction

Molecular dynamics simulations of the Caenorhabditis elegans transcription factor SKN-1 bound to its cognate DNA site show that the protein–DNA interface undergoes significant dynamics on the microsecond timescale. A detailed analysis of the simulation shows that movements of two key arginine side chains between the major groove and the backbone of DNA generate distinct conformational sub-states that each recognize only part of the consensus binding sequence of SKN-1, while the experimentally observed binding specificity results from a time-averaged view of the dynamic recognition occurring within this complex.

Molecular Origin of DNA Kinking by Transcription Factors

Binding of transcription factor (TF) proteins with DNA may cause severe kinks in the latter. Here, we investigate the molecular origin of the DNA kinks observed in the TF-DNA complexes using small molecule intercalation pathway, crystallographic analysis, and free energy calculations involving four different transcription factor (TF) protein-DNA complexes. We find that although protein binding may bend the DNA, bending alone is not sufficient to kink the DNA. We show that partial, not complete, intercalation is required to form the kink at a particular place in the DNA. It turns out that while amino acid alone can induce the desired kink through partial intercalation, protein provides thermodynamic stabilization of the kinked state in TF-DNA complexes.

SOXE transcription factors form selective dimers on non-compact DNA motifs through multifaceted interactions between dimerization and high-mobility group domains

The SOXE transcription factors SOX8, SOX9 and SOX10 are master regulators of mammalian development directing sex determination, gliogenesis, pancreas specification and neural crest development. We identified a set of palindromic SOX binding sites specifically enriched in regulatory regions of melanoma cells. SOXE proteins homodimerize on these sequences with high cooperativity. In contrast to other transcription factor dimers, which are typically rigidly spaced, SOXE group proteins can bind cooperatively at a wide range of dimer spacings. Using truncated forms of SOXE proteins, we show that a single dimerization (DIM) domain, that precedes the DNA binding high mobility group (HMG) domain, is sufficient for dimer formation, suggesting that DIM : HMG rather than DIM:DIM interactions mediate the dimerization. All SOXE members can also heterodimerize in this fashion, whereas SOXE heterodimers with SOX2, SOX4, SOX6 and SOX18 are not supported. We propose a structural model where SOXE-specific intramolecular DIM:HMG interactions are allosterically communicated to the HMG of juxtaposed molecules. Collectively, SOXE factors evolved a unique mode to combinatorially regulate their target genes that relies on a multifaceted interplay between the HMG and DIM domains. This property potentially extends further the diversity of target genes and cell-specific functions that are regulated by SOXE proteins.

Pioneer transcription factors target partial DNA motifs on nucleosomes to initiate reprogramming

Pioneer transcription factors (TFs) access silent chromatin and initiate cell-fate changes, using diverse types of DNA binding domains (DBDs). FoxA, the paradigm pioneer TF, has a winged helix DBD that resembles linker histone and thereby binds its target sites on nucleosomes and in compacted chromatin. Herein, we compare the nucleosome and chromatin targeting activities of Oct4 (POU DBD), Sox2 (HMG box DBD), Klf4 (zinc finger DBD), and c-Myc (bHLH DBD), which together reprogram somatic cells to pluripotency. Purified Oct4, Sox2, and Klf4 proteins can bind nucleosomes in vitro, and in vivo they preferentially target silent sites enriched for nucleosomes. Pioneer activity relates simply to the ability of a given DBD to target partial motifs displayed on the nucleosome surface. Such partial motif recognition can occur by coordinate binding between factors. Our findings provide insight into how pioneer factors can target naive chromatin sites.

A unique HMG-box domain of mouse Maelstrom binds structured RNA but not double stranded DNA

Piwi-interacting piRNAs are a major and essential class of small RNAs in the animal germ cells with a prominent role in transposon control. Efficient piRNA biogenesis and function require a cohort of proteins conserved throughout the animal kingdom. Here we studied Maelstrom (MAEL), which is essential for piRNA biogenesis and germ cell differentiation in flies and mice. MAEL contains a high mobility group (HMG)-box domain and a Maelstrom-specific domain with a presumptive RNase H-fold. We employed a combination of sequence analyses, structural and biochemical approaches to evaluate and compare nucleic acid binding of mouse MAEL HMG-box to that of canonical HMG-box domain proteins (SRY and HMGB1a). MAEL HMG-box failed to bind double-stranded (ds)DNA but bound to structured RNA. We also identified important roles of a novel cluster of arginine residues in MAEL HMG-box in these interactions. Cumulatively, our results suggest that the MAEL HMG-box domain may contribute to MAEL function in selective processing of retrotransposon RNA into piRNAs. In this regard, a cellular role of MAEL HMG-box domain is reminiscent of that of HMGB1 as a sentinel of immunogenic nucleic acids in the innate immune response.

Structure-specific nucleic acid recognition by L-motifs and their diverse roles in expression and regulation of the genome

The high-mobility group (HMG) domain containing proteins regulate transcription, DNA replication and recombination. They adopt L-shaped folds and are structure-specific DNA binding motifs. Here, I define the L-motif super-family that consists of DNA-binding HMG-box proteins and the L-motif of the histone mRNA binding domain of stem-loop binding protein (SLBP). The SLBP L-motif and HMG-box domains adopt similar L-shaped folds with three α-helices and two or three small hydrophobic cores that stabilize the overall fold, but have very different and distinct modes of nucleic acid recognition. A comparison of the structure, dynamics, protein-protein and nucleic acid interactions, and regulation by PTMs of the SLBP and the HMG-box L-motifs reveals the versatile and diverse modes by which L-motifs utilize their surfaces for structure-specific recognition of nucleic acids to regulate gene expression.

DNA-mediated cooperativity facilitates the co-selection of cryptic enhancer sequences by SOX2 and PAX6 transcription factors

Sox2 and Pax6 are transcription factors that direct cell fate decision during neurogenesis, yet the mechanism behind how they cooperate on enhancer DNA elements and regulate gene expression is unclear. By systematically interrogating Sox2 and Pax6 interaction on minimal enhancer elements, we found that cooperative DNA recognition relies on combinatorial nucleotide switches and precisely spaced, but cryptic composite DNA motifs. Surprisingly, all tested Sox and Pax paralogs have the capacity to cooperate on such enhancer elements. NMR and molecular modeling reveal very few direct protein-protein interactions between Sox2 and Pax6, suggesting that cooperative binding is mediated by allosteric interactions propagating through DNA structure. Furthermore, we detected and validated several novel sites in the human genome targeted cooperatively by Sox2 and Pax6. Collectively, we demonstrate that Sox-Pax partnerships have the potential to substantially alter DNA target specificities and likely enable the pleiotropic and context-specific action of these cell-lineage specifiers.

Structure-function relationships in human testis-determining factor SRY: an aromatic buttress underlies the specific DNA-bending surface of a high mobility group (HMG) box

Human testis determination is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in SRY cause 46 XY gonadal dysgenesis with female somatic phenotype (Swyer syndrome) and confer a high risk of malignancy (gonadoblastoma). Such mutations cluster in the SRY high mobility group (HMG) box, a conserved motif of specific DNA binding and bending. To explore structure-function relationships, we constructed all possible substitutions at a site of clinical mutation (W70L). Our studies thus focused on a core aromatic residue (position 15 of the consensus HMG box) that is invariant among SRY-related HMG box transcription factors (the SOX family) and conserved as aromatic (Phe or Tyr) among other sequence-specific boxes. In a yeast one-hybrid system sensitive to specific SRY-DNA binding, the variant domains exhibited reduced (Phe and Tyr) or absent activity (the remaining 17 substitutions). Representative nonpolar variants with partial or absent activity (Tyr, Phe, Leu, and Ala in order of decreasing side-chain volume) were chosen for study in vitro and in mammalian cell culture. The clinical mutation (Leu) was found to markedly impair multiple biochemical and cellular activities as respectively probed through the following: (i) in vitro assays of specific DNA binding and protein stability, and (ii) cell culture-based assays of proteosomal degradation, nuclear import, enhancer DNA occupancy, and SRY-dependent transcriptional activation. Surprisingly, however, DNA bending is robust to this or the related Ala substitution that profoundly impairs box stability. Together, our findings demonstrate that the folding, trafficking, and gene-regulatory function of SRY requires an invariant aromatic "buttress" beneath its specific DNA-bending surface.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

Model of a DNA-protein complex of the architectural monomeric protein MC1 from Euryarchaea

In Archaea the two major modes of DNA packaging are wrapping by histone proteins or bending by architectural non-histone proteins. To supplement our knowledge about the binding mode of the different DNA-bending proteins observed across the three domains of life, we present here the first model of a complex in which the monomeric Methanogen Chromosomal protein 1 (MC1) from Euryarchaea binds to the concave side of a strongly bent DNA. In laboratory growth conditions MC1 is the most abundant architectural protein present in Methanosarcina thermophila CHTI55. Like most proteins that strongly bend DNA, MC1 is known to bind in the minor groove. Interaction areas for MC1 and DNA were mapped by Nuclear Magnetic Resonance (NMR) data. The polarity of protein binding was determined using paramagnetic probes attached to the DNA. The first structural model of the DNA-MC1 complex we propose here was obtained by two complementary docking approaches and is in good agreement with the experimental data previously provided by electron microscopy and biochemistry. Residues essential to DNA-binding and -bending were highlighted and confirmed by site-directed mutagenesis. It was found that the Arg25 side-chain was essential to neutralize the negative charge of two phosphates that come very close in response to a dramatic curvature of the DNA.

Phosphorylation-coupled intramolecular dynamics of unstructured regions in chromatin remodeler FACT

The intrinsically disordered region (IDR) of a protein is an important topic in molecular biology. The functional significance of IDRs typically involves gene-regulation processes and is closely related to posttranslational modifications such as phosphorylation. We previously reported that the Drosophila facilitates chromatin transcription (FACT) protein involved in chromatin remodeling contains an acidic ID fragment (AID) whose phosphorylation modulates FACT binding to nucleosomes. Here, we performed dynamic atomic force microscopy and NMR analyses to clarify how the densely phosphorylated AID masks the DNA binding interface of the high-mobility-group domain (HMG). Dynamic atomic force microscopy of the nearly intact FACT revealed that a small globule temporally appears but quickly vanishes within each mobile tail-like image, corresponding to the HMG-containing IDR. The lifespan of the globule increases upon phosphorylation. NMR analysis indicated that phosphorylation induces no ordered structure but increases the number of binding sites in AID to HMG with an adjacent basic segment, thereby retaining the robust electrostatic intramolecular interaction within FACT even in the presence of DNA. These data lead to the conclusion that the inhibitory effect of nucleosome binding is ascribed to the increase in the probability of encounter between HMG and the phosphorylated IDR.

Analysis of Sry duplications on the Rattus norvegicus Y-chromosome

Background

Gene copy number variation plays a large role in the evolution of genomes. In Rattus norvegicus and other rodent species, the Y-chromosome has accumulated multiple copies of Sry loci. These copy number variations have been previously linked with changes in phenotype of animal models such as the spontaneously hypertensive rat (SHR). This study characterizes the Y-chromosome in the Sry region of Rattus norvegicus, while addressing functional variations seen in the Sry protein products.

Results

Eleven Sry loci have been identified in the SHR with one (nonHMG Sry) containing a frame shift mutation. The nonHMGSry is found and conserved in the related WKY and SD rat strains. Three new, previously unidentified, Sry loci were identified in this study (Sry3BII, Sry4 and Sry4A) in both SHR and WKY. Repetitive element analysis revealed numerous LINE-L1 elements at regions where conservation is lost among the Sry copies. In addition we have identified a retrotransposed copy of Med14 originating from spliced mRNA, two autosomal genes (Ccdc110 and HMGB1) and a normal mammalian Y-chromosome gene (Zfy) in the Sry region of the rat Y-chromosome. Translation of the sequences of each Sry gene reveals eight proteins with amino acid differences leading to changes in nuclear localization and promoter activation of a Sry-responsive gene. Sry-β (coded by the Sry2 locus) has an increased cytoplasmic fraction due to alterations at amino acid 21. Sry-γ has altered gene regulation of the Sry1 promoter due to changes at amino acid 76.

Conclusions

The duplication of Sry on the Rattus norvegicus Y-chromosome has led to proteins with altered functional ability that may have been selected for functions in addition to testis determination. Additionally, several other genes not normally found on the Y-chromosome have duplicated new copies into the region around the Sry genes. These suggest a role of active transposable elements in the evolution of the mammalian Y-chromosome in species such as Rattus norvegicus.

Conformational and thermodynamic hallmarks of DNA operator site specificity in the copper sensitive operon repressor from Streptomyces lividans

Metal ion homeostasis in bacteria relies on metalloregulatory proteins to upregulate metal resistance genes and enable the organism to preclude metal toxicity. The copper sensitive operon repressor (CsoR) family is widely distributed in bacteria and controls the expression of copper efflux systems. CsoR operator sites consist of G-tract containing pseudopalindromes of which the mechanism of operator binding is poorly understood. Here, we use a structurally characterized CsoR from Streptomyces lividans (CsoR^Sl) together with three specific operator targets to reveal the salient features pertaining to the mechanism of DNA binding. We reveal that CsoR^Sl binds to its operator site through a 2-fold axis of symmetry centred on a conserved 5′-TAC/GTA-3′ inverted repeat. Operator recognition is stringently dependent not only on electropositive residues but also on a conserved polar glutamine residue. Thermodynamic and circular dichroic signatures of the CsoR^Sl–DNA interaction suggest selectivity towards the A-DNA-like topology of the G-tracts at the operator site. Such properties are enhanced on protein binding thus enabling the symmetrical binding of two CsoR^Sl tetramers. Finally, differential binding modes may exist in operator sites having more than one 5′-TAC/GTA-3′ inverted repeat with implications in vivo for a mechanism of modular control.

Inherited human sex reversal due to impaired nucleocytoplasmic trafficking of SRY defines a male transcriptional threshold

Human testis determination is initiated by SRY (sex determining region on Y chromosome). Mutations in SRY cause gonadal dysgenesis with female somatic phenotype. Two subtle variants (V60L and I90M in the high-mobility group box) define inherited alleles shared by an XY sterile daughter and fertile father. Whereas specific DNA binding and bending are unaffected in a rat embryonic pre-Sertoli cell line, the variants exhibited selective defects in nucleocytoplasmic shuttling due to impaired nuclear import (V60L; mediated by Exportin-4) or export (I90M; mediated by chromosome region maintenance 1). Decreased shuttling limits nuclear accumulation of phosphorylated (activated) SRY, in turn reducing occupancy of DNA sites regulating Sertoli-cell differentiation [the testis-specific SRY-box 9 (Sox9) enhancer]. Despite distinct patterns of biochemical and cell-biological perturbations, V60L and I90M each attenuated Sox9 expression in transient transfection assays by twofold. Such attenuation was also observed in studies of V60A, a clinical variant associated with ovotestes and hence ambiguity between divergent cell fates. This shared twofold threshold is reminiscent of autosomal syndromes of transcription-factor haploinsufficiency, including XY sex reversal associated with mutations in SOX9. Our results demonstrate that nucleocytoplasmic shuttling of SRY is necessary for robust initiation of testicular development. Although also characteristic of ungulate orthologs, such shuttling is not conserved among rodents wherein impaired nuclear export of the high-mobility group box and import-dependent phosphorylation are compensated by a microsatellite-associated transcriptional activation domain. Human sex reversal due to subtle defects in the nucleocytoplasmic shuttling of SRY suggests that its transcriptional activity lies near the edge of developmental ambiguity.

Microsatellite-encoded domain in rodent Sry functions as a genetic capacitor to enable the rapid evolution of biological novelty

The male program of therian mammals is determined by Sry, a transcription factor encoded by the Y chromosome. Specific DNA binding is mediated by a high mobility group (HMG) box. Expression of Sry in the gonadal ridge activates a Sox9-dependent gene regulatory network leading to testis formation. A subset of Sry alleles in superfamily Muroidea (order Rodentia) is remarkable for insertion of an unstable DNA microsatellite, most commonly encoding (as in mice) a CAG repeat-associated glutamine-rich domain. We provide evidence, based on an embryonic pre-Sertoli cell line, that this domain functions at a threshold length as a genetic capacitor to facilitate accumulation of variation elsewhere in the protein, including the HMG box. The glutamine-rich domain compensates for otherwise deleterious substitutions in the box and absence of nonbox phosphorylation sites to ensure occupancy of DNA target sites. Such compensation enables activation of a male transcriptional program despite perturbations to the box. Whereas human SRY requires nucleocytoplasmic shuttling and coupled phosphorylation, mouse Sry contains a defective nuclear export signal analogous to a variant human SRY associated with inherited sex reversal. We propose that the rodent glutamine-rich domain has (i) fostered accumulation of cryptic intragenic variation and (ii) enabled unmasking of such variation due to DNA replicative slippage. This model highlights genomic contingency as a source of protein novelty at the edge of developmental ambiguity and may underlie emergence of non-Sry-dependent sex determination in the radiation of Muroidea.

CRITERIA FOR AN UPDATED CLASSIFICATION OF HUMAN TRANSCRIPTION FACTOR DNA-BINDING DOMAINS

By binding to cis-regulatory elements in a sequence-specific manner, transcription factors regulate the activity of nearby genes. Here, we discuss the criteria for a comprehensive classification of human TFs based on their DNA-binding domains. In particular, classification of basic leucine zipper (bZIP) and zinc finger factors is exemplarily discussed. The resulting classification can be used as a template for TFs of other biological species.

Impact of protein/protein interactions on global intermolecular translocation rates of the transcription factors Sox2 and Oct1 between DNA cognate sites analyzed by z-exchange NMR spectroscopy

Oct1 and Sox2 synergistically regulate developmental genes by binding to adjacent sites within promoters. We have investigated the kinetics of global intermolecular translocation of Sox2 and Oct1 between cognate sites located on different DNA molecules by z-exchange NMR spectroscopy. In the Hoxb1 promoter, the Sox2 and Oct1 sites are immediately adjacent to one another, and the intermolecular translocation rates are too slow to be measured by z-exchange spectroscopy. By introducing a 3-bp insertion between the Sox2 and Oct1 sites to mimic the spacing in the FGF4 enhancer, the interprotein contact surface is reduced, and the translocation rates are increased. Interaction between Sox2 and the POU-specific domain (POU(S)) of Oct1 does not affect the translocation mechanism but modulates the rates. Translocation involves only jumping (dissociation and reassociation) for Sox2, but both jumping and direct intersegment transfer (no dissociation into free solution) for Oct1. The dissociation (k(off) ∼1.5 s(-1)) and association (k(on) ∼5.1 × 10(9) m(-1)s(-1)) rate constants for Sox2 are reduced 4-fold and increased 5-fold, respectively, in the presence of Oct1. k(off) (∼3.5 s(-1)) for Oct1 is unaffected by Sox2, whereas k(on) (∼1.3 × 10(9) m(-1)s(-1)) is increased ∼13-fold. The direct intermolecular translocation rate (k(inter) ∼1.8 × 10(4) m(-1)s(-1)) for the POU(S) domain of Oct1 is reduced 2-fold by Sox2, whereas that for the POU homeodomain (POU(HD)) of Oct1 (k(inter) ∼ 1.7 × 10(4) m(-1)s(-1)) remains unaltered, consistent with the absence of contacts between Sox2 and POU(HD). The data suggest a model for the sequence of binding events involved in synergistic gene regulation by Sox2 and Oct1.

The crystal structure of the Sox4 HMG domain-DNA complex suggests a mechanism for positional interdependence in DNA recognition

It has recently been proposed that the sequence preferences of DNA-binding TFs (transcription factors) can be well described by models that include the positional interdependence of the nucleotides of the target sites. Such binding models allow for multiple motifs to be invoked, such as principal and secondary motifs differing at two or more nucleotide positions. However, the structural mechanisms underlying the accommodation of such variant motifs by TFs remain elusive. In the present study we examine the crystal structure of the HMG (high-mobility group) domain of Sox4 [Sry (sex-determining region on the Y chromosome)-related HMG box 4] bound to DNA. By comparing this structure with previously solved structures of Sox17 and Sox2, we observed subtle conformational differences at the DNA-binding interface. Furthermore, using quantitative electrophoretic mobility-shift assays we validated the positional interdependence of two nucleotides and the presence of a secondary Sox motif in the affinity landscape of Sox4. These results suggest that a concerted rearrangement of two interface amino acids enables Sox4 to accommodate primary and secondary motifs. The structural adaptations lead to altered dinucleotide preferences that mutually reinforce each other. These analyses underline the complexity of the DNA recognition by TFs and provide an experimental validation for the conceptual framework of positional interdependence and secondary binding motifs.

Interplay between minor and major groove-binding transcription factors Sox2 and Oct1 in translocation on DNA studied by paramagnetic and diamagnetic NMR

The pathways whereby Sox2 scans DNA to locate its specific binding site are investigated by NMR in specific and nonspecific Sox2·DNA complexes and in a specific ternary complex with Oct1 on the Hoxb1 regulatory element. Direct transfer of Sox2 between nonspecific sites on different DNA molecules occurs without dissociation into free solution at a rate of ∼10(6) M(-1) s(-1), whereas one-dimensional sliding proceeds with a diffusion constant of ≥0.1 μm(2)·s(-1). Translocation of Sox2 from one specific DNA site to another occurs via jumping, involving complete dissociation into free solution (k(d) ∼5-6 s(-1)) followed by reassociation (k(a) ∼5 × 10(8) M(-1) s(-1)). In the presence of Oct1 bound to an adjacent specific site, k(d) is reduced by more than 10-fold. Paramagnetic relaxation measurements, however, demonstrate that sparsely populated (<1%), transient states involving nonspecifically bound Sox2 in rapid exchange with specifically bound Sox2 are sampled in both binary Sox2·DNA- and ternary Oct1·Sox2·Hoxb1-DNA-specific complexes. Moreover, Sox2 modulates the mechanism of translocation of Oct1. Both Sox2 and the Oct1 POU(HD) domain are transiently released from the specific ternary complex by sliding to an adjacent nonspecific site, followed by direct transfer to another DNA molecule, whereas the Oct1 POU(S) domain is fixed to its specific site through direct interactions with Sox2. Intermolecular translocation of POU(HD) results in the formation of a bridged intermediate spanning two DNA molecules, enhancing the probability of complete intermolecular translocation of Oct1. By way of contrast, in the specific Oct1·DNA binary complex, POU(S) undergoes direct intermolecular translocation, whereas POU(HD) scans the DNA by sliding.