Energetic and structural considerations for the mechanism of protein sliding along DNA in the nonspecific BamHI-DNA complex

Noncanonical DNA Cleavage by BamHI Endonuclease in Laterally Confined DNA Monolayers Is a Step Function of DNA Density and Sequence

Cleavage of DNA at noncanonical recognition sequences by restriction endonucleases (star activity) in bulk solution can be promoted by global experimental parameters, including enzyme or substrate concentration, temperature, pH, or buffer composition. To study the effect of nanoscale confinement on the noncanonical behaviour of BamHI, which cleaves a single unique sequence of 6 bp, we used AFM nanografting to generate laterally confined DNA monolayers (LCDM) at different densities, either in the form of small patches, several microns in width, or complete monolayers of thiol-modified DNA on a gold surface. We focused on two 44-bp DNAs, each containing a noncanonical BamHI site differing by 2 bp from the cognate recognition sequence. Topographic AFM imaging was used to monitor end-point reactions by measuring the decrease in the LCDM height with respect to the surrounding reference surface. At low DNA densities, BamHI efficiently cleaves only its cognate sequence while at intermediate DNA densities, noncanonical sequence cleavage occurs, and can be controlled in a stepwise (on/off) fashion by varying the DNA density and restriction site sequence. This study shows that endonuclease action on noncanonical sites in confined nanoarchitectures can be modulated by varying local physical parameters, independent of global chemical parameters.

Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes

Synthesis and scission of phosphodiester bonds in DNA and RNA regulate vital processes within the cell. Enzymes that catalyze these reactions operate mostly via the recognized two-metal-ion mechanism. Our analysis reveals that basic amino acids and monovalent cations occupy structurally conserved positions nearby the active site of many two-metal-ion enzymes for which high-resolution (<3 Å) structures are known, including DNA and RNA polymerases, nucleases such as Cas9, and splicing ribozymes. Integrating multiple-sequence and structural alignments with molecular dynamics simulations, electrostatic potential maps, and mutational data, we found that these elements always interact with the substrates, suggesting that they may play an active role for catalysis, in addition to their electrostatic contribution. We discuss possible mechanistic implications of this expanded two-metal-ion architecture, including inferences on medium-resolution cryoelectron microscopy structures. Ultimately, our analysis may inspire future experiments and strategies for enzyme engineering or drug design to modulate nucleic acid processing.

•

Basic residues in the active site of two-metal-ion enzymes are structurally conserved

•

These residues are also conserved in evolution

•

Mutagenesis suggests these residues may exert an effect on DNA- and RNA processing

•

Our work offers insights into CRISPR/Cas9, spliceosome, and DNA/RNA polymerases

On the Possibility of Facilitated Diffusion of Dendrimers Along DNA

We investigate the electrostatics, energetics, and dynamics of dendrimer-DNA interactions that mimic protein-DNA complexes as a means to design facilitated mechanisms by which dendrimers can slide and search DNA for targets. By using all-atom molecular dynamics simulations, we calculated the free energy profiles of dendrimer-binding around the DNA via umbrella sampling. We also calculated electrostatic interaction maps in comparison to proteins, as well as the dynamical changes induced by DNA-dendrimer interactions via NMR-measurable order parameters. Our results show that for dendrimers to go around DNA, there is a free-energy barrier of 8.5 kcal/mol from the DNA major groove to DNA minor groove, with a minimum in the major groove. This barrier height makes it unlikely for an all-amine dendrimer to slide along DNA longitudinally, but following a helical path may be possible along the major groove. Comparison of the nonbonded interaction energy and the interaction free-energy profiles reveal a considerable entropic cost as the dendrimer binds to DNA. This is also supported by the mobility patterns obtained from NMR-measurable order parameter values, which show a decreased mobility of the dendrimer N-H bond vectors in the DNA-binding mode.

DNA recognition process of the lactose repressor protein studied via metadynamics and umbrella sampling simulations

The lactose repressor, LacI, finds its DNA target sites via a process that is faster than what it is expected from a diffusion-driven mechanism. This is possible thanks to nonspecific binding of LacI to DNA, followed by diffusion along the DNA molecule. The diffusion of the protein along DNA might lead to a fast-searching mechanism only if LacI binds with comparable strength to different nonspecific sequences and if, in addition, the value of the binding energy remarkably decreases in the presence of a binding site. The first condition would be favored by loose interactions with the base edges, while the second would take advantage from the opposite situation. In order to understand how the protein satisfies these two opposing requirements, the DNA recognition process was studied by a combination of umbrella sampling and metadynamics simulations. The simulations revealed that when aligned with a specific sequence, LacI establishes polar interactions with the base edges that require ∼4 kcal/mol to be disrupted. In contrast, these interactions are not stable when the protein is aligned with nonspecific sequences. These results confirm that LacI is able to efficiently recognize a specific sequence while sliding along DNA before any structural change of the protein-DNA complex occurs.

Transcription-factor binding and sliding on DNA studied using micro- and macroscopic models

Transcription factors search for specific operator sequences by alternating rounds of 3D diffusion with rounds of 1D diffusion (sliding) along the DNA. The details of such sliding have largely been beyond direct experimental observation. For this purpose we devised an analytical formulation of umbrella sampling along a helical coordinate, and from extensive and fully atomistic simulations we quantified the free-energy landscapes that underlie the sliding dynamics and dissociation kinetics for the LacI dimer. The resulting potential of mean force distributions show a fine structure with an amplitude of 1 k(B)T for sliding and 12 k(B)T for dissociation. Based on the free-energy calculations the repressor slides in close contact with DNA for 8 bp on average before making a microscopic dissociation. By combining the microscopic molecular-dynamics calculations with Brownian simulation including rotational diffusion from the microscopically dissociated state we estimate a macroscopic residence time of 48 ms at the same DNA segment and an in vitro sliding distance of 240 bp. The sliding distance is in agreement with previous in vitro sliding-length estimates. The in vitro prediction for the macroscopic residence time also compares favorably to what we measure by single-molecule imaging of nonspecifically bound fluorescently labeled LacI in living cells. The investigation adds to our understanding of transcription-factor search kinetics and connects the macro-/mesoscopic rate constants to the microscopic dynamics.

Weak frustration regulates sliding and binding kinetics on rugged protein-DNA landscapes

A fundamental step in gene-regulatory activities, such as repression, transcription, and recombination, is the binding of regulatory DNA-binding proteins (DBPs) to specific targets in the genome. To rapidly localize their regulatory genomic sites, DBPs reduce the dimensionality of the search space by combining three-dimensional (3D) diffusion in solution with one-dimensional (1D) sliding along DNA. However, the requirement to form a thermodynamically stable protein-DNA complex at the cognate genomic target sequence imposes a challenge on the protein because, as it navigates one-dimensionally along the genome, it may come in close contact with sites that share partial or even complete sequence similarity with the functional DNA sequence. This puzzling issue creates a conflict between two basic requirements: finding the cognate site quickly and stably binding it. Here, we structurally assessed the interface adopted by a variety of DBPs to bind DNA specifically and nonspecifically, and found that many DBPs utilize one interface to specifically recognize a DNA sequence and another to assist in propagating along the DNA through nonspecific associations. While these two interfaces overlap each other in some proteins, they present partial overlap in others and frustrate the protein-DNA interface. Using coarse-grained molecular dynamics simulations, we demonstrate that the existence of frustration in DBPs is a compromise between rapid 1D diffusion along other regions in the genome (high frustration smoothens the landscape for sliding) and rapid formation of a stable and essentially active protein-DNA complex (low frustration reduces the free energy barrier for switching between the two binding modes).

DNA-recognition process described by MD simulations of the lactose repressor protein on a specific and a non-specific DNA sequence

The lactose repressor protein may bind DNA in two possible configurations: a specific one, if the DNA sequence corresponds to a binding site, and a non-specific one otherwise. To find its target sequences, the lactose repressor first binds non-specifically to DNA, and subsequently, it rapidly searches for a binding site. Atomic structures of non-specific and specific complexes are available from crystallographic and nuclear magnetic resonance experiments. However, what remains unknown is a detailed description of the steps that transform the non-specific complex into the specific one. Here, how the protein first recognizes its binding site has been studied using molecular dynamics simulations. The picture that emerges is that of a protein that is as mobile when interacting with non-specific DNA sequences as when free in solution. This high degree of mobility allows the protein to rapidly sample different DNA sequences. In contrast, when the protein encounters a binding site, the configuration ensemble collapses, and the protein sliding movements along the DNA sequence become scarce. The binding energies in the specific and non-specific complexes were analysed using the Molecular Mechanics Poisson Boltzmann Surface Area approach. These results represent a first step towards a throughout characterization of the DNA-recognition process.

The roles of intrinsic disorder in orchestrating the Wnt-pathway

The canonical Wnt-pathway plays a number of crucial roles in the development of organism. Malfunctions of this pathway lead to various diseases including cancer. In the inactivated state, this pathway involves five proteins, Axin, CKI-α, GSK-3β, APC, and β-catenin. We analyzed these proteins by a number of computational tools, such as PONDR(r)VLXT, PONDR(r)VSL2, MoRF-II predictor and Hydrophobic Cluster Analysis (HCA) to show that each of the Wnt-pathway proteins contains several intrinsically disordered regions. Based on a comprehensive analysis of published data we conclude that these disordered regions facilitate protein-protein interactions, post-translational modifications, and signaling. The scaffold protein Axin and another large protein, APC, act as flexible concentrators in gathering together all other proteins involved in the Wnt-pathway, emphasizing the role of intrinsically disordered regions in orchestrating the complex protein-protein interactions. We further explore the intricate roles of highly disordered APC in regulation of β-catenin function. Intrinsically disordered APC helps the collection of β-catenin from cytoplasm, facilitates the b-catenin delivery to the binding sites on Axin, and controls the final detachment of β-catenin from Axin.

Frustration in protein-DNA binding influences conformational switching and target search kinetics

Rapid recognition of DNA target sites involves facilitated diffusion through which alternative sites are searched on genomic DNA. A key mechanism facilitating the localization of the target by a DNA-binding protein (DBP) is one-dimensional diffusion (sliding) in which electrostatic forces attract the protein to the DNA. As the protein reaches its target DNA site, it switches from purely electrostatic binding to a specific set of interactions with the DNA bases that also involves hydrogen bonding and van der Waals forces. High overlap between the DBP patches used for nonspecific and specific interactions with DNA may enable an immediate transition between the two binding modes following target site localization. By contrast, an imperfect overlap may result in greater frustration between the two potentially competing binding modes and consequently slower switching between them. A structural analysis of 125 DBPs indicates frustration between the two binding modes that results in a large difference between the orientations of the protein to the DNA when it slides compared to when it specifically interacts with DNA. Coarse-grained molecular dynamics simulations of in silico designed peptides comprising the full range of frustrations between the two interfaces show slower transition from nonspecific to specific DNA binding as the overlap between the patches involved in the two binding modes decreases. The complex search kinetics may regulate the search by eliminating trapping of the protein in semispecific sites while sliding.

Diffusion of MMPs on the surface of collagen fibrils: the mobile cell surface-collagen substratum interface

Remodeling of the extracellular matrix catalyzed by MMPs is central to morphogenetic phenomena during development and wound healing as well as in numerous pathologic conditions such as fibrosis and cancer. We have previously demonstrated that secreted MMP-2 is tethered to the cell surface and activated by MT1-MMP/TIMP-2-dependent mechanism. The resulting cell-surface collagenolytic complex (MT1-MMP)(2)/TIMP-2/MMP-2 can initiate (MT1-MMP) and complete (MMP-2) degradation of an underlying collagen fibril. The following question remained: What is the mechanism of substrate recognition involving the two structures of relatively restricted mobility, the cell surface enzymatic complex and a collagen fibril embedded in the ECM? Here we demonstrate that all the components of the complex are capable of processive movement on a surface of the collagen fibril. The mechanism of MT1-MMP movement is a biased diffusion with the bias component dependent on the proteolysis of its substrate, not adenosine triphosphate (ATP) hydrolysis. It is similar to that of the MMP-1 Brownian ratchet we described earlier. In addition, both MMP-2 and MMP-9 as well as their respective complexes with TIMP-1 and -2 are capable of Brownian diffusion on the surface of native collagen fibrils without noticeable dissociation while the dimerization of MMP-9 renders the enzyme immobile. Most instructive is the finding that the inactivation of the enzymatic activity of MT1-MMP has a detectable negative effect on the cell force developed in miniaturized 3D tissue constructs. We propose that the collagenolytic complex (MT1-MMP)(2)/TIMP-2/MMP-2 represents a Mobile Cell Surface-Collagen Substratum Interface. The biological implications of MT1-MMP acting as a molecular ratchet tethered to the cell surface in complex with MMP-2 suggest a new mechanism for the role of spatially regulated peri-cellular proteolysis in cell-matrix interactions.

Nonspecific DNA-protein interaction: why proteins can diffuse along DNA

Recent single molecule experiments have reported that DNA binding proteins (DNA-BPs) can diffuse along DNA. This suggests that interactions between proteins and DNA play a role during the target search even far from their specific site on DNA. Here we show by means of Monte Carlo simulations and analytical calculations that there is a counterintuitive repulsion between the two oppositely charged macromolecules at a nanometer range. For the concave shape of DNA-BPs, and for realistic protein charge densities, we find that the DNA-protein interaction free energy has a minimum at a finite surface-to-surface separation, in which proteins can easily slide. When a protein encounters its target, the free energy barrier is completely counterbalanced by the H-bond interaction, thus enabling the sequence recognition.

Human alkyladenine DNA glycosylase employs a processive search for DNA damage

DNA repair proteins conduct a genome-wide search to detect and repair sites of DNA damage wherever they occur. Human alkyladenine DNA glycosylase (AAG) is responsible for recognizing a variety of base lesions, including alkylated and deaminated purines, and initiating their repair via the base excision repair pathway. We have investigated the mechanism by which AAG locates sites of damage using an oligonucleotide substrate containing two sites of DNA damage. This substrate was designed so that AAG randomly binds to either of the two lesions. AAG-catalyzed base excision creates a repair intermediate, and the subsequent partitioning between dissociation and diffusion to the second site can be quantified from the rates of formation of the different products. Our results demonstrate that AAG has the ability to slide for short distances along DNA at physiological salt concentrations. The processivity of AAG decreases with increasing ionic strength to become fully distributive at high ionic strengths, suggesting that electrostatic interactions between the negatively charged DNA and the positively charged DNA binding surface are important for nonspecific DNA binding. Although the amino terminus of the protein is dispensable for glycosylase activity at a single site, we find that deletion of the 80 amino-terminal amino acids significantly decreases the processivity of AAG. These observations support the idea that diffusion on undamaged DNA contributes to the search for sites of DNA damage.

Effect of somatic mutation on DNA binding properties of anti-DNA autoantibodies

Autoantibodies that bind DNA are a hallmark of systemic lupus erythematosus. A subset of autoantibody*DNA complexes localize to kidney tissue and lead to damage and even death. 11F8, 9F11, and 15B10 are clonally related anti-DNA autoantibodies isolated from an autoimmune mouse. 11F8 binds ssDNA in a sequence-specific manner and causes tissue damage, while 9F11 and 15B10 bind ssDNA non-specifically and are benign. Among these antibodies, DNA binding properties are mediated by five amino acid differences in primary sequence. Thermodynamic and kinetic parameters associated with recognition of structurally different DNA sequences were determined for each antibody to provide insight toward recognition strategies, and to explore a link between binding properties and disease pathogenesis. A model of 11F8 bound to its high affinity consensus sequence provides a foundation for understanding the differences in thermodynamic and kinetic parameters between the three mAbs. Our data suggest that 11F8 utilizes the proposed ssDNA recognition motif including (Y32)V(L), a hydrogen bonding residue at (91)V(L), and an aromatic residue at the tip of the third heavy chain complementarity determining region. Interestingly, a somatic mutation to arginine at (31)V(H) in 11F8 may afford additional binding site contacts including (R31)V(H), (R96)V(H), and (R98)V(H) that could determine specificity.

Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study

Four 20 ns molecular dynamics simulations have been performed with two counterions, K⁺ or Na⁺, at two water contents, 15 or 20 H₂O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5′-ATGCAGTCAG) × d(5′-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K⁺ and Na⁺ binding to DNA are different. K⁺ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na⁺ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K⁺ and Na⁺ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the K–DNA systems is faster than in the Na–DNA systems: Na⁺ organizes and immobilizes water structure around itself and near DNA while for K⁺ water is less organized and more dynamic.

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.

RNA polymerase can track a DNA groove during promoter search

Many proteins select special DNA sequences to form functional complexes. In one possible mechanism, protein molecules would scan DNA sequences by tracking a groove without complete dissociation. Upon dragging single molecules of DNA over a surface carrying fixed Escherichia coli RNA polymerase holoenzyme, we detected rotation of individual DNA molecules, providing direct evidence that a DNA-binding protein can track a DNA groove. These results confirm our previous observations of longitudinal movement of RNA polymerase along fixed, extended DNA and, moreover, imply that groove tracking facilitates scanning of DNA sequences.

DNA Binding Provides a Molecular Strap Activating the Adenovirus Proteinase

Human adenovirus proteinase (AVP) requires two cofactors for maximal activity: pVIc, a peptide derived from the C terminus of adenovirus precursor protein pVI, and the viral DNA. Synchrotron protein footprinting was used to map the solvent accessible cofactor binding sites and to identify conformational changes associated with the binding of cofactors to AVP. The binding of pVIc alone or pVIc and DNA together to AVP triggered significant conformational changes adjacent to the active site cleft sandwiched between the two AVP subdomains. In addition, upon binding of DNA to AVP, it was observed that specific residues on each of the two major subdomains were significantly protected from hydroxyl radicals. Based on the locations of these protected side-chain residues and conserved aromatic and positively charged residues within AVP, a three-dimensional model of DNA binding was constructed. The model indicated that DNA binding can alter the relative orientation of the two AVP domains leading to the partial activation of AVP by DNA. In addition, both pVIc and DNA may independently alter the active site conformation as well as drive it cooperatively to fully activate AVP.

Obstacle Bypass in Protein Motion along DNA by Two-dimensional Rather than One-dimensional Sliding

Site-specific DNA-binding proteins locate their target sites by facilitated diffusion. Several proteins have been shown to slide along DNA in vitro. However, whereas sliding is often envisaged as one-dimensional tracking of the DNA major groove, such a mechanism would not allow linear diffusion over long distances in vivo, where short stretches of free DNA are delimited by bound proteins. I propose a two-dimensional sliding mechanism, in which the protein diffuses freely on the cylindrical DNA surface, and I present experiments that can distinguish between one- and higher-dimensional diffusion along the DNA contour length. At 100 mm NaCl, translocation of EcoRI restriction endonuclease between sites on two DNA helices connected by a Holliday junction is as efficient as between sites on the same helix, indicating a three-dimensional mechanism. At 25 mm NaCl, translocation between sites on the same DNA helix is more efficient, indicating a role for sliding at low ionic strength. Obstacles attached to the major groove of one face of the DNA helix did not interfere with sliding, regardless of their orientation relative to the cleavage sites. This result is compatible with two-dimensional but not one-dimensional sliding. As illustrated by Monte-Carlo simulation, two-dimensional sliding may not only allow proteins to move around nucleosomes in vivo but also reduce the redundancy of their search for the target site.