Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Abl2 repairs microtubules and phase separates with tubulin to promote microtubule nucleation

Abl family kinases are evolutionarily conserved regulators of cell migration and morphogenesis. Genetic experiments in Drosophila suggest that Abl family kinases interact functionally with microtubules to regulate axon guidance and neuronal morphogenesis. Vertebrate Abl2 binds to microtubules and promotes their plus-end elongation, both in vitro and in cells, but the molecular mechanisms by which Abl2 regulates microtubule (MT) dynamics are unclear. We report here that Abl2 regulates MT assembly via condensation and direct interactions with both the MT lattice and tubulin dimers. We find that Abl2 promotes MT nucleation, which is further facilitated by the ability of the Abl2 C-terminal half to undergo liquid-liquid phase separation (LLPS) and form co-condensates with tubulin. Abl2 binds to regions adjacent to MT damage, facilitates MT repair via fresh tubulin recruitment, and increases MT rescue frequency and lifetime. Cryo-EM analyses strongly support a model in which Abl2 engages tubulin C-terminal tails along an extended MT lattice conformation at damage sites to facilitate repair via fresh tubulin recruitment. Abl2Δ688-790, which closely mimics a naturally occurring splice isoform, retains binding to the MT lattice but does not bind tubulin, promote MT nucleation, or increase rescue frequency. In COS-7 cells, MT reassembly after nocodazole treatment is greatly slowed in Abl2 knockout COS-7 cells compared with wild-type cells, and these defects are rescued by re-expression of Abl2, but not Abl2Δ688-790. We propose that Abl2 locally concentrates tubulin to promote MT nucleation and recruits it to defects in the MT lattice to enable repair and rescue.

Negatively charged, intrinsically disordered regions can accelerate target search by DNA-binding proteins

In eukaryotes, many DNA/RNA-binding proteins possess intrinsically disordered regions (IDRs) with large negative charge, some of which involve a consecutive sequence of aspartate (D) or glutamate (E) residues. We refer to them as D/E repeats. The functional role of D/E repeats is not well understood, though some of them are known to cause autoinhibition through intramolecular electrostatic interaction with functional domains. In this work, we investigated the impacts of D/E repeats on the target DNA search kinetics for the high-mobility group box 1 (HMGB1) protein and the artificial protein constructs of the Antp homeodomain fused with D/E repeats of varied lengths. Our experimental data showed that D/E repeats of particular lengths can accelerate the target association in the overwhelming presence of non-functional high-affinity ligands ('decoys'). Our coarse-grained molecular dynamics (CGMD) simulations showed that the autoinhibited proteins can bind to DNA and transition into the uninhibited complex with DNA through an electrostatically driven induced-fit process. In conjunction with the CGMD simulations, our kinetic model can explain how D/E repeats can accelerate the target association process in the presence of decoys. This study illuminates an unprecedented role of the negatively charged IDRs in the target search process.

Competitive binding of TET1 and DNMT3A/B cooperates the DNA methylation pattern in human embryonic stem cells

DNMT3A/B and TET1 play indispensable roles in regulating DNA methylation that undergoes extensive reprogramming during mammalian embryogenesis. Yet the competitive and cooperative relationships between TET1 and DNMT3A/B remain largely unknown in the human embryonic stem cells. Here, we revealed that the main DNA-binding domain of TET1 contains more positive charges by using charge reduction of amino acid alphabet, followed by DNMT3A and DNMT3B. The genome-wide binding profiles showed that TET1 prefers binding to the proximal promoters and CpG islands compared with DNMT3A/B. Moreover, the binding regions of these three transcription factors can be divided into specific and co-binding regions. And a stronger inhibitory effect of DNMT3A on TET1 demethylation was observed in co-binding regions. Furthermore, we integrated TET1 knockout data to further discuss the competitive binding patterns of TET1 and DNMT3A/B. The lack of TET1 increased the occupation of DNMT3A/B at the specific binding regions of TET1 causing focal hypermethylation. The knockout of TET1 was also accompanied by a reduction of DNMT3A/B binding in the co-binding regions, further confirming the cooperative binding function between TET1 and DNMT3A/B. In conclusion, our studies found that the competitive binding of TET1 and DNMT3A/B cooperatively shapes the global DNA methylation pattern in human embryonic stem cells.

Negatively Charged Disordered Regions are Prevalent and Functionally Important Across Proteomes

Intrinsically disordered regions (IDRs) of proteins are often characterized by a high fraction of charged residues, but differ in their overall net charge and in the organization of the charged residues. The function-encoding information stored via IDR charge composition and organization remains elusive. Here, we aim to decipher the sequence-function relationship in IDRs by presenting a comprehensive bioinformatic analysis of the charge properties of IDRs in the human, mouse, and yeast proteomes. About 50% of the proteins comprise at least a single IDR, which is either positively or negatively charged. Highly negatively charged IDRs are longer and possess greater net charge per residue compared with highly positively charged IDRs. A striking difference between positively and negatively charged IDRs is the characteristics of the repeated units, specifically, of consecutive Lys or Arg residues (K/R repeats) and Asp or Glu (D/E repeats) residues. D/E repeats are found to be about five times longer than K/R repeats, with the longest found containing 49 residues. Long stretches of consecutive D and E are found to be more prevalent in nucleic acid-related proteins. They are less common in prokaryotes, and in eukaryotes their abundance increases with genome size. The functional role of D/E repeats and the profound differences between them and K/R repeats are discussed.

Fast Leaps between Millisecond Confinements Govern Ase1 Diffusion along Microtubules

Diffusion is the most fundamental mode of protein translocation within cells. Confined diffusion of proteins along the electrostatic potential constituted by the surface of microtubules, although modeled meticulously in molecular dynamics simulations, has not been experimentally observed in real-time. Here, interferometric scattering microscopy is used to directly visualize the movement of the microtubule-associated protein Ase1 along the microtubule surface at nanometer and microsecond resolution. Millisecond confinements of Ase1 and fast leaps between these positions of dwelling preferentially occurring along the microtubule protofilaments are resolved, revealing Ase1's mode of diffusive translocation along the microtubule's periodic surface. The derived interaction potential closely matches the tubulin-dimer periodicity and the distribution of the electrostatic potential on the microtubule lattice. It is anticipated that mapping the interaction landscapes for different proteins on microtubules, finding plausible energetic barriers of different positioning and heights, can provide valuable insights into regulating the dynamics of essential cytoskeletal processes, such as intracellular cargo trafficking, cell division, and morphogenesis, all of which rely on diffusive translocation of proteins along microtubules.

Dynamic Autoinhibition of the HMGB1 Protein via Electrostatic Fuzzy Interactions of Intrinsically Disordered Regions

Highly negatively charged segments containing only aspartate or glutamate residues ("D/E repeats") are found in many eukaryotic proteins. For example, the C-terminal 30 residues of the HMGB1 protein are entirely D/E repeats. Using nuclear magnetic resonance (NMR), fluorescence, and computational approaches, we investigated how the D/E repeats causes the autoinhibition of HMGB1 against its specific binding to cisplatin-modified DNA. By varying ionic strength in a wide range (40-900 mM), we were able to shift the conformational equilibrium between the autoinhibited and uninhibited states toward either of them to the full extent. This allowed us to determine the macroscopic and microscopic equilibrium constants for the HMGB1 autoinhibition at various ionic strengths. At a macroscopic level, a model involving the autoinhibited and uninhibited states can explain the salt concentration-dependent binding affinity data. Our data at a microscopic level show that the D/E repeats and other parts of HMGB1 undergo electrostatic fuzzy interactions, each of which is weaker than expected from the macroscopic autoinhibitory effect. This discrepancy suggests that the multivalent nature of the fuzzy interactions enables strong autoinhibition at a macroscopic level despite the relatively weak intramolecular interaction at each site. Both experimental and computational data suggest that the D/E repeats interact preferentially with other intrinsically disordered regions (IDRs) of HMGB1. We also found that mutations mimicking post-translational modifications relevant to nuclear export of HMGB1 can moderately modulate DNA-binding affinity, possibly by impacting the autoinhibition. This study illuminates a functional role of the fuzzy interactions of D/E repeats.

Modulating Microtubules: A Molecular Perspective on the Effects of Tail Modifications

Microtubules (MTs), an essential component of the eukaryotic cytoskeleton, are a lattice of polymerized tubulin dimers and are crucial for various cellular processes. The genetic and chemical diversity of tubulin and their disordered tails gives rise to a "tubulin code". The functional role of tubulin post-translational modifications (PTMs), which contribute to the chemical diversity of the tubulin code, is gradually being unraveled. However, variation in the length and spatial organization of tubulin poly-modifications leads to an enormous combinatorial PTM space, which is difficult to study experimentally. Hence, the impact of the combinatorial tubulin PTM space on the biophysical properties of tubulin tails and their interactions with other proteins remains elusive. Here, we combine all-atom and coarse-grained molecular dynamics simulations to elucidate the biophysical implications of the large combinatorial tubulin PTM space in the context of an MT lattice. We find that tail-body interactions are more dominant in the tubulin dimer than in an MT lattice, and are more significant for the tails of α compared with β tubulin. In addition, polyglutamylation, but not polyglycylation, expands the dimensions of the tubulin tails. Polyglutamylation also leads to a decrease in the diffusion rate of MT-associated protein EB1 on MTs, while polyglycylation often increases diffusion rate. These observations are generally not sensitive to the organization of the polymodifications. The effect of PTMs on MT charge density and tail dynamics are also discussed. Overall, this study presents a molecular quantification of the biophysical properties of tubulin tails and their polymodifications, and provides predictions on the functional importance of tubulin PTMs.

What Are the Molecular Requirements for Protein Sliding along DNA?

DNA-binding proteins rely on linear diffusion along the longitudinal DNA axis, supported by their nonspecific electrostatic affinity for DNA, to search for their target recognition sites. One may therefore expect that the ability to engage in linear diffusion along DNA is universal to all DNA-binding proteins, with the detailed biophysical characteristics of that diffusion differing between proteins depending on their structures and functions. One key question is whether the linear diffusion mechanism is defined by translation coupled with rotation, a mechanism that is often termed sliding. We conduct coarse-grained and atomistic molecular dynamics simulations to investigate the minimal requirements for protein sliding along DNA. We show that coupling, while widespread, is not universal. DNA-binding proteins that slide along DNA transition to uncoupled translation–rotation (i.e., hopping) at higher salt concentrations. Furthermore, and consistently with experimental reports, we find that the sliding mechanism is the less dominant mechanism for some DNA-binding proteins, even at low salt concentrations. In particular, the toroidal PCNA protein is shown to follow the hopping rather than the sliding mechanism.

Prerecognition Diffusion Mechanism of Human DNA Mismatch Repair Proteins along DNA: Msh2-Msh3 versus Msh2-Msh6

DNA mismatch repair (MMR) is an important postreplication process that eliminates mispaired or unpaired nucleotides to ensure genomic replication fidelity. In humans, Msh2-Msh6 and Msh2-Msh3 are the two mismatch repair initiation factors that recognize DNA lesions. While X-ray crystal structures exist for these proteins in complex with DNA lesions, little is known about their structures during the initial search along nonspecific double-stranded DNA, because they are short-lived and difficult to determine experimentally. In this study, various computational approaches were used to sidestep these difficulties. All-atom and coarse-grained simulations based on the crystal structures of Msh2-Msh3 and Msh2-Msh6 showed no translation along the DNA, suggesting that the initial search conformation differs from the lesion-bound crystal structure. We modeled probable search-mode structures of MSH proteins and showed, using coarse-grained molecular dynamics simulations, that they can perform rotation-coupled diffusion on DNA, which is a suitable and efficient search mechanism for their function and one predicted earlier by fluorescence resonance energy transfer and fluorescence microscopy studies. This search mechanism is implemented by electrostatic interactions among the mismatch-binding domain (MBD), the clamp domains, and the DNA backbone. During simulations, their diffusion rate did not change significantly with an increasing salt concentration, which is consistent with observations from experimental studies. When the gap between their DNA-binding clamps was increased, Msh2-Msh3 diffused mostly via the clamp domains while Msh2-Msh6 still diffused using the MBD, reproducing the experimentally measured lower diffusion coefficient of Msh2-Msh6. Interestingly, Msh2-Msh3 was capable of dissociating from the DNA, whereas Msh2-Msh6 always diffused on the DNA duplex. This is consistent with the experimental observation that Msh2-Msh3, unlike Msh2-Msh6, can overcome obstacles such as nucleosomes. Our models provide a molecular picture of the different mismatch search mechanisms undertaken by Msh2-Msh6 and Msh2-Msh3, despite the similarity of their structures.