Condensin extrudes DNA loops in steps up to hundreds of base pairs that are generated by ATP binding events

The condensin SMC protein complex organizes chromosomal structure by extruding loops of DNA. Its ATP-dependent motor mechanism remains unclear but likely involves steps associated with large conformational changes within the ∼50 nm protein complex. Here, using high-resolution magnetic tweezers, we resolve single steps in the loop extrusion process by individual yeast condensins. The measured median step sizes range between 20–40 nm at forces of 1.0–0.2 pN, respectively, comparable with the holocomplex size. These large steps show that, strikingly, condensin typically reels in DNA in very sizeable amounts with ∼200 bp on average per single extrusion step at low force, and occasionally even much larger, exceeding 500 bp per step. Using Molecular Dynamics simulations, we demonstrate that this is due to the structural flexibility of the DNA polymer at these low forces. Using ATP-binding-impaired and ATP-hydrolysis-deficient mutants, we find that ATP binding is the primary step-generating stage underlying DNA loop extrusion. We discuss our findings in terms of a scrunching model where a stepwise DNA loop extrusion is generated by an ATP-binding-induced engagement of the hinge and the globular domain of the SMC complex.     

Euler buckling and nonlinear kinking of double-stranded DNA

The bending stiffness of double-stranded DNA (dsDNA) at high curvatures is fundamental to its biological activity, yet this regime has been difficult to probe experimentally, and literature results have not been consistent. We created a ‘molecular vise’ in which base-pairing interactions generated a compressive force on sub-persistence length segments of dsDNA. Short dsDNA strands (<41 base pairs) resisted this force and remained straight; longer strands became bent, a phenomenon called ‘Euler buckling’. We monitored the buckling transition via Förster Resonance Energy Transfer (FRET) between appended fluorophores. For low-to-moderate concentrations of monovalent salt (up to ∼150 mM), our results are in quantitative agreement with the worm-like chain (WLC) model of DNA elasticity, without the need to invoke any ‘kinked’ states. Greater concentrations of monovalent salts or 1 mM Mg²⁺ induced an apparent softening of the dsDNA, which was best accounted for by a kink in the region of highest curvature. We tested the effects of all single-nucleotide mismatches on the DNA bending. Remarkably, the propensity to kink correlated with the thermodynamic destabilization of the mismatched DNA relative the perfectly complementary strand, suggesting that the kinked state is locally melted. The molecular vise is exquisitely sensitive to the sequence-dependent linear and nonlinear elastic properties of dsDNA.     

Double Strand Break Unwinding and Resection by the Mycobacterial Helicase-Nuclease AdnAB in the Presence of Single Strand DNA-binding Protein (SSB)

Mycobacterial AdnAB is a heterodimeric DNA helicase-nuclease and 3′ to 5′ DNA translocase implicated in the repair of double strand breaks (DSBs). The AdnA and AdnB subunits are each composed of an N-terminal motor domain and a C-terminal nuclease domain. Inclusion of mycobacterial single strand DNA-binding protein (SSB) in reactions containing linear plasmid dsDNA allowed us to study the AdnAB helicase under conditions in which the unwound single strands are coated by SSB and thereby prevented from reannealing or promoting ongoing ATP hydrolysis. We found that the AdnAB motor catalyzed processive unwinding of 2.7–11.2-kbp linear duplex DNAs at a rate of ∼250 bp s⁻¹, while hydrolyzing ∼5 ATPs per bp unwound. Crippling the AdnA phosphohydrolase active site did not affect the rate of unwinding but lowered energy consumption slightly, to ∼4.2 ATPs bp⁻¹. Mutation of the AdnB phosphohydrolase abolished duplex unwinding, consistent with a model in which the “leading” AdnB motor propagates a Y-fork by translocation along the 3′ DNA strand, ahead of the “lagging” AdnA motor domain. By tracking the resection of the 5′ and 3′ strands at the DSB ends, we illuminated a division of labor among the AdnA and AdnB nuclease modules during dsDNA unwinding, whereby the AdnA nuclease processes the unwound 5′ strand to liberate a short oligonucleotide product, and the AdnB nuclease incises the 3′ strand on which the motor translocates. These results extend our understanding of presynaptic DSB processing by AdnAB and engender instructive comparisons with the RecBCD and AddAB clades of bacterial helicase-nuclease machines.     

A novel mode of nuclease action is revealed by the bacterial Mre11/Rad50 complex

The Mre11/Rad50 complex is a central player in various genome maintenance pathways. Here, we report a novel mode of nuclease action found for the Escherichia coli Mre11/Rad50 complex, SbcC₂/D₂ complex (SbcCD). SbcCD cuts off the top of a cruciform DNA by making incisions on both strands and continues cleaving the dsDNA stem at ∼10-bp intervals. Using linear-shaped DNA substrates, we observed that SbcCD cleaved dsDNA using this activity when the substrate was 110 bp long, but that on shorter substrates the cutting pattern was changed to that predicted for the activity of a 3′-5′ exonuclease. Our results suggest that SbcCD processes hairpin and linear dsDNA ends with this novel DNA end-dependent binary endonuclease activity in response to substrate length rather than using previously reported activities. We propose a model for this mode of nuclease action, which provides new insight into SbcCD activity at a dsDNA end.     

Revising a Personal Genome by Comparing and Combining Data from Two Different Sequencing Platforms

For the robust practice of genomic medicine, sequencing results must be compatible, regardless of the sequencing technologies and algorithms used. Presently, genome sequencing is still an imprecise science and is complicated by differences in the chemistry, coverage, alignment, and variant-calling algorithms. We identified ∼3.33 million single nucleotide variants (SNVs) and ∼3.62 million SNVs in the SJK genome using SOLiD and Illumina data, respectively. Approximately 3 million SNVs were concordant between the two platforms while 68,532 SNVs were discordant; 219,616 SNVs were SOLiD-specific and 516,080 SNVs were Illumina-specific (i.e., platform-specific). Concordant, discordant, and platform-specific SNVs were further analyzed and characterized. Overall, a large portion of heterozygous SNVs that were discordant with genotyping calls of single nucleotide polymorphism chips were highly confident. Approximately 70% of the platform-specific SNVs were located in regions containing repetitive sequences. Such platform-specificity may arise from differences between platforms, with regard to read length (36 bp and 72 bp vs. 50 bp), insert size (∼100–300 bp vs. ∼1–2 kb), sequencing chemistry (sequencing-by-synthesis using single nucleotides vs. ligation-based sequencing using oligomers), and sequencing quality. When data from the two platforms were merged for variant calling, the proportion of callable regions of the reference genome increased to 99.66%, which was 1.43% higher than the average callability of the two platforms, representing ∼40 million bases. In this study, we compared the differences in sequencing results between two sequencing platforms. Approximately 90% of the SNVs were concordant between the two platforms, yet ∼10% of the SNVs were either discordant or platform-specific, indicating that each platform had its own strengths and weaknesses. When data from the two platforms were merged, both the overall callability of the reference genome and the overall accuracy of the SNVs improved, demonstrating the likelihood that a re-sequenced genome can be revised using complementary data.     

Identification of Allosteric Disulfides from Prestress Analysis

Disulfide bonds serve to form physical cross-links between residues in protein structures, thereby stabilizing the protein fold. Apart from this purely structural role, they can also be chemically active, participating in redox reactions, and they may even potentially act as allosteric switches controlling protein functions. Specific types of disulfide bonds have been identified in static protein structures from their distinctive pattern of dihedral bond angles, and the allosteric function of such bonds is purported to be related to the torsional strain they store. Using all-atom molecular-dynamics simulations for ∼700 disulfide bonded proteins, we analyzed the intramolecular mechanical forces in 20 classes of disulfide bonds. We found that two particular classes, the −RHStaple and the −/+RHHook disulfides, are indeed more stressed than other disulfide bonds, but the stress is carried primarily by stretching of the S-S bond and bending of the neighboring bond angles, rather than by dihedral torsion. This stress corresponds to a tension force of magnitude ∼200 pN, which is balanced by repulsive van der Waals interactions between the cysteine Cα atoms. We confirm stretching of the S-S bond to be a general feature of the −RHStaples and the −/+RHHooks by analyzing ∼20,000 static protein structures. Given that forced stretching of S-S bonds is known to accelerate their cleavage, we propose that prestress of allosteric disulfide bonds has the potential to alter the reactivity of a disulfide, thereby allowing us to readily switch between functional states.     

A Highlights from MBoC Selection: Quantitative Analysis of the Mechanism of Endocytic Actin Patch Assembly and Disassembly in Fission Yeast

We used quantitative confocal microscopy to measure the numbers of 16 proteins tagged with fluorescent proteins during assembly and disassembly of endocytic actin patches in fission yeast. The peak numbers of each molecule that accumulate in patches varied <30–50% between individual patches. The pathway begins with accumulation of 30–40 clathrin molecules, sufficient to build a hemisphere at the tip of a plasma membrane invagination. Thereafter precisely timed waves of proteins reach characteristic peak numbers: endocytic adaptor proteins (∼120 End4p and ∼230 Pan1p), activators of Arp2/3 complex (∼200 Wsp1p and ∼340 Myo1p) and ∼300 Arp2/3 complexes just ahead of a burst of actin assembly into short, capped and highly cross-linked filaments (∼7000 actins, ∼200 capping proteins, and ∼900 fimbrins). Coronin arrives last as all other components disperse upon patch internalization and movement over ∼10 s. Patch internalization occurs without recruitment of dynamins. Mathematical modeling, described in the accompanying paper (Berro et al., 2010, MBoC 21: 2803–2813), shows that the dendritic nucleation hypothesis can account for the time course of actin assembly into a branched network of several hundred filaments 100–200 nm long and that patch disassembly requires actin filament fragmentation in addition to depolymerization from the ends.     

Spontaneous sharp bending of DNA: role of melting bubbles

The role of centrally located and distributed base pair mismatches (‘melting bubbles’) on localized bending and stiffness of short dsDNA fragments is evaluated using time-dependent fluorescence lifetime measurements. Distributed melting bubbles are found to induce larger bending angles and decreased levels of stiffness in DNA than centrally located ones of comparable overall size. Our results indicate that spontaneous local opening-up of the DNA duplex could facilitate sharp bending of short DNA strands even in the absence of DNA binding proteins. We also find that the occurrence of two closely spaced melting bubbles will generally be favored when a large energetic barrier must be overcome in forming the desired bent DNA structure.     

Novel nucleosomal particles containing core histones and linker DNA but no histone H1

Eukaryotic chromosomal DNA is assembled into regularly spaced nucleosomes, which play a central role in gene regulation by determining accessibility of control regions. The nucleosome contains ∼147 bp of DNA wrapped ∼1.7 times around a central core histone octamer. The linker histone, H1, binds both to the nucleosome, sealing the DNA coils, and to the linker DNA between nucleosomes, directing chromatin folding. Micrococcal nuclease (MNase) digests the linker to yield the chromatosome, containing H1 and ∼160 bp, and then converts it to a core particle, containing ∼147 bp and no H1. Sequencing of nucleosomal DNA obtained after MNase digestion (MNase-seq) generates genome-wide nucleosome maps that are important for understanding gene regulation. We present an improved MNase-seq method involving simultaneous digestion with exonuclease III, which removes linker DNA. Remarkably, we discovered two novel intermediate particles containing 154 or 161 bp, corresponding to 7 bp protruding from one or both sides of the nucleosome core. These particles are detected in yeast lacking H1 and in H1-depleted mouse chromatin. They can be reconstituted in vitro using purified core histones and DNA. We propose that these ‘proto-chromatosomes’ are fundamental chromatin subunits, which include the H1 binding site and influence nucleosome spacing independently of H1.     

Structural mechanism of DNA recognition by the p204 HIN domain

The interferon gamma-inducible protein 16 (IFI16) and its murine homologous protein p204 function in non-sequence specific dsDNA sensing; however, the exact dsDNA recognition mechanisms of IFI16/p204, which harbour two HIN domains, remain unclear. In the present study, we determined crystal structures of p204 HINa and HINb domains, which are highly similar to those of other PYHIN family proteins. Moreover, we obtained the crystal structure of p204 HINab domain in complex with dsDNA and provided insights into the dsDNA binding mode. p204 HINab binds dsDNA mainly through α2 helix of HINa and HINb, and the linker between them, revealing a similar HIN:DNA binding mode. Both HINa and HINb are vital for HINab recognition of dsDNA, as confirmed by fluorescence polarization assays. Furthermore, a HINa dimerization interface was observed in structures of p204 HINa and HINab:dsDNA complex, which is involved in binding dsDNA. The linker between HINa and HINb reveals dynamic flexibility in solution and changes its direction at ∼90° angle in comparison with crystal structure of HINab:dsDNA complex. These structural information provide insights into the mechanism of DNA recognition by different HIN domains, and shed light on the unique roles of two HIN domains in activating the IFI16/p204 signaling pathway.     

Activation of Protein Kinase PKR Requires Dimerization-induced cis-Phosphorylation within the Activation Loop

Protein kinase R (PKR) functions in a plethora of cellular processes, including viral and cellular stress responses, by phosphorylating the translation initiation factor eIF2α. The minimum requirements for PKR function are homodimerization of its kinase and RNA-binding domains, and autophosphorylation at the residue Thr-446 in a flexible loop called the activation loop. We investigated the interdependence between dimerization and Thr-446 autophosphorylation using the yeast Saccharomyces cerevisiae model system. We showed that an engineered PKR that bypassed the need for Thr-446 autophosphorylation (PKR^T446∼P-bypass mutant) could function without a key residue (Asp-266 or Tyr-323) that is essential for PKR dimerization, suggesting that dimerization precedes and stimulates activation loop autophosphorylation. We also showed that the PKR^T446∼P-bypass mutant was able to phosphorylate eIF2α even without its RNA-binding domains. These two significant findings reveal that PKR dimerization and activation loop autophosphorylation are mutually exclusive yet interdependent processes. Also, we provide evidence that Thr-446 autophosphorylation during PKR activation occurs in a cis mechanism following dimerization.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Analysis of the Role of the C-Terminal Tail in the Regulation of the Epidermal Growth Factor Receptor

The ∼230-residue C-terminal tail of the epidermal growth factor receptor (EGFR) is phosphorylated upon activation. We examined whether this phosphorylation is affected by deletions within the tail and whether the two tails in the asymmetric active EGFR dimer are phosphorylated differently. We monitored autophosphorylation in cells using flow cytometry and found that the first ∼80 residues of the tail are inhibitory, as demonstrated previously. The entire ∼80-residue span is important for autoinhibition and needs to be released from both kinases that form the dimer. These results are interpreted in terms of crystal structures of the inactive kinase domain, including two new ones presented here. Deletions in the remaining portion of the tail do not affect autophosphorylation, except for a six-residue segment spanning Tyr 1086 that is critical for activation loop phosphorylation. Phosphorylation of the two tails in the dimer is asymmetric, with the activator tail being phosphorylated somewhat more strongly. Unexpectedly, we found that reconstitution of the transmembrane and cytoplasmic domains of EGFR in vesicles leads to a peculiar phenomenon in which kinase domains appear to be trapped between stacks of lipid bilayers. This artifactual trapping of kinases between membranes enhances an intrinsic functional asymmetry in the two tails in a dimer.     

Simple Elastic Network Models for Exhaustive Analysis of Long Double-Stranded DNA Dynamics with Sequence Geometry Dependence

Simple elastic network models of DNA were developed to reveal the structure-dynamics relationships for several nucleotide sequences. First, we propose a simple all-atom elastic network model of DNA that can explain the profiles of temperature factors for several crystal structures of DNA. Second, we propose a coarse-grained elastic network model of DNA, where each nucleotide is described only by one node. This model could effectively reproduce the detailed dynamics obtained with the all-atom elastic network model according to the sequence-dependent geometry. Through normal-mode analysis for the coarse-grained elastic network model, we exhaustively analyzed the dynamic features of a large number of long DNA sequences, approximately ∼150 bp in length. These analyses revealed positive correlations between the nucleosome-forming abilities and the inter-strand fluctuation strength of double-stranded DNA for several DNA sequences.     

Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms

Although satellite DNAs are well-explored components of heterochromatin and centromeres, little is known about emergence, dispersal and possible impact of comparably structured tandem repeats (TRs) on the genome-wide scale. Our bioinformatics analysis of assembled Tribolium castaneum genome disclosed significant contribution of TRs in euchromatic chromosomal arms and clear predominance of satellite DNA-typical 170 bp monomers in arrays of ≥5 repeats. By applying different experimental approaches, we revealed that the nine most prominent TR families Cast1–Cast9 extracted from the assembly comprise ∼4.3% of the entire genome and reside almost exclusively in euchromatic regions. Among them, seven families that build ∼3.9% of the genome are based on ∼170 and ∼340 bp long monomers. Results of phylogenetic analyses of 2500 monomers originating from these families show high-sequence dynamics, evident by extensive exchanges between arrays on non-homologous chromosomes. In addition, our analysis shows that concerted evolution acts more efficiently on longer than on shorter arrays. Efficient genome-wide distribution of nine TR families implies the role of transposition only in expansion of the most dispersed family, and involvement of other mechanisms is anticipated. Despite similarities in sequence features, FISH experiments indicate high-level compartmentalization of centromeric and euchromatic tandem repeats.