Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain–DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1–DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein–DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein–DNA complexes. The order of stability of protein–DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein–DNA complexes. Among specific protein–DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.     

Backbone dynamics of sequence specific recognition and binding by the yeast Pho4 bHLH domain probed by NMR

Backbone dynamics of the basic/helix-loop-helix domain of Pho4 from Saccharomyces cerevisae have been probed by NMR techniques, in the absence of DNA, nonspecifically bound to DNA and bound to cognate DNA. Alpha proton chemical shift indices and nuclear Overhauser effect patterns were used to elucidate the secondary structure in these states. These secondary structures are compared to the co-crystal complex of Pho4 bound to a cognate DNA sequence (Shimizu T. Toumoto A, Ihara K, Shimizu M, Kyogou Y, Ogawa N, Oshima Y, Hakoshima T, 1997, EMBO J 15: 4689-4697). The dynamic information provides insight into the nature of this DNA binding domain as it progresses from free in solution to a specifically bound DNA complex. Relative to the unbound form, we show that formation of either the nonspecific and cognate DNA bound complexes involves a large change in conformation and backbone dynamics of the basic region. The nonspecific and cognate complexes, however, have nearly identical secondary structure and backbone dynamics. We also present evidence for conformational flexibility at a highly conserved glutamate basic region residue. These results are discussed in relation to the mechanism of sequence specific recognition and binding.     

Mean field analysis of FKBP12 complexes with FK506 and rapamycin: Implications for a role of crystallographic water molecules in molecular recognition and specificity

Mean field analysis of FKBP12 complexes with FK506 and rapamycin has been performed by using structures obtained from molecular docking simulations on a simple, yet robust molecular recognition energy landscape. When crystallographic water molecules are included in the simulations as an extension of the FKBP12 protein surface, there is an appreciable stability gap between the energy of the native FKBP12–FK506 complex and energies of conformations with the “native-like” binding mode. By contrast, the energy spectrum of the FKBP12–rapamycin complex is dense regardless of the presence of the water molecules. The stability gap in the FKBP12–FK506 system is determined by two critical water molecules from the effector region that participate in a network of specific hydrogen bond interactions. This interaction pattern protects the integrity and precision of the composite ligand-protein effector surface in the binary FKBP12–FK506 complex and is preserved in the crystal structure of the FKBP12–FK506–calcineurin ternary complex. These features of the binding energy landscapes provide useful insights into specific and nonspecific aspects of FK506 and rapamycin recognition. Proteins 28:313–324, 1997. © 1997 Wiley-Liss, Inc.     

NMR studies of restriction enzyme-DNA interactions: role of conformation in sequence specificity

Sequence specific DNA binding proteins are thought to adopt distinct conformations when binding to target (cognate) and nontarget (noncognate) sequences. There is both biochemical and crystallographic evidence that this behavior is important in mediating sequence recognition by the Mg(II)-dependent type II restriction enzymes. Despite this, there are few systematic comparisons of the structural behavior of these enzymes in various complexes. Here, (1)H-(15)N HSQC NMR spectroscopy is applied to PvuII endonuclease (2 x 18 kDa) in an effort to better understand the relationship between sequence recognition and enzyme conformational behavior. Spectra of the free enzyme collected in the absence and presence of metal ions indicate that while there is a modest backbone conformational response upon binding Ca(II), this does not occur with Mg(II). Substrate binding itself is accompanied by very dramatic spectral changes consistent with a large-scale conformational response. HSQC spectra of the enzyme bound to cognate (specific) and noncognate (nonspecific) oligonucleotides in the presence of Ca(II) are dramatically distinct, revealing for the first time the structural uniqueness of a PvuII cognate complex in solution. The strong correlation between NMR spectral overlap and crystallographic data (C(alpha) rmsd) permits characterization of the nonspecific PvuII complex as being more similar to the free enzyme than to the specific complex. Collectively, these data support the notion that it is the DNA, not the metal ion, which promotes a unique conformational response by the enzyme. It therefore follows that the principle role of metal ions in complex formation is one of driving substrate affinity and stability rather than conformationally priming the enzyme for substrate binding and sequence recognition. These results not only provide valuable insights into the mechanism of protein-DNA interactions but also demonstrate the utility of NMR spectroscopy in structure-function studies of these representative nucleic acid systems.     

Differences between EcoRI nonspecific and "star" sequence complexes revealed by osmotic stress

The binding of the restriction endonuclease EcoRI to DNA is exceptionally specific. Even a single basepair change ("star" sequence) from the recognition sequence, GAATTC, decreases the binding free energy of EcoRI to values nearly indistinguishable from nonspecific binding. The difference in the number of waters sequestered by the protein-DNA complexes of the "star" sequences TAATTC and CAATTC and by the specific sequence complex determined from the dependence of binding free energy on water activity is also practically indistinguishable at low osmotic pressures from the 110 water molecules sequestered by nonspecific sequence complexes. Novel measurements of the dissociation rates of noncognate sequence complexes and competition equilibrium show that sequestered water can be removed from "star" sequence complexes by high osmotic pressure, but not from a nonspecific complex. By 5 Osm, the TAATTC "star" sequence complex has lost almost 90 of the approximately 110 waters initially present. It is more difficult to remove water from the CAATTC "star" sequence complex. The sequence dependence of water loss correlates with the known sequence dependence of "star" cleavage activity.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

Dynamics of repressor-operator recognition: the Tn10-encoded tetracycline resistance control

Binding of the Tet repressor to nonspecific and specific DNA leads to quenching of the Tet fluorescence by approximately 22% and approximately 35%, respectively. This effect is used for a direct, quantitative characterization of the binding equilibria and dynamics involved in the recognition of the operator by its repressor. From the dependence of the nonspecific binding constant on the ion concentration, it is concluded that nonspecific binding is almost completely driven by the entropy change resulting from the release of three to four Na+ ions from the double helix upon protein binding. Formation of the specific complex is driven by a higher entropy term resulting from the release of seven to eight Na+ ions and in addition by a free energy term of -33 kJ/mol from nonelectrostatic interactions, which are attributed to the specific contacts. The dynamics of the repressor-operator recognition are resolved by stopped-flow measurements at various salt concentrations and for different DNA chain lengths into two separate steps. The first step follows a second-order mechanism and results in an intermediate complex associated with formation of about three to four electrostatic contacts between protein and DNA; apparently, this complex is equivalent to the nonspecific complex. The existence of an intermediate is also indicated by experiments in mixed Na+-Mg2+ buffers, which can be described with high accuracy by competition of Mg2+ and protein. The intermediate complex is formed at a rate of 3 X 10(8) M-1 s-1 and is converted in the second reaction step to the specific complex with a rate constant of 6 X 10(4) s-1, which is almost independent of the salt concentration. Our interpretation and the parameters obtained from our model are confirmed by competition of nonspecific DNA with operator DNA for repressor binding. The observed maximal rate constant of 3 X 10(8) M-1 s-1 is very close to theoretical predictions for the association without a sliding mechanism. The very small dependence of the observed rate constants on the chain length shows that the Tet repressor is not able to slide over any substantial distance even at low salt concentrations. The question of a potential contribution from sliding under our experimental conditions is critically discussed. The absence of sliding in the case of the Tet repressor under physiological conditions is compared with the high sliding efficiency of the lac repressor and is discussed with respect to possible molecular mechanisms of sliding in relation to biological function.     

Selection of DNA binding sites by regulatory proteins. Functional specificity and pseudosite competition

The frequency of base-pair occurrence in a set of recognition sequences for a particular DNA-binding protein is strongly related to the contributions to the binding free energy from the individual base pairs. Thus from the statistics of base-pair choice, it is possible to estimate the relative binding strengths of any base-pair sequences and to predict the effect of point mutations in specific sites. On the same basis, one can describe the binding properties of random DNA sequences and thereby the expected competitive effects from all the nonspecific DNA sites in the genome of a living cell. The statistical selection theory [Berg & von Hippel.J. Mol. Biol. 193 (1987) 723-750] describing these relations is extended and tested with computer simulations. The theory is shown to hold up well also in the case when base pairs contribute cooperatively to the binding interaction. The simulations also demonstrate the effects of the statistical small-sample uncertainty that appears due to the limited size of all sets of recognition sites identified.     

The EcoRI restriction endonuclease with bacteriophage lambda DNA. Equilibrium binding studies.

The EcoRI restriction endonuclease was found by the filter binding technique to form stable complexes, in the absence of Mg2+, with the DNA from derivatives of bacteriophage lambda that either contain or lack EcoRI recognition sites. The amount of complex formed at different enzyme concentrations followed a hyperbolic equilibrium-binding curve with DNA molecules containing EcoRI recognition sites, but a sigmoidal equilibrium-binding curve was obtained with a DNA molecule lacking EcoRI recognition sites. The EcoRI enzyme displayed the same affinity for individual recognition sites on lambda DNA, even under conditions where it cleaves these sites at different rates. The binding of the enzyme to a DNA molecule lacking EcoRI sites was decreased by Mg2+. These observations indicate that (a) the EcoRI restriction enzyme binds preferentially to its recognition site on DNA, and that different reaction rates at different recognition sites are due to the rate of breakdown of this complex; (b) the enzyme also binds to other DNA sequences, but that two molecules of enzyme, in a different protein conformation, are involved in the formation of the complex at non-specific consequences; (c) the different affinities of the enzyme for the recognition site and for other sequences on DNA, coupled with the different protein conformations, account for the specificity of this enzyme for the cleavage of DNA at this recognition site; (d) the decrease in the affinity of the enzyme for DNA, caused by Mg2+, liberates binding energy from the DNA-protein complex that can be used in the catalytic reaction.     

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of Uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k_cat and reduced K_m relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K_m) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

Interaction of restriction and modification enzyme EcoRII with synthetic DNA fragments. VII. The study of complex-formation of endonuclease EcoRII with substrates containing natural and modified recognition sites

As shown by a nitrocellulose filter binding assay, in the absence of Mg2+ EcoRII restriction endonuclease binds specifically to a set of synthetic concatemer DNA duplexes of varying chain length, containing natural and modified recognition sites of this enzyme. The binding of the substrates with the central AT, TT or AA-pair in the recognition site decreases at AT greater than TT much greater than AA. Substitution of the pyrophosphate bond at the cleavage site for the phosphodiester or phosphoramide bond produces little influence on the stability of the complexes. The affinity of the enzyme for nonspecific sites is two orders of magnitude less than that for the specific EcoRII sequences. Equilibrium association constant for a substrate with one recognition site is 3.9 X 10(8) M-1. Addition of Mg2+ leads to the destabilization of the EcoRII endonuclease complex with DNA duplex, containing pyrophosphate bonds. The dissociation rate constants and the lifetime of the EcoRII endonuclease--synthetic substrates complexes have been determined.     

In vitro initial attachment of HIV-1 integrase to viral ends: control of the DNA specific interaction by the oligomerization state

HIV-1 integrase (IN) oligomerization and DNA recognition are crucial steps for the subsequent events of the integration reaction. Recent advances described the involvement of stable intermediary complexes including dimers and tetramers in the in vitro integration processes, but the initial attachment events and IN positioning on viral ends are not clearly understood. In order to determine the role of the different IN oligomeric complexes in these early steps, we performed in vitro functional analysis comparing IN preparations having different oligomerization properties. We demonstrate that in vitro IN concerted integration activity on a long DNA substrate containing both specific viral and nonspecific DNA sequences is highly dependent on binding of preformed dimers to viral ends. In addition, we show that IN monomers bound to nonspecific DNA can also fold into functionally different oligomeric complexes displaying nonspecific double-strand DNA break activity in contrast to the well known single strand cut catalyzed by associated IN. Our results imply that the efficient formation of the active integration complex highly requires the early correct positioning of monomeric integrase or the direct binding of preformed dimers on the viral ends. Taken together the data indicates that IN oligomerization controls both the enzyme specificity and activity.     

Role of substrate binding forces in exchange-only transport systems: I. Transition-state theory

Summary An analysis of transition-state models for exchange-only transport shows that substrate binding forces, carrier conformational changes, and coupled substrate flow are interrelated. For a system to catalyze exchange but not net transport, addition of the substrate must convert the carrier from an immobile to a mobile form. The reduction in the energy barrier to movement is necessarily paid for out of the intrinsic binding energy between the substrate and the transport site, and is dependent on the formation of two different types of complex: a loose complex initially and a tight complex in the transition state in carrier movement. Hence the site should at first be incompletely organized for optimal binding but, following a conformational change, complementary to the substrate structure in the transition state. The conformational change, which may involve the whole protein, would be induced by cooperative interactions between the substrate and several groups within the site, involving a chelate effect. The tightness of coupling, i.e., the ratio of exchange to net transport, is directly proportional to the increased binding energy in the transition state, a relationship which allows the virtual substrate dissociation constant in the transition state to be calculated from experimental rate and half-saturation constants. Because the transition state is present in minute amount, strong bonding here does not enhance the substrate's affinity, and specificity may, therefore, be expressed in maximum exchange rates alone. However, where substrates largely convert the carrier to a transport intermediate whose mobility is the same with all substrates, specificity is also expressed in affinity. Hence the expression of substrate specificity provides evidence on the translocation mechanism.     

Thermodynamic parameters governing interaction of EcoRI endonuclease with specific and nonspecific DNA sequences

Equilibrium binding of EcoRI endonuclease to DNA has been analyzed by nitrocellulose filter and preferential DNA cleavage methods. Association constants for pBR322 and a 34-base pair molecule containing the EcoRI site of this plasmid in a central position were determined to be 1.9 X 10(11) M-1 and 1.0 X 10(11) M-1 at 37 degrees C, respectively, with the stoichiometry of binding being 0.8 +/- 0.1 mol of endonuclease dimer per mol of DNA. In contrast, the affinity of the enzyme for a pBR322 derivative from which the EcoRI site has been deleted is 3.2 X 10(9) M-1 as judged by competitive binding experiments. If it is assumed that each base pair can define the beginning of a nonspecific binding site, this value corresponds to an affinity for nonspecific sites of 7.4 X 10(5) M-1. Furthermore, the affinity of the endonuclease for the EcoRI-methylated sequence is at least three orders of magnitude less than that for the unmodified recognition site. The dependence on temperature and ionic strength of the equilibrium constant governing specific interactions has also been examined. The temperature dependence of the reaction indicates that entropy increase accounts for 70% of the free energy of specific binding at 37 degrees C. Affinity of the endonuclease for the EcoRI site is highly dependent on NaCl concentration. Analysis of this dependence according to the theory of Record and colleagues (Record, T. M., Jr., Lohman, T. M., and deHaseth, P. (1976) J. Mol. Biol. 107, 145-158) has implicated 8 ion pairs in the stability of specific complexes, a value identical with the number of phosphate contacts determined by ethylation interference analysis (Lu, A. L., Jack, W. E., and Modrich, P. (1981) J. Biol. Chem. 256, 13200-13206). Extrapolation to 1 M NaCl suggests that nonelectrostatic interactions account for 40% of the free energy change associated with specific complex formation.     

Differences in hydration coupled to specific and nonspecific competitive binding and to specific DNA Binding of the restriction endonuclease BamHI

Using the osmotic stress technique together with a self-cleavage assay we measure directly differences in sequestered water between specific and nonspecific DNA-BamHI complexes as well as the numbers of water molecules released coupled to specific complex formation. The difference between specific and nonspecific binding free energy of the BamHI scales linearly with solute osmolal concentration for seven neutral solutes used to set water activity. The observed osmotic dependence indicates that the nonspecific DNA-BamHI complex sequesters some 120-150 more water molecules than the specific complex. The weak sensitivity of the difference in number of waters to the solute identity suggests that these waters are sterically inaccessible to solutes. This result is in close agreement with differences in the structures determined by x-ray crystallography. We demonstrate additionally that when the same solutes that were used in competition experiments are used to probe changes accompanying the binding of free BamHI to its specific DNA sequence, the measured number of water molecules released in the binding process is strikingly solute-dependent (with up to 10-fold difference between solutes). This result is expected for reactions resulting in a large change in a surface exposed area.     

Human adenovirus proteinase: DNA binding and stimulation of proteinase activity by DNA.

The interaction of the human adenovirus proteinase (AVP) with various DNAs was characterized. AVP requires two cofactors for maximal activity, the 11-amino acid residue peptide from the C-terminus of adenovirus precursor protein pVI (pVIc) and the viral DNA. DNA binding was monitored by changes in enzyme activity or by fluorescence anisotropy. The equilibrium dissociation constants for the binding of AVP and AVP-pVIc complexes to 12-mer double-stranded (ds) DNA were 63 and 2.9 nM, respectively. DNA binding was not sequence specific; the stoichiometry of binding was proportional to the length of the DNA. Three molecules of the AVP-pVIc complex bound to 18-mer dsDNA and six molecules to 36-mer dsDNA. When AVP-pVIc complexes bound to 12-mer dsDNA, two sodium ions were displaced from the DNA. A Delta of -4.6 kcal for the nonelectrostatic free energy of binding indicated that a substantial component of the binding free energy results from nonspecific interactions between the AVP-pVIc complex and DNA. The cofactors altered the interaction of the enzyme with the fluorogenic substrate (Leu-Arg-Gly-Gly-NH)2-rhodamine. In the absence of any cofactor, the Km was 94.8 microM and the kcat was 0.002 s(-1). In the presence of adenovirus DNA, the Km decreased 10-fold and the kcat increased 11-fold. In the presence of pVIc, the Km decreased 10-fold and the kcat increased 118-fold. With both cofactors present, the kcat/Km ratio increased 34000-fold, compared to that with AVP alone. Binding to DNA was coincident with stimulation of proteinase activity by DNA. Although other proteinases have been shown to bind to DNA, stimulation of proteinase activity by DNA is unprecedented. A model is presented suggesting that AVP moves along the viral DNA looking for precursor protein cleavage sites much like RNA polymerase moves along DNA looking for a promoter.     

Mg2+ confers DNA binding specificity to the EcoRV restriction endonuclease.

The EcoRV mutant D90A which carries an amino acid substitution in its active center does not cleave DNA. Therefore, it is possible to perform DNA binding experiments with the EcoRV-D90A mutant both in the absence and in the presence of Mg2+. Like wild-type EcoRV [Taylor et al. (1991) Biochemistry 30, 8743-8753], it does not show a pronounced specificity for binding to its recognition site in the absence of Mg2+ as judged by the appearance of multiple shifted bands in an electrophoretic mobility shift assay with a 377-bp DNA fragment carrying a single EcoRV recognition sequence. In the presence of Mg2+, however, only one band corresponding to a 1:1 complex appears even with a high excess of protein over DNA. This complex most likely is the specific one, because its formation is suppressed much more effectively by a 13-bp oligodeoxynucleotide with an EcoRV site than by a corresponding oligodeoxynucleotide without an EcoRV site. The preferential interaction of the EcoRV-D90A mutant with specific DNA in the presence of Mg2+ was also demonstrated directly: a 20-bp oligodeoxynucleotide with an EcoRV site is bound with KAss = 4 x 10(8) M-1, while a corresponding oligodeoxynucleotide without an EcoRV site is bound with KAss less than or equal to 1 x 10(5) M-1. From these data it appears that Mg2+ confers DNA binding specificity to this mutant by lowering the affinity to nonspecific sites and raising the affinity to specific sites as compared to binding in the absence of Mg2+. It is concluded that this is also true for wild-type EcoRV.     

DNA-protein complexes during attachment-site synapsis in Mu DNA transposition.

Initial events in Mu DNA transposition involve specific recognition of Mu DNA ends (att sites) and an internal enhancer site by the Mu transposase (A protein). This interaction between A protein and Mu DNA sequences present on a supercoiled DNA substrate leads to the formation of a stable synaptic complex in which the att ends are nicked, prior to DNA strand transfer. This study examines the properties of a synaptic complex proficient for DNA transposition. We show that the A protein binds as a monomer to its binding sites, and causes the DNA to bend through approximately 90 degrees at each site. All six att binding sites (three at each Mu end) are occupied by A within the synaptic complex. Three of these sites are loosely held and can be emptied of A upon challenge with heparin. A synaptic complex with only three sites occupied is stable and is fully competent in the subsequent strand-transfer step of transposition.