DNA length tunes the fluidity of DNA-based condensates

Living organisms typically store their genomic DNA in a condensed form. Mechanistically, DNA condensation can be driven by macromolecular crowding, multivalent cations, or positively charged proteins. At low DNA concentration, condensation triggers the conformational change of individual DNA molecules into a compacted state, with distinct morphologies. Above a critical DNA concentration, condensation goes along with phase separation into a DNA-dilute and a DNA-dense phase. The latter DNA-dense phase can have different material properties and has been reported to be rather liquid-like or solid-like depending on the characteristics of the DNA and the solvent composition. Here, we systematically assess the influence of DNA length on the properties of the resulting condensates. We show that short DNA molecules with sizes below 1 kb can form dynamic liquid-like assemblies when condensation is triggered by polyethylene glycol and magnesium ions, binding of linker histone H1, or nucleosome reconstitution in combination with linker histone H1. With increasing DNA length, molecules preferentially condense into less dynamic more solid-like assemblies, with phage λ-DNA with 48.5 kb forming mostly solid-like assemblies under the conditions assessed here. The transition from liquid-like to solid-like condensates appears to be gradual, with DNA molecules of roughly 1–10 kb forming condensates with intermediate properties. Titration experiments with linker histone H1 suggest that the fluidity of condensates depends on the net number of attractive interactions established by each DNA molecule. We conclude that DNA molecules that are much shorter than a typical human gene are able to undergo liquid-liquid phase separation, whereas longer DNA molecules phase separate by default into rather solid-like condensates. We speculate that the local distribution of condensing factors can modulate the effective length of chromosomal domains in the cell. We anticipate that the link between DNA length and fluidity established here will improve our understanding of biomolecular condensates involving DNA.     

Linker scanning mutagenesis of the internal ribosome entry site of poliovirus RNA.

The initiation of cap-independent translation of poliovirus mRNA occurs as a result of ribosome entry at an internal site(s) within the 5' noncoding region. A series of linker scanning mutations was constructed to define the genetic determinants of RNA-protein interactions that lead to high-fidelity translation of this unusual viral mRNA. The mutations are located within two distinct stem-loop structures in the 5' noncoding region of poliovirus RNA that constitute a major portion of a putative internal ribosome entry site. On the basis of our data derived from genetic and biochemical assays, the stability of one of the stem-loop structures appears to be essential for translation initiation via internal binding of ribosomes. However, the second stem-loop structure may function in a manner that requires base pairing and proper spacing between specific nucleotide sequences. By employing RNA electrophoretic mobility shift assays, an RNA-protein interaction was detected for this latter stem-loop structure that does not occur in RNAs containing mutations which perturb the predicted hairpin structure. Analysis of in vivo-selected virus revertants, in combination with mobility shift assays, suggests that extensive genetic rearrangement can lead to restoration of 5' noncoding region functions, possibly by the repositioning of specific RNA sequence or structure motifs.     

A 15N NMR mobility study on the dicalcium P43M calbindin D9k and its mono-La3+-substituted form

Calbindin D(9k) is a dicalcium binding protein consisting of two helix-loop-helix EF-hand motifs joined together by a flexible linker region where one metal ion can bind to each of the two loops. A proline residue at position 43 in the linker region displays cis-trans isomerism in the wild-type (WT) protein. Such isomerism appeared to be removed by substituting the proline with a glycine or a methionine in the P43G or P43M mutant. We have extended the available mobility studies on the P43M mutant through amide (15)N R(1), R(2), and R(1)(rho)() measurements. This has revealed unexpected conformational equilibria on the millisecond time scale involving residues 38, 42-44, and 46 in the linker region and residues 18 and 19 in calcium binding site I with similar energy barriers. These data are discussed in comparison with those available for the WT, as well as the apo-, mono-, and disubstituted P43G mutant. Quantification of water-amide proton exchange rates using saturation transfer and qualitative application of (15)N-(CLEANEX-PM)-FHSQC shows the values are in agreement with high mobility for the above-mentioned residues. Cross correlation between N-H dipole-dipole relaxation and (15)N CSA relaxation indicates that some of these mobility differences may extend to the sub-nanosecond time scale. Similar data were also obtained for the derivative where the calcium ion in the C-terminal loop was replaced with lanthanum. The results presented here show that, contrary to expectations, there are significant differences in dynamics between the dicalcium state of P43G and P43M and that these differences are not confined to the flexible linker region containing the point mutation. They also demonstrate that substitution of a lanthanide ion for calcium, which is a common procedure, does not significantly alter the mobility of the native protein.     

Linker Dependence of Avidity in Multivalent Interactions Between Disordered Proteins

Multidomain proteins often interact through several independent binding sites connected by disordered linkers. The architecture of such linkers affects avidity by modulating the effective concentration of intramolecular binding. The linker dependence of avidity has been estimated theoretically using simple physical models, but such models have not been tested experimentally because the effective concentrations could not be measured directly. We have developed a model system for bivalent protein interactions connected by disordered linkers, where the effective concentration can be measured using a competition experiment. We characterized the bivalent protein interactions kinetically and thermodynamically for a variety of linker lengths and interaction strengths. In total, this allowed us to critically assess the existing theoretical models of avidity in disordered, multivalent interactions. As expected, the onset of avidity occurs when the effective concentration reached the dissociation constant of the weakest interaction. Avidity decreased monotonously with linker length, but only by a third of what is predicted by theoretical models. We suggest that the length dependence of avidity is attenuated by compensating mechanisms such as linker interactions or entanglement. The direct role of linkers in avidity suggests they provide a generic mechanism for allosteric regulation of disordered, multivalent proteins.     

A specific role for the C-terminal region of the Poly(A)-binding protein in mRNA decay

mRNA poly(A) tails affect translation, mRNA export and mRNA stability, with translation initiation involving a direct interaction between eIF4G and the poly(A)-binding protein Pab1. The latter factor contains four RNA recognition motifs followed by a C-terminal region composed of a linker and a PABC domain. We show here that yeast mutants lacking the C-terminal domains of Pab1 display specific synthetic interactions with mutants in the 5′-3′ mRNA decay pathway. Moreover, these mutations impair mRNA decay in vivo without significantly affecting mRNA export or translation. Inhibition of mRNA decay occurs through slowed deadenylation. In vitro analyses demonstrate that removal of the Pab1 linker domain directly interferes with the ability of the Pop2–Ccr4 complex to deadenylate the Pab1-bound poly(A). Binding assays demonstrate that this results from a modulation of poly(A) packaging by the Pab1 linker region. Overall, our results demonstrate a direct involvement of Pab1 in mRNA decay and reveal the modular nature of this factor, with different domains affecting various cellular processes. These data suggest new models involving the modulation of poly(A) packaging by Pab1 to control mRNA decay.     

The influence of protein-lipid interactions on the order-disorder conformational transitions of the hydrocarbon chain.

The phases of simple systems involving one type of protein (lysozyme or cytochrome c) and one type of lipid (phosphatidic acid) have been characterized by X-ray crystallography, chemical analysis and spin-labeling technique as a function of temperature. They are of the lamellar type with alternative protein monolayers and lipid bilayers. According to the pH, two types of lamellar phases are obtained, one where the lipid-protein interactions are mainly hydrophobic, the other where they are electrostatic. In both cases, a phase transition occurs as temperature is lowered, between a high temperature phase, where all the lipids are in the liquid-like state, and another phase where some lipid chains are rigid. In the case of the phases with electrostatic interaction, it is shown that the onset of the order-disorder transition is shifted towards low temperature as compared with the homologous lipid-water phase and that the protein content of the phase decreases as the ratio of the liquid to rigid hydrocarbon chains decreases. This leads us to suggest that in the systems studied in this work the proteins interact only with lipid in the liquid-like state. In the case of the phases with hydrophobic interaction, it is shown that the extent of hydrophobic interaction between protein and lipid increases as the unsaturation of the hydrocarbon chains increases. The onset of the order-disorder transition shows a greater shift towards low temperature than the one observed in the case of the phase with electrostatic interaction.     

Long-range correlations improve understanding of the influence of network structure on contact dynamics

Models of infectious diseases are characterized by a phase transition between extinction and persistence. A challenge in contemporary epidemiology is to understand how the geometry of a host’s interaction network influences disease dynamics close to the critical point of such a transition. Here we address this challenge with the help of moment closures. Traditional moment closures, however, do not provide satisfactory predictions close to such critical points. We therefore introduce a new method for incorporating longer-range correlations into existing closures. Our method is technically simple, remains computationally tractable and significantly improves the approximation’s performance. Our extended closures thus provide an innovative tool for quantifying the influence of interaction networks on spatially or socially structured disease dynamics. In particular, we examine the effects of a network’s clustering coefficient, as well as of new geometrical measures, such as a network’s square clustering coefficients. We compare the relative performance of different closures from the literature, with or without our long-range extension. In this way, we demonstrate that the normalized version of the Bethe approximation-extended to incorporate long-range correlations according to our method-is an especially good candidate for studying influences of network structure. Our numerical results highlight the importance of the clustering coefficient and the square clustering coefficient for predicting disease dynamics at low and intermediate values of transmission rate, and demonstrate the significance of path redundancy for disease persistence.     

Modulation of assembly of TDP-43 low-complexity domain by heparin: From droplets to amyloid fibrils

TAR DNA-binding protein 43 (TDP-43) is an RNA-regulating protein that carries out many cellular functions through liquid-liquid phase separation (LLPS). The LLPS of TDP-43 is mediated by its C-terminal low-complexity domain (TDP43-LCD) corresponding to the region 267–414. In neurodegenerative disorders amyotrophic lateral sclerosis and frontotemporal dementia, pathological inclusions of the TDP-43 are found that are rich in the C-terminal fragments of ～25 and ～35 kDa, of which TDP43-LCD is a part. Thus, understanding the assembly process of TDP43-LCD is essential, given its involvement in the formation of both functional liquid-like assemblies and solid- or gel-like pathological aggregates. Here, we show that the solution pH and salt modulate TDP43-LCD LLPS. A gradual reduction in the pH below its isoelectric point of 9.8 results in a monotonic decrease of TDP43-LCD LLPS due to charge-charge repulsion between monomers, while at pH 6 and below no LLPS was observed. The addition of heparin to TDP43-LCD solution at pH 6, at a 1:2 heparin-to-TDP43-LCD molar ratio, promotes TDP43-LCD LLPS, while at higher concentration, it disrupts LLPS through a reentrant phase transition. Upon incubation at pH 6, TDP43-LCD undergoes gelation without phase separation. However, in the reentrant regime in the presence of a high heparin concentration, it forms thick amyloid aggregates that are significantly more SDS resistant than the gel. The results indicate that the material nature of the TDP43-LCD assembly products can be modulated by heparin which is significant in the context of liquid-to-solid phase transition observed in TDP-43 proteinopathies. Our findings are also crucial in relation to similar transitions that could occur due to alteration in the molecular level interactions among various multivalent biomolecules involving other LCDs and RNAs.     

The influence of protein-lipid interactions on the order-disorder conformational transitions of the hydrocarbon chain

The phases of simple systems involving one type of protein (lysozyme or cytochrome $\text{c}$) and one type of lipid (phosphatidic acid) have been characterized by X-ray crystallography, chemical analysis and spin-labeling technique as a function of temperature. They are of the lamellar type with alternative protein monolayers and lipid bilayers. According to the pH, two types of lamellar phases are obtained, one where the lipid-protein interactions are mainly hydrophobic, the other where they are electrostatic. In both cases, a phase transition occurs as temperature is lowered, between a high temperature phase, where all the lipids are in the liquid-like state, and another phase where some lipid chains are rigid. In the case of the phases with electrostatic interaction, it is shown that the onset of the order-disorder transition is shifted towards low temperature as compared with the homologous lipid-water phase and that the protein content of the phase decreases as the ratio of the liquid to rigid hydrocarbon chains decreases. This leads us to suggest that in the systems studied in this work the proteins interact only with lipid in the liquid-like state. In the case of the phases with hydrophobic interaction, it is shown that the extent of hydrophobic interaction between protein and lipid increases as the unsaturation of the hydrocarbon chains increases. The onset of the order-disorder transition shows a greater shift towards low temperture than the one observed in the case of the phase with electrostatic interaction.     

Different cation sensitivities and binding site domains of Na+-Ca2+-K+ and Na+-Ca2+ exchangers

We examined inhibitory effects of external multivalent cations Ni(2+), Co(2+), Cd(2+), La(3+), Mg(2+), and Mn(2+) on reverse-mode exchange of the K(+)-dependent Na(+)/Ca(2+) exchanger NCKX2 and the K(+)-independent exchanger NCX1 expressed in CCL-39 cells by measuring the rate of Ca(2+) uptake with radioisotope tracer and electrophysiological techniques. The apparent affinities for block of Ca(2+) uptake by multivalent cations was higher in NCKX2 than NCX1, and the rank order of inhibitory potencies among these cations was different. Additional experiments also showed that external Li(+) stimulated reverse-mode exchange by NCX1, but not NCKX2 in the presence of 5 mM K(+). Thus, both exchangers exhibited differential sensitivities to not only K(+) but also many other external cations. We attempted to locate the putative binding sites within the alpha motifs for multivalent cations by site-directed mutagenesis experiments. The cation affinities of NCKX2 were altered by mutations of amino acid residues in the alpha-1 motif, but not by mutations in the alpha-2 motif. These results contrast with those for NCX1 where mutations in both alpha-1 and alpha-2 motifs have been shown previously to affect cation affinities. Susceptibility tests with sulfhydryl alkylating agents suggested that the alpha-1 and alpha-2 motifs are situated extracellularly and intracellularly, respectively, in both exchangers. A topological model is proposed in which the extracellular-facing alpha-1 motif forms an external cation binding site that includes key residues N203, G207C, and I209 in NCKX2, while both alpha-1 and alpha-2 motifs together form the binding sites in NCX1.     

The molecular basis of selective DNA binding by the BRG1 AT-hook and bromodomain

The ATP-dependent BAF chromatin remodeling complex plays a critical role in gene regulation by modulating chromatin architecture, and is frequently mutated in cancer. Indeed, subunits of the BAF complex are found to be mutated in >20% of human tumors. The mechanism by which BAF properly navigates chromatin is not fully understood, but is thought to involve a multivalent network of histone and DNA contacts. We previously identified a composite domain in the BRG1 ATPase subunit that is capable of associating with both histones and DNA in a multivalent manner. Mapping the DNA binding pocket revealed that it contains several cancer mutations. Here, we utilize SELEX-seq to investigate the DNA specificity of this composite domain and NMR spectroscopy and molecular modelling to determine the structural basis of DNA binding. Finally, we demonstrate that cancer mutations in this domain alter the mode of DNA association.     

Alpha-helical linker of an artificial 6-zinc finger peptide contributes to selective DNA binding to a discontinuous recognition sequence

Although many zinc finger motifs have been developed to recognize specific DNA triplets, a rational way to selectively skip a particular non-recognized gap in the DNA sequence has never been established. We have now created a 6-zinc finger peptide with an alpha-helix linker, Sp1ZF6(EAAAR)4, which selectively binds to the discontinuous recognition sites in the same phase (10 bp gap) against the opposite phase (5 bp gap) of the DNA helix. The linker peptide (EAAAR)4 forms an alpha-helix structure stabilized by salt bridges, and the helical length is estimated to be about 30 A, corresponding to that of the 10 bp DNA. The gel shift assays demonstrate that Sp1ZF6(EAAAR)4 preferably binds to the 10 bp-gapped target rather than the 5 bp-gapped target. The CD spectra show that the alpha-helical content of the (EAAAR)4 linker is higher in the complex with the 10 bp-gapped target than in the complex with the 5 bp-gapped target. The present results indicate that the alpha-helical linker is suitable for binding to the recognition sites in the same phase and that the linker induces the loss of binding affinity to the recognition sites with the opposite phase. The engineering of a helix-structured linker in the 6-zinc finger peptides should be one of the most promising approaches for selectively targeting discontinuous recognition sites depending on their phase situations.     

GGA1 interacts with the adaptor protein AP-1 through a WNSF sequence in its hinge region

The Golgi-associated gamma-adaptin-related ADP-ribosylation factor-binding proteins (GGAs) are critical components of the transport machinery that mediates the trafficking of the mannose 6-phosphate receptors and associated cargo from the trans-Golgi network to the endosomes. The GGAs colocalize in vivo with the clathrin adaptor protein AP-1 and bind to AP-1 in vitro, suggesting that the two proteins may cooperate in packaging the mannose 6-phosphate receptors into clathrin-coated vesicles at the trans-Golgi network. Here, we demonstrate that the sequence, (382)WNSF(385), in the hinge region of GGA1 mediates its interaction with the AP-1 gamma-ear. The Trp and Phe constitute critical amino acids in this interaction. The binding of Rabaptin5 to the AP-1 gamma-ear, which occurs through a FXXPhi motif, is inhibited by a peptide encoding the GGA1 (382)WNSF(385) sequence. Moreover, mutations in the AP-1 gamma-ear that abolish its interaction with Rabaptin5 also preclude its association with GGA1. These results suggest that the GGA1 WXXF-type and Rabaptin5 FXXPhi-type motifs bind to the same or highly overlapping sites in the AP-1 gamma-ear. This binding is modulated by residues adjacent to the core motifs.     

Quantifying intramolecular binding in multivalent interactions: a structure-based synergistic study on Grb2-Sos1 complex

Numerous signaling proteins use multivalent binding to increase the specificity and affinity of their interactions within the cell. Enhancement arises because the effective binding constant for multivalent binding is larger than the binding constants for each individual interaction. We seek to gain both qualitative and quantitative understanding of the multivalent interactions of an adaptor protein, growth factor receptor bound protein-2 (Grb2), containing two SH3 domains interacting with the nucleotide exchange factor son-of-sevenless 1 (Sos1) containing multiple polyproline motifs separated by flexible unstructured regions. Grb2 mediates the recruitment of Sos1 from the cytosol to the plasma membrane where it activates Ras by inducing the exchange of GDP for GTP. First, using a combination of evolutionary information and binding energy calculations, we predict an additional polyproline motif in Sos1 that binds to the SH3 domains of Grb2. This gives rise to a total of five polyproline motifs in Sos1 that are capable of binding to the two SH3 domains of Grb2. Then, using a hybrid method combining molecular dynamics simulations and polymer models, we estimate the enhancement in local concentration of a polyproline motif on Sos1 near an unbound SH3 domain of Grb2 when its other SH3 domain is bound to a different polyproline motif on Sos1. We show that the local concentration of the Sos1 motifs that a Grb2 SH3 domain experiences is approximately 1000 times greater than the cellular concentration of Sos1. Finally, we calculate the intramolecular equilibrium constants for the crosslinking of Grb2 on Sos1 and use thermodynamic modeling to calculate the stoichiometry. With these equilibrium constants, we are able to predict the distribution of complexes that form at physiological concentrations. We believe this is the first systematic analysis that combines sequence, structure, and thermodynamic analyses to determine the stoichiometry of the complexes that are dominant in the cellular environment.     

Rac1 protein rescues neurite retraction caused by G2019S leucine-rich repeat kinase 2 (LRRK2)

Mutations in leucine-rich repeat kinase 2 (LRRK2) are currently the most common genetic cause of familial late-onset Parkinson disease, which is clinically indistinguishable from idiopathic disease. The most common pathological mutation in LRRK2, G2019S LRRK2, is known to cause neurite retraction. However, molecular mechanisms underlying regulation of neurite length by LRRK2 are unknown. Here, we demonstrate a novel interaction between LRRK2 and the Rho GTPase, Rac1, which plays a critical role in actin cytoskeleton remodeling necessary for the maintenance of neurite morphology. LRRK2 binds strongly to endogenous or expressed Rac1, while showing weak binding to Cdc42 and no binding to RhoA. Co-expression with LRRK2 increases Rac1 activity, as shown by increased binding to the p21-activated kinase, which modulates actin cytoskeletal dynamics. LRRK2 constructs carrying mutations that inactivate the kinase or GTPase activities do not activate Rac1. Interestingly, LRRK2 does not increase levels of membrane-bound Rac1 but dramatically changes the cellular localization of Rac1, causing polarization, which is augmented further when LRRK2 is co-expressed with constitutively active Rac1. Four different disease-related mutations in LRRK2 altered binding to Rac1, with the G2019S and R1441C LRRK2 mutations attenuating Rac1 binding and the Y1699C and I2020T LRRK2 mutations increasing binding. Co-expressing Rac1 in SH-SY5Y cells rescues the G2019S mutant phenotype of neurite retraction. We hypothesize that pathological mutations in LRRK2 attenuates activation of Rac1, causing disassembly of actin filaments, leading to neurite retraction. The interactions between LRRK2 and Rho GTPases provide a novel pathway through which LRRK2 might modulate cellular dynamics and contribute to the pathophysiology of Parkinson disease.     

Mapping the interaction surface of linker histone H1(0) with the nucleosome of native chromatin in vivo

H1 linker histones stabilize the nucleosome, limit nucleosome mobility and facilitate the condensation of metazoan chromatin. Here, we have combined systematic mutagenesis, measurement of in vivo binding by photobleaching microscopy, and structural modeling to determine the binding geometry of the globular domain of the H1(0) linker histone variant within the nucleosome in unperturbed, native chromatin in vivo. We demonstrate the existence of two distinct DNA-binding sites within the globular domain that are formed by spatial clustering of multiple residues. The globular domain is positioned via interaction of one binding site with the major groove near the nucleosome dyad. The second site interacts with linker DNA adjacent to the nucleosome core. Multiple residues bind cooperatively to form a highly specific chromatosome structure that provides a mechanism by which individual domains of linker histones interact to facilitate chromatin condensation.