Structure of the repressor-operator complex of bacteriophage 434

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Principles of Amino Acid and Nucleotide Revealed by Binding Affinities between Homogeneous Oligopeptides and Single-stranded DNA Molecule s

We have determined the binding strengths between nucleotides of adenine, thymine, guanine and cytosine in homogeneous single stranded DNAs and homo-octapeptides consisting of 20 common amino acids. We use a bead-based fluorescence assay for these measurements in which octapeptides are immobilized on the bead surface and ssDNAs are in solutions. The results provide a molecular basis for analyzing selectivity, specificity and polymorphisms of amino-acid-nucleotide interactions. Comparative analyses of the distribution of the binding energies reveal unique binding strengths patterns assignable to each pair of DNA nucleotide and amino acid originating from the chemical structures. Pronounced favorable (such as Arg-G , etc.) and unfavorable (such as Ile-T , etc.) binding interactions can be identified in selected groups of amino acid and nucleotide pairs that could provide basis to elucidate energetics of amino-acid-nucleotide interactions. Such interaction selectivity, specificity and polymorphism manifest the contributions from DNA backbone, DNA bases, as well as main chain and side chain of the amino acids.

Reading the chromatinized genome

Eukaryotic DNA-binding proteins operate in the context of chromatin, where nucleosomes are the elementary building blocks. Nucleosomal DNA is wrapped around a histone core, thereby rendering a large fraction of the DNA surface inaccessible to DNA-binding proteins. Nevertheless, first responders in DNA repair and sequence-specific transcription factors bind DNA target sites obstructed by chromatin. While early studies examined protein binding to histone-free DNA, it is only now beginning to emerge how DNA sequences are interrogated on nucleosomes. These readout strategies range from the release of nucleosomal DNA from histones, to rotational/translation register shifts of the DNA motif, and nucleosome-specific DNA binding modes that differ from those observed on naked DNA. Since DNA motif engagement on nucleosomes strongly depends on position and orientation, we argue that motif location and nucleosome positioning co-determine protein access to DNA in transcription and DNA repair.

How structural biology transformed studies of transcription regulation

The past 4 decades have seen remarkable advances in our understanding of the structural basis of gene regulation. Technological advances in protein expression, nucleic acid synthesis, and structural biology made it possible to study the proteins that regulate transcription in the context of ever larger complexes containing proteins bound to DNA. This review, written on the occasion of the 50th anniversary of the founding of the Protein Data Bank focuses on the insights gained from structural studies of protein-DNA complexes and the role the PDB has played in driving this research. I cover highlights in the field, beginning with X-ray crystal structures of the first DNA-binding domains to be studied, through recent cryo-EM structures of transcription factor binding to nucleosomal DNA.

Motif orientation matters: Structural characterization of TEAD1 recognition of genomic DNA

TEAD transcription factors regulate gene expression through interactions with DNA and other proteins. They are crucial for the development of eukaryotic organisms and to control the expression of genes involved mostly in cell proliferation and differentiation; however, their deregulation can lead to tumorigenesis. To study the interactions of TEAD1 with M-CAT motifs and their inverted versions, the K_D of each complex was determined, and H/D exchange, quantitative chemical cross-linking, molecular docking, and smFRET were utilized for structural characterization. ChIP-qPCR was employed to correlate the results with a cell line model. The results obtained showed that although the inverted motif has 10× higher K_D, the same residues were affected by the presence of M-CAT in both orientations. Molecular docking and smFRET revealed that TEAD1 binds the inverted motif rotated 180°. In addition, the inverted motif was proven to be occupied by TEAD1 in Jurkat cells, suggesting that the low-affinity binding sites present in the human genome may possess biological relevance.

Extreme divergence between one-to-one orthologs: the structure of N15 Cro bound to operator DNA and its relationship to the λ Cro complex

The gene cro promotes lytic growth of phages through binding of Cro protein dimers to regulatory DNA sites. Most Cro proteins are one-to-one orthologs, yet their sequence, structure and binding site sequences are quite divergent across lambdoid phages. We report the cocrystal structure of bacteriophage N15 Cro with a symmetric consensus site. We contrast this complex with an orthologous structure from phage λ, which has a dissimilar binding site sequence and a Cro protein that is highly divergent in sequence, dimerization interface and protein fold. The N15 Cro complex has less DNA bending and smaller DNA-induced changes in protein structure. N15 Cro makes fewer direct contacts and hydrogen bonds to bases, relying mostly on water-mediated and Van der Waals contacts to recognize the sequence. The recognition helices of N15 Cro and λ Cro make mostly nonhomologous and nonanalogous contacts. Interface alignment scores show that half-site binding geometries of N15 Cro and λ Cro are less similar to each other than to distantly related CI repressors. Despite this divergence, the Cro family shows several code-like protein–DNA sequence covariations. In some cases, orthologous genes can achieve a similar biological function using very different specific molecular interactions.

Engineering altered protein-DNA recognition specificity

Protein engineering is used to generate novel protein folds and assemblages, to impart new properties and functions onto existing proteins, and to enhance our understanding of principles that govern protein structure. While such approaches can be employed to reprogram protein-protein interactions, modifying protein-DNA interactions is more difficult. This may be related to the structural features of protein-DNA interfaces, which display more charged groups, directional hydrogen bonds, ordered solvent molecules and counterions than comparable protein interfaces. Nevertheless, progress has been made in the redesign of protein-DNA specificity, much of it driven by the development of engineered enzymes for genome modification. Here, we summarize the creation of novel DNA specificities for zinc finger proteins, meganucleases, TAL effectors, recombinases and restriction endonucleases. The ease of re-engineering each system is related both to the modularity of the protein and the extent to which the proteins have evolved to be capable of readily modifying their recognition specificities in response to natural selection. The development of engineered DNA binding proteins that display an ideal combination of activity, specificity, deliverability, and outcomes is not a fully solved problem, however each of the current platforms offers unique advantages, offset by behaviors and properties requiring further study and development.

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Genes are sections of DNA that carry the instructions needed to build other molecules including all the proteins that the cell needs to fulfill its role. The information in the DNA is stored as a code consisting of four chemical bases, often referred to simply as “A”, “C”, “G” and “T”. The order or sequence of these bases determines the role of a protein. Many organisms – including humans – are built of many different types of cells that perform unique roles. Almost all cells carry the same genetic information, but proteins called transcription factors can regulate the activity of genes so that only a relevant subset of genes is switched on at a particular time.

Transcription factors glide along DNA and bind to short DNA sequences by attaching to the DNA bases directly or through bridges made up of water molecules. Two physical concepts known as enthalpy and entropy determine the strength of the connection. Enthalpy relates to how strong the chemical bonds that form between the transcription factors and the DNA bases are, compared to a situation where the transcription factor and DNA do not form a complex and bind to water molecules around them. Entropy measures the disorder of the system – the more disordered the solvent and protein-DNA complex are compared to solvent-containing free DNA and protein, the stronger the binding. A water molecule that bridges a DNA base with an amino-acid of a protein contributes to enthalpy, but results in loss of entropy, because the system becomes more ordered since the water molecule can no longer move freely.

Most transcription factors can only bind to DNA sequences that are very similar to each other, but some transcription factors can recognize several different kinds of sequences, and until now it was not clear how they could do this. Morgunova et al. studied four different human transcription factors that can each bind to two distinct DNA sequences. The results showed that the transcription factors bound to both DNA sequences with similar strength, but via different mechanisms. For one DNA sequence, an enthalpy-based mechanism essentially ‘froze’ the transcription factor to the DNA through rigid water bridges. The other DNA sequence was bound equally strongly but through moving water molecules, because this increased the entropy of the system. It is possible that these mechanisms could also apply to many other molecules that interact with each other through water-molecule bridges.

A better knowledge of the chemical bonds between transcription factors and DNA bases may in future help efforts to develop new treatments that depend on molecules being able to bind to other molecules. In addition, these findings may one day help scientists to predict how strongly two molecules will interact simply by knowing the structures of the molecules involved.

Structures of an all-α protein running along the DNA major groove

Despite over 3300 protein–DNA complex structures have been reported in the past decades, there remain some unknown recognition patterns between protein and target DNA. The silkgland-specific transcription factor FMBP-1 from the silkworm Bombyx mori contains a unique DNA-binding domain of four tandem STPRs, namely the score and three amino acid peptide repeats. Here we report three structures of this STPR domain (termed BmSTPR) in complex with DNA of various lengths. In the presence of target DNA, BmSTPR adopts a zig-zag structure of three or four tandem α-helices that run along the major groove of DNA. Structural analyses combined with binding assays indicate BmSTPR prefers the AT-rich sequences, with each α-helix covering a DNA sequence of 4 bp. The successive AT-rich DNAs adopt a wider major groove, which is in complementary in shape and size to the tandem α-helices of BmSTPR. Substitutions of DNA sequences and affinity comparison further prove that BmSTPR recognizes the major groove mainly via shape readout. Multiple-sequence alignment suggests this unique DNA-binding pattern should be highly conserved for the STPR domain containing proteins which are widespread in animals. Together, our findings provide structural insights into the specific interactions between a novel DNA-binding protein and a unique deformed B-DNA.

A glimpse of structural biology through X-ray crystallography

Since determination of the myoglobin structure in 1957, X-ray crystallography, as the anchoring tool of structural biology, has played an instrumental role in deciphering the secrets of life. Knowledge gained through X-ray crystallography has fundamentally advanced our views on cellular processes and greatly facilitated development of modern medicine. In this brief narrative, I describe my personal understanding of the evolution of structural biology through X-ray crystallography-using as examples mechanistic understanding of protein kinases and integral membrane proteins-and comment on the impact of technological development and outlook of X-ray crystallography.

In Vitro Reconstitution and Crystallization of Cas9 Endonuclease Bound to a Guide RNA and a DNA Target

The programmable RNA-guided DNA cleavage activity of the bacterial CRISPR-associated endonuclease Cas9 is the basis of genome editing applications in numerous model organisms and cell types. In a binary complex with a dual crRNA:tracrRNA guide or single-molecule guide RNA, Cas9 targets double-stranded DNAs harboring sequences complementary to a 20-nucleotide segment in the guide RNA. Recent structural studies of the enzyme have uncovered the molecular mechanism of RNA-guided DNA recognition. Here, we provide protocols for electrophoretic mobility shift and fluorescence-detection size exclusion chromatography assays used to probe DNA binding by Cas9 that allowed us to reconstitute and crystallize the enzyme in a ternary complex with a guide RNA and a bona fide target DNA. The procedures can be used for further mechanistic investigations of the Cas9 endonuclease family and are potentially applicable to other multicomponent protein-nucleic acid complexes.

Veritas per structuram

When I entered graduate school in 1963, the golden age of molecular biology had just begun, and myoglobin was the only protein with a known high-resolution structure. The romance of working out the structure of a virus by X-ray crystallography nonetheless captured both my imagination and the ensuing 15 years of my scientific life, during which "protein crystallography" began to morph into "structural biology." The course of the research recounted here follows the broader, 50-year trajectory of structural biology, as I could rarely resist opportunities to capitalize on new technologies when they opened some interesting part of biology to three-dimensional rigor. That fascination shows no sign of subsiding.

Promoter activation by CII, a potent transcriptional activator from bacteriophage 186

The lysogeny promoting protein CII from bacteriophage 186 is a potent transcriptional activator, capable of mediating at least a 400-fold increase in transcription over basal activity. Despite being functionally similar to its counterpart in phage λ, it shows no homology at the level of protein sequence and does not belong to any known family of transcriptional activators. It also has the unusual property of binding DNA half-sites that are separated by 20 base pairs, center to center. Here we investigate the structural and functional properties of CII using a combination of genetics, in vitro assays, and mutational analysis. We find that 186 CII possesses two functional domains, with an independent activation epitope in each. 186 CII owes its potent activity to activation mechanisms that are dependent on both the σ(70) and α C-terminal domain (αCTD) components of RNA polymerase, contacting different functional domains. We also present evidence that like λ CII, 186 CII is proteolytically degraded in vivo, but unlike λ CII, 186 CII proteolysis results in a specific, transcriptionally inactive, degradation product with altered self-association properties.

The N-terminal domain of the repressor of Staphylococcus aureus phage Φ11 possesses an unusual dimerization ability and DNA binding affinity

Bacteriophage Φ11 uses Staphylococcus aureus as its host and, like lambdoid phages, harbors three homologous operators in between its two divergently oriented repressor genes. None of the repressors of Φ11, however, showed binding to all three operators, even at high concentrations. To understand why the DNA binding mechanism of Φ11 repressors does not match that of lambdoid phage repressors, we studied the N-terminal domain of the Φ11 lysogenic repressor, as it harbors a putative helix-turn-helix motif. Our data revealed that the secondary and tertiary structures of the N-terminal domain were different from those of the full-length repressor. Nonetheless, the N-terminal domain was able to dimerize and bind to the operators similar to the intact repressor. In addition, the operator base specificity, binding stoichiometry, and binding mechanism of this domain were nearly identical to those of the whole repressor. The binding affinities of the repressor and its N-terminal domain were reduced to a similar extent when the temperature was increased to 42°C. Both proteins also adequately dislodged a RNA polymerase from a Φ11 DNA fragment carrying two operators and a promoter. Unlike the intact repressor, the binding of the N-terminal domain to two adjacent operator sites was not cooperative in nature. Taken together, we suggest that the dimerization and DNA binding abilities of the N-terminal domain of the Φ11 repressor are distinct from those of the DNA binding domains of other phage repressors.

PROTEIN STRUCTURE: INSIGHTS FROM GRAPH THEORY

The sequence and structure of a large body of proteins are becoming increasingly available. It is desirable to explore mathematical tools for efficient extraction of information from such sources. The principles of graph theory, which was earlier applied in fields such as electrical engineering and computer networks are now being adopted to investigate protein structure, folding, stability, function and dynamics. This review deals with a brief account of relevant graphs and graph theoretic concepts. The concepts of protein graph construction are discussed. The manner in which graphs are analyzed and parameters relevant to protein structure are extracted, are explained. The structural and biological information derived from protein structures using these methods is presented.

Applications of graph theory in protein structure identification

There is a growing interest in the identification of proteins on the proteome wide scale. Among different kinds of protein structure identification methods, graph-theoretic methods are very sharp ones. Due to their lower costs, higher effectiveness and many other advantages, they have drawn more and more researchers' attention nowadays. Specifically, graph-theoretic methods have been widely used in homology identification, side-chain cluster identification, peptide sequencing and so on. This paper reviews several methods in solving protein structure identification problems using graph theory. We mainly introduce classical methods and mathematical models including homology modeling based on clique finding, identification of side-chain clusters in protein structures upon graph spectrum, and de novo peptide sequencing via tandem mass spectrometry using the spectrum graph model. In addition, concluding remarks and future priorities of each method are given.

Functional analysis of archaeal MBF1 by complementation studies in yeast

Background

Multiprotein-bridging factor 1 (MBF1) is a transcriptional co-activator that bridges a sequence-specific activator (basic-leucine zipper (bZIP) like proteins (e.g. Gcn4 in yeast) or steroid/nuclear-hormone receptor family (e.g. FTZ-F1 in insect)) and the TATA-box binding protein (TBP) in Eukaryotes. MBF1 is absent in Bacteria, but is well- conserved in Eukaryotes and Archaea and harbors a C-terminal Cro-like Helix Turn Helix (HTH) domain, which is the only highly conserved, classical HTH domain that is vertically inherited in all Eukaryotes and Archaea. The main structural difference between archaeal MBF1 (aMBF1) and eukaryotic MBF1 is the presence of a Zn ribbon motif in aMBF1. In addition MBF1 interacting activators are absent in the archaeal domain. To study the function and therefore the evolutionary conservation of MBF1 and its single domains complementation studies in yeast (mbf1Δ) as well as domain swap experiments between aMBF1 and yMbf1 were performed.

Results

In contrast to previous reports for eukaryotic MBF1 (i.e. Arabidopsis thaliana, insect and human) the two archaeal MBF1 orthologs, TMBF1 from the hyperthermophile Thermoproteus tenax and MMBF1 from the mesophile Methanosarcina mazei were not functional for complementation of an Saccharomyces cerevisiae mutant lacking Mbf1 (mbf1Δ). Of twelve chimeric proteins representing different combinations of the N-terminal, core domain, and the C-terminal extension from yeast and aMBF1, only the chimeric MBF1 comprising the yeast N-terminal and core domain fused to the archaeal C-terminal part was able to restore full wild-type activity of MBF1.

However, as reported previously for Bombyx mori, the C-terminal part of yeast Mbf1 was shown to be not essential for function. In addition phylogenetic analyses revealed a common distribution of MBF1 in all Archaea with available genome sequence, except of two of the three Thaumarchaeota; Cenarchaeum symbiosum A and Nitrosopumilus maritimus SCM1.

Conclusions

The absence of MBF1-interacting activators in the archaeal domain, the presence of a Zn ribbon motif in the divergent N-terminal domain of aMBF1 and the complementation experiments using archaeal- yeast chimeric proteins presented here suggests that archaeal MBF1 is not able to functionally interact with the transcription machinery and/or Gcn4 of S. cerevisiae. Based on modeling and structural prediction it is tempting to speculate that aMBF1 might act as a single regulator or non-essential transcription factor, which directly interacts with DNA via the positive charged linker or the basal transcription machinery via its Zn ribbon motif and the HTH domain. However, also alternative functions in ribosome biosynthesis and/or functionality have been discussed and therefore further experiments are required to unravel the function of MBF1 in Archaea.

Reviewers

This article was reviewed by William Martin, Patrick Forterre, John van der Oost and Fabian Blombach (nominated by Eugene V Koonin (United States)). For the full reviews, please go to the Reviewer's Reports section.

Modeling metal cation-phosphate interactions in nucleic acids: activated dissociation of Mg+, Al+, Cu+, and Zn+ complexes of triethyl phosphate

Threshold collision-induced dissociation techniques are employed to determine the activation energies (AEs) and bond dissociation energies (BDEs) of metal cation-triethyl phosphate complexes, M(+)(TEP), where M(+) = Mg(+), Al(+), Cu(+), and Zn(+). Activated dissociation resulting in loss of ethene, C(2)H(4), corresponds to the primary and lowest energy pathway for all four systems examined. Sequential loss of additional C(2)H(4) molecules and loss of the intact TEP ligand is also observed at elevated energies. Theoretical calculations at the B3LYP/6-31G* level of theory are used to determine the structures, vibrational frequencies, and rotational constants of neutral TEP and the M(+)(TEP) complexes, transition states, intermediates, and products of the activated dissociation of these complexes. Theoretical AEs and BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level using the B3LYP/6-31G* optimized geometries. The agreement between the calculated and measured AEs for elimination of C(2)H(4) is excellent for all four systems. In contrast, less satisfactory agreement between theory and experiment is found for the M(+)-TEP BDEs and may indicate limitations in the competitive model used to analyze these high energy dissociation pathways. The influence of the valence orbital occupation of the metal cation on the binding and activation propensities for elimination of ethene from TEP is examined. The binding of metal cations to TEP is compared to that of the nucleobases to assess the binding preferences of metal cations to nucleic acids.

Recognition of AT-rich DNA binding sites by the MogR repressor

The MogR transcriptional repressor of the intracellular pathogen Listeria monocytogenes recognizes AT-rich binding sites in promoters of flagellar genes to downregulate flagellar gene expression during infection. We describe here the 1.8 A resolution crystal structure of MogR bound to the recognition sequence 5' ATTTTTTAAAAAAAT 3' present within the flaA promoter region. Our structure shows that MogR binds as a dimer. Each half-site is recognized in the major groove by a helix-turn-helix motif and in the minor groove by a loop from the symmetry-related molecule, resulting in a "crossover" binding mode. This oversampling through minor groove interactions is important for specificity. The MogR binding site has structural features of A-tract DNA and is bent by approximately 52 degrees away from the dimer. The structure explains how MogR achieves binding specificity in the AT-rich genome of L. monocytogenes and explains the evolutionary conservation of A-tract sequence elements within promoter regions of MogR-regulated flagellar genes.

Studies on drug-DNA complexes, adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2), by phosphorus-31 nuclear magnetic resonance spectroscopy

The complexes of adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2) are studied by one- and two-dimensional (31)P nuclear magnetic resonance spectroscopy (NMR) at 500 MHz in the temperature range 275-328 K and as a function of drug to DNA ratio (0.0-2.0). The binding of drug to DNA is clearly evident in (31)P-(31)P exchange NOESY spectra that shows two sets of resonances in slow chemical exchange. The phosphate resonances at the intercalating steps, T1pG2/C1pG2 and C5pA6/C5pG6, shift downfield up to 1.7 ppm and that at the adjacent step shift downfield up to 0.7 ppm, whereas the central phosphate A3pT4 is relatively unaffected. The variations of chemical shift with drug to DNA ratio and temperature as well as linewidths are different in each of the two complexes. These observations reflect change in population of B(I)/B(II) conformation, stretching of backbone torsional angle zeta, and distortions in O-P-O bond angles that occur on binding of drug to DNA. To the best of our knowledge, there are no solution studies on 4'-epiadriamycin, a better tolerated drug, and binding of daunomycin or its analogue to d-(TGATCA)(2) hexamer sequence. The studies report the use of (31)P NMR as a tool to differentiate various complexes. The specific differences may well be the reasons that are responsible for different antitumor action of these drugs due to different binding ability and distortions in DNA.

Mutational analysis of an extracytoplasmic-function sigma factor to investigate its interactions with RNA polymerase and DNA

The extracytoplasmic-function (ECF) family of sigma factors comprises a large group of proteins required for synthesis of a wide variety of extracytoplasmic products by bacteria. Residues important for core RNA polymerase (RNAP) binding, DNA melting, and promoter recognition have been identified in conserved regions 2 and 4.2 of primary sigma factors. Seventeen residues in region 2 and eight residues in region 4.2 of an ECF sigma factor, PvdS from Pseudomonas aeruginosa, were selected for alanine-scanning mutagenesis on the basis of sequence alignments with other sigma factors. Fourteen of the mutations in region 2 had a significant effect on protein function in an in vivo assay. Four proteins with alterations in regions 2.1 and 2.2 were purified as His-tagged fusions, and all showed a reduced affinity for core RNAP in vitro, consistent with a role in core binding. Region 2.3 and 2.4 mutant proteins retained the ability to bind core RNAP, but four mutants had reduced or no ability to cause core RNA polymerase to bind promoter DNA in a band-shift assay, identifying residues important for DNA binding. All mutations in region 4.2 reduced the activity of PvdS in vivo. Two of the region 4.2 mutant proteins were purified, and each showed a reduced ability to cause core RNA polymerase to bind to promoter DNA. The results show that some residues in PvdS have functions equivalent to those of corresponding residues in primary sigma factors; however, they also show that several residues not shared with primary sigma factors contribute to protein function.

High-resolution crystal structure of the restriction-modification controller protein C.AhdI from Aeromonas hydrophila

Restriction-modification (R-M) systems serve to protect the host bacterium from invading bacteriophage. The multi-component system includes a methyltransferase, which recognizes and methylates a specific DNA sequence, and an endonuclease which recognises the same sequence and cleaves within or close to this site. The endonuclease will only cleave DNA that is unmethylated at the specific site, thus host DNA is protected while non-host DNA is cleaved. However, following DNA replication, expression of the endonuclease must be delayed until the host DNA is appropriately methylated. In many R-M systems, this regulation is achieved at the transcriptional level via the controller protein, or C-protein. We have solved the first X-ray structure of an R-M controller protein, C.AhdI, to 1.69 A resolution using selenomethionine MAD. C.AhdI is part of a Type IIH R-M system from the pathogen Aeromonas hydrophila. The structure reveals an all-alpha protein that contains a classical helix-turn-helix (HTH) domain and can be assigned to the Xre family of transcriptional regulators. Unlike its monomeric structural homologues, an extended helix generates an interface that results in dimerisation of the free protein. The dimer is electrostatically polarised and a positively charged surface corresponds to the position of the DNA recognition helices of the HTH domain. Comparison with the structure of the lambda cI ternary complex suggests that C.AhdI activates transcription through direct contact with the sigma70 subunit of RNA polymerase.

Monovalent Cations Regulate DNA Sequence Recognition by 434 Repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.

Recognition of 5'-YpG-3' sequences by coupled stacking/hydrogen bonding interactions with amino acid residues

The combined biochemical and structural study of hundreds of protein-DNA complexes has indicated that sequence-specific interactions are mediated by two mechanisms termed direct and indirect readout. Direct readout involves direct interactions between the protein and base-specific atoms exposed in the major and minor grooves of DNA. For indirect readout, the protein recognizes DNA by sensing conformational variations in the structure dependent on nucleotide sequence, typically through interactions with the phosphodiester backbone. Based on our recent structure of Ndt80 bound to DNA in conjunction with a search of the existing PDB database, we propose a new method of sequence-specific recognition that utilizes both direct and indirect readout. In this mode, a single amino acid side-chain recognizes two consecutive base-pairs. The 3'-base is recognized by canonical direct readout, while the 5'-base is recognized through a variation of indirect readout, whereby the conformational flexibility of the particular dinucleotide step, namely a 5'-pyrimidine-purine-3' step, facilitates its recognition by the amino acid via cation-pi interactions. In most cases, this mode of DNA recognition helps explain the sequence specificity of the protein for its target DNA.

Copper homeostasis in Enterococcus hirae

Copper is an essential component of life because of its convenient redox potential of 200-800 mV when bound to protein. Extensive insight into copper homeostasis has only emerged in the last decade and Enterococcus hirae has served as a paradigm for many aspects of the process. The cop operon of E. hirae regulates copper uptake, availability, and export. It consists of four genes that encode a repressor, CopY, a copper chaperone, CopZ, and two CPx-type copper ATPases, CopA and CopB. Most of these components have been conserved across the three evolutionary kingdoms. The four Cop proteins have been studied in vivo as well as in vitro and their function is understood in some detail.

The role of the minor groove substituents in indirect readout of DNA sequence by 434 repressor

The sequence of non-contacted bases at the center of the 434 repressor binding site affects the strength of the repressor-DNA complex by influencing the structure and flexibility of DNA (Koudelka, G. B., and Carlson, P. (1992) Nature 355, 89-91). We synthesized 434 repressor binding sites that differ in their central sequence base composition to test the importance of minor groove substituents and/or the number of base pair hydrogen bonds between these base pairs on DNA structure and strength of the repressor-DNA complex. We show here that the number of base pair H-bonds between the central bases apparently has no role in determining the relative affinity of a DNA site for repressor. Instead we find that the affinity of DNA for repressor depends on the absence or presence the N2-NH(2) group on the purine bases at the binding site center. The N2-NH(2) group on bases at the center of the 434 binding site appears to destabilize 434 repressor-DNA complexes by decreasing the intimacy of the specific repressor-DNA contacts, while increasing the reliance on protein contacts to the DNA phosphate backbone. Thus, the presence of an N2-NH(2) group on the purines at the center of a binding site globally alters the precise conformation of the protein-DNA interface.

Flexibility of alpha-helices: results of a statistical analysis of database protein structures

Alpha-helices stand out as common and relatively invariant secondary structural elements of proteins. However, alpha-helices are not rigid bodies and their deformations can be significant in protein function (e.g. coiled coils). To quantify the flexibility of alpha-helices we have performed a structural principal-component analysis of helices of different lengths from a representative set of protein folds in the Protein Data Bank. We find three dominant modes of flexibility: two degenerate bend modes and one twist mode. The data are consistent with independent Gaussian distributions for each mode. The mode eigenvalues, which measure flexibility, follow simple scaling forms as a function of helix length. The dominant bend and twist modes and their harmonics are reproduced by a simple spring model, which incorporates hydrogen-bonding and excluded volume. As an application, we examine the amount of bend and twist in helices making up all coiled-coil proteins in SCOP. Incorporation of alpha-helix flexibility into structure refinement and design is discussed.

Retroevolution of lambda Cro toward a stable monomer

The Cro protein from bacteriophage lambda has a dimeric alpha+beta fold that evolved from an ancestral all-alpha monomer. The sequence mutations responsible for this dramatic structural evolution are unknown. Here we use analysis of sequence alignments to show that Ala-33, a small side chain in the hydrophobic "ball-and-socket" dimer interface of lambda Cro, was a much larger tryptophan side chain at a previous point in evolution. The retroevolutionary lambda Cro-A33W mutant shows a 10-fold reduction in dimerization affinity relative to the wild type as well as a large increase in monomer thermal stability (Delta T(m) > 10 degrees C), apparently due to partial filling of the hydrophobic socket from within the same monomer. An additional mutation in the dimer interface, F58D, almost completely abolishes detectable dimerization while maintaining the high monomer stability. The secondary structure content of the monomerized versions of lambda Cro is similar to that of the wild-type protein, and the tertiary structure of the monomer appears relatively well defined. These results (i) support a model in which the ball-and-socket dimer interface of lambda Cro was created by altered volume mutations within a limited branch of the Cro lineage and (ii) suggest the possibility that the evolution of the alpha+beta dimer from an all-alpha monomer proceeded through an alpha+beta monomer intermediate.

Operator recognition by the phage 434 cI repressor: MD simulations of free and bound 50-bp DNA reveal important differences between the OR1 and OR2 sites

Using molecular dynamics simulations in explicit solvent, we investigated the behavior of a 50-bp DNA sequence containing the 434 bacteriophage operators OR1 and OR2 separated by an 8-bp spacer. Two simulations of 1 ns each were carried out, with DNA alone and with DNA complexed to dimers of the R1-69 DNA binding domain of the phage 434 cI repressor protein at the OR1 and OR2 sites. Strong correlations among average structural parameters are observed between our simulations and available experimental data for the bound OR1/OR2 subsites. In the free state, some differences appear between the three relevant fragments (OR1, the spacer, and OR2). Unbound OR1 exhibits a large, shallow major groove into which the base atoms protrude and is also bent toward the major groove. This structure is maintained because structural fluctuations are weak. Unbound OR2 resembles canonical B-DNA although the structural parameters show greater fluctuations, essentially due to a malleable step (the innermost CpA/TpG), absent in OR1. Complexation with the proteins slightly alters the base positions but strongly modifies the sugar and backbone motions. The most crucial repressor effects are changes in the flexibility of the OR1/OR2 sites. Structural fluctuations are enhanced for OR1, conferring a favorable energetic contribution to the OR1 binding, whereas they are reduced for OR2. Therefore, both structural and dynamic properties of DNA suggest OR1 is the most attractive site for the repressor, which may explain the different binding association constants observed for the OR1 and OR2 sites. Finally, we also investigated the impact of the protein on the DNA backbone dynamics and find that direct or indirect interactions facilitate the DNA structural variations required for achieving complementarity with the protein.

Binding sites of different geometries for the 16-3 phage repressor

Prokaryotic repressor-operator systems provide exemplars for the sequence-specific interactions between DNA and protein. The crucial atomic contacts of the two macromolecules are attained in a compact, geometrically defined structure of the DNA-protein complex. The pitch of the DNA interface seems an especially sensitive part of this architecture because changes in its length introduce new spacing and rotational relations in one step. We discovered a natural system that may serve as a model for investigating this problem: the repressor of the 16-3 phage of Rhizobium meliloti (helix-turn-helix class protein) possesses inherent ability to accommodate to various DNA twistings. It binds the cognate operators, which are 5'-ACAA-4 bp-TTGT-3' (O(L)) and 5'-ACAA-6 bp-TTGT-3' (O(R)) and thus differ 2 bp in length, and consequently the two half-sites will be rotated with respect to each other by 72 degrees in the idealized B-DNA (64 degrees by dinucleotide steps calculations). Furthermore, a synthetic intermediate (DNA sequence) 5'-ACAA-5 bp-TTGT-3' (O5) also binds specifically the repressor. The natural operators and bound repressors can form higher order DNA-protein complexes and perform efficient repression, whereas the synthetic operator-repressor complex cannot do either. The natural operators are bent when complexed with the repressor, whereas the O5 operator does not show bending in electrophoretic mobility assay. Possible structures of the complexes are discussed.

BphS, a key transcriptional regulator of bph genes involved in polychlorinated biphenyl/biphenyl degradation in Pseudomonas sp. KKS102

The bph genes in Pseudomonas sp. KKS102, which are involved in the degradation of polychlorinated biphenyl/biphenyl, are induced in the presence of biphenyl. In this study our goal was to understand the regulatory mechanisms involved in the inducible expression. The bph genes (bphEGF(orf4)A1A2A3BCD(orf1)A4R) constitute an operon, and its expression is strongly dependent on the pE promoter located upstream of the bphE gene. A bphS gene, whose deduced amino acid sequence showed homology with the GntR family transcriptional repressors, was identified at the upstream region of the bphE gene. Disruption of the bphS gene resulted in constitutive expression of bph genes, suggesting that the bphS gene product negatively regulated the pE promoter. The gel retardation and DNase footprinting analyses demonstrated specific binding of BphS to the pE promoter region and identified four BphS binding sites that were located within and immediately downstream of the -10 box of the pE promoter. The four binding sites were functional in repression because their respective elimination resulted in derepression of the pE promoter. The binding of BphS was abolished in the presence of 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid, an intermediate compound in the biphenyl degradation pathway. We concluded that the negative regulator BphS plays a central role in the regulation of bph gene expression through its action at the pE promoter.

Selection and design of high affinity DNA ligands for mutant single-chain derivatives of the bacteriophage 434 repressor

Single-chain repressor RR(Tres) is a derivative of bacteriophage 434 repressor, which contains covalently dimerized DNA-binding domains (amino acids 1-69) of the phage 434 repressor. In this single-chain molecule, the wild type domain R is connected to the mutant domain R(TRES) by a recombinant linker in a head-to-tail arrangement. The DNA-contacting amino acids of R(TRES) at the -1, 1, 2, and 5 positions of the a3 helix are T, R, E, S respectively. By using a randomized DNA pool containing the central sequence -CATACAAGAAAGNNNNNNTTT-, a cyclic,in vitro DNA-binding site selection was performed. The selected population was cloned and the individual members were characterized by determining their binding affinities to RR(Tres) The results showed that the optimal operators contained the TTAC or TTCC sequences in the underlined positions as above, and that the Kd values were in the 1 x 10(-12) mol/L-1 x 10(11) mol/L concentration range. Since the affinity of the natural 434 repressor to its natural operator sites is in the 1 x 10(-9) mol/L range, the observed binding affinity increase is remarkable. It was also found that binding affinity was strongly affected by the flanking bases of the optimal tetramer binding sites, especially by the base at the 5' position. We constructed a new homodimeric single-chain repressor R(TRES)R(TRES) and its DNA-binding specificity was tested by using a series of new operators designed according to the recognition properties previously determined for the R(TREs) domain. These operators containing the consensus sequenceGTAAGAAARNTTACN orGGAAGAAARNTTCCN (R is A or G) were recognized by R(TRES)R(TRES) specifically, and with high binding affinity. Thus, by using a combination of random selection and rational design principles, we have discovered novel, high affinity protein-DNA interactions with new specificity. This method can potentially be used to obtain new binding specificity for other DNA-binding proteins.

The new architectonics: an invitation to structural biology

The philosophy of art might offer an epistemological basis for talking about the complexity of biological molecules in a meaningful way. The analysis of artistic compositions requires the resolution of intrinsic tensions between disparate sensory categories-color, line and form-not unlike those encountered in looking at the surfaces of protein molecules, where charge, polarity, hydrophobicity, and shape compete for our attentions. Complex living systems exhibit behaviors such as contraction waves moving along muscle fibers, or shivers passing through the growth cones of migrating neurons, that are easy to describe with common words, but difficult to explain in terms of the language of chemistry. The problem follows from a lack of everyday experience with processes that move towards equilibrium by switching between crystalline order and chain-like disorder, a commonplace occurrence in the submicroscopic world of proteins. Since most of what is understood about protein function comes from studies of isolated macromolecules in solution, a serious gap exists between what we know and what we would like to know about organized biological systems. Closing this gap can be achieved by recognizing that protein molecules reside in gradients of Gibbs free energy, where local forces and movements can be large compared with Brownian motion. Architectonics, a term borrowed from the philosophical literature, symbolizes the eventual union of the structure of theories-how our minds construct the world-with the theory of structures-or how stability is maintained in the chaotic world of microsystems.

DNA sequence requirements for the activation of 434 P(RM) transcription by 434 repressor

A dimer of the 434 repressor bound at O(R)2 activated transcription initiation from P(RM) by contacting RNA polymerase. Although DNA-binding site mutations at either end of O(R)2 decreased the ability of the repressor to activate P(RM) transcription, mutations proximal to the promoter had a greater effect on transcription activation. Orienting a repressor subunit bearing the altered specificity Gln-28 --> Ala mutation to the halfsite of O(R)2 proximal to the P(RM) promoter decreased the repressor's ability to activate transcription initiation at 434 P(RM) to a much greater extent than if this subunit was placed in the O(R)2 half-site distal to P(RM). In addition to showing that the downstream (promoter proximal) subunit of the O(R)2-bound 434 repressor functions in activating 434 P(RM), the results indicated that DNA sequence-dependent conformational changes alter the efficiency with which the repressor activates P(RM) transcription. These unexpected findings highlight the importance of the structure of the repressor-DNA interface in activating transcription from P(RM).

Inhibition of HIV-1 enhancer-controlled transcription by artificial enhancer-binding peptides derived from bacteriophage 434 repressor

An artificial HIV-1 enhancer-binding 42-residue peptide (R42) that had been derived from bacteriophage 434 repressor inhibited the cell-free in vitro transcription of HIV-1 enhancer-containing plasmids [Hehlgans, T., Stolz, M., Klauser, S., Cui, T., Salgam, P., Brenz Verca, S., Widmann, M., Leiser, A., Städler, K. & Gutte, B. (1993) FEBS Lett. 315, 51-55; Caderas, G. (1997) PhD Thesis, University of Zürich]. Here we show that, after N-terminal extension of R42 with a viral nuclear localization signal, the resulting nucR42 peptide was active in intact cells. NucR42 could be detected immunologically in nuclear extracts and produced a 60-70% reduction of the rate of transcription of an HIV-1 enhancer-carrying plasmid in COS-1 cells that had been cotransfected with the HIV enhancer plasmid, an expression plasmid for nucR42, and a control. NucR42 was also synthesized chemically and the synthetic product characterized by HPLC, mass spectrometry, and quantitative amino acid analysis. Band shift, footprint, and in vitro transcription assays in the presence of exogenous NF-kappaBp50 indicated that the binding sites of nucR42 and NF-kappaB on the HIV enhancers overlapped and that a relatively small excess of nucR42 sufficed to displace NF-kappaBp50. Band shift and in vitro transcription experiments showed also that exchange of the 434 repressor-derived nine-residue recognition helix of nucR42 for four glycines abolished the HIV enhancer binding specificity whereas leucine zipper- or retro-leucine zipper-mediated dimerization of R42 analogues increased it suggesting the potential application of such dimeric HIV enhancer-binding peptides as intracellular inhibitors of HIV replication.

Identification of side-chain clusters in protein structures by a graph spectral method

This paper presents a novel method to detect side-chain clusters in protein three-dimensional structures using a graph spectral approach. Protein side-chain interactions are represented by a labeled graph in which the nodes of the graph represent the Cbeta atoms and the edges represent the distance between the Cbeta atoms. The distance information and the non-bonded connectivity of the residues are represented in the form of a matrix called the Laplacian matrix. The constructed matrix is diagonalized and clustering information is obtained from the vector components associated with the second lowest eigenvalue and cluster centers are obtained from the vector components associated with the top eigenvalues. The method uses global information for clustering and a single numeric computation is required to detect clusters of interest. The approach has been adopted here to detect a variety of side-chain clusters and identify the residue which makes the largest number of interactions among the residues forming the cluster (cluster centers). Detecting such clusters and cluster centers are important from a protein structure and folding point of view. The crucial residues which are important in the folding pathway as determined by PhiF values (which is a measure of the effect of a mutation on the stability of the transition state of folding) as obtained from protein engineering methods, can be identified from the vector components corresponding to the top eigenvalues. Expanded clusters are detected near the active and binding site of the protein, supporting the nucleation condensation hypothesis for folding. The method is also shown to detect domains in protein structures and conserved side-chain clusters in topologically similar proteins.

Molecular dynamics simulation in solvent of the bacteriophage 434 cI repressor protein DNA binding domain amino acids (R1-69) in complex with its cognate operator (OR1) DNA sequence

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the 434 cI Repressor protein DNA binding domain (DBD) amino acids (R1-69) and DNA of operator (OR1) and its flanks consisting of 28 nucleotide base pairs. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the 434 cI repressor DNA recognition helix 3 formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the OR1 DNA major groove halfsites and flanking regions. In addition, hydrophilic amino acids within the loop between helix 3 and helix 4 have strong electrostatic attraction to codon-anticodon nucleotides located within the central nucleotides of the minor groove between the OR1 major groove halfsites. These interactions together induced significant structural changes in the operator DNA manifested by overtwisting of the central nucleotide base pairs and narrowing of the minor groove between the DNA major groove halfsites. Finally, these findings offer a code for site specific DNA recognition by the 434 cI repressor protein.