DNA curvature and phosphate neutralization: an important aspect of specific protein binding

A theoretical study is presented of the influence of salt bridges between protein cationic side chains and DNA phosphates on DNA conformation and flexibility. Two DNA sequences are studied containing respectively the HNF3 and CAP binding sites. The effect of salt bridges is modelled by the neutralisation of net phosphate charges for the groups involved in such interactions in the complex. Energy optimised conformations are obtained by molecular mechanics calculations using the JUMNA program. Base sequence dependence is studied by moving the phosphate neutralisation pattern along the sequence, while normal mode analysis is used to evaluate DNA flexibility. The results show that phosphate neutralisation has a strong influence on DNA conformation. For the HNF3 binding sequence, the free oligomer is bent in direction very different from that observed in the protein complex. Phosphate neutralisation changes this direction by 45 degrees to within only 4 degrees of the direction in the complex, without changing the bending angle. For the CAP binding sequence, the free oligomer is already intrinsically curved in the direction favoured by the protein, but phosphate neutralisation increases the bending angle. For both oligomers studied these effects are strongly sequence dependent. It is also shown that oligomer flexibility cannot be explained by a simple superposition of the properties of successive dinucleotide steps. Important long range coupling effects are observed. However, for both sequence studied, phosphate neutralisation however leads to a reduction in oligomer flexibility.     

Crystal structure of the specific DNA-binding domain of Tc3 transposase of C.elegans in complex with transposon DNA.

The crystal structure of the complex between the N-terminal DNA-binding domain of Tc3 transposase and an oligomer of transposon DNA has been determined. The specific DNA-binding domain contains three alpha-helices, of which two form a helix-turn-helix (HTH) motif. The recognition of transposon DNA by the transposase is mediated through base-specific contacts and complementarity between protein and sequence-dependent deformations of the DNA. The HTH motif makes four base-specific contacts with the major groove, and the N-terminus makes three base-specific contacts with the minor groove. The DNA oligomer adopts a non-linear B-DNA conformation, made possible by a stretch of seven G:C base pairs at one end and a TATA sequence towards the other end. Extensive contacts (seven salt bridges and 16 hydrogen bonds) of the protein with the DNA backbone allow the protein to probe and recognize the sequence-dependent DNA deformation. The DNA-binding domain forms a dimer in the crystals. Each monomer binds a separate transposon end, implying that the dimer plays a role in synapsis, necessary for the simultaneous cleavage of both transposon termini.     

Charge neutralization and DNA bending by the Escherichia coli catabolite activator protein

We are interested in the role of asymmetric phosphate neutralization in DNA bending induced by proteins. We describe an experimental estimate of the actual electrostatic contribution of asymmetric phosphate neutralization to the bending of DNA by the Escherichia coli catabolite activator protein (CAP), a prototypical DNA-bending protein. Following assignment of putative electrostatic interactions between CAP and DNA phosphates based on X-ray crystal structures, appropriate phosphates in the CAP half-site DNA were chemically neutralized by methylphosphonate substitution. DNA shape was then evaluated using a semi-synthetic DNA electrophoretic phasing assay. Our results confirm that the unmodified CAP DNA half-site sequence is intrinsically curved by 26° in the direction enhanced in the complex with protein. In the absence of protein, neutralization of five appropriate phosphates increases DNA curvature to 32° (~23% increase), in the predicted direction. Shifting the placement of the neutralized phosphates changes the DNA shape, suggesting that sequence-directed DNA curvature can be modified by the asymmetry of phosphate neutralization. We suggest that asymmetric phosphate neutralization contributes favorably to DNA bending by CAP, but cannot account for the full DNA deformation.     

Modeling DNA deformations induced by minor groove binding proteins

Molecular modeling is used to demonstrate that the major structural deformations of DNA caused by four different minor groove binding proteins, TBP, SRY, LEF-1, and PurR, can all be mimicked by stretching the double helix between two 3'-phosphate groups flanking the binding region. This deformation reproduces the widening of the minor groove and the overall bending and unwinding of DNA caused by protein binding. It also reproduces the principal kinks associated with partially intercalated amino acid side chains, observed with such interactions. In addition, when protein binding involves a local transition to an A-like conformation, phosphate neutralization, via the formation of protein-DNA salt bridges, appears to favor the resulting deformation.     

Surface salt bridges modulate DNA wrapping by the type II DNA-binding protein TF1

The histone-like protein HU is involved in compaction of the bacterial genome. Up to 37 bp of DNA may be wrapped about some HU homologues in a process that has been proposed to depend on a linked disruption of surface salt bridges that liberates cationic side chains for interaction with the DNA. Despite significant sequence conservation between HU homologues, binding sites from 9 to 37 bp have been reported. TF1, an HU homologue that is encoded by Bacillus subtilis bacteriophage SPO1, has nM affinity for 37 bp preferred sites in DNA with 5-hydroxymethyluracil (hmU) in place of thymine. On the basis of electrophoretic mobility shift assays, we show that TF1-DNA complex formation is associated with a net release of only approximately 0.5 cations. The structure of TF1 suggests that Asp13 can form a dehydrated surface salt bridge with Lys23; substitution of Asp13 with Ala increases the net release of cations to approximately 1. These data are consistent with complex formation linked to disruption of surface salt bridges. Substitution of Glu90 with Ala, which would expose Lys87 predicted to contact DNA immediately distal to a proline-mediated DNA kink, causes an increase in affinity and an abrogation of the preference for hmU-containing DNA. We propose that hmU preference is due to finely tuned interactions at the sites of kinking that expose a differential flexibility of hmU- and T-containing DNA. Our data further suggest that the difference in binding site size for HU homologues is based on a differential ability to stabilize the DNA kinks.     

Stability and rigidity/flexibility—Two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile–mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

The salt dependence of the interferon regulatory factor 1 DNA binding domain binding to DNA reveals ions are localized around protein and DNA

The equilibrium dissociation constant of the DNA binding domain of interferon regulatory factor 1 (IRF1 DBD) for its DNA binding site depends strongly on salt concentration and salt type. These dependencies are consistent with IRF1 DBD binding to DNA, resulting in the release of cations from the DNA and both release of anions from the protein and uptake of a cation by the protein. We demonstrated this by utilizing the fact that the release of fluoride from protein upon complex formation does not contribute to the salt concentration dependence of binding and by studying mutants in which charged residues in IRF1 DBD that form salt bridges with DNA phosphates are changed to alanine. The salt concentration dependencies of the dissociation constants of wild-type IRF1 DBD and the mutants R64A, D73A, K75A, and D73A/K75A were measured in buffer containing NaF, NaCl, or NaBr. The salt concentration and type dependencies of the mutants relative to wild-type IRF1 DBD provide evidence of charge neutralization by solution ions for R64 and by a salt bridge between D73 and K75 in buffer containing chloride or bromide salts. These data also allowed us to determine the number, type, and localization of condensed ions around both IRF1 DBD and its DNA binding site.     

Molecular beacon-equilibrium cyclization detection of DNA-protein complexes

Molecular beacon detection of equilibrium cyclization (MBEC) is a novel, high sensitivity technique that can allow DNA-protein complex formation to be studied under diverse conditions in a cost effective and rapid manner that can be adapted to high throughput screening. To demonstrate the ease and utility of applying MBEC to the investigation of the K(D) values of protein-DNA complexes, the sequence-specific Escherichia coli integration host factor (IHF) protein has been used as a test system. Competition between a labeled MBEC DNA construct and unlabeled duplex DNA for IHF binding allows the determination of K(D) values as a function of the DNA duplex sequence. This allows sequence specificity to be monitored while using only a single molecular beacon-labeled DNA. The robustness of MBEC for monitoring protein-DNA complex formation has been further demonstrated by determining the K(D) values as a function of salt concentration to investigate the net number of salt bridges formed in sequence-specific and -nonspecific IHF-DNA complexes. These MBEC results have been compared with those from other approaches.     

Interaction of human SRY protein with DNA: A molecular dynamics study

Molecular dynamics simulations have been conducted to study the interaction of human sex-determining region Y (hSRY) protein with DNA. For this purpose, simulations of the hSRY high mobility group (HMG) domain (hSRY-HMG) with and without its DNA target site, a DNA octamer, and the DNA octamer alone have been carried out, employing the NMR solution structure of hSRY-HMG–DNA complex as a starting model. Analyses of the simulation results demonstrated that the interaction between hSRY and DNA was hydrophobic, just a few hydrogen bonds and only one water molecule as hydrogen-bonding bridge were observed at the protein–DNA interface. These two hydrophobic cores in the hSRY-HMG domain were the physical basis of hSRY-HMG–DNA specific interaction. They not only maintained the stability of the complex, but also primarily caused the DNA deformation. The salt bridges formed between the positive-charged residues of hSRY and phosphate groups of DNA made the phosphate electroneutral, which was advantageous for the deformation of DNA and the formation of a stable complex. We predicted the structure of hSRY-HMG domain in the free state and found that both hSRY and DNA changed their conformations to achieve greater complementarity of geometries and properties during the binding process; that is, the protein increased the angle between its long and short arms to accommodate the DNA, and the DNA became bent severely to adapt to the protein, although the conformational change of DNA was more severe than that of the hSRY-HMG domain. The sequence specificity and the role of residue Met9 are also discussed. Proteins 31:417–433, 1998. © 1998 Wiley-Liss, Inc.     

Role of hydrophobic core on the thermal stability of proteins-molecular dynamics simulations on a single point mutant of sso7d

Abstract The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Salt bridges: geometrically specific, designable interactions

Salt bridges occur frequently in proteins, providing conformational specificity and contributing to molecular recognition and catalysis. We present a comprehensive analysis of these interactions in protein structures by surveying a large database of protein structures. Salt bridges between Asp or Glu and His, Arg, or Lys display extremely well-defined geometric preferences. Several previously observed preferences are confirmed, and others that were previously unrecognized are discovered. Salt bridges are explored for their preferences for different separations in sequence and in space, geometric preferences within proteins and at protein-protein interfaces, co-operativity in networked salt bridges, inclusion within metal-binding sites, preference for acidic electrons, apparent conformational side chain entropy reduction on formation, and degree of burial. Salt bridges occur far more frequently between residues at close than distant sequence separations, but, at close distances, there remain strong preferences for salt bridges at specific separations. Specific types of complex salt bridges, involving three or more members, are also discovered. As we observe a strong relationship between the propensity to form a salt bridge and the placement of salt-bridging residues in protein sequences, we discuss the role that salt bridges might play in kinetically influencing protein folding and thermodynamically stabilizing the native conformation. We also develop a quantitative method to select appropriate crystal structure resolution and B-factor cutoffs. Detailed knowledge of these geometric and sequence dependences should aid de novo design and prediction algorithms.     

Role of hydrophobic core on the thermal stability of proteins - molecular dynamics simulations on a single point mutant of Sso7d abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 s long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Composition-dependent energetic contribution of complex salt bridges to collagen stability

Complex salt bridges, on which three or more charged residues interplay simultaneously, cannot be considered as addition of individual salt bridges. This is still an intriguing problem in protein folding and stability. Here, we used an obligated ABC-type collagen heterotrimer as a platform to study the relationship between energetic contributions and conformational details of three-body complex salt bridges anchored by positively charged residues, K and R. Eight complex salt bridges were constructed by engineering point mutations in the heterotrimer. The circular dichroism measurements showed that the K-anchored complex salt bridges were stronger than the R-anchored ones. The molecular dynamics simulation revealed that both types of salt bridges had distinct dynamic features. The energetic contribution of K-anchored salt bridges was mainly determined by strong single bridges. In the R-anchored complex salt bridges, both side-chain electrostatic interactions and side-chain-backbone hydrogen bonding were involved. An empirical equation was proposed to predict the energetic contributions with high accuracy (R² = 0.93). This work could help us take insights into the mechanisms of composition-dependent behaviors of the complex salt bridges on protein surface.     

The activity of the EcoRV restriction endonuclease is influenced by flanking DNA sequences both inside and outside the DNA-protein complex.

The EcoRV restriction endonuclease cleaves DNA not only at its recognition sequence but also at most other sequences that differ from the recognition site by one base pair. Compared to the reaction at the recognition site, the reactions at noncognate sites are slow but 1 out of the 12 noncognate sites on the plasmid pAT153 is cleaved more than 50 times faster than any other. The increase in the reaction rate at the preferred noncognate site, relative to other sites, was caused by the DNA sequences in the 4 base pairs from either side of the site. For enhanced activity by EcoRV, particular bases were needed immediately adjacent to the site, inside the DNA-protein complex. At these loci, the protein interacts with the phosphate groups in the DNA and the flanking sequence may control the activity of the enzyme by determining the conformation of the DNA, thus aligning the phosphate contacts. But the preferential cleavage also depended on sequences further away from the site, at loci outside the complex. At external positions, beyond the reach of the protein, the EcoRV enzyme required flanking sequences that give rise to flexibility in DNA conformation. These may facilitate the distortion of the DNA required for catalysis by EcoRV.     

Role of Hydrophobic Core on the Thermal Stability of Proteins—Molecular Dynamics Simulations on a Single Point Mutant of Sso7d

Abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Solution structure of the Pseudomonas putida protein PpPutA45 and its DNA complex

Proline utilization A (PutA) is a membrane-associated multifunctional enzyme that catalyzes the oxidation of proline to glutamate in a two-step process. In certain, gram-negative bacteria such as Pseudomonas putida, PutA also acts as an auto repressor in the cytoplasm, when an insufficient concentration of proline is available. Here, the N-terminal residues 1-45 of PutA from P. putida (PpPutA45) are shown to be responsible for DNA binding and dimerization. The solution structure of PpPutA45 was determined using NMR methods, where the protein is shown to be a symmetrical homodimer (12 kDa) consisting of two ribbon-helix-helix (RHH) structures. DNA sequence recognition by PpPutA45 was determined using DNA gel mobility shift assays and NMR chemical shift perturbations (CSPs). PpPutA45 was shown to bind a 14 base-pair DNA oligomer (5'-GCGGTTGCACCTTT-3'). A model of the PpPutA45-DNA oligomer complex was generated using Haddock 2.1. The antiparallel beta-sheet that results from PpPutA45 dimerization serves as the DNA recognition binding site by inserting into the DNA major groove. The dimeric core of four alpha-helices provides a structural scaffold for the beta-sheet from which residues Thr5, Gly7, and Lys9 make sequence-specific contacts with the DNA. The structural model implies flexibility of Lys9 which can make hydrogen bond contacts with either guanine or thymine. The high sequence and structure conservation of the PutA RHH domain suggest interdomain interactions play an important role in the evolution of the protein.     

Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contacts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue–residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue–residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy. Proteins 28:494–514, 1997. © 1997 Wiley-Liss, Inc.     

The B- to A-DNA transition and the reorganization of solvent at the DNA surface

DNA geometry depends on relative humidity. Using the CHARMM22 force field to push B-DNA to A-DNA, a molecular dynamics simulation of a mixed-sequence 24-basepair DNA double-stranded oligomer, starting from B-DNA, was carried out to explore both the mechanism of the transition and the evolution of hydration patterns on the surface of DNA. Over the 11-ns trajectory, the transition recapitulates the slide-first, roll-later mechanism, is opposed by DNA electrostatics, and is favored by an increasing amount of condensed sodium ions. Hydration was characterized by counting the hydrogen bonds between water and DNA, and by the number of water bridges linking two DNA atoms. The number of hydrogen bonds between water and DNA remains constant during the transition, but there is a 40% increase in the number of water bridges, in agreement with the principle of economy of hydration. Water bridges emerge as delicate sensors of both structure and dynamics of DNA. Both local flexibility and the frustration of the water network on the surface of DNA probably account for the low populations and short residence times of the bridges, and for the lubricant role of water in ligand-DNA interactions.     

Surface salt bridges stabilize the GCN4 leucine zipper.

We present a study of the role of salt bridges in stabilizing a simplified tertiary structural motif, the coiled-coil. Changes in GCN4 sequence have been engineered that introduce trial patterns of single and multiple salt bridges at solvent exposed sites. At the same sites, a set of alanine mutants was generated to provide a reference for thermodynamic analysis of the salt bridges. Introduction of three alanines stabilizes the dimer by 1.1 kcal/mol relative to the wild-type. An arrangement corresponding to a complex type of salt bridge involving three groups stabilizes the dimer by 1.7 kcal/ mol, an apparent elevation of the melting temperature relative to wild type of about 22 degrees C. While identifying local from nonlocal contributions to protein stability is difficult, stabilizing interactions can be identified by use of cycles. Introduction of alanines for side chains of lower helix propensity and complex salt bridges both stabilize the coiled-coil, so that combining the two should yield melting temperatures substantially higher than the starting species, approaching those of thermophilic sequences.