Cooperative DNA-binding and sequence-recognition mechanism of aristaless and clawless

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Cooperative Gsx2-DNA Binding Requires DNA Bending and a Novel Gsx2 Homeodomain Interface

The conserved Gsx homeodomain (HD) transcription factors specify neural cell fates in animals from flies to mammals. Like many HD proteins, Gsx factors bind A/T-rich DNA sequences prompting the question - how do HD factors that bind similar DNA sequences in vitro regulate specific target genes in vivo? Prior studies revealed that Gsx factors bind DNA both as a monomer on individual A/T-rich sites and as a cooperative homodimer to two sites spaced precisely seven base pairs apart. However, the mechanistic basis for Gsx DNA binding and cooperativity are poorly understood. Here, we used biochemical, biophysical, structural, and modeling approaches to (1) show that Gsx factors are monomers in solution and require DNA for cooperative complex formation; (2) define the affinity and thermodynamic binding parameters of Gsx2/DNA interactions; (3) solve a high-resolution monomer/DNA structure that reveals Gsx2 induces a 20° bend in DNA; (4) identify a Gsx2 protein-protein interface required for cooperative DNA binding; and (5) determine that flexible spacer DNA sequences enhance Gsx2 cooperativity on dimer sites. Altogether, our results provide a mechanistic basis for understanding the protein and DNA structural determinants that underlie cooperative DNA binding by Gsx factors, thereby providing a deeper understanding of HD specificity.

Coevolution of the Tlx homeobox gene with medusa development (Cnidaria: Medusozoa)

Cnidarians display a wide diversity of life cycles. Among the main cnidarian clades, only Medusozoa possesses a swimming life cycle stage called the medusa, alternating with a benthic polyp stage. The medusa stage was repeatedly lost during medusozoan evolution, notably in the most diverse medusozoan class, Hydrozoa. Here, we show that the presence of the homeobox gene Tlx in Cnidaria is correlated with the presence of the medusa stage, the gene having been lost in clades that ancestrally lack a medusa (anthozoans, endocnidozoans) and in medusozoans that secondarily lost the medusa stage. Our characterization of Tlx expression indicate an upregulation of Tlx during medusa development in three distantly related medusozoans, and spatially restricted expression patterns in developing medusae in two distantly related species, the hydrozoan Podocoryna carnea and the scyphozoan Pelagia noctiluca. These results suggest that Tlx plays a key role in medusa development and that the loss of this gene is likely linked to the repeated loss of the medusa life cycle stage in the evolution of Hydrozoa.

Structural and functional analysis of the human cone‐rod homeobox transcription factor

The cone‐rod homeobox (CRX) protein is a critical K50 homeodomain transcription factor responsible for the differentiation and maintenance of photoreceptor neurons in the vertebrate retina. Mutant alleles in the human gene encoding CRX result in a variety of distinct blinding retinopathies, including retinitis pigmentosa, cone‐rod dystrophy, and Leber congenital amaurosis. Despite the success of using in vitro biochemistry, animal models, and genomics approaches to study this clinically relevant transcription factor over the past 25 years since its initial characterization, there are no high‐resolution structures in the published literature for the CRX protein. In this study, we use bioinformatic approaches and small‐angle X‐ray scattering (SAXS) structural analysis to further understand the biochemical complexity of the human CRX homeodomain (CRX‐HD). We find that the CRX‐HD is a compact, globular monomer in solution that can specifically bind functional cis‐regulatory elements encoded upstream of retina‐specific genes. This study presents the first structural analysis of CRX, paving the way for a new approach to studying the biochemistry of this protein and its disease‐causing mutant protein variants.

Unheeded SARS-CoV-2 proteins? A deep look into negative-sense RNA

SARS-CoV-2 is a novel positive-sense single-stranded RNA virus from the Coronaviridae family (genus Betacoronavirus), which has been established as causing the COVID-19 pandemic. The genome of SARS-CoV-2 is one of the largest among known RNA viruses, comprising of at least 26 known protein-coding loci. Studies thus far have outlined the coding capacity of the positive-sense strand of the SARS-CoV-2 genome, which can be used directly for protein translation. However, it has been recently shown that transcribed negative-sense viral RNA intermediates that arise during viral genome replication from positive-sense viruses can also code for proteins. No studies have yet explored the potential for negative-sense SARS-CoV-2 RNA intermediates to contain protein-coding loci. Thus, using sequence and structure-based bioinformatics methodologies, we have investigated the presence and validity of putative negative-sense ORFs (nsORFs) in the SARS-CoV-2 genome. Nine nsORFs were discovered to contain strong eukaryotic translation initiation signals and high codon adaptability scores, and several of the nsORFs were predicted to interact with RNA-binding proteins. Evolutionary conservation analyses indicated that some of the nsORFs are deeply conserved among related coronaviruses. Three-dimensional protein modeling revealed the presence of higher order folding among all putative SARS-CoV-2 nsORFs, and subsequent structural mimicry analyses suggest similarity of the nsORFs to DNA/RNA-binding proteins and proteins involved in immune signaling pathways. Altogether, these results suggest the potential existence of still undescribed SARS-CoV-2 proteins, which may play an important role in the viral lifecycle and COVID-19 pathogenesis.

Structure of the unique tetrameric STENOFOLIA homeodomain bound with target promoter DNA

Homeobox transcription factors are key regulators of morphogenesis and development in both animals and plants. In plants, the WUSCHEL-related homeobox (WOX) family of transcription factors function as central organizers of several developmental programs ranging from embryo patterning to meristematic stem-cell maintenance through transcriptional activation and repression mechanisms. The Medicago truncatula STENOFOLIA (STF) gene is a master regulator of leaf-blade lateral development. Here, the crystal structure of the homeodomain (HD) of STF (STF-HD) in complex with its promoter DNA is reported at 2.1 Å resolution. STF-HD binds DNA as a tetramer, enclosing nearly the entire bound DNA surface. The STF-HD tetramer is partially stabilized by docking of the C-terminal tail of one protomer onto a conserved hydrophobic surface on the head of another protomer in a head-to-tail manner. STF-HD specifically binds TGA motifs, although the promoter sequence also contains TAAT motifs. Helix α3 not only serves a canonical role as a base reader in the major groove, but also provides DNA binding in the minor groove through basic residues located at its C-terminus. The structural and functional data in planta reported here provide new insights into the DNA-binding mechanisms of plant-specific HDs from the WOX family of transcription factors.

Dynamic Molecular Evolution of Mammalian Homeobox Genes: Duplication, Loss, Divergence and Gene Conversion Sculpt PRD Class Repertoires

The majority of homeobox genes are highly conserved across animals, but the eutherian-specific ETCHbox genes, embryonically expressed and highly divergent duplicates of CRX, are a notable exception. Here we compare the ETCHbox genes of 34 mammalian species, uncovering dynamic patterns of gene loss and tandem duplication, including the presence of a large tandem array of LEUTX loci in the genome of the European rabbit (Oryctolagus cuniculus). Despite extensive gene gain and loss, all sampled species possess at least two ETCHbox genes, suggesting their collective role is indispensable. We find evidence for positive selection and show that TPRX1 and TPRX2 have been the subject of repeated gene conversion across the Boreoeutheria, homogenising their sequences and preventing divergence, especially in the homeobox region. Together, these results are consistent with a model where mammalian ETCHbox genes are dynamic in evolution due to functional overlap, yet have collective indispensable roles.

Molecular Cloning, Purification and Characterization of Mce1R of Mycobacterium tuberculosis

The mce1 operon of Mycobacterium tuberculosis, important for lipid metabolism/transport, host cell invasion, modulation of host immune response and pathogenicity, is under the transcriptional control of Mce1R. Hence characterizing Mce1R is an important step for novel anti-tuberculosis drug discovery. The present study reports functional and in silico characterization of Mce1R. In this work, we have computationally modeled the structure of Mce1R and have validated the structure by computational and experimental methods. Mce1R has been shown to harbor the canonical VanR-like structure with a flexible N-terminal 'arm', carrying conserved positively charged residues, most likely involved in the operator DNA binding. The mce1R gene has been cloned, expressed, purified and its DNA-binding activity has been measured in vitro. The K_d value for Mce1R-operator DNA interaction has been determined to be 0.35 ± 0.02 µM which implies that Mce1R binds to DNA with moderate affinity compared to the other FCD family of regulators. So far, this is the first report for measuring the DNA-binding affinity of any VanR-type protein. Despite significant sequence similarity at the N-terminal domain, the wHTH motif of Mce1R exhibits poor conservancy of amino acid residues, critical for DNA-binding, thus results in moderate DNA-binding affinity. The N-terminal DNA-binding domain is structurally dynamic while the C-terminal domain showed significant stability and such profile of structural dynamics is most likely to be preserved in the structural orthologs of Mce1R. In addition to this, a cavity has been detected in the C-terminal domain of Mce1R which contains a few conserved residues. Comparison with other FCD family of regulators suggests that most of the conserved residues might be critical for binding to specific ligand. The max pKd value and drug score for the cavity are estimated to be 9.04 and 109 respectively suggesting that the cavity represents a suitable target site for novel anti-tuberculosis drug discovery approaches.

Investigating Developmental and Epileptic Encephalopathy Using Drosophila melanogaster

Developmental and epileptic encephalopathies (DEEs) are the spectrum of severe epilepsies characterized by early-onset, refractory seizures occurring in the context of developmental regression or plateauing. Early infantile epileptic encephalopathy (EIEE) is one of the earliest forms of DEE, manifesting as frequent epileptic spasms and characteristic electroencephalogram findings in early infancy. In recent years, next-generation sequencing approaches have identified a number of monogenic determinants underlying DEE. In the case of EIEE, 85 genes have been registered in Online Mendelian Inheritance in Man as causative genes. Model organisms are indispensable tools for understanding the in vivo roles of the newly identified causative genes. In this review, we first present an overview of epilepsy and its genetic etiology, especially focusing on EIEE and then briefly summarize epilepsy research using animal and patient-derived induced pluripotent stem cell (iPSC) models. The Drosophila model, which is characterized by easy gene manipulation, a short generation time, low cost and fewer ethical restrictions when designing experiments, is optimal for understanding the genetics of DEE. We therefore highlight studies with Drosophila models for EIEE and discuss the future development of their practical use.

Distortion of double-stranded DNA structure by the binding of the restriction DNA glycosylase R.PabI

R.PabI is a restriction DNA glycosylase that recognizes the sequence 5′-GTAC-3′ and hydrolyses the N-glycosidic bond of adenine in the recognition sequence. R.PabI drastically bends and unwinds the recognition sequence of double-stranded DNA (dsDNA) and flips the adenine and guanine bases in the recognition sequence into the catalytic and recognition sites on the protein surface. In this study, we determined the crystal structure of the R.PabI-dsDNA complex in which the dsDNA is drastically bent by the binding of R.PabI but the base pairs are not unwound. This structure is predicted to be important for the indirect readout of the recognition sequence by R.PabI. In the complex structure, wedge loops of the R.PabI dimer are inserted into the minor groove of dsDNA to stabilize the deformed dsDNA structure. A base stacking is distorted between the two wedge-inserted regions. R.PabI is predicted to utilize the distorted base stacking for the detection of the recognition sequence.

Structural changes in DNA-binding proteins on complexation

Characterization and prediction of the DNA-biding regions in proteins are essential for our understanding of how proteins recognize/bind DNA. We analyze the unbound (U) and the bound (B) forms of proteins from the protein–DNA docking benchmark that contains 66 binary protein–DNA complexes along with their unbound counterparts. Proteins binding DNA undergo greater structural changes on complexation (in particular, those in the enzyme category) than those involved in protein–protein interactions (PPI). While interface atoms involved in PPI exhibit an increase in their solvent-accessible surface area (ASA) in the bound form in the majority of the cases compared to the unbound interface, protein–DNA interactions indicate increase and decrease in equal measure. In 25% structures, the U form has missing residues which are located in the interface in the B form. The missing atoms contribute more toward the buried surface area compared to other interface atoms. Lys, Gly and Arg are prominent in disordered segments that get ordered in the interface on complexation. In going from U to B, there may be an increase in coil and helical content at the expense of turns and strands. Consideration of flexibility cannot distinguish the interface residues from the surface residues in the U form.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease

Purpose

RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease.

Methods

Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members.

Results

Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry.

Conclusion

This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.

Arginine side chain stacking with peptide plane stabilizes the protein helix conformation in a cooperative way

A combined experimental and computational study is performed for arginine side chain stacking with the protein α-helix. Theremostability measurements of Aristaless homeodomain, a helical protein, suggest that mutating the arginine residue R106, R137 or R141, which has the guanidino side chain stacking with the peptide plane, to alanine, destabilizes the protein. The R-PP stacking has an energy of ∼0.2-0.4 kcal/mol. This stacking interaction mainly comes from dispersion and electrostatics, based on MP2 calculations with the energy decomposition analysis. The calculations also suggest that the stacking stabilizes 2 backbone-backbone h-bonds (i→i-4 and i-3→i-7) in a cooperative way. Desolvation and electrostatic polarization are responsible for cooperativity with the i→i-4 and i-3→i-7 h-bonds, respectively. This cooperativity is supported by a protein α-helices h-bond survey in the pdb databank where stacking shortens the corresponding h-bond distances.

Biallelic mutations in the homeodomain of NKX6-2 underlie a severe hypomyelinating leukodystrophy

Hypomyelinating leukodystrophies are genetically heterogeneous disorders with overlapping clinical and neuroimaging features reflecting variable abnormalities in myelin formation. We report on the identification of biallelic inactivating mutations in NKX6-2, a gene encoding a transcription factor regulating multiple developmental processes with a main role in oligodendrocyte differentiation and regulation of myelin-specific gene expression, as the cause underlying a previously unrecognized severe variant of hypomyelinating leukodystrophy. Five affected subjects (three unrelated families) were documented to share biallelic inactivating mutations affecting the NKX6-2 homeobox domain. A trio-based whole exome sequencing analysis in the first family detected a homozygous frameshift change [c.606delinsTA; p.(Lys202Asnfs*?)]. In the second family, homozygosity mapping coupled to whole exome sequencing identified a homozygous nucleotide substitution (c.565G>T) introducing a premature stop codon (p.Glu189*). In the third family, whole exome sequencing established compound heterozygosity for a non-conservative missense change affecting a key residue participating in DNA binding (c.599G>A; p.Arg200Gln) and a nonsense substitution (c.589C>T; p.Gln197*), in both affected siblings. The clinical presentation was homogeneous, with four subjects having severe motor delays, nystagmus and absent head control, and one individual showing gross motor delay at the age of 6 months. All exhibited neuroimaging that was consistent with hypomyelination. These findings define a novel, severe form of leukodystrophy caused by impaired NKX6-2 function.

Tissue-specific enhancer repression through molecular integration of cell signaling inputs

Drosophila leg morphogenesis occurs under the control of a relatively well-known genetic cascade, which mobilizes both cell signaling pathways and tissue-specific transcription factors. However, their cross-regulatory interactions, deployed to refine leg patterning, remain poorly characterized at the gene expression level. Within the genetically interacting landscape that governs limb development, the bric-à-brac2 (bab2) gene is required for distal leg segmentation. We have previously shown that the Distal-less (Dll) homeodomain and Rotund (Rn) zinc-finger activating transcription factors control limb-specific bab2 expression by binding directly a single critical leg/antennal enhancer (LAE) within the bric-à-brac locus. By genetic and molecular analyses, we show here that the EGFR-responsive C15 homeodomain and the Notch-regulated Bowl zinc-finger transcription factors also interact directly with the LAE enhancer as a repressive duo. The appendage patterning gene bab2 is the first identified direct target of the Bowl repressor, an Odd-skipped/Osr family member. Moreover, we show that C15 acts on LAE activity independently of its regular partner, the Aristaless homeoprotein. Instead, we find that C15 interacts physically with the Dll activator through contacts between their homeodomain and binds competitively with Dll to adjacent cognate sites on LAE, adding potential new layers of regulation by C15. Lastly, we show that C15 and Bowl activities regulate also rn expression. Our findings shed light on how the concerted action of two transcriptional repressors, in response to cell signaling inputs, shapes and refines gene expression along the limb proximo-distal axis in a timely manner.

Limb morphogenesis is controlled by a well-known genetic cascade, mobilizing both cell signaling and tissue-specific transcription factors (TFs). However, how their concerted action refines gene expression remains to be deciphered. It is thus crucial to understand how cell signaling inputs are integrated by transcriptional “enhancers”. The Drosophila leg provides a good paradigm to dissect the molecular mechanisms underlying gene regulation. Here, we used the bric-a-brac2 (bab2) gene as a model to study the integrated regulation of patterning genes implicated in tarsal segmentation. bab2 expression in the leg primordium is dynamic and complex, going from initial broad distal expression to precisely positioned tarsal rings. By genetic and molecular analyses, we show here that the cell signaling-responding TFs C15 and Bowl interact directly with the limb-specific bab2 enhancer as a repressive duo. Moreover, C15 acts independently of its partner Aristaless through physical interaction with the Dll activator. We propose that Dll induces early circular bab2 expression pattern, then EGFR signaling-induced C15 in the distalmost cells competes with Dll for LAE binding and resolves bab2 pattern as a ring. Taken together our data shed light on how the concerted action of a quartet of transcription factors reshapes gene expression during limb proximo-distal axis development.

Tetrameric structure of the restriction DNA glycosylase R.PabI in complex with nonspecific double-stranded DNA

R.PabI is a type II restriction enzyme that recognizes the 5′-GTAC-3′ sequence and belongs to the HALFPIPE superfamily. Although most restriction enzymes cleave phosphodiester bonds at specific sites by hydrolysis, R.PabI flips the guanine and adenine bases of the recognition sequence out of the DNA helix and hydrolyzes the N-glycosidic bond of the flipped adenine in a similar manner to DNA glycosylases. In this study, we determined the structure of R.PabI in complex with double-stranded DNA without the R.PabI recognition sequence by X-ray crystallography. The 1.9 Å resolution structure of the complex showed that R.PabI forms a tetrameric structure to sandwich the double-stranded DNA and the tetrameric structure is stabilized by four salt bridges. DNA binding and DNA glycosylase assays of the R.PabI mutants showed that the residues that form the salt bridges (R70 and D71) are essential for R.PabI to find the recognition sequence from the sea of nonspecific sequences. R.PabI is predicted to utilize the tetrameric structure to bind nonspecific double-stranded DNA weakly and slide along it to find the recognition sequence.

A Novel Approach to Predict Core Residues on Cancer-Related DNA-Binding Domains

Protein-DNA interactions are involved in different cancer pathways. In particular, the DNA-binding domains of proteins can determine where and how gene regulatory regions are bound in different cell lines at different stages. Therefore, it is essential to develop a method to predict and locate the core residues on cancer-related DNA-binding domains. In this study, we propose a computational method to predict and locate core residues on DNA-binding domains. In particular, we have selected the cancer-related DNA-binding domains for in-depth studies, namely, winged Helix Turn Helix family, homeodomain family, and basic Helix-Loop-Helix family. The results demonstrate that the proposed method can predict the core residues involved in protein-DNA interactions, as verified by the existing structural data. Given its good performance, various aspects of the method are discussed and explored: for instance, different uses of prediction algorithm, different protein domains, and hotspot threshold setting.

Molecular dynamics simulations show altered secondary structure of clawless in binary complex with DNA providing insights into aristaless-clawless-DNA ternary complex formation

Aristaless (Al) and clawless (Cll) homeodomains that are involved in leg development in Drosophila melanogaster are known to bind cooperatively to 5'-(T/C)TAATTAA(T/A)(T/A)G-3' DNA sequence, but the mechanism of their binding to DNA is unknown. Molecular dynamics (MD) studies have been carried out on binary, ternary, and reconstructed protein-DNA complexes involving Al, Cll, and DNA along with binding free energy analysis of these complexes. Analysis of MD trajectories of Cll-3A01, binary complex reveals that C-terminal end of helixIII of Cll, unwind in the absence of Al and remains so in reconstructed ternary complex, Cll-3A01-Al. In addition, this change in secondary structure of Cll does not allow it to form protein-protein interactions with Al in the ternary reconstructed complex. However, secondary structure of Cll and its interactions are maintained in other reconstructed ternary complex, Al-3A01-Cll where Cll binds to Al-3A01, binary complex to form ternary complex. These interactions as observed during MD simulations compare well with those observed in ternary crystal structure. Thus, this study highlights the role of helixIII of Cll and protein-protein interactions while proposing likely mechanism of recognition in ternary complex, Al-Cll-DNA.

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain-DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1-DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein-DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein-DNA complexes. The order of stability of protein-DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein-DNA complexes. Among specific protein-DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.

Single-molecule imaging and kinetic analysis of cooperative cofilin-actin filament interactions

The actin filament-severing protein actin depolymerizing factor (ADF)/cofilin is ubiquitously distributed among eukaryotes and modulates actin dynamics. The cooperative binding of cofilin to actin filaments is crucial for the concentration-dependent unconventional modulation of actin dynamics by cofilin. In this study, the kinetic parameters associated with the cooperative binding of cofilin to actin filaments were directly evaluated using a single-molecule imaging technique. The on-rate of cofilin binding to the actin filament was estimated to be 0.06 µM(-1)⋅s(-1) when the cofilin concentration was in the range of 30 nM to 1 µM. A dwell time histogram of cofilin bindings decays exponentially to give an off-rate of 0.6 s(-1). During long-term cofilin binding events (>0.4 s), additional cofilin bindings were observed in the vicinity of the initial binding site. The on-rate for these events was 2.3-fold higher than that for noncontiguous bindings. Super-high-resolution image analysis of the cofilin binding location showed that the on-rate enhancement occurred within 65 nm of the original binding event. By contrast, the cofilin off-rate was not affected by the presence of prebound cofilin. Neither decreasing the temperature nor increasing the viscosity of the test solution altered the on-rates, off-rates, or the cooperative parameter (ω) of the binding. These results indicate that cofilin binding enhances additional cofilin binding in the vicinity of the initial binding site (ca. 24 subunits), but it does not affect the off-rate, which could be the molecular mechanism of the cooperative binding of cofilin to actin filaments.

Solution NMR structures of homeodomains from human proteins ALX4, ZHX1, and CASP8AP2 contribute to the structural coverage of the Human Cancer Protein Interaction Network

High-quality solution NMR structures of three homeodomains from human proteins ALX4, ZHX1 and CASP8AP2 were solved. These domains were chosen as targets of a biomedical theme project pursued by the Northeast Structural Genomics Consortium. This project focuses on increasing the structural coverage of human proteins associated with cancer.

Cephalopod eye evolution was modulated by the acquisition of Pax-6 splicing variants

Previous studies have reported that the developmental processes of vertebrate eyes are controlled by four Pax-6 splicing variants, each modulating different downstream genes, whereas those of insect eyes are controlled by duplicated Pax-6 genes. Cephalopods belong to the Protostomes but possess a camera-type eye similar to those in vertebrates. We examined Pax-6 variations in the squid and found five types of Pax-6 splicing variants but no duplication of the Pax-6 gene. In the five splicing variants, the splicing patterns were produced by the combination of two additional exons to the ortholog and one jettisoned exon containing most of the Homeobox domain (HD). These five variants show spatio-temporal patterns of gene expression during development in the squid. Our study suggests that cephalopods acquired Pax-6 splicing variants independent of those in vertebrates and that these variants were similarly utilized in the development of the squid eye.

Solution structure of CEH-37 homeodomain of the nematode Caenorhabditis elegans

The nematode Caenorhabditis elegans protein CEH-37 belongs to the paired OTD/OTX family of homeobox-containing homeodomain proteins. CEH-37 shares sequence similarity with homeodomain proteins, although it specifically binds to double-stranded C. elegans telomeric DNA, which is unusual to homeodomain proteins. Here, we report the solution structure of CEH-37 homeodomain and molecular interaction with double-stranded C. elegans telomeric DNA using nuclear magnetic resonance (NMR) spectroscopy. NMR structure shows that CEH-37 homeodomain is composed of a flexible N-terminal region and three α-helices with a helix-turn-helix (HTH) DNA binding motif. Data from size-exclusion chromatography and fluorescence spectroscopy reveal that CEH-37 homeodomain interacts strongly with double-stranded C. elegans telomeric DNA. NMR titration experiments identified residues responsible for specific binding to nematode double-stranded telomeric DNA. These results suggest that C. elegans homeodomain protein, CEH-37 could play an important role in telomere function via DNA binding.

The solution structures of two prophage homologues of the bacteriophage λ Ea8.5 protein reveal a newly discovered hybrid homeodomain/zinc-finger fold

A cluster of genes in the exoxis region of bacteriophage λ are capable of inhibiting the initiation of DNA synthesis in Escherichia coli. The most indispensible gene in this region is ea8.5. Here, we report the nuclear magnetic resonance structures of two ea8.5 orthologs from enteropathogenic E. coli and Pseudomonas putida prophages. Both proteins are characterized by a fused homeodomain/zinc-finger fold that escaped detection by primary sequence search methods. While these folds are both associated with a nucleic acid binding function, the amino acid composition suggests otherwise, leading to the possibility that Ea8.5 associates with other viral and host proteins.

Beamline AR-NW12A: high-throughput beamline for macromolecular crystallography at the Photon Factory

AR-NW12A is an in-vacuum undulator beamline optimized for high-throughput macromolecular crystallography experiments as one of the five macromolecular crystallography (MX) beamlines at the Photon Factory. This report provides details of the beamline design, covering its optical specifications, hardware set-up, control software, and the latest developments for MX experiments. The experimental environment presents state-of-the-art instrumentation for high-throughput projects with a high-precision goniometer with an adaptable goniometer head, and a UV-light sample visualization system. Combined with an efficient automounting robot modified from the SSRL SAM system, a remote control system enables fully automated and remote-access X-ray diffraction experiments.

Investigating dynamic and energetic determinants of protein nucleic acid recognition: analysis of the zinc finger zif268-DNA complexes

BACKGROUND

Protein-DNA recognition underlies fundamental biological processes ranging from transcription to replication and modification. Herein, we present a computational study of the sequence modulation of internal dynamic properties and of intraprotein networks of aminoacid interactions that determine the stability and specificity of protein-DNA complexes.

RESULTS

To this aim, we apply novel theoretical approaches to analyze the dynamics and energetics of biological systems starting from MD trajectories. As model system, we chose different sequences of Zinc Fingers (ZF) of the Zif268 family bound with different sequences of DNA. The complexes differ for their experimental stability properties, but share the same overall 3 D structure and do not undergo structural modifications during the simulations. The results of our analysis suggest that the energy landscape for DNA binding may be populated by dynamically different states, even in the absence of major conformational changes. Energetic couplings between residues change in response to protein and/or DNA sequence variations thus modulating the selectivity of recognition and the relative importance of different regions for binding.

CONCLUSIONS

The results show differences in the organization of the intra-protein energy-networks responsible for the stabilization of the protein conformations recognizing and binding DNA. These, in turn, are reflected into different modulation of the ZF's internal dynamics. The results also show a correlation between energetic and dynamic properties of the different proteins and their specificity/selectivity for DNA sequences. Finally, a dynamic and energetic model for the recognition of DNA by Zinc Fingers is proposed.