Evaluation of free energy landscape for base-amino acid interactions using ab initio force field and extensive sampling

Structural data of protein-DNA complex show redundancy and flexibility in base-amino acid interactions. To understand the origin of the specificity in protein-DNA recognition, we calculated the interaction free energy, enthalpy, entropy, and minimum energy maps for AT-Asn, GC-Asn, AT-Ser, and GC-Ser by means of a set of ab initio force field with extensive conformational sampling. We found that the most preferable interactions in these pairs are stabilized by hydrogen bonding, and are mainly enthalpy driven. However, minima in the free energy maps are not necessarily the same as those in the minimum energy map or enthalpy maps, due to the entropic effect. The effect of entropy is particularly important in the case of GC-Asn. Experimentally observed structures of base-amino acid interactions are within preferable regions in the calculated free energy maps, where there are many different interaction configurations with similar energy. The full geometry optimization procedure using ab initio molecular orbital method was applied to get the optimal interaction geometries for AT-Asn, GC-Asn, AT-Ser, and GC-Ser. We found that there are various base-amino acid combinations with similar interaction energies. These results suggest that the redundancy and conformational flexibility in the base-amino acid interactions play an important role in the protein-DNA recognition.     

The delta-endotoxin proteins accumulate in Escherichia coli as a protein-DNA complex that can be dissociated by hydrophobic interaction chromatography

The insecticidal protein CryIAc accumulated to form inclusion bodies in Escherichia coli upon overexpression of the cloned gene. The solubilized inclusion bodies contained the delta-endotoxin in association with DNA fragments of about 25 kb. The protein-DNA complex could be dissociated and the delta-endotoxin purified by hydrophobic interaction chromatography on phenyl-Sepharose. The DNA was washed out in the high-salt buffer while the delta-endotoxin was bound to the matrix and was eluted at 4 degrees C by a stepwise decreasing potassium chloride gradient. The DNA-protein complex also contained plasmids harbored by the host strain. The plasmid DNA associated with the complex became competent to transform E. coli only after it was dissociated from the delta-endotoxin. The hydrophobic interaction chromatography provides an efficient method for the purification of DNA-free activated toxin.     

Probing protein-DNA interactions by unzipping a single DNA double helix

We present unzipping force analysis of protein association (UFAPA) as a novel and versatile method for detection of the position and dynamic nature of protein-DNA interactions. A single DNA double helix was unzipped in the presence of DNA-binding proteins using a feedback-enhanced optical trap. When the unzipping fork in a DNA reached a bound protein molecule we observed a dramatic increase in the tension in the DNA, followed by a sudden tension reduction. Analysis of the unzipping force throughout an unbinding "event" revealed information about the spatial location and dynamic nature of the protein-DNA complex. The capacity of UFAPA to spatially locate protein-DNA interactions is demonstrated by noncatalytic restriction mapping on a 4-kb DNA with three restriction enzymes (BsoBI, XhoI, and EcoRI). A restriction map for a given restriction enzyme was generated with an accuracy of approximately 25 bp. UFAPA also allows direct determination of the site-specific equilibrium association constant (K(A)) for a DNA-binding protein. This capability is demonstrated by measuring the cation concentration dependence of K(A) for EcoRI binding. The measured values are in good agreement with previous measurements of K(A) over an intermediate range of cation concentration. These results demonstrate the potential utility of UFAPA for future studies of site-specific protein-DNA interactions.     

An amino acid map of inter-residue contact energies using metric multi-dimensional scaling

We present an amino map based on their inter-residue contact energies using the Miyazawa-Jernigan matrix. This work is based on the method of metric multi-dimensional scaling (MMDS). The MMDS map shows, among other things, that the MJ contact energies imply the hydrophobic-hydrophilic nature of the amino acid residues. With the help of the map we are able to compare and draw inferences from uncorrelated data sets such as BLOSUM and PAM with MJ methods. We also use a hierarchical clustering method on our MMDS distance matrix to group the amino acids and arrive at an optimum number of groups for simplifying the amino acid set.     

Analyzing protein-DNA recognition mechanisms

We present a computational algorithm that can be used to analyze the generic mechanisms involved in protein-DNA recognition. Our approach is based on energy calculations for the full set of base sequences that can be threaded onto the DNA within a protein-DNA complex. It is able to reproduce experimental consensus binding sequences for a variety of DNA binding proteins and also correlates well with the order of measured binding free energies. These results suggest that the crystal structure of a protein-DNA complex can be used to identify all potential binding sequences. By analyzing the energy contributions that lead to base sequence selectivity, it is possible to quantify the concept of direct versus indirect recognition and to identify a new concept describing whether the protein-DNA interaction and DNA deformation terms select optimal binding sites by acting in accord or in disaccord.     

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.     

Geometric analysis and comparison of protein-DNA interfaces: why is there no simple code for recognition?

Structural studies of protein-DNA complexes have shown that there are many distinct families of DNA-binding proteins, and have shown that there is no simple "code" describing side-chain/base interactions. However, systematic analysis and comparison of protein-DNA complexes has been complicated by the diversity of observed contacts, the sheer number of complexes currently available and the absence of any consistent method of comparison that retains detailed structural information about the protein-DNA interface. To address these problems, we have developed geometric methods for characterizing the local structural environment in which particular side-chain/base interactions are observed. In particular, we develop methods for analyzing and comparing spatial relationships at the protein-DNA interface. Our method involves attaching local coordinate systems to the DNA bases and to the C(alpha) atoms of the peptide backbone (these are relatively rigid structural units). We use these tools to consider how the position and orientation of the polypeptide backbone (with respect to the DNA) helps to determine what contacts are possible at any given position in a protein-DNA complex. Here, we focus on base contacts that are made in the major groove, and we use spatial relationships in analyzing: (i) the observed patterns of side-chain/base interactions; (ii) observed helix docking orientations; (iii) family/subfamily relationships among DNA-binding proteins; and (iv) broader questions about evolution, altered specificity mutants and the limits for the design of new DNA-binding proteins. Our analysis, which highlights differences in spatial relationships in different complexes and at different positions in a complex, helps explain why there is no simple, general code for protein-DNA recognition.     

An effective approach for generating a three-Cys<Subscript>2</Subscript>His<Subscript>2</Subscript> zinc-finger-DNA complex model by docking

Background  

Determination of protein-DNA complex structures with both NMR and X-ray crystallography remains challenging in many cases. High Ambiguity-Driven DOCKing (HADDOCK) is an information-driven docking program that has been used to successfully model many protein-DNA complexes. However, a protein-DNA complex model whereby the protein wraps around DNA has not been reported. Defining the ambiguous interaction restraints for the classical three-Cys₂His₂ zinc-finger proteins that wrap around DNA is critical because of the complicated binding geometry. In this study, we generated a Zif268-DNA complex model using three different sets of ambiguous interaction restraints (AIRs) to study the effect of the geometric distribution on the docking and used this approach to generate a newly reported Sp1-DNA complex model.     

Effector-stimulated single molecule protein-DNA interactions of a quorum-sensing system in Sinorhizobium meliloti

Intercellular communication by means of small signal molecules coordinates gene expression among bacteria. This population density-dependent regulation is known as quorum sensing. The symbiotic nitrogen-fixing bacterium Sinorhizobium meliloti Rm1021 possesses the Sin quorum sensing system based on N-acyl homoserine lactones (AHL) as signal molecules. Here, we demonstrate that the LuxR-type regulator ExpR binds specifically to a target sequence in the sinRI locus in the presence of different AHLs with acyl side chains from 8 to 20 carbons. Dynamic force spectroscopy based on the atomic force microscope provided detailed information about the molecular mechanism of binding upon activation by six different AHLs. These single molecule experiments revealed that the mean lifetime of the bound protein-DNA complex varies depending on the specific effector molecule. The small differences between individual AHLs also had a pronounced influence on the structure of protein-DNA interaction: The reaction length of dissociation varied from 2.6 to 5.8 A. In addition, dynamic force spectroscopy experiments indicate that N-heptanoyl-DL-homoserine lactone binds to ExpR but is not able to stimulate protein-DNA interaction.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Rational design of DNA sequence-specific zinc fingers

We developed a rational scheme for designing DNA binding proteins. The scheme was applied for a zinc finger protein and the designed sequences were experimentally characterized with high DNA sequence specificity. Starting with the backbone of a known finger structure, we initially calculated amino acid sequences compatible with the expected structure and the secondary structures of the designed fingers were then experimentally confirmed. The DNA-binding function was added to the designed finger by reconsidering a section of the amino acid sequence and computationally selecting amino acids to have the lowest protein-DNA interaction energy for the target DNA sequences. Among the designed proteins, one had a gap between the lowest and second lowest protein-DNA interaction energies that was sufficient to give DNA sequence-specificity.     

Potential application of FoldX force field based protein modeling in zinc finger nucleases design

Engineered sequence-specific zinc finger nucleases (ZFNs) make the highly efficient modification of eukaryotic genomes possible. However, most current strategies for developing zinc finger nucleases with customized sequence specificities require the construction of numerous tandem arrays of zinc finger proteins (ZFPs), and subsequent largescale in vitro validation of their DNA binding affinities and specificities via bacterial selection. The labor and expertise required in this complex process limits the broad adoption of ZFN technology. An effective computational assisted design strategy will lower the complexity of the production of a pair of functional ZFNs. Here we used the FoldX force field to build 3D models of 420 ZFP-DNA complexes based on zinc finger arrays developed by the Zinc Finger Consortium using OPEN (oligomerized pool engineering). Using nonlinear and linear regression analysis, we found that the calculated protein-DNA binding energy in a modeled ZFP-DNA complex strongly correlates to the failure rate of the zinc finger array to show significant ZFN activity in human cells. In our models, less than 5% of the three-finger arrays with calculated protein-DNA binding energies lower than −13.132 kcal mol⁻¹ fail to form active ZFNs in human cells. By contrast, for arrays with calculated protein-DNA binding energies higher than −5 kcal mol⁻¹, as many as 40% lacked ZFN activity in human cells. Therefore, we suggest that the FoldX force field can be useful in reducing the failure rate and increasing efficiency in the design of ZFNs.     

Site-specific enthalpic regulation of DNA transcription at bacteriophage lambda OR.

Binding of cI repressor to DNA fragments containing the three specific binding sites of the right operator (OR) of bacteriophage lambda was studied in vitro over the temperature range 5-37 degrees C by quantitative footprint titration. The individual-site isotherms, obtained for binding repressor dimers to each site of wild-type OR and to appropriate mutant operator templates, were analyzed for the Gibbs energies of intrinsic binding and pairwise cooperative interactions. It is found that dimer affinity for each of the three sites varies inversely with temperature, i.e., the binding reactions are enthalpy driven, unlike many protein-DNA reactions. By contrast, the magnitude of the pairwise cooperativity terms describing interaction between adjacently site-bound repressor dimers is quite small. This result in combination with the recent finding that repressor monomer-dimer assembly is highly enthalpy driven (with delta H degrees = -16 kcal mol-1) [Koblan, K. S., & Ackers, G. K. (1991) Biochemistry 30, 7817-7821] indicates that the associative contacts between site-bound repressors that mediate cooperativity are unlikely to be the same as those responsible for dimerization. The intrinsic binding enthalpies for all three sites are negative (exothermic) and nearly temperature-invariant, indicating no heat capacity changes on the scale of those inferred in other protein-DNA systems. However, the three operator sites are affected differentially by temperature: the intrinsic binding free energies for sites OR1 and OR3 change in parallel over the entire range, delta H0OR1 = -23.3 +/- 4.0 kcal mol-1 and delta H0OR3 = -22.7 +/- 1.2 kcal mol-1.(ABSTRACT TRUNCATED AT 250 WORDS)