What Are the Molecular Requirements for Protein Sliding along DNA?

Probing the modulation in facilitated diffusion guided by DNA-protein interactions in target search processes

Many fundamental biophysical processes involving gene regulation and gene editing rely, at the molecular level, on an intricate methodology of searching and locating the precise target base pair sequence on the genome by specific binding proteins. A unique mechanism, known as 'facilitated diffusion', which is a combination of 1D sliding along with 3D movement, is considered to be the key step for such events. This also explains the relatively much shorter timescale of the target searching process, compared to other diffusion-controlled biophysical processes. In this work, we aim to probe the modulation of target search dynamics of a protein moiety by estimating the rate of the target search process, and the statistics of the search rounds and timescales accomplished by the 1D and 3D motions, based on first passage time (FPT) calculations. This is studied with its characteristics getting influenced by various given conditions such as, when the DNA is rigid or flexible, and when the target is placed at different locations on the DNA. The current theoretical framework includes a Brownian dynamics simulation setup adopting a straightforward coarse-grained model for a diffusing protein on DNA. Moreover, this theoretical analysis provides insights into the complex target search dynamics by highlighting the significance of the chain dynamics in the mechanistic details of the facilitated diffusion process.

Interactions of proteins with heparan sulfate

Heparan sulfate (HS) is a glycosaminoglycan, polysaccharides that are considered to have arisen in the last common unicellular ancestor of multicellular animals. In this light, the large interactome of HS and its myriad functions in relation to the regulation of cell communication are not surprising. The binding of proteins to HS determines their localisation and diffusion, essential for embryonic development and homeostasis. Following the biosynthesis of the initial heparosan polymer, the subsequent modifications comprise an established canonical pathway and a minor pathway. The more frequent former starts with N-deacetylation and N-sulfation of GlcNAc residues, the latter with C-5 epimerisation of a GlcA residue adjacent to a GlcNAc. The binding of proteins to HS is driven by ionic interactions. The multivalent effect arising from the many individual ionic bonds between a single protein and a polysaccharide chain results in a far stronger interaction than would be expected from an ion-exchange process. In many instances, upon binding, both parties undergo substantial conformational change, the resulting hydrogen and van der Waal bonds contributing significant free energy to the binding reaction. Nevertheless, ionic bonds dominate the protein-polysaccharide interaction kinetically. Together with the multivalent effect, this provides an explanation for the observed trapping of HS-binding proteins in extracellular matrix. Importantly, individual ionic bonds have been observed to be dynamic; breaking and reforming, while the protein remains bound to the polysaccharide. These considerations lead to a model for 1D diffusion of proteins in extracellular matrix on HS, involving mechanisms such as sliding, chain switching and rolling.

Three-State Diffusion Model of DNA Glycosylase Translocation along Stretched DNA as Revealed by Free Energy Landscapes at the All-Atom Level

DNA glycosylases play key roles in the maintenance of genomic integrity. These enzymes effectively find rare damaged sites in DNA and participate in subsequent base excision repair. Single-molecule and ensemble experiments have revealed key aspects of this damage-site searching mechanism and the involvement of facilitated diffusion. In this study, we describe free energy landscapes of enzyme translocation along nonspecific DNA obtained using a fully atomistic molecular dynamics (MD) simulation of a well-known DNA glycosylase, human 8-oxoguanine DNA glycosylase 1 (hOGG1). Based on an analysis of simulated free energy profiles, we propose a three-state model for the damage-site searching mechanism. In the three states, named the L1, L2, and L3 states, the L1 state is a helical sliding mode in close contact with DNA, whereas the L2 state is a major- or minor-groove tracking mode in loose contact with DNA and the L3 state is a two-dimensional freely diffusing mode during which hOGG1 is somewhat removed from the DNA surface (∼24 Å away from the surface). This three-state model well describes key experimental findings obtained from single-molecule and ensemble experiments and provides a unified molecular picture of the DNA lesion-searching mechanism of hOGG1.

Nonspecific vs. specific DNA binding free energetics of a transcription factor domain protein

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling.

Single-molecule analysis reveals TDG exhibits multiple modes of linear diffusion to process 5-formylcytosine

DNA methylation plays a key role in epigenetics, with 60-80% of CpG sites containing 5-methylcytosine. Base excision repair (BER) is suggested to be the main pathway involved in active DNA demethylation. 5-formylctyosine (5fC), an oxidized moiety of methylated cytosine, is recognized and removed by thymine DNA glycosylase (TDG) to generate an abasic site. TDG binds avidly to abasic sites and is product inhibited. Using single molecule fluorescence experiments, we saw TDG interact with DNA containing 5fC specifically and non-specifically with lifetimes of 72.9 and 7.5 seconds, respectively. These results indicate that TDG cleaves the 5fC and stays bound for an extended time at the generated abasic site. Mean squared displacement analysis and a two color TDG experiment indicate that TDG exhibits multiple modes of linear diffusion, including hopping and sliding, in search of a lesion. The catalytically crippled variants, N140A and R275A/L, have a reduced binding lifetime compared to wild type and Mean Squared Displacement (MSD) analysis indicates that R275L/A moves on the DNA with a faster diffusivity. These results indicate that mutating R275, but not N140 interferes with damage recognition by TDG. Our findings give insight into how TDG searches for its lesions in long stretches of undamaged DNA.

Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions

The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.

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.

Conformational Analysis of Charged Homo-Polypeptides

Many proteins have intrinsically disordered regions (IDRs), which are often characterized by a high fraction of charged residues with polyampholytic (i.e., mixed charge) or polyelectrolytic (i.e., uniform charge) characteristics. Polyelectrolytic IDRs include consecutive positively charged Lys or Arg residues (K/R repeats) or consecutive negatively charged Asp or Glu residues (D/E repeats). In previous research, D/E repeats were found to be about five times longer than K/R repeats and to be much more common in eukaryotes. Within these repeats, a preference is often observed for E over D and for K over R. To understand the greater prevalence of D/E over K/R repeats and the higher abundance of E and K, we simulated the conformational ensemble of charged homo-polypeptides (polyK, polyR, polyD, and polyE) using molecular dynamics simulations. The conformational preferences and dynamics of these polyelectrolytic polypeptides change with changes in salt concentration. In particular, polyD and polyE are more sensitive to salt than polyK and polyR, as polyD and polyE tend to adsorb more divalent cations, which leads to their having more compact conformations. We conclude with a discussion of biophysical explanations for the relative abundance of charged amino acids and particularly for the greater abundance of D/E repeats over K/R repeats.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

Theoretical insights into the full description of DNA target search by subdiffusing proteins

DNA binding proteins (DBPs) diffuse in the cytoplasm to recognise and bind with their respective target sites on the DNA to initiate several biologically important processes. The first passage time distributions (FPTDs) of DBPs are useful in quantifying the timescales of the most-probable search paths in addition to the mean value of the distribution which, strikingly, are decades of order apart in time. However, extremely crowded in vivo conditions or the viscoelasticity of the cellular medium among other factors causes biomolecules to exhibit anomalous diffusion which is usually overlooked in most theoretical studies. We have obtained approximate analytical expressions of a general FPTD and the two characteristic timescales that are valid for any single subdiffusing protein searching for its target in vivo. Our results can be applied to single-particle tracking experiments of target search.

Salt Dependence of DNA Binding Activity of Human Transcription Factor Dlx3

Distal-less 3 (Dlx3) is a homeobox-containing transcription factor and plays a crucial role in the development and differentiation process. Human Dlx3 consists of two transactivation domains and a homeobox domain (HD) that selectively binds to the consensus site (5'-TAATT-3') of the DNA duplex. Here, we performed chemical shift perturbation experiments on Dlx3-HD in a complex with a 10-base-paired (10-bp) DNA duplex under various salt conditions. We also acquired the imino proton spectra of the 10-bp DNA to monitor the changes in base-pair stabilities during titration with Dlx3-HD. Our study demonstrates that Dlx3-HD selectively recognizes its consensus DNA sequences through the α3 helix and L1 loop regions with a unique dynamic feature. The dynamic properties of the binding of Dlx3-HD to its consensus DNA sequence can be modulated by varying the salt concentrations. Our study suggested that this unique structural and dynamic feature of Dlx3-HD plays an important role in target DNA recognition, which might be associated with tricho-dento-osseous syndrome.

Dissociation Pathways of the p53 DNA Binding Domain from DNA and Critical Roles of Key Residues Elucidated by dPaCS-MD/MSM

p53 is a transcriptional factor that regulates cell response to a variety of stresses. About a half of all human tumors contain p53 mutations, and the accumulation of mutations in the DNA binding domain of p53 (p53-DBD) can cause destabilization of p53 and its complex with DNA. To identify the key residues of the p53-DBD/DNA binding and to understand the dissociation mechanisms of the p53-DBD/DNA complex, the dissociation process of p53-DBD from a DNA duplex that contains the consensus sequence (the specific target of p53-DBD) was investigated by a combination of dissociation parallel cascade selection molecular dynamics (dPaCS-MD) and the Markov state model (MSM). This combination (dPaCS-MD/MSM) enabled us to simulate dissociation of the two large molecules based on an all-atom model with a short simulation time (11.2 ± 2.2 ns per trial) and to analyze dissociation pathways, free energy landscape (FEL), and binding free energy. Among 75 trials of dPaCS-MD, p53-DBD dissociated first from the major groove and then detached from the minor groove in 93% of the cases, while 7% of the cases unbinding from the minor groove occurred first. Minor groove binding is mainly stabilized by R248, identified as the most important residue that tightly binds deep inside the minor groove. The standard binding free energy calculated from the FEL was -10.9 ± 0.4 kcal/mol, which agrees with an experimental value of -11.1 kcal/mol. These results indicate that the dPaCS-MD/MSM combination can be a powerful tool to investigate dissociation mechanisms of two large molecules. Analysis of the p53 key residues for DNA binding indicates high correlations with cancer-related mutations, confirming that impairment of the interactions between p53-DBD and DNA can be frequently related to cancer.

Structure and receptor recognition by the Lassa virus spike complex

Lassa virus (LASV) is a human pathogen, causing substantial morbidity and mortality^1,2. Similar to other Arenaviridae, it presents a class-I spike complex on its surface that facilitates cell entry. The virus's cellular receptor is matriglycan, a linear carbohydrate that is present on α-dystroglycan^3,4, but the molecular mechanism that LASV uses to recognize this glycan is unknown. In addition, LASV and other arenaviruses have a unique signal peptide that forms an integral and functionally important part of the mature spike^5-8; yet the structure, function and topology of the signal peptide in the membrane remain uncertain^9-11. Here we solve the structure of a complete native LASV spike complex, finding that the signal peptide crosses the membrane once and that its amino terminus is located in the extracellular region. Together with a double-sided domain-switching mechanism, the signal peptide helps to stabilize the spike complex in its native conformation. This structure reveals that the LASV spike complex is preloaded with matriglycan, suggesting the mechanism of binding and rationalizing receptor recognition by α-dystroglycan-tropic arenaviruses. This discovery further informs us about the mechanism of viral egress and may facilitate the rational design of novel therapeutics that exploit this binding site.