Principles of protein-DNA recognition revealed in the structural analysis of Ndt80-MSE DNA complexes

Transcription Factor Mavib-1 Negatively Regulates Conidiation by Affecting Utilization of Carbon and Nitrogen Source in Metarhizium acridum

Conidium is the main infection unit and reproductive unit of pathogenic fungi. Exploring the mechanism of conidiation and its regulation contributes to understanding the pathogenicity of pathogenic fungi. Vib-1, a transcription factor, was reported to participate in the conidiation process. However, the regulation mechanism of Vib-1 in conidiation is still unclear. In this study, we analyzed the function of Vib-1 and its regulation mechanism in conidiation through knocking out and overexpression of Vib-1 in entomopathogenic fungus Metarhizium acridum. Results showed that the colonial growth of Mavib-1 disruption mutant (ΔMavib-1) was significantly decreased, and conidiation was earlier compared to wild type (WT), while overexpression of Mavib-1 led to a delayed conidiation especially when carbon or nitrogen sources were insufficient. Overexpression of Mavib-1 resulted in a conidiation pattern shift from microcycle conidiation to normal conidiation on nutrient-limited medium. These results indicated that Mavib-1 acted as a positive regulator in vegetative growth and a negative regulator in conidiation by affecting utilization of carbon and nitrogen sources in M. acridum. Transcription profile analysis demonstrated that many genes related to carbon and nitrogen source metabolisms were differentially expressed in ΔMavib-1 and OE strains compared to WT. Moreover, Mavib-1 affects the conidial germination, tolerance to UV-B and heat stresses, cell wall integrity, conidial surface morphology and conidial hydrophobicity in M. acridum. These findings unravel the regulatory mechanism of Mavib-1 in fungal growth and conidiation, and enrich the knowledge to conidiation pattern shift of filamentous fungi.

Family-wide Characterization of Methylated DNA Binding Ability of Arabidopsis MBDs

13 MBD-containing genes (AtMBD1-13) have been identified in Arabidopsis thaliana so far, however, their DNA binding ability is still controversial. Here, we systematically measured the DNA binding affinities of these MBDs by ITC and EMSA binding assays, except for those of pseudogenes AtMBD3 and AtMBD13, and found that only AtMBD6 and AtMBD7 function as methylated DNA readers. We also found that the MBD of AtMBD5 exhibits very weak binding to methylated DNA compared to that of AtMBD6. To further investigate the structural basis of AtMBDs in binding to methylated DNA, we determined the complex structure of the AtMBD6 MBD with a 12mer mCG DNA and the apo structure of the AtMBD5 MBD. Structural analysis coupled with mutagenesis studies indicated that, in addition to the conserved arginine fingers contributing to the DNA binding specificity, the residues located in the loop1 and α1 are also essential for the methylated DNA binding of these MBDs in Arabidopsis thaliana, which explains why AtMBD5 MBD and the other AtMBDs display very weak or no binding to methylated DNA. Thus, our study here systematically demonstrates the DNA binding ability of the MBDs in Arabidopsis thaliana, which also provides a general guideline in understanding the DNA binding ability of the MBDs in other plants as a whole.

Structural basis for the ability of MBD domains to bind methyl-CG and TG sites in DNA

Cytosine methylation is a well-characterized epigenetic mark and occurs at both CG and non-CG sites in DNA. Both methylated CG (mCG)- and mCH (H = A, C, or T)-containing DNAs, especially mCAC-containing DNAs, are recognized by methyl-CpG-binding protein 2 (MeCP2) to regulate gene expression in neuron development. However, the molecular mechanism involved in the binding of methyl-CpG-binding domain (MBD) of MeCP2 to these different DNA motifs is unclear. Here, we systematically characterized the DNA-binding selectivities of the MBD domains in MeCP2 and MBD1-4 with isothermal titration calorimetry-based binding assays, mutagenesis studies, and X-ray crystallography. We found that the MBD domains of MeCP2 and MBD1-4 bind mCG-containing DNAs independently of the sequence identity outside the mCG dinucleotide. Moreover, some MBD domains bound to both methylated and unmethylated CA dinucleotide-containing DNAs, with a preference for the CAC sequence motif. We also found that the MBD domains bind to mCA or nonmethylated CA DNA by recognizing the complementary TG dinucleotide, which is consistent with an overlooked ligand of MeCP2, i.e. the matrix/scaffold attachment regions (MARs/SARs) with a consensus sequence of 5'-GGTGT-3' that was identified in early 1990s. Our results also explain why MeCP2 exhibits similar binding affinity to both mCA- and hmCA-containing dsDNAs. In summary, our results suggest that in addition to mCG sites, unmethylated CA or TG sites also serve as DNA-binding sites for MeCP2 and other MBD-containing proteins. This discovery expands the genome-wide activity of MBD-containing proteins in gene regulation.

The N-acetylglucosamine catabolic gene cluster in Trichoderma reesei is controlled by the Ndt80-like transcription factor RON1

Chitin is an important structural constituent of fungal cell walls composed of N-acetylglucosamine (GlcNAc) monosaccharides, but catabolism of GlcNAc has not been studied in filamentous fungi so far. In the yeast Candida albicans, the genes encoding the three enzymes responsible for stepwise conversion of GlcNAc to fructose-6-phosphate are clustered. In this work, we analysed GlcNAc catabolism in ascomycete filamentous fungi and found that the respective genes are also clustered in these fungi. In contrast to C. albicans, the cluster often contains a gene for an Ndt80-like transcription factor, which we named RON1 (regulator of N-acetylglucosamine catabolism 1). Further, a gene for a glycoside hydrolase 3 protein related to bacterial N-acetylglucosaminidases can be found in the GlcNAc gene cluster in filamentous fungi. Functional analysis in Trichoderma reesei showed that the transcription factor RON1 is a key activator of the GlcNAc gene cluster and essential for GlcNAc catabolism. Furthermore, we present an evolutionary analysis of Ndt80-like proteins in Ascomycota. All GlcNAc cluster genes, as well as the GlcNAc transporter gene ngt1, and an additional transcriptional regulator gene, csp2, encoding the homolog of Neurospora crassa CSP2/GRHL, were functionally characterised by gene expression analysis and phenotypic characterisation of knockout strains in T. reesei.

Structures of complexes comprised of Fischerella transcription factor HetR with Anabaena DNA targets

HetR is an essential regulator of heterocyst development in cyanobacteria. Many mutations in HetR render Anabaena incapable of nitrogen fixation. The protein binds to a DNA palindrome upstream of hetP and other genes. We have determined the crystal structures of HetR complexed with palindromic DNA targets, 21, 23, and 29 bp at 2.50-, 3.00-, and 3.25-Å resolution, respectively. The highest-resolution structure shows fine details of specific protein-DNA interactions. The lower-resolution structures with longer DNA duplexes have similar interaction patterns and show how the flap domains interact with DNA in a sequence nonspecific fashion. Fifteen of 15 protein-DNA contacts predicted on the basis of the structure were confirmed by single amino acid mutations that abolished binding in vitro and complementation in vivo. A striking feature of the structure is the association of glutamate 71 from each subunit of the HetR dimer with three successive cytosines in each arm of the palindromic target, a feature that is conserved among all known heterocyst-forming cyanobacteria sequenced to date.

A common mode of recognition for methylated CpG

Much is known about vertebrate DNA methylation, however it is not known how methylated CpG within particular sequences is recognized. Two recent structures of C2H2 zinc finger (ZnF) proteins in complex with methylated DNA reveal a common recognition mode for 5-methylcytosine (5mC) that involves a 5mC-Arg-G triad. In the two ZnF proteins, an arginine that precedes the first Zn-binding histidine (RH motif) can interact with a 5mCpG or TpG dinucleotide. Among a family of >300 human Krüppel-associated box (KRAB) domain containing ZnF proteins examined, two-thirds contained at least one ZnF that included an RH motif. We propose that the RH-ZnF motifs provide specificity for 5mCpG, whereas the neighboring Zn fingers recognize the surrounding DNA sequence context.

The Sum1/Ndt80 transcriptional switch and commitment to meiosis in Saccharomyces cerevisiae

Cells encounter numerous signals during the development of an organism that induce division, differentiation, and apoptosis. These signals need to be present for defined intervals in order to induce stable changes in the cellular phenotype. The point after which an inducing signal is no longer needed for completion of a differentiation program can be termed the "commitment point." Meiotic development in the yeast Saccharomyces cerevisiae (sporulation) provides a model system to study commitment. Similar to differentiation programs in multicellular organisms, the sporulation program in yeast is regulated by a transcriptional cascade that produces early, middle, and late sets of sporulation-specific transcripts. Although critical meiosis-specific events occur as early genes are expressed, commitment does not take place until middle genes are induced. Middle promoters are activated by the Ndt80 transcription factor, which is produced and activated shortly before most middle genes are expressed. In this article, I discuss the connection between Ndt80 and meiotic commitment. A transcriptional regulatory pathway makes NDT80 transcription contingent on the prior expression of early genes. Once Ndt80 is produced, the recombination (pachytene) checkpoint prevents activation of the Ndt80 protein. Upon activation, Ndt80 triggers a positive autoregulatory loop that leads to the induction of genes that promote exit from prophase, the meiotic divisions, and spore formation. The pathway is controlled by multiple feed-forward loops that give switch-like properties to the commitment transition. The conservation of regulatory components of the meiotic commitment pathway and the recently reported ability of Ndt80 to increase replicative life span are discussed.

Recognition of methylated DNA through methyl-CpG binding domain proteins

DNA methylation is a key regulatory control route in epigenetics, involving gene silencing and chromosome inactivation. It has been recognized that methyl-CpG binding domain (MBD) proteins play an important role in interpreting the genetic information encoded by methylated DNA (mDNA). Although the function of MBD proteins has attracted considerable attention and is well characterized, the mechanism underlying mDNA recognition by MBD proteins is still poorly understood. In this article, we demonstrate that the methyl-CpG dinucleotides are recognized at the MBD–mDNA interface by two MBD arginines through an interplay of hydrogen bonding and cation-π interaction. Through molecular dynamics and quantum-chemistry calculations we investigate the methyl-cytosine recognition process and demonstrate that methylation enhances MBD–mDNA binding by increasing the hydrophobic interfacial area and by strengthening the interaction between mDNA and MBD proteins. Free-energy perturbation calculations also show that methylation yields favorable contribution to the binding free energy for MBD–mDNA complex.

From nonspecific DNA-protein encounter complexes to the prediction of DNA-protein interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.

Many essential biological activities require interactions between DNA and proteins. These proteins usually use certain amino acids, called DNA-binding sites, to recognize their specific DNA targets. To facilitate the search of its specific DNA targets, a DNA-binding protein often associates with nonspecific DNA and then diffuses along the DNA. Due to the weak interactions between nonspecific DNA and the protein, structural characterization of nonspecific DNA–protein complexes is experimentally challenging. This paper describes a computational modeling study on nonspecific DNA–protein complexes and comparative analysis with respect to specific DNA–protein complexes. The study found that the specific DNA-binding sites on a protein are typically favorable for nonspecific DNA and that nonspecific and specific DNA–protein interaction modes are quite similar. This similarity may reflect an important sampling step in the search for the specific DNA target sequence by a DNA-binding protein. On the basis of these observations, a novel method was proposed for predicting DNA-binding sites and binding modes of a DNA-binding protein without knowing its specific DNA target sequence. Ultimately, the combination of this method and protein structure prediction may lead the way to high throughput modeling of DNA–protein interactions.

Investigation of transcription factor Ndt80 affinity differences for wild type and mutant DNA: a molecular dynamics study

Molecular dynamics simulations and free energy calculations have been performed on the transcription factor Ndt80 either in complex with the native DNA sequence or with a mutant DNA with a switched central base pair, C5-G5' to G5-C5'. This mutant has been shown to have a 100-fold decrease in binding affinity of Ndt80, and in this study we explain this both structurally and energetically. The major interactions between the protein and the DNA were maintained in the simulations, apart from around the mutation site. The crystal structure of the Ndt80-DNA complex revealed that R177 makes a base specific bidentate interaction with G5' which is part of a conserved 5'-YpG-3' step. In the simulation with the mutant DNA, the side chain of R177 changes conformation and makes three new stable hydrogen bonds to the DNA backbone. This in turn induces a conformational change in the DNA backbone of the T6'-G5' step from the unusual BII state to the canonical BI state. The affinity difference for the protein-DNA complex with the native DNA compared with the mutant DNA is only about 3 kcal/mol. The free energy calculations of the base pair switch indicated a larger difference than what was found experimentally, about 7.7 kcal/mol, but this is explained in structural terms using the 10 ns simulations of the solvated complexes and the rearrangement of the R177 side chain.

Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site

The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.