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

cAMP Activation of the cAMP Receptor Protein, a Model Bacterial Transcription Factor

The active and inactive structures of the Escherichia coli cAMP receptor protein (CRP), a model bacterial transcription factor, are compared to generate a paradigm in the cAMP-induced activation of CRP. The resulting paradigm is shown to be consistent with numerous biochemical studies of CRP and CRP*, a group of CRP mutants displaying cAMP-free activity. The cAMP affinity of CRP is dictated by two factors: (i) the effectiveness of the cAMP pocket and (ii) the protein equilibrium of apo-CRP. How these two factors interplay in determining the cAMP affinity and cAMP specificity of CRP and CRP* mutants are discussed. Both the current understanding and knowledge gaps of CRP-DNA interactions are also described. This review ends with a list of several important CRP issues that need to be addressed in the future.

Study of DNA binding and bending by Bacillus subtilis GabR, a PLP-dependent transcription factor

BACKGROUND

GabR is a transcriptional regulator belonging to the MocR/GabR family, characterized by a N-terminal wHTH DNA-binding domain and a C-terminal effector binding and/or oligomerization domain, structurally homologous to aminotransferases (ATs). In the presence of γ-aminobutyrate (GABA) and pyridoxal 5'-phosphate (PLP), GabR activates the transcription of gabT and gabD genes involved in GABA metabolism.

METHODS

Here we report a biochemical and atomic force microscopy characterization of Bacillus subtilis GabR in complex with DNA. Complexes were assembled in vitro to study their stoichiometry, stability and conformation.

RESULTS

The fractional occupancy of the GabR cognate site suggests that GabR binds as a dimer with Kd of 10nM. Upon binding GabR bends the DNA by 80° as measured by anomalous electrophoretic mobility. With GABA we observed a decrease in affinity and conformational rearrangements compatible with a less compact nucleo-protein complex but no changes of the DNA bending angle. By employing promoter and GabR mutants we found that basic residues of the positively charged groove on the surface of the AT domain affect DNA affinity.

CONCLUSIONS

The present data extend current understanding of the GabR-DNA interaction and the effect of GABA and PLP. A model for the GabR-DNA complex, corroborated by a docking simulation, is proposed.

GENERAL SIGNIFICANCE

Characterization of the GabR DNA binding mode highlights the key role of DNA bending and interactions with bases outside the canonical direct repeats, and might be of general relevance for the action mechanism of MocR transcription factors.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

DNA phosphate crowding correlates with protein cationic side chain density and helical curvature in protein/DNA crystal structures

Sequence-specific binding of proteins to their DNA targets involves a complex spectrum of processes that often induce DNA conformational variation in the bound complex. The forces imposed by protein binding that cause the helical deformations are intimately interrelated and difficult to parse or rank in importance. To investigate the role of electrostatics in helical deformation, we quantified the relationship between protein cationic residue density (Cpc) and DNA phosphate crowding (Cpp). The correlation between Cpc and Cpp was then calculated for a subset of 58 high resolution protein-DNA crystal structures. Those structures containing strong Cpc/Cpp correlation (>±0.25) were likely to contain DNA helical curvature. Further, the correlation factor sign predicted the direction of helical curvature with positive (16 structures) and negative (seven structures) correlation containing concave (DNA curved toward protein) and convex (DNA curved away from protein) curvature, respectively. Protein-DNA complexes without significant Cpc/Cpp (36 structures) correlation (-0.25<0<0.25) tended to contain DNA without significant curvature. Interestingly, concave and convex complexes also include more arginine and lysine phosphate contacts, respectively, whereas linear complexes included essentially equivalent numbers of Lys/Arg phosphate contacts. Together, these findings suggest an important role for electrostatic interactions in protein-DNA complexes involving helical curvature.

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes -2.3 kcal mol(-1) to -0.9 kcal mol(-1) (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of +1.3 kcal mol(-1). Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.

Examining the ribonuclease H primer grip of HIV-1 reverse transcriptase by charge neutralization of RNA/DNA hybrids

The crystal structure of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) bound to an RNA/DNA hybrid reveals an extensive network of contacts with the phosphate backbone of the DNA strand ∼4–9 bp downstream from the ribonuclease H (RNase H) catalytic center. Collectively designated as ‘the RNase H primer grip’, this motif contains a phosphate binding pocket analogous to the human and Bacillus halodurans RNases H. The notion that the RNase H primer grip mediates the trajectory of RNA/DNA hybrids accessing the RNase H active site suggests that locally neutralizing the phosphate backbone may be exploited to manipulate nucleic acid flexibility. To examine this, we introduced single and tandem methylphosphonate substitutions through the region of the DNA primer contacted by the RNase H primer grip and into the RNase H catalytic center. The ability of mutant hybrids to support RNase H and DNA polymerase activity was thereafter examined. In addition, site-specific chemical footprinting was used to evaluate movement of the DNA polymerase and RNase H domains. We show here that minor alteration to the RNase H primer can have a dramatic effect on enzyme positioning, and discuss these findings in light of recent crystallography of human RNase H containing an RNA/DNA hybrid.

Water-mediated interactions in the CRP-cAMP-DNA complex: does water mediate sequence-specific binding at the DNA primary-kink site?

The cyclic AMP receptor protein (CRP) of Escherichia coli binds preferentially to DNA sequences possessing a T:A base pair at position 6 (at which the DNA becomes kinked), but with which it does not form any direct interactions. It has been proposed that indirect readout is involved in CRP-DNA binding, in which specificity for this base pair is primarily related to sequence effects on the energetic susceptibility of the DNA to kink formation. In the current study, the possibility of contributions to indirect readout by water-mediated hydrogen bonding of CRP with the T:A base pair was investigated. A 1.0 ns molecular dynamics simulation of the CRP-cAMP-DNA complex in explicit solvent was performed, and assessed for water-mediated CRP-DNA hydrogen bonds; results were compared to several X-ray crystal structures of comparable complexes. While several water-mediated CRP-DNA hydrogen bonds were identified, none of these involved the T:A base pair at position 6. Therefore, the sequence specificity for this base pair is not likely enhanced by water-mediated hydrogen bonding with the CRP.

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

Positively charged C-terminal subdomains of EcoRV endonuclease: contributions to DNA binding, bending, and cleavage

The carboxy-terminal subdomains of the homodimeric EcoRV restriction endonuclease each bear a net charge of +4 and are positioned on the inner concave surface of the 50 degree DNA bend that is induced by the enzyme. A complete kinetic and structural analysis of a truncated EcoRV mutant lacking these domains was performed to assess the importance of this diffuse charge in facilitating DNA binding, bending, and cleavage. At the level of formation of an enzyme-DNA complex, the association rate for the dimeric mutant enzyme was sharply decreased by 10(3)-fold, while the equilibrium dissociation constant was weakened by nearly 10(6)-fold compared with that of wild-type EcoRV. Thus, the C-terminal subdomains strongly stabilize the enzyme-DNA ground-state complex in which the DNA is known to be bent. Further, the extent of DNA bending as observed by fluorescence resonance energy transfer was also significantly decreased. The crystal structure of the truncated enzyme bound to DNA and calcium ions at 2.4 A resolution reveals that the global fold is preserved and suggests that a divalent metal ion crucial to catalysis is destabilized in the active site. This may explain the 100-fold decrease in the rate of metal-dependent phosphoryl transfer observed for the mutant. These results show that diffuse positive charge associated with the C-terminal subdomains of EcoRV plays a key role in DNA association, bending, and cleavage.

Changes in DNA bending and flexing due to tethered cations detected by fluorescence resonance energy transfer

Local DNA deformation arises from an interplay among sequence-related base stacking, intrastrand phosphate repulsion, and counterion and water distribution, which is further complicated by the approach and binding of a protein. The role of electrostatics in this complex chemistry was investigated using tethered cationic groups that mimic proximate side chains. A DNA duplex was modified with one or two centrally located deoxyuracils substituted at the 5-position with either a flexible 3-aminopropyl group or a rigid 3-aminopropyn-1-yl group. End-to-end helical distances and duplex flexibility were obtained from measurements of the time-resolved Förster resonance energy transfer between 5′- and 3′-linked dye pairs. A novel analysis utilized the first and second moments of the G(t) function, which encompasses only the energy transfer process. Duplex flexibility is altered by the presence of even a single positive charge. In contrast, the mean 5′–3′ distance is significantly altered by the introduction of two adjacently tethered cations into the double helix but not by a single cation: two adjacent aminopropyl groups decrease the 5′–3′ distance while neighboring aminopropynyl groups lengthen the helix.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ∼30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Catabolite activator protein: DNA binding and transcription activation

Recently determined structures of the Escherichia coli catabolite activator protein (CAP) in complex with DNA, and in complex with the RNA polymerase alpha subunit C-terminal domain (alphaCTD) and DNA, have yielded insights into how CAP binds DNA and activates transcription. Comparison of multiple structures of CAP-DNA complexes has revealed the contributions of direct and indirect readout to DNA binding by CAP. The structure of the CAP-alphaCTD-DNA complex has provided the first structural description of interactions between a transcription activator and its functional target within the general transcription machinery. Using the structure of the CAP-alphaCTD-DNA complex, the structure of an RNA polymerase-DNA complex, and restraints from biophysical, biochemical and genetic experiments, it has been possible to construct detailed three-dimensional models of intact class I and class II transcription activation complexes.

Reflections on apparent DNA bending by charge variants of bZIP proteins

Basic-leucine zipper (bZIP) proteins have been studied intensely as transcription factors. It has been proposed that the bZIP domain might modulate transcription activation through the induction of conformational changes in the DNA binding site. We have been interested in using bZIP peptides as convenient models with which to study the role of asymmetric phosphate neutralization in DNA bending. DNA bending experiments have yielded discordant results for bZIP peptides studied by electrophoretic- vs solution-based assays. We review the history of DNA bending assays involving bZIP peptides and introduce the reader to examples of discordant results. Our recent published experiments designed to clarify this field of study will then be reviewed. The engineering of protein fusions has established that electrophoretic phasing assays are relatively insensitive to precise protein structure/conformation and instead appear to report DNA bending, as influenced by protein charge. New applications of time-resolved fluorescence resonance energy transfer (FRET) have allowed for the first time corroboration of electrophoretic phasing assays with solution-based FRET measurements. We report that two conventional DNA bending assays that rely on DNA ligation cannot be applied to analysis of the bZIP peptides we studied due to ligation inhibition.