Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: DNA binding specificity based on energetics of DNA kinking

Molecular stripping underpins derepression of a toxin-antitoxin system

Transcription factors control gene expression; among these, transcriptional repressors must liberate the promoter for derepression to occur. Toxin-antitoxin (TA) modules are bacterial elements that autoregulate their transcription by binding the promoter in a T:A ratio-dependent manner, known as conditional cooperativity. The molecular basis of how excess toxin triggers derepression has remained elusive, largely because monitoring the rearrangement of promoter-repressor complexes, which underpin derepression, is challenging. Here, we dissect the autoregulation of the Salmonella enterica tacAT3 module. Using a combination of assays targeting DNA binding and promoter activity, as well as structural characterization, we determine the essential TA and DNA elements required to control transcription, and we reconstitute a repression-to-derepression path. We demonstrate that excess toxin triggers molecular stripping of the repressor complex off the DNA through multiple allosteric changes causing DNA distortion and ultimately leading to derepression. Thus, our work provides important insight into the mechanisms underlying conditional cooperativity.

Allosteric Response of DNA Recognition Helices of Catabolite Activator Protein to cAMP and DNA Binding

The homodimeric catabolite activator protein (CAP) regulates the transcription of several bacterial genes based on the cellular concentration of cyclic adenosine monophosphate (cAMP). The binding of cAMP to CAP triggers allosteric communication between the cAMP binding domains (CBD) and DNA binding domains (DBD) of CAP, which entails repositioning of DNA recognition helices (F-helices) in the DBD to dock favorably to the target DNA. Despite considerable progress, much remains to be understood about the mechanistic details of DNA recognition by CAP and about the map of allosteric pathways involved in CAP-mediated gene transcription. The present study uses molecular dynamics and umbrella sampling simulations to investigate the mechanism of cAMP- and DNA-induced changes in the conformation and energetics of F-helices observed during the allosteric regulation of CAP by cAMP and the subsequent binding to the DNA promoter region. Using novel collective variables, the free energy profiles associated with the orientation and dynamics of F-helices in the unliganded, cAMP-bound, and cAMP-DNA-bound states of CAP are calculated and compared. The binding-induced alterations in the resultant free energy profiles reveal important flexibility constraints imposed on DBD upon cAMP and DNA binding. A comprehensive analysis of residue-wise interaction maps reveals potential allosteric pathways between CBD and DBD that facilitate the allosteric transduction of regulatory signals in CAP. The revelation that the predicted allosteric pathways crisscross the intersubunit interface offers important clues on the microscopic origin of the intersubunit cooperativity and dimer stability of CAP.

Singular value decomposition for the correlation of atomic fluctuations with arbitrary angle

Many proteins exhibit a critical property called allostery, which enables intra-molecular transmission of information between distal sites. Microscopically, allosteric response is closely related to correlated atomic fluctuations. Conventional correlation analysis correlates the atomic fluctuations at two sites by taking the dot product (DP) between the fluctuations, which accounts only for the parallel and antiparallel components. Here, we present a singular value decomposition (SVD) method that analyzes the correlation coefficient of fluctuation dynamics with an arbitrary angle between the correlated directions. In a model allosteric system, the second PDZ domain (PDZ2) in the human PTP1E protein, approximately one third of the strong correlations have near-perpendicular directions, which are underestimated in the conventional method. The discrimination becomes more prominent for residue pairs with larger separation. The results of the proposed SVD method are more consistent with the experimentally determined PDZ2 dynamics than those of conventional method. In addition, the SVD method improved the prediction accuracy of the allosteric sites in a dataset of 23 known allosteric monomer proteins. The proposed method may inspire extended investigation not only into allostery, but also into protein dynamics and drug design.

The Arginine Pairs and C-Termini of the Sso7c4 from Sulfolobus solfataricus Participate in Binding and Bending DNA

The Sso7c4 from Sulfolobus solfataricus forms a dimer, which is believed to function as a chromosomal protein involved in genomic DNA compaction and gene regulation. Here, we present the crystal structure of wild-type Sso7c4 at a high resolution of 1.63 Å, showing that the two basic C-termini are disordered. Based on the fluorescence polarization (FP) binding assay, two arginine pairs, R11/R22' and R11'/R22, on the top surface participate in binding DNA. As shown in electron microscopy (EM) images, wild-type Sso7c4 compacts DNA through bridging and bending interactions, whereas the binding of C-terminally truncated proteins rigidifies and opens DNA molecules, and no compaction of the DNA occurs. Moreover, the FP, EM and fluorescence resonance energy transfer (FRET) data indicated that the two basic and flexible C-terminal arms of the Sso7c4 dimer play a crucial role in binding and bending DNA. Sso7c4 has been classified as a repressor-like protein because of its similarity to Escherichia coli Ecrep 6.8 and Ecrep 7.3 as well as Agrobacterium tumefaciens ACCR in amino acid sequence. Based on these data, we proposed a model of the Sso7c4-DNA complex using a curved DNA molecule in the catabolite activator protein-DNA complex. The DNA end-to-end distance measured with FRET upon wild-type Sso7c4 binding is almost equal to the distance measured in the model, which supports the fidelity of the proposed model. The FRET data also confirm the EM observation showing that the binding of wild-type Sso7c4 reduces the DNA length while the C-terminal truncation does not. A functional role for Sso7c4 in the organization of chromosomal DNA and/or the regulation of gene expression through bridging and bending interactions is suggested.

Dynamic Allostery of the Catabolite Activator Protein Revealed by Interatomic Forces

The Catabolite Activator Protein (CAP) is a showcase example for entropic allostery. For full activation and DNA binding, the homodimeric protein requires the binding of two cyclic AMP (cAMP) molecules in an anti-cooperative manner, the source of which appears to be largely of entropic nature according to previous experimental studies. We here study at atomic detail the allosteric regulation of CAP with Molecular dynamics (MD) simulations. We recover the experimentally observed entropic penalty for the second cAMP binding event with our recently developed force covariance entropy estimator and reveal allosteric communication pathways with Force Distribution Analyses (FDA). Our observations show that CAP binding results in characteristic changes in the interaction pathways connecting the two cAMP allosteric binding sites with each other, as well as with the DNA binding domains. We identified crucial relays in the mostly symmetric allosteric activation network, and suggest point mutants to test this mechanism. Our study suggests inter-residue forces, as opposed to coordinates, as a highly sensitive measure for structural adaptations that, even though minute, can very effectively propagate allosteric signals.

The Catabolite Activator Protein (CAP) is a well-studied example for how cellular catabolite levels are integrated into the gene regulation. Its affinity for a specific stretch of DNA can be switched on by the binding of two nucleotide molecules termed cAMP to its two protomers. Even though the nucleotides occupy structurally identical binding pockets, the second cAMP binding occurs at an affinity orders of magnitude lower than the first cAMP binding. The question arises how, in the absence of structural changes, the first binding can affect the second. An answer from experiments has been that the communication is largely of entropic nature, i.e. the second cAMP binding would lead to a pronounced reduction in atomic fluctuations of the protein without affecting the atomic mean positions. We here revisited this question by performing Molecular Dynamics simulations. By measuring correlations of forces, a newly derived method outperforming the more common coordinate-based approach, we could recover the previously determined entropic penalty. In addition, however, we observed unobtrusive structural changes of side-chain interactions leading to the occlusion of the second binding pocket that add a critical ‘enthalpic’ component hitherto overlooked. Our study provides a mechanistic view onto the intriguing anti-cooperativity of CAP.

A novel indirect sequence readout component in the E. coli cyclic AMP receptor protein operator

The cyclic AMP receptor protein (CRP) from Escherichia coli has been extensively studied for several decades. In particular, a detailed characterization of CRP interaction with DNA has been obtained. The CRP dimer recognizes a consensus sequence AANTGTGANNNNNNTCACANTT through direct amino acid nucleobase interactions in the major groove of the two operator half-sites. Crystal structure analyses have revealed that the interaction results in two strong kinks at the TG/CA steps closest to the 6-base-pair spacer (N6). This spacer exhibits high sequence variability among the more than 100 natural binding sites in the E. coli genome, but the exact role of the N6 region in CRP interaction has not previously been systematic examined. Here we employ an in vitro selection system based on a randomized N6 spacer region to demonstrate that CRP binding to the lacP1 site may be enhanced up to 14-fold or abolished by varying the N6 spacer sequences. Furthermore, on the basis of sequence analysis and uranyl (UO2(2+)) probing data, we propose that the underlying mechanism relies on N6 deformability.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Comprehensive insights into Mycobacterium tuberculosis DevR (DosR) regulon activation switch

DevR regulon function is believed to be crucial for the survival of Mycobacterium tuberculosis during dormancy. In this study, we undertook a comprehensive analysis of the DevR regulon. All the regulon promoters were assigned to four classes based on the number of DevR binding sites (Dev boxes). A minimum of two boxes are essential for complete interaction and their tandem arrangement is an architectural hallmark at all promoters. Initial interaction of DevR with the conserved box is essential for its cooperative binding to adjacent sites bearing low to very poor sequence conservation and is the universal mechanism underlying DevR-mediated transcriptional induction. The functional importance of tandem arrangement was established by analyzing promoter variants harboring Dev boxes with altered spacing. Conserved sequence logos were generated from 47 binding sequences which included 24 newly discovered Dev boxes. In each half site of an 18-bp binding motif, G₅ and C₇ are essential for DevR binding. Finally, we show that DevR regulon induction occurs in a temporal manner and genes that are induced early are also usually powerfully induced. The information theory-based approach along with binding and temporal expression studies provide us with comprehensive insights into the complex pattern of DevR regulon activation.

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression

BACKGROUND

A nucleosome is the fundamental repeating unit of the eukaryotic chromosome. It has been shown that the positioning of a majority of nucleosomes is primarily controlled by factors other than the intrinsic preference of the DNA sequence. One of the key questions in this context is the role, if any, that can be played by the variability of nucleosomal DNA structure.

RESULTS

In this study, we have addressed this question by analysing the variability at the dinucleotide and trinucleotide as well as longer length scales in a dataset of nucleosome X-ray crystal structures. We observe that the nucleosome structure displays remarkable local level structural versatility within the B-DNA family. The nucleosomal DNA also incorporates a large number of kinks.

CONCLUSIONS

Based on our results, we propose that the local and global level versatility of B-DNA structure may be a significant factor modulating the formation of nucleosomes in the vicinity of high-plasticity genes, and in varying the probability of binding by regulatory proteins. Hence, these factors should be incorporated in the prediction algorithms and there may not be a unique 'template' for predicting putative nucleosome sequences. In addition, the multimodal distribution of dinucleotide parameters for some steps and the presence of a large number of kinks in the nucleosomal DNA structure indicate that the linear elastic model, used by several algorithms to predict the energetic cost of nucleosome formation, may lead to incorrect results.

First-principles determination of molecular conformations of cyclic adenosine 3',5'-monophosphate in gas phase and aqueous solution

Extensive first-principles electronic structure calculations were performed in this study to explore the possible molecular structures and their concentration distribution of an intracellular second messenger, that is, cyclic adenosine 3',5'-monophosphate (cAMP), and its protonated form (cAMPH) in the gas phase and aqueous solution. The calculations resulted in prediction of four different stable conformers of cAMP and eight different stable conformers of cAMPH and their relative Gibbs free energies in the gas phase and aqueous solution. All of the computational results consistently demonstrate that the predominant conformers of cAMP and cAMPH are always the cAMP-chair-anti and cAMPH-chair2-syn conformers, respectively, in both the gas phase and aqueous solution. It has been demonstrated that the free energy barriers calculated for the intertransformation reactions between different conformers are very low (below approximately 6 kcal/mol) such that the intertransformation reactions between different conformers are very fast so that the concentration distribution of the system can quickly reach the thermodynamic equilibration during the process of binding with a protein. The calculated phenomenological pKa of 3.66 is in good agreement with the experimental pKa of 3.9 reported in literature, suggesting that the computational predictions resulted from this study are reasonable.

Altered oligomerization properties of N316 mutants of Escherichia coli TyrR

The transcriptional regulator TyrR is known to undergo a dimer-to-hexamer conformational change in response to aromatic amino acids, through which it controls gene expression. In this study, we identified N316D as the second-site suppressor of Escherichia coli TyrR(E274Q), a mutant protein deficient in hexamer formation. N316 variants exhibited altered in vivo regulatory properties, and the most drastic changes were observed for TyrR(N316D) and TyrR(N316R) mutants. Gel filtration analyses revealed that the ligand-mediated oligomer formation was enhanced and diminished for TyrR(N316D) and TyrR(N316R), respectively, compared with the wild-type TyrR. ADP was substituted for ATP in the oligomer formation of TyrR(N316D).

CRP binding and transcription activation at CRP-S sites

In Haemophilus influenzae, as in Escherichia coli, the cAMP receptor protein (CRP) activates transcription from hundreds of promoters by binding symmetrical DNA sites with the consensus half-site 5'-A(1)A(2)A(3)T(4)G(5)T(6)G(7)A(8)T(9)C(10)T(11). We have previously identified 13 H. influenzae CRP sites that differ from canonical (CRP-N) sites in the following features: (1) Both half-sites of these noncanonical (CRP-S) sites have C(6) instead of T(6), although they otherwise have an unusually high level of identity with the binding site consensus. (2) Only promoters with CRP-S sites require both the CRP and Sxy proteins for transcription activation. To study the functional significance of CRP-S site sequences, we purified H. influenzae (Hi)CRP and compared its DNA binding properties to those of the well-characterized E. coli (Ec)CRP. All EcCRP residues that contact DNA are conserved in HiCRP, and both proteins demonstrated a similar high affinity for the CRP-N consensus sequence. However, whereas EcCRP bound specifically to CRP-S sites in vitro, HiCRP did not. By systematically substituting base pairs in native promoters and in the CRP-N consensus sequence, we confirmed that HiCRP is highly specific for the perfect core sequence T(4)G(5)T(6)G(7)A(8) and is more selective than EcCRP at other positions in CRP sites. Even though converting C(6)-->T(6) greatly enhanced HiCRP binding to a CRP-S site, this had the unexpected effect of nearly abolishing promoter activity. A+T-rich sequences upstream of CRP-S sites were also found to be required for promoter activation, raising the possibility that Sxy binds these A+T sequences to simultaneously enable CRP-DNA binding and assist in RNA polymerase recruitment.

Syn, anti, and finally both conformations of cyclic AMP are involved in the CRP-dependent transcription initiation mechanism in E. coli lac operon

The cyclic AMP receptor protein (CRP) of Escherichia coli regulates the activity of more than 150 genes. Allosteric changes in CRP structure accompanied by cAMP binding, initiate transcription through protein binding to specific DNA sequences. Initially, researchers proposed a two-site cAMP-binding model for CRP-dependent transcription activation since biophysical methods showed two transitions during titration experiments. Three conformational states were considered; apo-CRP, CRP:(cAMP)(1) and CRP:(cAMP)(2), and CRP:(cAMP)(1) was proposed as the active form in this initial model. X-ray data indicated an anti conformation and in contrast NMR experiments suggested a syn conformation for bound cAMPs. For years this paradigm about ligand conformation has been ambiguous. When CRP was crystallized with four bound cAMP in the last decade, two cAMPs were assigned to syn and the other two to anti conformations. Again three conformational states were suggested; apo-CRP, CRP:(cAMP)(2), and CRP:(cAMP)(4). This new structure changed the view of CRP allosteric activation from a two-site model to a four-site model in the literature and the new model has been supported by biochemical and genetic data so far. According to the accepted model, binding of the first two cAMP molecules displays positive cooperativity, however, binding of the last two cAMP molecules shows negative cooperativity. This resolved the conflict between dynamic and static experimental observations. However, this new model cannot explain the initiation mechanism as previously proposed because functionally active CRP has only one cAMP equivalent. Gene regulation and transcription factors are involved in regulating both prokaryotic and eukaryotic metabolism. Although gene regulation and expression are much more complex in eukaryotes, CRP-mediated transcription initiation is a model of general interest to life sciences and medicine. Therefore, the aim of this review is to summarize recent works and developments on the cAMP-dependent CRP activation mechanism in E. coli.

Dissecting direct and indirect readout of cAMP receptor protein DNA binding using an inosine and 2,6-diaminopurine in vitro selection system

The DNA interaction of the Escherichia coli cyclic AMP receptor protein (CRP) represents a typical example of a dual recognition mechanism exhibiting both direct and indirect readout. We have dissected the direct and indirect components of DNA recognition by CRP employing in vitro selection of a random library of DNA-binding sites containing inosine (I) and 2,6-diaminopurine (D) instead of guanine and adenine, respectively. Accordingly, the DNA helix minor groove is structurally altered due to the ‘transfer’ of the 2-amino group of guanine (now I) to adenine (now D), whereas the major groove is functionally intact. The majority of the selected sites contain the natural consensus sequence TGTGAN₆TCACA (i.e. TITIDN₆TCDCD). Thus, direct readout of the consensus sequence is independent of minor groove conformation. Consequently, the indirect readout known to occur in the TG/CA base pair step (primary kink site) in the consensus sequence is not affected by I–D substitutions. In contrast, the flanking regions are selected as I/C rich sequences (mostly I-tracts) instead of A/T rich sequences which are known to strongly increase CRP binding, thereby demonstrating almost exclusive indirect readout of helix structure/flexibility in this region through (anisotropic) flexibility of I-tracts.

Using the tools and resources of the RCSB protein data bank

The Protein Data Bank (PDB; http://www.pdb.org) is the world-wide repository for three-dimensional structural data determined using various experimental methods. The options and procedures for searching and downloading structural data from the Research Collaboratory for Structural Bioinformatics (RCSB) PDB are described here, along with tools for assessing the quality of structures. Several types of information are associated with each structure deposition, including atomic coordinates of the structure, experimental data used to solve it, sequences of all macromolecules in the structures, details about the structure solution method, images showing different views of the structure, derived geometric data, and a variety of links to other resources. These data and resources may be used for understanding the function and stability of the molecule and for designing biochemical, genetic, or other experiments. They can also be used for molecular modeling and drug design.

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.

3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures

We present a set of protocols showing how to use the 3DNA suite of programs to analyze, rebuild and visualize three-dimensional nucleic-acid structures. The software determines a wide range of conformational parameters, including the identities and rigid-body parameters of interacting bases and base-pair steps, the nucleotides comprising helical fragments, the area of overlap of stacked bases and so on. The reconstruction of three-dimensional structure takes advantage of rigorously defined rigid-body parameters, producing rectangular block representations of the nucleic-acid bases and base pairs and all-atom models with approximate sugar-phosphate backbones. The visualization components create vector-based drawings and scenes that can be rendered as raster-graphics images, allowing for easy generation of publication-quality figures. The utility programs use geometric variables to control the view and scale of an object, for comparison of related structures. The commands run in seconds even for large structures. The software and related information are available at http://3dna.rutgers.edu/.

Visualization analysis of inter-fragment interaction energies of CRP-cAMP-DNA complex based on the fragment molecular orbital method

A visualization method for inter-fragment interaction energies (IFIEs) of biopolymers is presented on the basis of the fragment molecular orbital (FMO) method. The IFIEs appropriately illustrate the information about the interaction energies between the fragments consisting of amino acids, nucleotides and other molecules. The IFIEs are usually analyzed in a matrix form called an IFIE matrix. Analyzing the IFIE matrix, we detect important fragments for the function of biomolecular systems and quantify the strength of interaction energies based on quantum chemistry, including the effects of charge transfer, electronic polarization and dispersion force. In this study, by analyzing a protein-DNA complex, we report a visual representation of the IFIE matrix, a so-called IFIE map. We comprehensively examine what information the IFIE map contains concerning structures and stabilities of the protein-DNA complex.

Control of the respiratory metabolism of Thermus thermophilus by the nitrate respiration conjugative element NCE

The strains of Thermus thermophilus that contain the nitrate respiration conjugative element (NCE) replace their aerobic respiratory chain by an anaerobic counterpart made of the Nrc-NADH dehydrogenase and the Nar-nitrate reductase in response to nitrate and oxygen depletion. This replacement depends on DnrS and DnrT, two homologues to sensory transcription factors encoded in a bicistronic operon by the NCE. DnrS is an oxygen-sensitive protein required in vivo to activate transcription on its own dnr promoter and on that of the nar operon, but not required for the expression of the nrc operon. In contrast, DnrT is required for the transcription of these three operons and also for the repression of nqo, the operon that encodes the major respiratory NADH dehydrogenase expressed during aerobic growth. Thermophilic in vitro assays revealed that low DnrT concentrations allows the recruitment of the T. thermophilus RNA polymerase sigma(A) holoenzyme to the nrc promoter and its transcription, whereas higher DnrT concentrations are required to repress transcription on the nqo promoter. In conclusion, our data show a complex autoinducible mechanism by which DnrT functions as the transcriptional switch that allows the NCE to take the control of the respiratory metabolism of its host during adaptation to anaerobic growth.

cAMP activation of PKA defines an ancient signaling mechanism

cAMP and the cAMP binding domain (CBD) constitute a ubiquitous regulatory switch that translates an extracellular signal into a biological response. The CBD contains alpha- and beta-subdomains with cAMP binding to a phosphate binding cassette (PBC) in the beta-sandwich. The major receptors for cAMP in mammalian cells are the regulatory subunits (R-subunits) of PKA where cAMP and the catalytic subunit compete for the same CBD. The R-subunits inhibit kinase activity, whereas cAMP releases that inhibition. Here, we use NMR to map at residue resolution the cAMP-dependent interaction network of the CBD-A domain of isoform Ialpha of the R-subunit of PKA. Based on H/D, H/H, and N(z) exchange data, we propose a molecular model for the allosteric regulation of PKA by cAMP. According to our model, cAMP binding causes long-range perturbations that propagate well beyond the immediate surroundings of the PBC and involve two key relay sites located at the C terminus of beta(2) (I163) and N terminus of beta(3) (D170). The I163 site functions as one of the key triggers of global unfolding, whereas the D170 locus acts as an electrostatic switch that mediates the communication between the PBC and the B-helix. Removal of cAMP not only disrupts the cap for the B' helix within the PBC, but also breaks the circuitry of cooperative interactions stemming from the PBC, thereby uncoupling the alpha- and beta-subdomains. The proposed model defines a signaling mechanism, conserved in every genome, where allosteric binding of a small ligand disrupts a large protein-protein interface.

A model for agonism and antagonism in an ancient and ubiquitous cAMP-binding domain

The cAMP-binding domain (CBD) is an ancient and conserved regulatory motif that allosterically modulates the function of a group of diverse proteins, thereby translating the cAMP signal into a controlled biological response. The main receptor for cAMP in mammals is the ubiquitous regulatory (R) subunit of protein kinase A. Despite the recognized significant potential for pharmacological applications of CBDs, currently only one group of competitive inhibitor antagonists is known: the (R(p))-cAMPS family of phosphorothioate cAMP analogs, in which the equatorial exocyclic oxygen of cAMP is replaced by sulfur. It is also known that the diastereoisomer (S(p))-cAMPS with opposite phosphorous chirality is a cAMP agonist, but the molecular mechanism of action of these analogs is currently not fully understood. Previous crystallographic and unfolding investigations point to the enhanced CBD dynamics as a key determinant of antagonism. Here, we investigate the (R(p))- and (S(p))-cAMPS-bound states of R(CBD-A) using a comparative NMR approach that reveals a clear chemical shift and dynamic NMR signature, differentiating the (S(p))-cAMPS agonist from the (R(p))-cAMPS antagonist. Based on these data, we have proposed a model for the (R(p)/S(p))-cAMPS antagonism and agonism in terms of steric and electronic effects on two main allosteric relay sites, Ile(163) and Asp(170), respectively, affecting the stability of a ternary inhibitory complex formed by the effector ligand, the regulatory and the catalytic subunits of protein kinase A. The proposed model not only rationalizes the existing data on the phosphorothioate analogs, but it will also facilitate the design of novel cAMP antagonists and agonists.

Non-canonical CRP sites control competence regulons in Escherichia coli and many other gamma-proteobacteria

Escherichia coli's cAMP receptor protein (CRP), the archetypal bacterial transcription factor, regulates over a hundred promoters by binding 22 bp symmetrical sites with the consensus core half-site TGTGA. However, Haemophilus influenzae has two types of CRP sites, one like E.coli's and one with the core sequence TGCGA that regulates genes required for DNA uptake (natural competence). Only the latter 'CRP-S' sites require both CRP and the coregulator Sxy for activation. To our knowledge, the TGTGA and TGCGA motifs are the first example of one transcription factor having two distinct binding-site motifs. Here we show that CRP-S promoters are widespread in the gamma-proteobacteria and demonstrate their Sxy-dependence in E.coli. Orthologs of most H.influenzae CRP-S-regulated genes are ubiquitous in the five best-studied gamma-proteobacteria families, Enterobacteriaceae, Pasteurellaceae, Pseudomonadaceae, Vibrionaceae and Xanthomonadaceae. Phylogenetic footprinting identified CRP-S sites in the promoter regions of the Enterobacteriaceae, Pasteurellaceae and Vibrionaceae orthologs, and canonical CRP sites in orthologs of genes known to be Sxy-independent in H.influenzae. Bandshift experiments confirmed that E.coli CRP-S sequences are low affinity binding sites for CRP, and mRNA analysis showed that they require CRP, cAMP (CRP's allosteric effector) and Sxy for gene induction. This work suggests not only that the gamma-proteobacteria share a common DNA uptake mechanism, but also that, in the three best studied families, their competence regulons share both CRP-S specificity and Sxy dependence.

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: recognition of pyrimidine-purine and purine-purine steps

The catabolite activator protein (CAP) bends DNA in the CAP-DNA complex, typically introducing a sharp DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A1A2A3T4G5T6G7A8T9C10T11 -3' ("primary kink"). In previous work, we showed that CAP recognizes the nucleotide immediately 5' to the primary-kink site, T6, through an "indirect-readout" mechanism involving sequence effects on energetics of primary-kink formation. Here, to understand further this example of indirect readout, we have determined crystal structures of CAP-DNA complexes containing each possible nucleotide at position 6. The structures show that CAP can introduce a DNA kink at the primary-kink site with any nucleotide at position 6. The DNA kink is sharp with the consensus pyrimidine-purine step T6G7 and the non-consensus pyrimidine-purine step C6G7 (roll angles of approximately 42 degrees, twist angles of approximately 16 degrees ), but is much less sharp with the non-consensus purine-purine steps A6G7 and G6G7 (roll angles of approximately 20 degrees, twist angles of approximately 17 degrees). We infer that CAP discriminates between consensus and non-consensus pyrimidine-purine steps at positions 6-7 solely based on differences in the energetics of DNA deformation, but that CAP discriminates between the consensus pyrimidine-purine step and non-consensus purine-purine steps at positions 6-7 both based on differences in the energetics of DNA deformation and based on qualitative differences in DNA deformation. The structures further show that CAP can achieve a similar, approximately 46 degrees per DNA half-site, overall DNA bend through a sharp DNA kink, a less sharp DNA kink, or a smooth DNA bend. Analysis of these and other crystal structures of CAP-DNA complexes indicates that there is a large, approximately 28 degrees per DNA half-site, out-of-plane component of CAP-induced DNA bending in structures not constrained by end-to-end DNA lattice interactions and that lattice contacts involving CAP tend to involve residues in or near biologically functional surfaces.

Intra- and intermolecular interactions between cyclic-AMP receptor protein and DNA:Ab initio fragment molecular orbital study

The ab initio fragment molecular orbital (FMO) calculations were performed for the cAMP receptor protein (CRP) complexed with a cAMP and DNA duplex to elucidate their sequence-specific binding and the stability of the DNA duplex, as determined by analysis of their inter- and intramolecular interactions. Calculations were performed with the AMBER94 force field and at the HF and MP2 levels with several basis sets. The interfragment interaction energies (IFIEs) were analyzed for interactions of CRP-cAMP with each base pair, DNA duplex with each amino acid residue, and each base pair with each residue. In addition, base-base interactions were analyzed including hydrogen bonding and stacking of DNA. In the interaction between DNA and CRP-cAMP, there was a significant charge transfer (CT) from the DNA to CRP, and this CT interaction played an important role as well as the electrostatic interactions. It is necessary to apply a quantum mechanical approach beyond the "classical" force-field approach to describe the sequence specificity. In the DNA intramolecular interaction, the dispersion interactions dominated the stabilization of the base-pair stacking interactions. Strong, attractive 1,2-stacking interactions and weak, repulsive 1,3-stacking interactions were observed. Comparison of the intramolecular interactions of free and complexed DNA revealed that the base-pairing interactions were stronger, and the stacking interactions were weaker, in the complexed structure. Therefore, the DNA duplex stability appears to change due to both the electrostatic and the CT interactions that take place under conditions of DNA-CRP binding.

Correlated evolutionary pressure at interacting transcription factors and DNA response elements can guide the rational engineering of DNA binding specificity

Understanding the molecular mechanisms of the specific interaction between transcription factor proteins and DNA is key to comprehend the regulation of gene expression and to develop technologies to engineer transcription factors. Thus far, although there have been several attempts to elucidate protein-DNA interaction through amino acid-base recognition codes, sequence based profiles, or physical models of interaction, the greatest successes in engineering DNA binding specificity remain experimental. Here we present the first systematic evidence of correlated evolutionary pressure at interacting amino acid residues and DNA base-pairs in transcription factors, and show that it can be used to rationally engineer DNA binding specificity. The correlation is between the relative evolutionary importance of protein residues and DNA bases, measured, respectively, in terms of the Evolutionary Trace (ET) rank and information entropy. The evolutionarily most important residues interact with the most conserved base-pairs within the response element while residues of least importance interact with the most variable base-pairs. The correlation averages 0.74 over 12 unrelated families of transcriptional regulators, including nuclear hormone receptors, basic helix-loop-helix, ETS- and homeo-domain family. To test the predictive power of this correlation, we targeted a mutational swap of top-ranked ET residues in a transcription factor, LRH-1. This redirects LRH-1 binding as predicted and showed that, in this case, evolutionary importance and binding specificity are coupled sufficiently strongly for the Evolutionary Trace to guide the computational design of DNA binding specificity. This establishes the existence of evolutionary importance correlation at protein-DNA interfaces, and demonstrates that it is a useful principle for the rational engineering of binding specificity.

Fluorescence quenching and kinetic studies of conformational changes induced by DNA and cAMP binding to cAMP receptor protein from Escherichia coli

Cyclic AMP receptor protein (CRP) regulates the expression of more then 100 genes in Escherichia coli. It is known that the allosteric activation of CRP by cAMP involves a long-distance signal transmission from the N-terminal cAMP-binding domain to the C-terminal domain of CRP responsible for the interactions with specific sequences of DNA. In this report we have used a CRP mutant containing a single Trp13 located in the N-terminal domain of the protein. We applied the iodide and acrylamide fluorescence quenching method in order to study how different DNA sequences and cAMP binding induce the conformational changes in the CRP molecule. The results presented provide evidence for the occurrence of a long-distance conformational signal transduction within the protein from the C-terminal DNA-binding domain to the N-terminal domain of CRP. This conformational signal transmission depends on the promoter sequence. We also used the stopped-flow and Forster resonance energy transfer between labeled Cys178 of CRP and fluorescently labeled DNA sequences to study the kinetics of DNA-CRP interactions. The results thus obtained lead to the conclusion that CRP can exist in several conformational states and that their distribution is affected by binding of both the cAMP and of specific DNA sequences.

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

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.

Two cAMP receptor proteins with different biochemical properties in the filamentous cyanobacterium Anabaena sp. PCC 7120

Two open reading frames (ORFs), alr0295 and alr2325, are found to encode putative cAMP receptor proteins (CRPs) in the genome of the filamentous cyanobacterium Anabaena sp. PCC 7120. These ORFs were named cAMP receptor protein-like gene A in Anabaena sp. PCC 7120 (ancrpA) and cAMP receptor protein-like gene B in Anabaena sp. PCC 7120 (ancrpB), respectively, and those translated products were investigated. The equilibrium dialysis measurements revealed that AnCrpA bound with cAMP specifically, while AnCrpB bound with both cAMP and cGMP, though the affinity for cGMP was weak. The binding affinity for cAMP of AnCrpA showed the lowest dissociation constant, approximately 0.8 microM, among bacterial CRPs. A gel mobility shift assay elucidated that AnCrpA and AnCrpB formed a complex with the consensus DNA sequence in the presence of cAMP, although AnCrpB did not have ordinary DNA-binding motifs.

Downstream DNA selectively affects a paused conformation of human RNA polymerase II

Transcriptional pausing by human RNA polymerase II (RNAPII) in the HIV-1 LTR is caused principally by a weak RNA:DNA hybrid that allows rearrangement of reactive or catalytic groups in the enzyme's active site. This rearrangement creates a transiently paused state called the unactivated intermediate that can backtrack into a more long-lived paused species. We report that three different regions of the not-yet-transcribed DNA just downstream of the pause site affect the duration of the HIV-1 pause, and also can influence pause formation. Downstream DNA in at least one region, a T-tract from +5 to +8, increases pause duration by specifically affecting the unactivated intermediate, without corresponding effects on the active or backtracked states. We suggest this effect depends on RNAPII-modulated DNA plasticity and speculate it is mediated by the "trigger loop" thought to participate in RNAP's catalytic cycle. These findings provide a new framework for understanding downstream DNA effects on RNAP.

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

From sequence to structure and back again: approaches for predicting protein-DNA binding

Gene regulation in higher organisms is achieved by a complex network of transcription factors (TFs). Modulating gene expression and exploring gene function are major aims in molecular biology. Furthermore, the identification of putative target genes for a certain TF serve as powerful tools for specific targeting of rational drugs.

Detecting the short and variable transcription factor binding sites (TFBSs) in genomic DNA is an intriguing challenge for computational and structural biologists. Fast and reliable computational methods for predicting TFBSs on a whole-genome scale offer several advantages compared to the current experimental methods that are rather laborious and slow. Two main approaches are being explored, advanced sequence-based algorithms and structure-based methods.

The aim of this review is to outline the computational and experimental methods currently being applied in the field of protein-DNA interactions. With a focus on the former, the current state of the art in modeling these interactions is discussed. Surveying sequence and structure-based methods for predicting TFBSs, we conclude that in order to achieve a sound and specific method applicable on genomic sequences it is desirable and important to bring these two approaches together.

Sequence specific DNA binding of Ets-1 transcription factor: molecular dynamics study on the Ets domain--DNA complexes

Molecular dynamics (MD) simulations for Ets-1 ETS domain-DNA complexes were performed to investigate the mechanism of sequence-specific recognition of the GGAA DNA core by the ETS domain. Employing the crystal structure of the Ets-1 ETS domain-DNA complex as a starting structure we carried out MD simulations of: (i). the complex between Ets-1 ETS domain and a 14 base-pair DNA containing GGAA core sequence (ETS-GGAA); (ii). the complex between the ETS domain and a DNA having single base-pair mutation, GGAG sequence (ETS-GGAG); and (iii). the 14 base-pair DNA alone (GGAA). Comparative analyses of the MD structures of ETS-GGAA and ETS-GGAG reveal that the DNA bending angles and the ETS domain-DNA phosphate interactions are similar in these complexes. These results support that the GGAA core sequence is distinguished from the mutated GGAG sequence by a direct readout mechanism in the Ets-1 ETS domain-DNA complex. Further analyses of the direct contacts in the interface between the helix-3 region of Ets-1 and the major groove of the core DNA sequence clearly show that the highly conserved arginine residues, Arg391 and Arg394, play a critical role in binding to the GGAA core sequence. These arginine residues make bidentate contacts with the nucleobases of GG dinucleotides in GGAA core sequence. In ETS-GGAA, the hydroxyl group of Tyr395 is hydrogen bonded to N7 nitrogen of A(3) (the third adenosine in the GGAA core), while the hydroxyl group makes a contact with N4 nitrogen of C(4') (the complementary nucleotide of the fourth guanosine G(4) in the GGAG sequence) in the ETS-GGAG complex. We have found that this difference in behavior of Tyr395 results in the relatively large motion of helix-3 in the ETS-GGAG complex, causing the collapse of bidentate contacts between Arg391/Arg394 and the GG dinucleotides in the GGAG sequence.

Base coupling in sequence-specific site recognition by the ETS domain of murine PU.1

The ETS domain of murine PU.1 tolerates a large number of DNA cognates bearing a central consensus 5'-GGAA-3' that is flanked by a diverse combination of bases on both sides. Previous attempts to define the sequence selectivity of this DNA binding domain by combinatorial methods have not successfully predicted observed patterns among in vivo promoter sequences in the genome, and have led to the hypothesis that energetic coupling occurs among the bases in the flanking sequences. To test this hypothesis, we determined, using thermodynamic cycles, the complex stabilities and base coupling energies of the PU.1 ETS domain for a set of 26 cognate variants (based on the lambdaB site of the Ig(lambda)2-4 enhancer, 5'-AATAAAAGGAAGTGAAACCAA-3') in which flanking sequences up to three bases upstream and/or two bases downstream of the core consensus are substituted. We observed that both cooperative and anticooperative coupling occurs commonly among the flanking sequences at all the positions investigated. This phenomenon extends at least three bases in the 5' side and is, at least on our experimental data, due exclusively to pairwise interactions between the flanking bases, and not changes in the local environment of the DNA groove floor. Energetic coupling also occurs between the flanking sides across the core consensus, suggesting long-range conformational effects along the DNA target and/or in the protein. Our data provide an energetic explanation for the pattern of flanking bases observed among in vivo promoter sequences and reconcile the apparent discrepancies raised by the combinatorial experiments. We also discuss the significance of base coupling in light of an indirect readout mechanism in ETS/DNA site recognition.

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.

Hidden Markov models from molecular dynamics simulations on DNA

An enhanced bioinformatics tool incorporating the participation of molecular structure as well as sequence in protein DNA recognition is proposed and tested. Boltzmann probability models of sequence-dependent DNA structure from all-atom molecular dynamics simulations were obtained and incorporated into hidden Markov models (HMMs) that can recognize molecular structural signals as well as sequence in protein-DNA binding sites on a genome. The binding of catabolite activator protein (CAP) to cognate DNA sequences was used as a prototype case for implementation and testing of the method. The results indicate that even HMMs based on probabilistic roll/tilt dinucleotide models of sequence-dependent DNA structure have some capability to discriminate between known CAP binding and nonbinding sites and to predict putative CAP binding sites in unknowns. Restricting HMMs to sequence only in regions of strong consensus in which the protein makes base specific contacts with the cognate DNA further improved the discriminatory capabilities of the HMMs. Comparison of results with controls based on sequence only indicates that extending the definition of consensus from sequence to structure improves the transferability of the HMMs, and provides further supportive evidence of a role for dynamical molecular structure as well as sequence in genomic regulatory mechanisms.

Solution structure of dAATAA and dAAUAA DNA bulges

The NMR structure analysis is described for two DNA molecules of identical stem sequences with a five base loop containing a pyrimidine, thymin or uracil, in between purines. These five unpaired nucleotides are bulged out and are known to induce a kink in the duplex structure. The dAATAA bulge DNA is kinked between the third and the fourth nucleotide. This contrasts with the previously studied dAAAAA bulge DNA where we found a kink between the fourth and fifth nucleotide. The total kinking angle is approximately 104 degrees for the dAATAA bulge. The findings were supported by electrophoretic data and fluorescence resonance energy transfer measurements of a similar DNA molecule end-labeled by suitable fluorescent dyes. For the dAAUAA bulge the NMR data result in a similar structure as reported for the dAATAA bulge with a kinking angle of approximately 87 degrees. The results are discussed in comparison with a rAAUAA RNA bulge found in a group I intron. Generally, the sequence-dependent structure of bulges is important to understand the role of DNA bulges in protein recognition.

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 degrees in the direction enhanced in the complex with protein. In the absence of protein, neutralization of five appropriate phosphates increases DNA curvature to 32 degrees (approximately 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.

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: alteration of DNA binding specificity through alteration of DNA kinking

The catabolite activator protein (CAP) sharply bends DNA in the CAP-DNA complex, introducing a DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A(1)A(2)A(3)T(4)G(5)T(6)G(7)A(8)T(9)C(10)T(11)-3' ("primary kink"). CAP recognizes the base-pair immediately 5' to the primary-kink site, T:A(6), through an "indirect-readout" mechanism involving sequence effects on the energetics of primary-kink formation. CAP recognizes the base-pair immediately 3' to the primary-kink site, G:C(7), through a "direct-readout" mechanism involving formation of a hydrogen bond between Glu181 of CAP and G:C(7). Here, we report that substitution of the carboxylate side-chain of Glu181 of CAP by the one-methylene-group-shorter carboxylate side-chain of Asp changes DNA binding specificity at position 6 of the DNA half site, changing specificity for T:A(6) to specificity for C:G(6), and we report a crystallographic analysis defining the structural basis of the change in specificity. The Glu181-->Asp substitution eliminates the primary kink and thus eliminates indirect-readout-based specificity for T:A(6). The Glu181-->Asp substitution does not eliminate hydrogen-bond formation with G:C(7), and thus does not eliminate direct-readout-based specificity for G:C(7).