Probing the physical basis for trp repressor-operator recognition

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Optimal performance of the tryptophan operon of E. coli: a stochastic, dynamical, mathematical-modeling approach

In this work, we develop a detailed, stochastic, dynamical model for the tryptophan operon of E. coli, and estimate all of the model parameters from reported experimental data. We further employ the model to study the system performance, considering the amount of biochemical noise in the trp level, the system rise time after a nutritional shift, and the amount of repressor molecules necessary to maintain an adequate level of repression, as indicators of the system performance regime. We demonstrate that the level of cooperativity between repressor molecules bound to the first two operators in the trp promoter affects all of the above enlisted performance characteristics. Moreover, the cooperativity level found in the wild-type bacterial strain optimizes a cost-benefit function involving low biochemical noise in the tryptophan level, short rise time after a nutritional shift, and low number of regulatory molecules.

Protein-DNA docking with a coarse-grained force field

Background

Protein-DNA interactions are important for many cellular processes, however structural knowledge for a large fraction of known and putative complexes is still lacking. Computational docking methods aim at the prediction of complex architecture given detailed structures of its constituents. They are becoming an increasingly important tool in the field of macromolecular assemblies, complementing particularly demanding protein-nucleic acids X ray crystallography and providing means for the refinement and integration of low resolution data coming from rapidly advancing methods such as cryoelectron microscopy.

Results

We present a new coarse-grained force field suitable for protein-DNA docking. The force field is an extension of previously developed parameter sets for protein-RNA and protein-protein interactions. The docking is based on potential energy minimization in translational and orientational degrees of freedom of the binding partners. It allows for fast and efficient systematic search for native-like complex geometry without any prior knowledge regarding binding site location.

Conclusions

We find that the force field gives very good results for bound docking. The quality of predictions in the case of unbound docking varies, depending on the level of structural deviation from bound geometries. We analyze the role of specific protein-DNA interactions on force field performance, both with respect to complex structure prediction, and the reproduction of experimental binding affinities. We find that such direct, specific interactions only partially contribute to protein-DNA recognition, indicating an important role of shape complementarity and sequence-dependent DNA internal energy, in line with the concept of indirect protein-DNA readout mechanism.

A new method for evaluating the specificity of indirect readout in protein-DNA recognition

Proteins recognize a specific DNA sequence not only through direct contact (direct readout) with base pairs but also through sequence-dependent conformation and/or flexibility of DNA (indirect readout). However, it is difficult to assess the contribution of indirect readout to the sequence specificity. What is needed is a straightforward method for quantifying its contributions to specificity. Using Bayesian statistics, we derived the probability of a particular sequence for a given DNA structure from the trajectories of molecular dynamics (MD) simulations of DNAs containing all possible tetramer sequences. Then, we quantified the specificity of indirect readout based on the information entropy associated with the probability. We tested this method with known structures of protein–DNA complexes. This method enabled us to correctly predict those regions where experiments suggested the involvement of indirect readout. The results also indicated new regions where the indirect readout mechanism makes major contributions to the recognition. The present method can be used to estimate the contribution of indirect readout without approximations to the distributions in the conformational ensembles of DNA, and would serve as a powerful tool to study the mechanism of protein–DNA recognition.

Numerical solution of the chemical master equation uniqueness and stability of the stationary distribution for chemical networks, and mRNA bursting in a gene network with negative feedback regulation

In this work, we introduce a couple of algorithms to compute the stationary probability distribution for the chemical master equation (CME) of arbitrary chemical networks. We further find the conditions guaranteeing the algorithms' convergence and the unity and stability of the stationary distribution. Next, we employ these algorithms to study the mRNA and protein probability distributions in a gene regulatory network subject to negative feedback regulation. In particular, we analyze the influence of the promoter activation/deactivation speed on the shape of such distributions. We find that a reduction of the promoter activation/deactivation speed modifies the shape of those distributions in a way consistent with the phenomenon known as mRNA (or transcription) bursting.

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.

Physical basis of the inducer-dependent cooperativity of the Central glycolytic genes Repressor/DNA complex

The Central glycolytic genes Repressor (CggR) from Bacillus subtilis belongs to the SorC family of transcription factors that control major carbohydrate metabolic pathways. Recent studies have shown that CggR binds as a tetramer to its tandem operator DNA sequences and that the inducer metabolite, fructose 1,6-bisphosphate (FBP), reduces the binding cooperativity of the CggR/DNA complex. Here, we have determined the effect of FBP on the size, shape and stoichiometry of CggR complexes with full-length and half-site operator sequence by small-angle X-ray scattering, size-exclusion chromatography, fluorescence cross-correlation spectroscopy and noncovalent mass spectrometry (MS). Our results show that CggR forms a compact tetrameric assembly upon binding to either the full-length operator or two half-site DNAs and that FBP triggers a tetramer-dimer transition that leaves a single dimer on the half-site or two physically independent dimers on the full-length target. Although the binding of other phospho-sugars was evidenced by MS, only FBP was found to completely disrupt dimer-dimer contacts. We conclude that inducer-dependent dimer-dimer bridging interactions constitute the physical basis for CggR cooperative binding to DNA and the underlying repression mechanism. This work provides experimental evidences for a cooperativity-based regulation model that should apply to other SorC family members.

Cycling expression and cooperative operator interaction in the trp operon of Escherichia coli

Oscillatory behaviour in the tryptophan operon of an Escherichia coli mutant strain lacking the enzyme-inhibition regulatory mechanism has been observed by Bliss et al. but not confirmed by others. This behaviour could be important from the standpoint of synthetic biology, whose goals include the engineering of intracellular genetic oscillators. This work is devoted to investigating, from a mathematical modelling point of view, the possibility that the trp operon of the E. coli inhibition-free strain expresses cyclically. For that we extend a previously introduced model for the regulatory pathway of the tryptophan operon in Escherichia coli to account for the observed multiplicity and cooperativity of repressor binding sites. Thereafter we investigate the model dynamics using deterministic numeric solutions, stochastic simulations, and analytic studies. Our results suggest that a quasi-periodic behaviour could be observed in the trp operon expression level of single bacteria.

Biosynthesis of the Aromatic Amino Acids

This chapter describes in detail the genes and proteins of Escherichia coli involved in the biosynthesis and transport of the three aromatic amino acids tyrosine, phenylalanine, and tryptophan. It provides a historical perspective on the elaboration of the various reactions of the common pathway converting erythrose-4-phosphate and phosphoenolpyruvate to chorismate and those of the three terminal pathways converting chorismate to phenylalanine, tyrosine, and tryptophan. The regulation of key reactions by feedback inhibition, attenuation, repression, and activation are also discussed. Two regulatory proteins, TrpR (108 amino acids) and TyrR (513 amino acids), play a major role in transcriptional regulation. The TrpR protein functions only as a dimer which, in the presence of tryptophan, represses the expression of trp operon plus four other genes (the TrpR regulon). The TyrR protein, which can function both as a dimer and as a hexamer, regulates the expression of nine genes constituting the TyrR regulon. TyrR can bind each of the three aromatic amino acids and ATP and under their influence can act as a repressor or activator of gene expression. The various domains of this protein involved in binding the aromatic amino acids and ATP, recognizing DNA binding sites, interacting with the alpha subunit of RNA polymerase, and changing from a monomer to a dimer or a hexamer are all described. There is also an analysis of the various strategies which allow TyrR in conjunction with particular amino acids to differentially affect the expression of individual genes of the TyrR regulon.

In vitro assays to characterize inhibitors of the activation of small G proteins by their guanine nucleotide exchange factors

Guanine nucleotide exchange factors (GEFs) are essential regulators of the spatiotemporal conditions of small GTP-binding protein (SMG) activation. Their cellular activities combine the biochemical stimulation of GDP/GTP exchange, which leads to the active conformation of the SMG, to the detection of upstream signals and, in some cases, interaction with downstream effectors. Inhibition of GEF activities by small molecules has become recently a very active field, both for understanding biology with the tools of chemistry and because GEFs are emerging as therapeutic targets. The natural compound brefeldin A (BFA) was the first inhibitor of a GEF to be characterized, and several inhibitors of SMG activation have since been discovered using a variety of screening methods. An essential step toward their use in basic research or as leads in therapeutics is the characterization of their mechanism of inhibition. GEFs function according to a multistep mechanism, involving transient ternary (nucleotide-bound) and binary (nucleotide-free) intermediates. This mechanism thereby offers many opportunities for blockage, but a thorough analysis is necessary to define the inhibition mechanism and the steps of the reaction that are affected by the inhibitor. Here, based on the case study of how BFA inhibits the activation of Arf activation by Sec7 domains, we describe a flowchart of assays to decipher the mechanism of inhibitors of the activation of SMGs by their GEFs.

Structure-based discovery of an inhibitor of Arf activation by Sec7 domains through targeting of protein-protein complexes

Small molecules that produce nonfunctional protein-protein complexes are an alternative to competitive inhibitors for the inhibition of protein functions. Here we target the activation of the small GTP-binding protein Arf1, a major regulator of membrane traffic, by the Sec7 catalytic domain of its guanine nucleotide exchange factor ARNO. The crystal structure of the Arf1-GDP/ARNO complex, which initiates the exchange reaction, was used to discover an inhibitor, LM11, using in silico screening of a flexible pocket near the Arf1/ARNO interface. Using fluorescence kinetics and anisotropy, NMR spectroscopy and mutagenesis, we show that LM11 acts following a noncompetitive mechanism in which the inhibitor targets both Arf1-GDP and the Arf1-GDP/ARNO complex and produces a nonfunctional Arf-GDP/ARNO complex whose affinity is similar to that of the native complex. In addition, LM11 recognizes features of both Arf and ARNO near the Arf/Sec7 interface, a characteristic reminiscent of the paradigm interfacial inhibitor Brefeldin A. We then show that LM11 is a cell-active inhibitor that impairs Arf-dependent trafficking structures at the Golgi. Furthermore, LM11 inhibits ARNO-dependent migration of Madin-Darby canine kidney (MDCK) cells, demonstrating that ARNO is a target of LM11 in cells. Remarkably, LM11 inhibits the activation of Arf1 but not Arf6 in vitro, pointing to a possible synergy between Arf1 and Arf6 activation by ARNO in cell migration. Our design method shows that flexible regions in protein-protein complexes provide drugable sites with the potential to develop novel tools for investigating and inhibiting signaling pathways.

Relationship between the oligomeric status of HIV-1 integrase on DNA and enzymatic activity

The 3'-processing of the extremities of viral DNA is the first of two reactions catalyzed by HIV-1 integrase (IN). High order IN multimers (tetramers) are required for complete integration, but it remains unclear which oligomer is responsible for the 3'-processing reaction. Moreover, IN tends to aggregate, and it is unknown whether the polymerization or aggregation of this enzyme on DNA is detrimental or beneficial for activity. We have developed a fluorescence assay based on anisotropy for monitoring release of the terminal dinucleotide product in real-time. Because the initial anisotropy value obtained after DNA binding and before catalysis depends on the fractional saturation of DNA sites and the size of IN.DNA complexes, this approach can be used to study the relationship between activity and binding/multimerization parameters in the same assay. By increasing the IN:DNA ratio, we found that the anisotropy increased but the 3'-processing activity displayed a characteristic bell-shaped behavior. The anisotropy values obtained in the first phase were predictive of subsequent activity and accounted for the number of complexes. Interestingly, activity peaked and then decreased in the second phase, whereas anisotropy continued to increase. Time-resolved fluorescence anisotropy studies showed that the most competent form for catalysis corresponds to a dimer bound to one viral DNA end, whereas higher order complexes such as aggregates predominate during the second phase when activity drops off. We conclude that a single IN dimer at each extremity of viral DNA molecules is required for 3'-processing, with a dimer of dimers responsible for the subsequent full integration.

Protein-DNA binding specificity predictions with structural models

Protein-DNA interactions play a central role in transcriptional regulation and other biological processes. Investigating the mechanism of binding affinity and specificity in protein-DNA complexes is thus an important goal. Here we develop a simple physical energy function, which uses electrostatics, solvation, hydrogen bonds and atom-packing terms to model direct readout and sequence-specific DNA conformational energy to model indirect readout of DNA sequence by the bound protein. The predictive capability of the model is tested against another model based only on the knowledge of the consensus sequence and the number of contacts between amino acids and DNA bases. Both models are used to carry out predictions of protein-DNA binding affinities which are then compared with experimental measurements. The nearly additive nature of protein-DNA interaction energies in our model allows us to construct position-specific weight matrices by computing base pair probabilities independently for each position in the binding site. Our approach is less data intensive than knowledge-based models of protein-DNA interactions, and is not limited to any specific family of transcription factors. However, native structures of protein-DNA complexes or their close homologs are required as input to the model. Use of homology modeling can significantly increase the extent of our approach, making it a useful tool for studying regulatory pathways in many organisms and cell types.

Target search of N sliding proteins on a DNA

At low to moderate ambient salt concentrations, DNA-binding proteins bind relatively tightly to DNA, and only very rarely detach. Intersegmental transfer due to DNA-looping can be excluded by applying an external pulling force to the DNA molecule. Under such conditions, we explore the targeting dynamics of N proteins sliding diffusively along DNA in search of their specific target sequence. At lower densities of binding proteins, we find a reduction of the characteristic search time proportional to N(-2), with corrections at higher concentrations. Rates for detachment and attachment of binding proteins are incorporated in the model. Our findings are in agreement with recent single molecule studies in the presence of bacteriophage T4 gene 32 protein for which the unbinding rate is much lower than the specific binding rate.

A knowledge-based energy function for protein-ligand, protein-protein, and protein-DNA complexes

We developed a knowledge-based statistical energy function for protein-ligand, protein-protein, and protein-DNA complexes by using 19 atom types and a distance-scale finite ideal-gas reference (DFIRE) state. The correlation coefficients between experimentally measured protein-ligand binding affinities and those predicted by the DFIRE energy function are around 0.63 for one training set and two testing sets. The energy function also makes highly accurate predictions of binding affinities of protein-protein and protein-DNA complexes. Correlation coefficients between theoretical and experimental results are 0.73 for 82 protein-protein (peptide) complexes and 0.83 for 45 protein-DNA complexes, despite the fact that the structures of protein-protein (peptide) and protein-DNA complexes were not used in training the energy function. The results of the DFIRE energy function on protein-ligand complexes are compared to the published results of 12 other scoring functions generated from either physical-based, knowledge-based, or empirical methods. They include AutoDock, X-Score, DrugScore, four scoring functions in Cerius 2 (LigScore, PLP, PMF, and LUDI), four scoring functions in SYBYL (F-Score, G-Score, D-Score, and ChemScore), and BLEEP. While the DFIRE energy function is only moderately successful in ranking native or near native conformations, it yields the strongest correlation between theoretical and experimental binding affinities of the testing sets and between rmsd values and energy scores of docking decoys in a benchmark of 100 protein-ligand complexes. The parameters and the program of the all-atom DFIRE energy function are freely available for academic users at http://theory.med.buffalo.edu.

Influence of DNA stiffness in protein–DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.

Kinetics of protein-DNA interaction: facilitated target location in sequence-dependent potential

Recognition and binding of specific sites on DNA by proteins is central for many cellular functions such as transcription, replication, and recombination. In the process of recognition, a protein rapidly searches for its specific site on a long DNA molecule and then strongly binds this site. Here we aim to find a mechanism that can provide both a fast search (1-10 s) and high stability of the specific protein-DNA complex (Kd=10(-15)-10(-8) M). Earlier studies have suggested that rapid search involves sliding of the protein along the DNA. Here we consider sliding as a one-dimensional diffusion in a sequence-dependent rough energy landscape. We demonstrate that, despite the landscape's roughness, rapid search can be achieved if one-dimensional sliding is accompanied by three-dimensional diffusion. We estimate the range of the specific and nonspecific DNA-binding energy required for rapid search and suggest experiments that can test our mechanism. We show that optimal search requires a protein to spend half of its time sliding along the DNA and the other half diffusing in three dimensions. We also establish that, paradoxically, realistic energy functions cannot provide both rapid search and strong binding of a rigid protein. To reconcile these two fundamental requirements we propose a search-and-fold mechanism that involves the coupling of protein binding and partial protein folding. The proposed mechanism has several important biological implications for search in the presence of other proteins and nucleosomes, simultaneous search by several proteins, etc. The proposed mechanism also provides a new framework for interpretation of experimental and structural data on protein-DNA interactions.

Intermolecular and Intramolecular Readout Mechanisms in Protein–DNA Recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.

RAR antagonists diminish the level of DNA binding by the RAR/RXR heterodimer

A purified RAR/RXR-DeltaAB heterodimer was obtained by production of His-tagged RAR and untagged RXR in Escherichia coli, followed by combined purification on a Ni(2+) affinity column using excess RXR extract, and finally a gel filtration chromatography step to isolate a pure heterodimer. The purified heterodimer preparation bound 9-cisRA at a level of 0.85-0.95 mol of binding sites per mole of protein monomer. Titration of a 26 kDa fluorescent labeled fragment of the SRC-1 coactivator protein with the purified heterodimer in the presence of the agonist 9-cisRA yielded a binding affinity near 300 nM, whereas no binding was observed in the absence of agonist. Binding of the purified heterodimer to a DR5 target was identical in the absence of ligand and in the presence of 9-cisRA. Competition by unlabeled specific and nonspecific DNA allowed us to demonstrate that the binding curve was bimodal. The first phase of binding was highly specific and of high affinity. This phase also exhibited a high degree of cooperativity in the binding profile. Nonspecific DNA efficiently competes for the second phase. Thus, the first phase of binding likely corresponds to the formation of the specific heterodimer complex in which heterodimerization is energetically coupled to DNA binding. While agonist binding had no effect on the apparent affinity of the heterodimer for DR5, a series of antagonists significantly destabilized the heterodimer-DR5 complex, either through a direct decrease in the affinity of the protein for the DNA or through destabilization of the heterodimer itself. Impeding the interaction between the heterodimer and DNA appears as an additional mechanism of antagonist action of varying efficiency, depending upon the chemical structure of the antagonist.

Trp repressor-operator binding: NMR and electrophoretic mobility shift studies of the effect of DNA sequence and corepressor binding on two Trp repressor-operator complexes

In Trp repressor-DNA complexes, most interactions either occur with phosphate groups or are water-mediated hydrogen bonds to bases. To examine the factors involved in DNA selectivity, we have studied Trp repressor binding to two operator sequences, trpR(S)() and trpO(M)(), with L-tryptophan or 5-methyltryptophan as corepressor. These operators contain all the consensus bases but differ at base pairs contacted by their phosphate groups. In electrophoretic mobility shift assays (EMSAs) the trpR(S)() sequence gives solely 1:1 protein-DNA complexes with either corepressor. The trpO(M )()sequence binds more weakly than trpR(S)(). It gives dissociating 2:1 complexes in EMSAs with L-tryptophan, but both 1:1 and 2:1 complexes are observed with 5-methyltryptophan or if glycerol is present in the gel. The backbone resonances of the TrpR-L-tryptophan-DNA complexes were assigned using triple-resonance experiments and selectively (15)N labeled protein. On changing the DNA sequence, the largest differences in the NMR spectra are at residues 78-81, at the turn of the helix-turn-helix motif and the tip of the recognition helix. I79 and A80 interact with the conserved bases of the operators, while G78 and T81 interact with phosphate groups at bases that differ between the two sequences. Changing the corepressor from L-tryptophan to 5-methyltryptophan causes effects at residues 52, 60, 61, and 85, which do not interact with the DNA. The spectra suggest that there is mutual induced fit between protein and DNA so that sequence changes at bases contacted only by the phosphate groups affect the environment of the protein at residues that bind to conserved bases elsewhere in the DNA.

Equilibrium binding assays reveal the elevated stoichiometry and salt dependence of the interaction between full-length human sex-determining region on the Y chromosome (SRY) and DNA

In an effort to better define the molecular mechanism of the functional specificity of human sex-determining region on the Y chromosome (SRY), we have carried out equilibrium binding assays to study the interaction of the full-length bacterial-expressed protein with a DNA response element derived from the CD3epsilon gene enhancer. These assays are based on the observation of the fluorescence anisotropy of a fluorescein moiety covalently bound to the target oligonucleotide. The low anisotropy value due to the fast tumbling of the free oligonucleotide in solution increases substantially upon binding the protein to the labeled target DNA. Our results indicate that the full-length human wild-type SRY (SRY(WT)) forms a complex of high stoichiometry with its target DNA. Moreover, we have demonstrated a strong salt dependence of both the affinity and specificity of the interaction. We have also addressed the DNA bending properties of full-length human SRY(WT) in solution by fluorescence resonance energy transfer and revealed that maximal bending is achieved with a protein to DNA ratio significantly higher than the classical 1:1. Oligomerization thus appears, at least in vitro, to be tightly coupled to SRY-DNA interactions. Alteration of protein-protein interactions observed for the mutant protein SRY(Y129N), identified in a patient presenting with 46,XY sex reversal, suggests that oligomerization may play an important role in vivo as well.

Structural analysis of conserved base pairs in protein-DNA complexes

Understanding of protein-DNA interactions is crucial for prediction of DNA-binding specificity of transcription factors and design of novel DNA-binding proteins. In this paper we develop a novel approach to analysis of protein-DNA interactions. We bring together two sources of information: (i) structures of protein-DNA complexes (PDB/NDB database) and (ii) experimentally obtained sites recognized by DNA-binding proteins. Sites are used to compute conservation (information content) of each base pair, which indicates relative importance of the base pair in specific recognition. The main result of this study is that conservation of base pairs in a site exhibits significant correlation with the number of contacts the base pairs have with the protein. In particular, base pairs that have more contacts with the protein are more conserved in evolution. As natural as it is, this result has never been reported before. We also observe that for most of the studied proteins, hydrogen bonds and hydrophobic interactions alone cannot explain the pattern of evolutionary conservation in the binding site suggesting cumulative contribution of different types of interactions to specific recognition. Implications for prediction of the DNA-binding specificity are discussed.

Indirect readout of the trp-repressor-operator complex by B-DNA's backbone conformation transitions

Although the trp-repressor-operator complex is one of the best studied transcriptional controlling systems, some questions regarding the specific recognition of the operator by the repressor remain. We performed a 2.35 ns long molecular dynamics simulation to clarify the influence of the two B-DNA backbone conformational substates B(I) and B(II) on complexation. The trp-repressor-operator is an ideal biological system for this study because experimental results have already figured out that the interaction between the internucleotide phosphates and the protein is essential for the formation of the high affinity complex. Our simulation supports these results, but more important it shows a strong correlation between the B(I)/B(II) phosphate substate and the number of interactions with this phosphate. In particular the B(I) <==> B(II) transitions occur synchronous to hydrogen bond breaking or formation. To the best of our knowledge, this was observed for the first time. Thus, we conclude that the sequence specific B(I)/B(II) behavior contributes via indirect readout to sequence specific recognition. These results have implication for the design of transcription-controlling drugs in view of the recently published influence of minor groove binders on the B(I)/B(II) pattern. The simulation also agrees with crystallographically observed hydration sites. This is consistent with experimental results and indicates the correctness of the model used.

Exocyclic groups in the minor groove influence the backbone conformation of DNA

Exocyclic groups in the minor groove of DNA modulate the affinity and positioning of nucleic acids to the histone protein. The addition of exocyclic groups decreases the formation of this protein-DNA complex, while their removal increases nucleosome formation. On the other hand, recent theoretical results show a strong correlation between the B(I)/B(II) phosphate backbone conformation and the hydration of the grooves of the DNA. We performed a simulation of the d(CGCGAATTCGCG)2 Drew Dickerson dodecamer and one simulation of the d(CGCIAATTCGCG)2 dodecamer in order to investigate the influence of the exocyclic amino group of guanine. The removal of the amino group introduces a higher intrinsic flexibility to DNA supporting the suggestions that make the enhanced flexibility responsible for the enlarged histone complexation affinity. This effect is attributed to changes in the destacking interactions of both strands of the DNA. The differences in the hydration of the minor groove could be the explanation of this flexibility. The changed hydration of the minor groove also leads to a different B(I)/B(II) substate pattern. Due to the fact that the histone preferentially builds contacts with the backbone of the DNA, we propose an influence of these B(I)/B(II) changes on the nucleosome formation process. Thus, we provide an additional explanation for the enhanced affinity to the histone due to removal of exocyclic groups. In terms of B(I)/B(II) we are also able to explain how minor groove binding ligands could affect the nucleosome assembly without disrupting the structure of DNA.

Phased A-tracts bind to the alpha subunit of RNA polymerase with increased affinity at low temperature

Previously we showed that the expression of a Clostridium perfringens phospholipase C gene (plc) is activated by promoter upstream phased A-tracts in a low temperature-dependent manner. In this paper we characterize the interaction between the alpha subunit of C. perfringens RNA polymerase and the phased A-tracts. Hydroxyl radical footprinting and fluorescence polarization assaying revealed that the alpha subunit binds to the minor grooves of the phased A-tracts through its C-terminal domain with increased affinity at low temperature. The result provides a molecular mechanism underlying the activation of the plc promoter by the phased A-tracts.

Characterization of variants altered at the N-terminal proline, a novel heme-axial ligand in CooA, the CO-sensing transcriptional activator

CooA, the carbon monoxide-sensing transcription factor from Rhodospirillum rubrum, binds CO through a heme moiety resulting in conformational changes that promote DNA binding. The crystal structure shows that the N-terminal Pro(2) of one subunit (Met(1) is removed post-translationally) provides one ligand to the heme of the other subunit in the CooA homodimer. To determine the importance of this novel ligand and the contiguous residues to CooA function, we have altered the N terminus through two approaches: site-directed mutagenesis and regional randomization, and characterized the resulting CooA variants. While Pro(2) appears to be optimal for CooA function, it is not essential and a variety of studied variants at this position have substantial CO-sensing function. Surprisingly, even alterations that add a residue (where Pro(2) is replaced by Met(1)-Tyr(2), for example) accumulate heme-containing CooA with functional properties that are similar to those of wild-type CooA. Other nearby residues, such as Phe(5) and Asn(6) appear to be important for either the structural integrity or the function of CooA. These results are contrasted with those previously reported for alteration of the His(77) ligand on the opposite side of the heme.

Distinctive cognate sequence discrimination, bound DNA conformation, and binding modes in the E2 C-terminal domains from prototype human and bovine papillomaviruses

The C-terminal DNA binding domain of the E2 protein is involved in transcriptional regulation and DNA replication in papillomaviruses. At low ionic strength, the domain has a tendency to form aggregates, a process readily reversible by the addition of salt. While fluorescence anisotropy measurements show a 1:1 stoichiometry at pH 5.5, we observed that a second HPV-16 E2 C-terminal dimer can bind per DNA site at pH 7.0. This was confirmed by displacement of bis-ANS binding, tryptophan fluorescence, native electrophoresis, and circular dichroism. The two binding events are nonequivalent, with a high-affinity binding involving one E2C dimer per DNA molecule with a K(D) of 0.18 +/- 0.02 nM and a lower affinity binding mode of 2.0 +/- 0.2 nM. The bovine (BPV-1) E2 C-terminal domain binds to an HPV-16 E2 site with 350-fold lower affinity than the human cognate domain and binds 7-fold less tightly even to a bovine-derived DNA site. The ability to discriminate between cognate and noncognate sequences is 50-fold higher for the human domain, and the latter is 180-fold better than the bovine at discriminating specific from nonspecific DNA. A substantial conformational change in bound DNA is observed by near-UV circular dichroism. The bovine domain imposes a different DNA conformation than that caused by the human counterpart, which could be explained by a more pronounced bent. Structure-function differences and biochemical properties of the complexes depend on the protein domain rather than on the DNA, in line with crystallographic evidence. Despite the strong sequence homology and overall folding topology, the differences observed may explain the distinctive transcriptional regulation in bovine and human viruses.

LexA repressor forms stable dimers in solution. The role of specific dna in tightening protein-protein interactions

Cooperativity in the interactions among proteins subunits and DNA is crucial for DNA recognition. LexA repressor was originally thought to bind DNA as a monomer, with cooperativity leading to tighter binding of the second monomer. The main support for this model was a high value of the dissociation constant for the LexA dimer (micromolar range). Here we show that the protein is a dimer at nanomolar concentrations under different conditions. The reversible dissociation of LexA dimer was investigated by the effects of hydrostatic pressure or urea, using fluorescence emission and polarization to monitor the dissociation process. The dissociation constant lies in the picomolar range (lower than 20 pM). LexA monomers associate with an unusual large volume change (340 ml/mol), indicating the burial of a large surface area upon dimerization. Whereas nonspecific DNA has no stabilizing effect, specific DNA induces tightening of the dimer and a 750-fold decrease in the K(d). In contrast to the previous model, a tight dimer rather than a monomer is the functional repressor. Accordingly, the LexA dimer only loses its ability to recognize a specific DNA sequence by RecA-induced autoproteolysis. Our work provides insights into the linkage between protein-protein interactions, DNA recognition, and DNA repair.

The basis for the super-repressor phenotypes of the AV77 and EK18 mutants of trp repressor

The DNA-binding properties of two super-repressor mutants of the Escherichia coli trp repressor, EK18 and AV77, have been investigated using steady-state fluorescence anisotropy measurements, in order to further elucidate the basis for their super-repressor phenotypes. Several suggestions have been previously proposed as the basis for the super-repressor phenotype of EK18 and AV77. For the negative to positive charge change EK18 mutant, increased electrostatic interactions between the EK18 mutant and the operator and increased protein-protein interactions between EK18 dimers have been suggested as contributing to the super-repressor phenotype of this mutant. We show that EK18 dimers actually bind to wild-type and variant operator sequences with a decrease in apparent cooperativity and an increase in affinity, compared to WTTR dimers. Thus, the EK18 super-repressor phenotype is not due to increased cooperative binding between EK18 dimers. These results support the hypothesis that the super-repressor phenotype of EK18 arises from increased electrostatic interactions between the mutant and DNA. In the case of the AV77 mutant, weaker binding affinity of apo-AV77 to non-specific DNA, increased selectivity of binding of AV77 for the operator, and a higher population of folded functional AV77 dimers available to bind the operator under limiting L-Trp conditions in vivo, have been proposed for the super-repressor phenotype of this mutant. We show that like the EK18 mutant, apoAV77 binds with higher affinity to non-specific DNA compared to apo-WTTR and that the holo-AV77 mutant does not bind with higher selectivity to the operator, has had been previously proposed. We therefore conclude that the super-repressor phenotype of the AV77 mutant is due to an increase in the population of folded, functional AV77 dimers, under limiting L-Trp conditions in vivo.

Studies of the Escherichia coli Trp repressor binding to its five operators and to variant operator sequences

The Escherichia coli Trp repressor binds to promoters of very different sequence and intrinsic activity. Its mode of binding to trp operator DNA has been studied extensively yet remains highly controversial. In order to examine the selectivity of the protein for DNA, we have used electromobility shift assays (EMSAs) to study its binding to synthetic DNA containing the core sequences of each of its five operators and of operator variants. Our results for DNA containing sequences of two of the operators, trpEDCBA and aroH are similar to those of previous studies. Up to three bands of lower mobility than the free DNA are obtained which are assigned to complexes of stoichiometry 1 : 1, 2 : 1 and 3 : 1 Trp repressor dimer to DNA. The mtr and aroL operators have not been studied previously in vitro. For DNA containing these sequences, we observe predominantly one retarded band in EMSA with mobility corresponding to 2 : 1 complexes. We have also obtained retardation of DNA containing the trpR operator sequence, which has only been previously obtained with super-repressor Trp mutants. This gives bands with mobilities corresponding to 1 : 1 and 2 : 1 complexes. In contrast, DNA containing containing a symmetrized trpR operator sequence, trpRs, gives a single retarded band with mobility corresponding solely to a 1 : 1 protein dimer-DNA complex. Using trpR operator variants, we show that a change in a single base pair in the core 20 base pairs can alter the number of retarded DNA bands in EMSA and the length of the DNase I footprint observed. This shows that the binding of the second dimer is sequence selective. We propose that the broad selectivity of Trp repressor coupled to tandem 2 : 1 binding, which we have observed with all five operator sequences, enables the Trp repressor to bind to a limited number of sites with diverse sequences. This allows it to co-ordinately control promoters of different intrinsic strength. This mechanism may be of importance in a number of promoters that bind multiple effector molecules.

Probing the role of water in the tryptophan repressor-operator complex

The Escherichia coli tryptophan repressor protein (TR) represses the transcription of several genes in response to the concentration of tryptophan in the environment. In the co-crystal structure of TR bound to a DNA fragment containing its target very few direct contacts between TR and the DNA were observed. In contrast, a number of solvent mediated contacts were apparent. NMR solution structures, however, did not resolve any solvent mediated bonds at the complex interface. To probe for the role of water in TR operator recognition, the effect of osmolytes on the interactions between TR and a target oligonucleotide bearing the operator site was examined. In the absence of specific solvent mediated hydrogen bonding interactions between the protein and the DNA, increasing osmolyte concentration is expected to strongly stabilize the TR operator interaction due to the large amount of macromolecular surface area buried upon complexation. The results of our studies indicate that xylose did not alter the binding affinity significantly, while glycerol and PEG had a small stabilizing effect. A study of binding as a function of betaine concentration revealed that this osmolyte at low concentration results in a stabilization of the 1:1 TR/operator complex, but at higher concentrations leads to a switching between binding modes to favor tandem binding. Analysis of the effects of betaine on the 1:1 complex suggest that this osmolyte has about 78% of the expected effect. If one accepts the analysis in terms of the number of water molecules excluded upon complexation, these results suggest that about 75 water molecules remain at the interface of the 1:1 dimer/DNA complex. This value is consistent with the number of water molecules found at the interface in the crystallographically determined structure and supports the notion that interfacial waters play an important thermodynamic role in the specific complexation of one TR dimer with its target DNA. However, the complexity of the effects of betaine and the small or negligible effects of the other osmolytes could also arise from osmolyte induced competition between antagonistic coupled reactions.