Protein-DNA docking with a coarse-grained force field

ABSTRACT: BACKGROUND: Protein-DNA interactions are important for many cellular processes, however structural knowledge for a large fraction of known and putative complexes is still lacking. Computational docking methods aim at the prediction of complex architecture given detailed structures of its constituents. They are becoming an increasingly important tool in the field of macromolecular assemblies, complementing particularly demanding protein-nucleic acids X ray crystallography and providing means for the refinement and integration of low resolution data coming from rapidly advancing methods such as cryoelectron microscopy. RESULTS: We present a new coarse-grained force field suitable for protein-DNA docking. The force field is an extension of previously developed parameter sets for protein-RNA and protein-protein interactions. The docking is based on potential energy minimization in translational and orientational degrees of freedom of the binding partners. It allows for fast and efficient systematic search for native-like complex geometry without any prior knowledge regarding binding site location. CONCLUSIONS: We find that the force field gives very good results for bound docking. The quality of predictions in the case of unbound docking varies, depending on the level of structural deviation from bound geometries. We analyze the role of specific protein-DNA interactions on force field performance, both with respect to complex structure prediction, and the reproduction of experimental binding affinities. We find that such direct, specific interactions only partially contribute to protein-DNA recognition, indicating an important role of shape complementarity and sequence-dependent DNA internal energy, in line with the concept of indirect protein-DNA readout mechanism.     

Non-specific protein-DNA interactions control I-CreI target binding and cleavage

Homing endonucleases represent protein scaffolds that provide powerful tools for genome manipulation, as these enzymes possess a very low frequency of DNA cleavage in eukaryotic genomes due to their high specificity. The basis of protein-DNA recognition must be understood to generate tailored enzymes that target the DNA at sites of interest. Protein-DNA interaction engineering of homing endonucleases has demonstrated the potential of these approaches to create new specific instruments to target genes for inactivation or repair. Protein-DNA interface studies have been focused mostly on specific contacts between amino acid side chains and bases to redesign the binding interface. However, it has been shown that 4 bp in the central DNA sequence of the 22-bp substrate of a homing endonuclease (I-CreI), which do not show specific protein-DNA interactions, is not devoid of content information. Here, we analyze the mechanism of target discrimination in this substrate region by the I-CreI protein, determining how it can occur independently of the specific protein-DNA interactions. Our data suggest the important role of indirect readout in this substrate region, opening the possibility for a fully rational search of new target sequences, thus improving the development of redesigned enzymes for therapeutic and biotechnological applications.     

Self-assembled monolayer-assisted silicon nanowire biosensor for detection of protein-DNA interactions in nuclear extracts from breast cancer cell

The large number of estrogen receptor (ER) binding sites of various sequence patterns requires a sensitive detection to differentiate between subtle differences in ER-DNA binding affinities. A self-assembled monolayer (SAM)-assisted silicon nanowire (SiNW) biosensor for specific and highly sensitive detection of protein-DNA interactions, remarkably in nuclear extracts prepared from breast cancer cells, is presented. As a typical model, estrogen receptor element (ERE, dsDNA) and estrogen receptor alpha (ERα, protein) binding was adopted in the work. The SiNW surface was coated with a vinyl-terminated SAM, and the termination of the surface was changed to carboxylic acid via oxidation. DNA modified with amine group was subsequently immobilized on the SiNW surface. Protein-DNA binding was finally investigated by the functionalized SiNW biosensor. X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM) were employed to characterize the stepwise functionalization of the SAM and DNA on bare silicon surface, and to visualize protein-DNA binding on the SiNW surface, respectively. We observed that ERα had high sequence specificity to the SiNW biosensor which was functionalized with three different EREs including wild-type, mutant and scrambled DNA sequences. We also demonstrate that the specific DNA-functionalized SiNW biosensor was capable of detecting ERα as low as 10 fM. Impressively, the developed SiNW biosensor was able to detect ERα-DNA interactions in nuclear extracts from breast cancer cells. The SAM-assisted SiNW biosensor, as a label-free and highly sensitive tool, shows a potential in studying protein-DNA interactions.     

Characterization of functional domains in human Claspin

Claspin is a mediator of the ATR-dependent DNA replication checkpoint in human cells and also promotes DNA replication fork progression and stability. Though Claspin has been shown to bind DNA and co-immunoprecipitate with other replication fork-associated proteins, the specific protein-protein and protein-DNA interactions that are important for Claspin function are not known. We therefore purified several domains of human Claspin and then tested for direct interactions of these fragments with several replication fork-associated proteins and with DNA. Our data show that the N terminus of Claspin binds to the replicative helicase co-factor Cdc45, the Timeless protein and a branched, replication fork-like DNA structure. In contrast, the C terminus of Claspin associates with DNA polymerase epsilon and Rad17-Replication Factor C (RFC). We conclude that multiple protein-DNA and protein-protein interactions may be important for Claspin function during DNA replication and DNA replication checkpoint signaling.Key words: Claspin, DNA replication, checkpoint, DNA damage, Cdc45, DNA polymerase, Rad17     

Combining chromatin immunoprecipitation and DNA footprinting: a novel method to analyze protein-DNA interactions in vivo

A variety of methods are available to analyze protein-DNA interactions in vivo. Two of the most prominent of these methods are chromatin immunoprecipitation (ChIP) and in vivo footprinting. Both of these procedures have specific limitations. For example, the ChIP assay fails to document where exactly a protein binds in vivo. The precipitation of a specific segment of DNA with antibodies directed against DNA-binding proteins does not necessarily indicate that the protein directly interacts with a sequence in the precipitate but could rather reflect protein-protein interactions. Furthermore, the results of in vivo footprinting studies are inconclusive if a DNA sequence is analyzed that is bound by a specific protein in only a certain fraction of cells. Finally, in vivo footprinting does not indicate which protein is bound at a specific site. We have developed a new procedure that combines the ChIP assay and DMS footprinting techniques. Using this method we show here that antibodies specific for USF1 and NF-E2 precipitate the murine beta-globin promoter in MEL cells. DMS footprinting analysis of the DNA precipitated with NF-E2 antibodies revealed a protection over a partial NF-E2-binding site in the beta-globin downstream promoter region. We believe that this novel method will generally benefit investigators interested in analyzing protein-DNA interactions in vivo.     

A simple physical model for the prediction and design of protein-DNA interactions

Protein-DNA interactions are crucial for many biological processes. Attempts to model these interactions have generally taken the form of amino acid-base recognition codes or purely sequence-based profile methods, which depend on the availability of extensive sequence and structural information for specific structural families, neglect side-chain conformational variability, and lack generality beyond the structural family used to train the model. Here, we take advantage of recent advances in rotamer-based protein design and the large number of structurally characterized protein-DNA complexes to develop and parameterize a simple physical model for protein-DNA interactions. The model shows considerable promise for redesigning amino acids at protein-DNA interfaces, as design calculations recover the amino acid residue identities and conformations at these interfaces with accuracies comparable to sequence recovery in globular proteins. The model shows promise also for predicting DNA-binding specificity for fixed protein sequences: native DNA sequences are selected correctly from pools of competing DNA substrates; however, incorporation of backbone movement will likely be required to improve performance in homology modeling applications. Interestingly, optimization of zinc finger protein amino acid sequences for high-affinity binding to specific DNA sequences results in proteins with little or no predicted specificity, suggesting that naturally occurring DNA-binding proteins are optimized for specificity rather than affinity. When combined with algorithms that optimize specificity directly, the simple computational model developed here should be useful for the engineering of proteins with novel DNA-binding specificities.     

Structure-based prediction of DNA-binding sites on proteins using the empirical preference of electrostatic potential and the shape of molecular surfaces

Protein-DNA interactions play an essential role in the genetic activities of life. Many structures of protein-DNA complexes are already known, but the common rules on how and where proteins bind to DNA have not emerged. Many attempts have been made to predict protein-DNA interactions using structural information, but the success rate is still about 80%. We analyzed 63 protein-DNA complexes by focusing our attention on the shape of the molecular surface of the protein and DNA, along with the electrostatic potential on the surface, and constructed a new statistical evaluation function to make predictions of DNA interaction sites on protein molecular surfaces. The shape of the molecular surface was described by a combination of local and global average curvature, which are intended to describe the small convex and concave and the large-scale concave curvatures of the protein surface preferentially appearing at DNA-binding sites. Using these structural features, along with the electrostatic potential obtained by solving the Poisson-Boltzmann equation numerically, we have developed prediction schemes with 86% and 96% accuracy for DNA-binding and non-DNA-binding proteins, respectively.     

A new analytical scale DNA affinity binding assay for analyses of specific protein-DNA interactions.

We describe a rapid analytical assay for identification of proteins binding to specific DNA sequences. The DAPSTER assay (DNA affinity preincubation specificity test of recognition assay) is a DNA affinity chromatography-based microassay that can discriminate between specific and nonspecific protein-DNA interactions. The assay is sensitive and can detect protein-DNA interactions and larger multicomponent complexes that can be missed by other analytical methods. Here we describe in detail the optimization and utilization of the DAPSTER assay to isolate AP-1 complexes and associated proteins in multimeric complexes bound to the AP-1 DNA element.     

Insights into protein-DNA interactions through structure network analysis

Protein-DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein-DNA complexes, gaining deeper insights into the nature of protein-DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein-DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein-DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein-DNA interface. A predominance of deoxyribose-amino acid clusters in beta-sheet proteins, distinction of the interface clusters in helix-turn-helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein-DNA complexes on the basis of the protein-DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein-DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.     

Detection of protein-protein interactions using nonimmune IgG and BirA-mediated biotinylation

Detection of protein-protein interactions in cells is crucial for understanding the biological functions of proteins, including their roles in signal transduction. However, current methods require specific antibodies both for immunoprecipitation and detection, making them expensive and sometimes unreliable. Here we describe protocols for protein-protein interaction assays that use nonimmune IgG-conjugated Sepharose to precipitate the IgG binding domain (ZZ) fused to the bait protein; the interaction partner is fused to Avitag and biotinylated by BirA so that it can be detected by a one-step blot with Dylight 680 streptavidin to detect the Avitag fusion protein. Since this method does not require specific antibodies and is inexpensive, sensitive, and reliable, it should be useful for detecting protein-protein interactions in cells.     

ChIPOTle: a user-friendly tool for the analysis of ChIP-chip data

ChIPOTle (Chromatin ImmunoPrecipitation On Tiled arrays) takes advantage of two unique properties of ChIP-chip data: the single-tailed nature of the data, caused by specific enrichment but not specific depletion of genomic fragments; and the predictable enrichment of DNA fragments adjacent to sites of direct protein-DNA interaction. Implemented as a Microsoft Excel macro written in Visual Basic, ChIPOTle uses a sliding window approach that yields improvements in the identification of bona fide sites of protein-DNA interaction.     

Cooperativity mutants of the gamma delta resolvase identify an essential interdimer interaction

gamma delta resolvase, a transposon-encoded site-specific recombinase, catalyzes the resolution of the cointegrate intermediate of gamma delta transposition. The recombination reaction involves the formation of a catalytic nucleoprotein complex whose structure is determined by specific protein-DNA and protein-protein interactions. We have isolated many resolvase mutants and have identified four that are unable to mediate a subclass of higher order protein-protein interactions necessary for recombination. This mutant phenotype is characterized by an inability to catalyze recombination, a loss of cooperative binding to res DNA, and a failure to induce looping out of the DNA between two resolvase binding sites within res. The amino acid side chains identified by the cooperativity mutants cluster on a surface of the protein that mediates an interaction between resolvase dimers in a crystallographic tetramer. We have therefore identified a region of resolvase that mediates an interdimer protein-protein interaction necessary for the formation of the recombinogenic synaptic intermediate.     

Tethered-Hopping Model for Protein-DNA Binding and Unbinding Based on Sox2-Oct1-Hoxb1 Ternary Complex Simulations

The sliding and hopping models encapsulate the essential protein-DNA binding process for binary complex formation and dissociation. However, the effects of a cofactor protein on the protein-DNA binding process that leads to the formation of a ternary complex remain largely unknown. Here we investigate the effect of the cofactor Sox2 on the binding and unbinding of Oct1 with the Hoxb1 control element. We simulate the association of Oct1 with Sox2-Hoxb1 using molecular dynamics simulations, and the dissociation of Oct1 from Sox2-Hoxb1 using steered molecular dynamics simulations, in analogy to a hopping event of Oct1. We compare the kinetic and thermodynamic properties of three model complexes (the wild-type and two mutants) in which the Oct1-DNA base-specific interactions or the Sox2-Oct1 protein-protein interactions are largely abolished. We find that Oct1-DNA base-specific interactions contribute significantly to the total interaction energy of the ternary complex, and that nonspecific Oct1-DNA interactions are sufficient for driving the formation of the protein-DNA interface. The Sox2-Oct1 protein-protein binding interface is largely hydrophobic, with remarkable shape complementarity. This interface promotes the formation of the ternary complex and slows the dissociation of Oct1 from its DNA-binding site. We propose a simple two-step reaction model of protein-DNA binding, called the tethered-hopping model, that explains the importance of the cofactor Sox2 and may apply to similar ternary protein-DNA complexes.     

Characterization of functional domains in human Claspin

Claspin is a mediator of the ATR-dependent DNA replication checkpoint in human cells and also promotes DNA replication fork progression and stability. Though Claspin has been shown to bind DNA and co-immunoprecipitate with other replication fork-associated proteins, the specific protein-protein and protein-DNA interactions that are important for Claspin function are not known. We therefore purified several domains of human Claspin and then tested for direct interactions of these fragments with several replication fork-associated proteins and with DNA. Our data show that the N terminus of Claspin binds to the replicative helicase co-factor Cdc45, the Timeless protein and a branched, replication fork-like DNA structure. In contrast, the C terminus of Claspin associates with DNA polymerase epsilon and Rad17-Replication Factor C (RFC). We conclude that multiple protein-DNA and protein-protein interactions may be important for Claspin function during DNA replication and DNA replication checkpoint signaling.     

Mapping of a higher order protein-DNA complex: two kinds of long-range interactions in lambda attL

To map the protein-protein and protein-DNA interactions involved in lambda site-specific recombination, Int cleavage assays with suicide substrates, nuclease protection patterns, gel retardation experiments, and quantitative Western blotting were applied to wild-type attL and attL mutants. The results lead to a model in which one IHF molecule bends the attL DNA and forms a higher order complex with the three bivalent Int molecules required for excisive recombination. It is proposed that each of the Int molecules binds in a unique manner: one bridges two DNA binding sites in cis, one is held via its high affinity amino-terminal DNA binding domain, and the third depends upon protein-protein interactions in addition to its low affinity carboxy-terminal DNA binding domain. This protein-DNA complex contains two unsatisfied DNA binding domains, each with a different sequence specificity, and is well suited to specific interactions with an appropriate recombination partner.