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

Probing conformational dynamics of DNA binding by CO-sensing transcription factor, CooA

The transcription factor CooA is a CRP/FNR (cAMP receptor protein/ fumarate and nitrate reductase) superfamily protein that uses heme to sense carbon monoxide (CO). Allosteric activation of CooA in response to CO binding is currently described as a series of discrete structural changes, without much consideration for the potential role of protein dynamics in the process of DNA binding. This work uses site-directed spin-label electron paramagnetic resonance spectroscopy (SDSL-EPR) to probe slow timescale (μs-ms) conformational dynamics of CooA with a redox-stable nitroxide spin label, and IR spectroscopy to probe the environment at the CO-bound heme. A series of cysteine substitution variants were created to selectively label CooA in key functional regions, the heme-binding domain, the 4/5-loop, the hinge region, and the DNA binding domain. The EPR spectra of labeled CooA variants are compared across three functional states: Fe(III) "locked off", Fe(II)-CO "on", and Fe(II)-CO bound to DNA. We observe changes in the multicomponent EPR spectra at each location; most notably in the hinge region and DNA binding domain, broadening the description of the CooA allosteric mechanism to include the role of protein dynamics in DNA binding. DNA-dependent changes in IR vibrational frequency and band broadening further suggest that there is conformational heterogeneity in the active WT protein and that DNA binding alters the environment of the heme-bound CO.

Neuromodulators as Interdomain Signaling Molecules Capable of Occupying Effector Binding Sites in Bacterial Transcription Factors

Hormones and neurotransmitters are important components of inter-kingdom signaling systems that ensure the coexistence of eukaryotes with their microbial community. Their ability to affect bacterial physiology, metabolism, and gene expression was evidenced by various experimental approaches, but direct penetration into bacteria has only recently been reported. This opened the possibility of considering neuromodulators as potential effectors of bacterial ligand-dependent regulatory proteins. Here, we assessed the validity of this assumption for the neurotransmitters epinephrine, dopamine, and norepinephrine and two hormones (melatonin and serotonin). Using flexible molecular docking for transcription factors with ligand-dependent activity, we assessed the ability of neuromodulators to occupy their effector binding sites. For many transcription factors, including the global regulator of carbohydrate metabolism, CRP, and the key regulator of lactose assimilation, LacI, this ability was predicted based on the analysis of several 3D models. By occupying the ligand binding site, neuromodulators can sterically hinder the interaction of the target proteins with the natural effectors or even replace them. The data obtained suggest that the direct modulation of the activity of at least some bacterial transcriptional factors by neuromodulators is possible. Therefore, the natural hormonal background may be a factor that preadapts bacteria to the habitat through direct perception of host signaling molecules.

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.

Structural basis for the inhibitory mechanism of auranofin and gold(I) analogues against Pseudomonas aeruginosa global virulence factor regulator Vfr

Pseudomonas aeruginosa is a leading cause of hospital-acquired infections. Treatment of P. aeruginosa infections is difficult given its multiple virulence mechanisms, intrinsic antibiotic resistance mechanisms, and biofilm-forming ability. Auranofin, an approved oral gold compound for rheumatoid arthritis treatment, was recently reported to inhibit the growth of multiple bacterial species. Here, we identify P. aeruginosa's global virulence factor regulator Vfr as one target of auranofin. We report the mechanistic insights into the inhibitory mechanism of auranofin and gold(I) analogues to Vfr through structural, biophysical, and phenotypic inhibition studies. This work suggests that auranofin and gold(I) analogues have potential to be developed as anti-virulence drugs against P. aeruginosa.

General Designs Reveal a Purine-Pyrimidine Structural Code in Human DNA

The human genome carries a vast amount of information within its DNA sequences. The chemical bases A, T, C, and G are the basic units of information content, that are arranged into patterns and codes. Expansive areas of the genome contain codes that are not yet well understood. To decipher these, mathematical and computational tools are applied here to study genomic signatures or general designs of sequences. A novel binary components analysis is devised and utilized. This seeks to isolate the physical and chemical properties of DNA bases, which reveals sequence design and function. Here, information theory tools break down the information content within DNA bases, in order to study them in isolation for their genomic signatures and non-random properties. In this way, the RY (purine/pyrimidine), WS (weak/strong), and KM (keto/amino) general designs are observed in the sequences. The results show that RY, KM, and WS components have a similar and stable overall profile across all human chromosomes. It reveals that the RY property of a sequence is most distant from randomness in the human genome with respect to the genomic signatures. This is true across all human chromosomes. It is concluded that there exists a widespread potential RY code, and furthermore, that this is likely a structural code. Ascertaining this feature of general design, and potential RY structural code has far-reaching implications. This is because it aids in the understanding of cell biology, growth, and development, as well as downstream in the study of human disease and potential drug design.

A conformational study of the 10-23 DNAzyme via programmed DNA self-assembly

This communication measures the inter-helical angle of the 10-23 DNAzyme-substrate complex by atomic force microscopy (AFM). Herein, we have devised a strategy to assemble the DNAzyme-substrate complex into a periodic DNA 2D array, which allows reliable study of the conformation of the 10-23 DNAzyme by AFM imaging and fast Fourier transform (FFT). Specifically, the angle between the two flanking helical domains of the catalytic core has been determined via the repeating distance of the 2D array. We expect that the same strategy can generally be applicable for studying other nucleic acid structures.

Dual-function Spot 42 RNA encodes a 15-amino acid protein that regulates the CRP transcription factor

SignificanceDual-function RNAs base pair with target messenger RNAs as small regulatory RNAs and encode small protein regulators. However, only a limited number of these dual-function regulators have been identified. In this study, we show that a well-characterized base-pairing small RNA surprisingly also encodes a 15-amino acid protein. The very small protein binds the cyclic adenosine monophosphate receptor protein transcription factor to block activation of some promoters, raising the question of how many other transcription factors are modulated by unidentified small proteins.

Signal Transmission in Escherichia coli Cyclic AMP Receptor Protein for Survival in Extreme Acidic Conditions

During the life cycle of enteric bacterium Escherichia coli, it encounters a wide spectrum of pH changes. The asymmetric dimer of the cAMP receptor protein, CRP, plays a key role in regulating the expressions of genes and the survival of E. coli. To elucidate the pH effects on the mechanism of signal transmission, we present a combination of results derived from ITC, crystallography, and computation. CRP responds to a pH change by inducing a differential effect on the affinity for the binding events to the two cAMP molecules, ensuing in a reversible conversion between positive and negative cooperativity at high and low pH, respectively. The structures of four crystals at pH ranging from 7.8 to 6.5 show that CRP responds by inducing a differential effect on the structures of the two subunits, particularly in the DNA binding domain. Employing the COREX/BEST algorithm, computational analysis shows the change in the stability of residues at each pH. The change in residue stability alters the connectivity between residues including those in cAMP and DNA binding sites. Consequently, the differential impact on the topology of the connectivity surface among residues in adjacent subunits is the main reason for differential change in affinity; that is, the pH-induced differential change in residue stability is the biothermodynamic basis for the change in allosteric behavior. Furthermore, the structural asymmetry of this homodimer amplifies the differential impact of any perturbations. Hence, these results demonstrate that the combination of these approaches can provide insights into the underlying mechanism of an apparent complex allostery signal and transmission in CRP.

The RNA Polymerase α Subunit Recognizes the DNA Shape of the Upstream Promoter Element

We demonstrate here that the α subunit C-terminal domain of Escherichia coli RNA polymerase (αCTD) recognizes the upstream promoter (UP) DNA element via its characteristic minor groove shape and electrostatic potential. In two compositionally distinct crystallized assemblies, a pair of αCTD subunits bind in tandem to the UP element consensus A-tract that is 6 bp in length (A₆-tract), each with their arginine 265 guanidinium group inserted into the minor groove. The A₆-tract minor groove is significantly narrowed in these crystal structures, as well as in computationally predicted structures of free and bound DNA duplexes derived by Monte Carlo and molecular dynamics simulations, respectively. The negative electrostatic potential of free A₆-tract DNA is substantially enhanced compared to that of generic DNA. Shortening the A-tract by 1 bp is shown to "knock out" binding of the second αCTD through widening of the minor groove. Furthermore, in computationally derived structures with arginine 265 mutated to alanine in either αCTD, either with or without the "knockout" DNA mutation, contact with the DNA is perturbed, highlighting the importance of arginine 265 in achieving αCTD-DNA binding. These results demonstrate that the importance of the DNA shape in sequence-dependent recognition of DNA by RNA polymerase is comparable to that of certain transcription factors.

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.

Change in binding states between catabolite activating protein and DNA induced by ligand-binding: molecular dynamics and ab initio fragment molecular orbital calculations

The transcription mechanism of genetic information from DNA to RNA is efficiently controlled by regulatory proteins, such as catabolite activator protein (CAP), and their ligands. When cyclic AMP (cAMP) binds to CAP, the complex forms a dimer and binds specifically to DNA to activate the transcription mechanism. On the other hand, when cyclic GMP (cGMP) binds to CAP, the complex has no marked effect on the mechanism. In our previous study, based on molecular dynamics (MD) and ab initio fragment molecular orbital (FMO) methods, we elucidated which residues of CAP are important for the specific interactions between CAP and DNA in the CAP-monomer+DNA + cAMP complex. However, this monomer model for CAP cannot describe real interactions between the CAP-dimer and DNA because CAPs form a dimer before binding to DNA. Accordingly, here, we investigated stable structures and their electronic states for the CAP-dimer+DNA complex with cAMP or cGMP ligand, to clarify the influence of ligand-binding on the interactions between CAP-dimer and DNA. The MD simulations elucidated that the DNA-binding domains of CAP-dimer behave differently depending on the ligand bound to the CAP-dimer. In addition, FMO calculations revealed that the binding energy between CAP-dimer and DNA for the CAP-dimer+DNA + cAMP complex is larger than that for the CAP-dimer+DNA + cGMP complex, being consistent with experiments. It was also highlighted that the Arg185 and Lys188 residues of CAP-dimer are important for the binding between CAP-dimer and DNA. These results provide useful information for proposing new compounds that efficiently control the transcription mechanism.

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.

Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks determine the levels of transcription of a given gene at a given time. Using three complementary computational techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are constrained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any requirement for external signalling factors. Our models combine atomistic molecular dynamics (MD) with coarse-grained simulations and statistical mechanical calculations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denaturation in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informational transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles.

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.

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.

Structure of catabolite activator protein with cobalt(II) and sulfate

The crystal structure of cyclic AMP-catabolite activator protein (CAP) from Escherichia coli containing cobalt(II) chloride and ammonium sulfate is reported at 1.97 Å resolution. Each of the two CAP subunits in the asymmetric unit binds one cobalt(II) ion, in each case coordinated by N-terminal domain residues His19, His21 and Glu96 plus an additional acidic residue contributed via a crystal contact. The three identified N-terminal domain cobalt-binding residues are part of a region of CAP that is important for transcription activation at class II CAP-dependent promoters. Sulfate anions mediate additional crystal lattice contacts and occupy sites corresponding to DNA backbone phosphate positions in CAP-DNA complex structures.

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.

Crystallization and preliminary X-ray analysis of the CRP-cAMP-DNA (full length) complex

The Escherichia coli cyclic AMP receptor protein (CRP) is a well known transcription activator protein. In this study, CRP was overexpressed, purified and cocrystallized with cAMP and a 38 bp full-length double-stranded DNA fragment. The full-length segment differed from the half-site fragments used in previous crystallization experiments and is more similar to the environment in vivo. CRP-cAMP-DNA crystals were obtained and diffracted to 2.9 Å resolution. The crystals belonged to space group P3121, with unit-cell parameters a = b = 76.03, c = 144.00 Å. The asymmetric unit was found to contain one protein molecule and half a 38 bp full-length double-stranded DNA fragment, with a Matthews coefficient of 2.62 Å(3) Da(-1) and a solvent content of 53.14%.

Functional specificity of a protein-DNA complex mediated by two arginines bound to the minor groove

A bacteriophage Ø29 transcriptional regulator, protein p4, interacts with its DNA target by employing two mechanisms: by direct readout of the chemical signatures of only one DNA base and by inducing local modification on the topology of short A tracts (indirect readout). p4 binds as a dimer to targets consisting of imperfect inverted repeats. Here we used molecular dynamic simulation to define interactions of a cluster of 12 positively charged amino acids of p4 with DNA and biochemical assays with modified DNA targets and mutated proteins to quantify the contribution of residues in the nucleoprotein complex. Our results show the implication of Arg54, with non-base-specific interaction in the central A tract, in p4 binding affinity. Despite being chemically equivalent and in identical protein monomers, the two Arg54 residues differed in their interactions with DNA. We discuss an indirect-readout mechanism for p4-DNA recognition mediated by dissimilar interaction of arginines penetrating the minor groove and the inherent properties of the A tract. Our findings extend the current understanding of protein-DNA recognition and contribute to the relevance of the sequence-dependent conformational malleability of the DNA, shedding light on the role of arginines in binding affinity. Characterization of mutant p4R54A shows that the residue is required for the activity of the protein as a transcriptional regulator.

Crystal structure of the Pseudomonas aeruginosa virulence factor regulator

Virulence factor regulator (Vfr) enhances Pseudomonas aeruginosa pathogenicity through its role as a global transcriptional regulator. The crystal structure of Vfr shows that it is a winged-helix DNA-binding protein like its homologue cyclic AMP receptor protein (CRP). In addition to an expected primary cyclic AMP-binding site, a second ligand-binding site is nestled between the N-terminal domain and the C-terminal helix-turn-helix domain. Unlike CRP, Vfr is a symmetric dimer in the absence of DNA. Removal of seven disordered N-terminal residues of Vfr prevents the growth of P. aeruginosa.

Comparing binding site information to binding affinity reveals that Crp/DNA complexes have several distinct binding conformers

We show that the cAMP receptor protein (Crp) binds to DNA as several different conformers. This situation has precluded discovering a high correlation between any sequence property and binding affinity for proteins that bend DNA. Experimentally quantified affinities of Synechocystis sp. PCC 6803 cAMP receptor protein (SyCrp1), the Escherichia coli Crp (EcCrp, also CAP) and DNA were analyzed to mathematically describe, and make human-readable, the relationship of DNA sequence and binding affinity in a given system. Here, sequence logos and weight matrices were built to model SyCrp1 binding sequences. Comparing the weight matrix model to binding affinity revealed several distinct binding conformations. These Crp/DNA conformations were asymmetrical (non-palindromic).

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Molecular interactions and protein-induced DNA hairpin in the transcriptional control of bacteriophage ø29 DNA

Studies on the regulation of phage Ø29 gene expression revealed a new mechanism to accomplish simultaneous activation and repression of transcription leading to orderly gene expression. Two phage-encoded early proteins, p4 and p6, bind synergistically to DNA, modifying the topology of the sequences encompassing early promoters A2c and A2b and late promoter A3 in a hairpin that allows the switch from early to late transcription. Protein p6 is a nucleoid-like protein that binds DNA in a non-sequence specific manner. Protein p4 is a sequence-specific DNA binding protein with multifaceted sequence-readout properties. The protein recognizes the chemical signature of only one DNA base on the inverted repeat of its target sequence through a direct-readout mechanism. In addition, p4 specific binding depends on the recognition of three A-tracts by indirect-readout mechanisms. The biological importance of those three A-tracts resides in their individual properties rather than in the global curvature that they may induce.

Monitoring multiple distances within a single molecule using switchable FRET

The analysis of structure and dynamics of biomolecules is important for understanding their function. Toward this aim, we introduce a method called 'switchable FRET', which combines single-molecule fluorescence resonance energy transfer (FRET) with reversible photoswitching of fluorophores. Typically, single-molecule FRET is measured within a single donor-acceptor pair and reports on only one distance. Although multipair FRET approaches that monitor multiple distances have been developed, they are technically challenging and difficult to extend, mainly because of their reliance on spectrally distinct acceptors. In contrast, switchable FRET sequentially probes FRET between a single donor and spectrally identical photoswitchable acceptors, dramatically reducing the experimental and analytical complexity and enabling direct monitoring of multiple distances. Our experiments on DNA molecules, a protein-DNA complex and dynamic Holliday junctions demonstrate the potential of switchable FRET for studying dynamic, multicomponent biomolecules.

DNA bending and looping in the transcriptional control of bacteriophage phi29

Recent studies on the regulation of phage phi29 gene expression reveal new ways to accomplish the processes required for the orderly gene expression in prokaryotic systems. These studies revealed a novel DNA-binding domain in the phage main transcriptional regulator and the nature and dynamics of the multimeric DNA-protein complex responsible for the switch from early to late gene expression. This review describes the features of the regulatory mechanism that leads to the simultaneous activation and repression of transcription, and discusses it in the context of the role of the topological modification of the DNA carried out by two phage-encoded proteins working synergistically with the DNA.

Profound asymmetry in the structure of the cAMP-free cAMP Receptor Protein (CRP) from Mycobacterium tuberculosis

The cyclic AMP receptor protein (CRP, also called catabolite gene activator protein or CAP) plays a key role in metabolic regulation in bacteria and has become a widely studied model allosteric transcription factor. On binding its effector cAMP in the N-terminal domain, CRP undergoes a structural transition to a conformation capable of specific DNA binding in the C-terminal domain and transcription initiation. The crystal structures of Escherichia coli CRP (EcCRP) in the cAMP-bound state, both with and without DNA, are known, although its structure in the off state (cAMP-free, apoCRP) remains unknown. We describe the crystal structure at 2.0A resolution of the cAMP-free CRP homodimer from Mycobacterium tuberculosis H(37)R(v) (MtbCRP), whose sequence is 30% identical with EcCRP, as the first reported structure of an off-state CRP. The overall structure is similar to that seen for the cAMP-bound EcCRP, but the apo MtbCRP homodimer displays a unique level of asymmetry, with a root mean square deviation of 3.5A between all Calpha positions in the two subunits. Unlike structures of on-state EcCRP and other homologs in which the C-domains are asymmetrically positioned but possess the same internal conformation, the two C-domains of apo MtbCRP differ both in hinge structure and in internal arrangement, with numerous residues that have completely different local environments and hydrogen bond interactions, especially in the hinge and DNA-binding regions. Comparison of the structures of apo MtbCRP and DNA-bound EcCRP shows how DNA binding would be inhibited in the absence of cAMP and supports a mechanism involving functional asymmetry in apoCRP.

DNA sequence-directed organization of chromatin: structure-based computational analysis of nucleosome-binding sequences

The folding of DNA on the nucleosome core particle governs many fundamental issues in eukaryotic molecular biology. In this study, an updated set of sequence-dependent empirical "energy" functions, derived from the structures of other protein-bound DNA molecules, is used to investigate the extent to which the architecture of nucleosomal DNA is dictated by its underlying sequence. The potentials are used to estimate the cost of deforming a collection of sequences known to bind or resist uptake in nucleosomes along various left-handed superhelical pathways and to deduce the features of sequence contributing to a particular structural form. The deformation scores reflect the choice of template, the deviations of structural parameters at each step of the nucleosome-bound DNA from their intrinsic values, and the sequence-dependent "deformability" of a given dimer. The correspondence between the computed scores and binding propensities points to a subtle interplay between DNA sequence and nucleosomal folding, e.g., sequences with periodically spaced pyrimidine-purine steps deform at low cost along a kinked template whereas sequences that resist deformation prefer a smoother spatial pathway. Successful prediction of the known settings of some of the best-resolved nucleosome-positioning sequences, however, requires a template with "kink-and-slide" steps like those found in high-resolution nucleosome structures.

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.

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.

Murine leukemia virus reverse transcriptase: structural comparison with HIV-1 reverse transcriptase

Recent X-ray crystal structure determinations of Moloney murine leukemia virus reverse transcriptase (MoMLV RT) have allowed for more accurate structure/function comparisons to HIV-1 RT than were formerly possible. Previous biochemical studies of MoMLV RT in conjunction with knowledge of sequence homologies to HIV-1 RT and overall fold similarities to RTs in general, provided a foundation upon which to build. In addition, numerous crystal structures of the MoMLV RT fingers/palm subdomain had also shed light on one of the critical functions of the enzyme, specifically polymerization. Now in the advent of new structural information, more intricate examination of MoMLV RT in its entirety can be realized, and thus the comparisons with HIV-1 RT may be more critically elucidated. Here, we will review the similarities and differences between MoMLV RT and HIV-1 RT via structural analysis, and propose working models for the MoMLV RT based upon that information.

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.

A computational investigation of allostery in the catabolite activator protein

The catabolite activator protein is a dimer that consists of two cAMP-binding subunits, each containing a C-terminus DNA-binding module and a N-terminus ligand binding domain. The system is well-known to exhibit negative cooperativity, whereby the binding of one cAMP molecule reduces the binding affinity of the other cAMP molecule by 2 orders of magnitude, despite the large separation between the cAMP binding pockets. Here we use extensive explicit-solvent molecular dynamics simulations (135 ns) to investigate the allosteric mechanism of CAP. Six trajectories were carried out for apo, singly liganded, and doubly liganded CAP, both in the presence and absence of DNA. Thorough analyses of the dynamics through the construction of dynamical cross-correlated maps, as well as essential dynamics analyses, indicated that the system experienced a switch in motion as a result of cAMP binding, in accordance with recent NMR experiments carried out on a truncated form of the protein. Analyses of conformer structures collected from the simulations revealed a remarkable event: the DNA-binding module was found to dissociate from the N-terminus ligand binding domain. An interesting aspect of this structural change is that it only occurred in unoccupied subunits, suggesting that the binding of cAMP provides additional stability to the system, consistent with the increase in entropy that was observed in our calculations and from isothermal titration calorimetry. Analysis of the distribution of intrinsic disorder propensities in CAP amino acid sequence using PONDR VLXT and VSL1 predictors revealed that the region connecting ligand-binding and DNA-binding domains of CAP have the potential to exhibit increased flexibility. We complemented these trajectories with free energy calculations following the MM-PBSA approach on more than 2000 snapshots that included 880 normal mode analysis. The resulting free energy differences between the singly liganded and doubly liganded states were in excellent agreement with isothermal titration calorimetry data. When the free energy calculations were carried out in the presence of DNA, we discovered that a switch in cooperativity occurred, so that the binding of the first cAMP promoted the binding of the other cAMP. The components of the free energy reveal that this effect is mainly entropic in nature, whereby the DNA reduces the degree of tightening that is observed in its absence, thereby promoting binding of the second cAMP. This finding prompted us to propose a new mechanism by which CAP triggers the transcription activation that is based on an order to disorder transition mediated by cAMP binding as well as DNA.

Understanding cAMP-dependent allostery by NMR spectroscopy: comparative analysis of the EPAC1 cAMP-binding domain in its apo and cAMP-bound states

cAMP (adenosine 3',5'-cyclic monophosphate) is a ubiquitous second messenger that activates a multitude of essential cellular responses. Two key receptors for cAMP in eukaryotes are protein kinase A (PKA) and the exchange protein directly activated by cAMP (EPAC), which is a recently discovered guanine nucleotide exchange factor (GEF) for the small GTPases Rap1 and Rap2. Previous attempts to investigate the mechanism of allosteric activation of eukaryotic cAMP-binding domains (CBDs) at atomic or residue resolution have been hampered by the instability of the apo form, which requires the use of mixed apo/holo systems, that have provided only a partial picture of the CBD apo state and of the allosteric networks controlled by cAMP. Here, we show that, unlike other eukaryotic CBDs, both apo and cAMP-bound states of the EPAC1 CBD are stable under our experimental conditions, providing a unique opportunity to define at an unprecedented level of detail the allosteric interactions linking two critical functional sites of this CBD. These are the phosphate binding cassette (PBC), where cAMP binds, and the N-terminal helical bundle (NTHB), which is the site of the inhibitory interactions between the regulatory and catalytic regions of EPAC. Specifically, the combined analysis of the cAMP-dependent changes in chemical shifts, 2 degrees structure probabilities, hydrogen/hydrogen exchange (H/H) and hydrogen/deuterium exchange (H/D) protection factors reveals that the long-range communication between the PBC and the NTHB is implemented by two distinct intramolecular cAMP-signaling pathways, respectively, mediated by the beta2-beta3 loop and the alpha6 helix. Docking of cAMP into the PBC perturbs the NTHB inner core packing and the helical probabilities of selected NTHB residues. The proposed model is consistent with the allosteric role previously hypothesized for L273 and F300 based on site-directed mutagenesis; however, our data show that such a contact is part of a significantly more extended allosteric network that, unlike PKA, involves a tight coupling between the alpha- and beta-subdomains of the EPAC CBD. The proposed mechanism of allosteric activation will serve as a basis to understand agonism and antagonism in the EPAC system and provides also a general paradigm for how small ligands control protein-protein interfaces.

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein-DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein-DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.

DNA sequence-specific recognition by a transcriptional regulator requires indirect readout of A-tracts

The bacteriophage Ø29 transcriptional regulator p4 binds to promoters of different intrinsic activities. The p4–DNA complex contains two identical protomers that make similar interactions with the target sequence 5′-AACTTTTT-15 bp-AAAATGTT-3′. To define how the various elements in the target sequence contribute to p4's affinity, we studied p4 binding to a series of mutated binding sites. The binding specificity depends critically on base pairs of the target sequence through both direct as well as indirect readout. There is only one specific contact between a base and an amino acid residue; other contacts take place with the phosphate backbone. Alteration of direct amino acid–base contacts, or mutation of non-contacted A·T base pairs at A-tracts abolished binding. We generated three 5 ns molecular dynamics (MD) simulations to investigate the basis for the p4–DNA complex specificity. Recognition is controlled by the protein and depends on DNA dynamic properties. MD results on protein–DNA contacts and the divergence of p4 affinity to modified binding sites reveal an inherent asymmetry, which is required for p4-specific binding and may be crucial for transcription regulation.

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.