Tandem binding in crystals of a trp repressor/operator half-site complex

Developing Community Resources for Nucleic Acid Structures

In this review, we describe the creation of the Nucleic Acid Database (NDB) at Rutgers University and how it became a testbed for the current infrastructure of the RCSB Protein Data Bank. We describe some of the special features of the NDB and how it has been used to enable research. Plans for the next phase as the Nucleic Acid Knowledgebase (NAKB) are summarized.

Synthetic biology approaches to dissecting linear motor protein function: towards the design and synthesis of artificial autonomous protein walkers

Molecular motors and machines are essential for all cellular processes that together enable life. Built from proteins with a wide range of properties, functionalities and performance characteristics, biological motors perform complex tasks and can transduce chemical energy into mechanical work more efficiently than human-made combustion engines. Sophisticated studies of biological protein motors have provided many structural and biophysical insights and enabled the development of models for motor function. However, from the study of highly evolved, biological motors, it remains difficult to discern detailed mechanisms, for example, about the relative role of different force generation mechanisms, or how information is communicated across a protein to achieve the necessary coordination. A promising, complementary approach to answering these questions is to build synthetic protein motors from the bottom up. Indeed, much effort has been invested in functional protein design, but so far, the "holy grail" of designing and building a functional synthetic protein motor has not been realized. Here, we review the progress made to date, and we put forward a roadmap for achieving the aim of constructing the first artificial, autonomously running protein motor. Specifically, we propose to break down the task into (i) enzymatic control of track binding, (ii) the engineering of asymmetry and (iii) the engineering of allosteric control for internal communication. We also propose specific approaches for solving each of these challenges.

Conserved Dynamic Mechanism of Allosteric Response to L-arg in Divergent Bacterial Arginine Repressors

Hexameric arginine repressor, ArgR, is the feedback regulator of bacterial L-arginine regulons, and sensor of L-arg that controls transcription of genes for its synthesis and catabolism. Although ArgR function, as well as its secondary, tertiary, and quaternary structures, is essentially the same in E. coli and B. subtilis, the two proteins differ significantly in sequence, including residues implicated in the response to L-arg. Molecular dynamics simulations are used here to evaluate the behavior of intact B. subtilis ArgR with and without L-arg, and are compared with prior MD results for a domain fragment of E. coli ArgR. Relative to its crystal structure, B. subtilis ArgR in absence of L-arg undergoes a large-scale rotational shift of its trimeric subassemblies that is very similar to that observed in the E. coli protein, but the residues driving rotation have distinct secondary and tertiary structural locations, and a key residue that drives rotation in E. coli is missing in B. subtilis. The similarity of trimer rotation despite different driving residues suggests that a rotational shift between trimers is integral to ArgR function. This conclusion is supported by phylogenetic analysis of distant ArgR homologs reported here that indicates at least three major groups characterized by distinct sequence motifs but predicted to undergo a common rotational transition. The dynamic consequences of L-arg binding for transcriptional activation of intact ArgR are evaluated here for the first time in two-microsecond simulations of B. subtilis ArgR. L-arg binding to intact B. subtilis ArgR causes a significant further shift in the angle of rotation between trimers that causes the N-terminal DNA-binding domains lose their interactions with the C-terminal domains, and is likely the first step toward adopting DNA-binding-competent conformations. The results aid interpretation of crystal structures of ArgR and ArgR-DNA complexes.

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.

Mechanistic insights into metal ion activation and operator recognition by the ferric uptake regulator

Ferric uptake regulator (Fur) plays a key role in the iron homeostasis of prokaryotes, such as bacterial pathogens, but the molecular mechanisms and structural basis of Fur-DNA binding remain incompletely understood. Here, we report high-resolution structures of Magnetospirillum gryphiswaldense MSR-1 Fur in four different states: apo-Fur, holo-Fur, the Fur-feoAB1 operator complex and the Fur-Pseudomonas aeruginosa Fur box complex. Apo-Fur is a transition metal ion-independent dimer whose binding induces profound conformational changes and confers DNA-binding ability. Structural characterization, mutagenesis, biochemistry and in vivo data reveal that Fur recognizes DNA by using a combination of base readout through direct contacts in the major groove and shape readout through recognition of the minor-groove electrostatic potential by lysine. The resulting conformational plasticity enables Fur binding to diverse substrates. Our results provide insights into metal ion activation and substrate recognition by Fur that suggest pathways to engineer magnetotactic bacteria and antipathogenic drugs.

Concentration-dependent effects on fully hydrated DNA at terahertz frequencies

Using terahertz time-domain spectroscopy (THz-TDS), the frequency-dependent dielectric constant of deoxyribonucleic acid (DNA) in solution was measured. The response of the buffer solution is dominated by two Debye modes in this frequency range, and, from an analysis of the concentration dependence, the presence of the DNA increases the main relaxation time and dielectric constant. This reflects the fact that the water in the hydration layer is more tightly bound under the influence of the DNA molecule in comparison to bulk water. This dynamical slowing down with increasing DNA concentration is similar to what is observed with purine nucleotides, but opposite to the behavior of pyrimidine nucleotides. In addition, a suspension model was used with the concentration-dependent data to isolate the dielectric response of the hydrated DNA molecule. The data for the hydrated DNA molecule is still dominated by a Debye response. It is also possible to determine the thickness of the hydration layer, and the DNA molecule influences the surrounding water out to 16 or 17 Å, which corresponds to about six effective hydration layers.

Regulatory mechanism of the light-activable allosteric switch LOV-TAP for the control of DNA binding: a computer simulation study

The spatio-temporal control of gene expression is fundamental to elucidate cell proliferation and deregulation phenomena in living systems. Novel approaches based on light-sensitive multiprotein complexes have recently been devised, showing promising perspectives for the noninvasive and reversible modulation of the DNA-transcriptional activity in vivo. This has lately been demonstrated in a striking way through the generation of the artificial protein construct light-oxygen-voltage (LOV)-tryptophan-activated protein (TAP), in which the LOV-2-Jα photoswitch of phototropin1 from Avena sativa (AsLOV2-Jα) has been ligated to the tryptophan-repressor (TrpR) protein from Escherichia coli. Although tremendous progress has been achieved on the generation of such protein constructs, a detailed understanding of their functioning as opto-genetical tools is still in its infancy. Here, we elucidate the early stages of the light-induced regulatory mechanism of LOV-TAP at the molecular level, using the noninvasive molecular dynamics simulation technique. More specifically, we find that Cys450-FMN-adduct formation in the AsLOV2-Jα-binding pocket after photoexcitation induces the cleavage of the peripheral Jα-helix from the LOV core, causing a change of its polarity and electrostatic attraction of the photoswitch onto the DNA surface. This goes along with the flexibilization through unfolding of a hairpin-like helix-loop-helix region interlinking the AsLOV2-Jα- and TrpR-domains, ultimately enabling the condensation of LOV-TAP onto the DNA surface. By contrast, in the dark state the AsLOV2-Jα photoswitch remains inactive and exerts a repulsive electrostatic force on the DNA surface. This leads to a distortion of the hairpin region, which finally relieves its tension by causing the disruption of LOV-TAP from the DNA.

Environment-dependent long-range structural distortion in a temperature-sensitive point mutant

Extensive environment-dependent rearrangement of the helix-turn-helix DNA recognition region and adjacent L-tryptophan binding pocket is reported in the crystal structure of dimeric E. coli trp aporepressor with point mutation Leu75Phe. In one of two subunits, the eight residues immediately C-terminal to the mutation are shifted forward in helical register by three positions, and the five following residues form an extrahelical loop accommodating the register shift. In contrast, the second subunit has wildtype-like conformation, as do both subunits in an isomorphous wildtype control structure. Treated together as an ensemble pair, the distorted and wildtype-like conformations of the mutant apoprotein agree more fully than either conformation alone with previously reported NOE measurements, and account more completely for its diverse biochemical and biophysical properties. The register-shifted segment Ile79-Ala80-Thr81-Ile82-Thr83 is helical in both conformations despite low helical propensity, suggesting an important structural role for the steric constraints imposed by β-branched residues in helical conformation.

Bacillus subtilis CodY operators contain overlapping CodY binding sites

CodY is a global transcriptional regulator that is activated by branched-chain amino acids. A palindromic 15-bp sequence motif, AATTTTCNGAAAATT, is associated with CodY DNA binding. A gel mobility shift assay was used to examine the effect of pH on the binding of Bacillus subtilis CodY to the hutPp and ureAp(3) promoters. CodY at pH 6.0 has higher affinity for DNA, more enhanced activation by isoleucine, and a lower propensity for nonspecific DNA binding than CodY at pH 8.0. DNase I footprinting was used to identify the CodY-protected regions in the hutPp and ureAp(3) promoters. The CodY-protected sequences for both promoters were found to contain multiple copies of the 15-bp motif with 6-bp overlaps. Mutational analysis of the hutPp regulatory region revealed that two overlapping sequence motifs were required for CodY-mediated regulation. The presence of overlapping sequence motifs in the regulatory regions of many B. subtilis CodY-regulated genes suggests that CodY binds to native operators that contain overlapping binding sites.

Dielectric properties of fully hydrated nucleotides in the terahertz frequency range

We use terahertz time domain spectroscopy (THz-TDS) to determine the complex frequency-dependent dielectric response of all four nucleotides at different dilute concentrations. In addition, the suspension model's ability to extract the dielectric response of just the nucleotide with the hydration shell epsilon(b) excluding the dielectric information relating to the bulk will be verified. The suspension model enables us to make the determination that the nucleotides have influences on the water molecules out to the fourth hydration shell. We use a two Debye relaxation fit model for water, all concentrations and all epsilon(b) values. We observed how the nucleotides affect the relaxation parameters in relation to that of pure bulk water. With this information, we notice a transition between purines and pyrimidines, where one is a hydrogen-bond network structure building type material with a low concentration increment and the other is a structure breaking type material with a low concentration decrement. Due to conductivity measurements, we determine that kinetic depolarization is a negligible affect compared to that of dielectric saturation, which we find to dominate where a decrement is found.

Functional rules for lac repressor-operator associations and implications for protein-DNA interactions

The Lac repressor has been used as a tool to understand protein-DNA recognition for many years. Recent experiments have demonstrated the ability of the Lac repressor to control gene expression in various eukaryotic systems, making the quest for an arsenal of protein-DNA binding partners desirable for potential therapeutic applications. Here, we present the results of the most exhaustive screen of Lac repressor-DNA binding partners to date, resulting in the elucidation of functional rules for Lac-DNA binding. Even within the confines of a single protein-DNA scaffold, modes of binding of different protein-DNA partners are sufficiently diverse so as to prevent elucidation of generalized rules for recognition for a single protein, much less an entire protein family.

A structural-alphabet-based strategy for finding structural motifs across protein families

Proteins with insignificant sequence and overall structure similarity may still share locally conserved contiguous structural segments; i.e. structural/3D motifs. Most methods for finding 3D motifs require a known motif to search for other similar structures or functionally/structurally crucial residues. Here, without requiring a query motif or essential residues, a fully automated method for discovering 3D motifs of various sizes across protein families with different folds based on a 16-letter structural alphabet is presented. It was applied to structurally non-redundant proteins bound to DNA, RNA, obligate/non-obligate proteins as well as free DNA-binding proteins (DBPs) and proteins with known structures but unknown function. Its usefulness was illustrated by analyzing the 3D motifs found in DBPs. A non-specific motif was found with a ‘corner’ architecture that confers a stable scaffold and enables diverse interactions, making it suitable for binding not only DNA but also RNA and proteins. Furthermore, DNA-specific motifs present ‘only’ in DBPs were discovered. The motifs found can provide useful guidelines in detecting binding sites and computational protein redesign.

Overlapping repressor binding sites regulate expression of the Methanococcus maripaludis glnK(1) operon

The euryarchaeal transcriptional repressor NrpR regulates a variety of nitrogen assimilation genes by 2-oxoglutarate-reversible binding to conserved palindromic operators. The number and positioning of these operators varies among promoter regions of regulated genes, suggesting NrpR can bind in different patterns. Particularly intriguing is the contrast between the nif and glnK(1) promoter regions of Methanococcus maripaludis, where two operators are present but with different configurations. Here we study NrpR binding and regulation at the glnK(1) promoter, where the two operator sequences overlap and occur on opposite faces of the double helix. We find that both operators function in binding, with a dimer of NrpR binding simultaneously to each overlapping operator. We show in vivo that the first operator plays a primary role in regulation and the second operator plays an enhancing role. This is the first demonstration of overlapping operators functioning in Archaea.

The Tumbleweed: towards a synthetic proteinmotor

Biomolecular motors have inspired the design and construction of artificial nanoscale motors and machines based on nucleic acids, small molecules, and inorganic nanostructures. However, the high degree of sophistication and efficiency of biomolecular motors, as well as their specific biological function, derives from the complexity afforded by protein building blocks. Here, we discuss a novel bottom-up approach to understanding biological motors by considering the construction of synthetic protein motors. Specifically, we present a design for a synthetic protein motor that moves along a linear track, dubbed the "Tumbleweed." This concept uses three discrete ligand-dependent DNA-binding domains to perform cyclically ligand-gated, rectified diffusion along a synthesized DNA molecule. Here we describe how de novo peptide design and molecular biology could be used to produce the Tumbleweed, and we explore the fundamental motor operation of such a design using numerical simulations. The construction of this and more sophisticated protein motors is an exciting challenge that is likely to enhance our understanding of the structure-function relationship in biological motors.

Signatures of protein-DNA recognition in free DNA binding sites

One obstacle to achieving complete understanding of the principles underlying sequence-dependent recognition of DNA is the paucity of structural data for DNA recognition sequences in their free (unbound) state. Here, we carried out crystallization screening of 50 DNA duplexes containing cognate protein binding sites and obtained new crystal structures of free DNA binding sites for three distinct modes of DNA recognition: anti-parallel beta strands (MetR), helix-turn-helix motif + hinge helices (PurR), and zinc fingers (Zif268). Structural changes between free and protein-bound DNA are manifested differently in each case. The new DNA structures reveal that distinctive sequence-dependent DNA geometry dominates recognition by MetR, protein-induced bending of DNA dictates recognition by PurR, and deformability of DNA along the A-B continuum is important in recognition by Zif268. Together, our findings show that crystal structures of free DNA binding sites provide new information about the nature of protein-DNA interactions and thus lend insights towards a structural code for DNA recognition.

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.

Light-activated DNA binding in a designed allosteric protein

An understanding of how allostery, the conformational coupling of distant functional sites, arises in highly evolvable systems is of considerable interest in areas ranging from cell biology to protein design and signaling networks. We reasoned that the rigidity and defined geometry of an alpha-helical domain linker would make it effective as a conduit for allosteric signals. To test this idea, we rationally designed 12 fusions between the naturally photoactive LOV2 domain from Avena sativa phototropin 1 and the Escherichia coli trp repressor. When illuminated, one of the fusions selectively binds operator DNA and protects it from nuclease digestion. The ready success of our rational design strategy suggests that the helical "allosteric lever arm" is a general scheme for coupling the function of two proteins.

Transcript analysis reveals an extended regulon and the importance of protein-protein co-operativity for the Escherichia coli methionine repressor

We have used DNA arrays to investigate the effects of knocking out the methionine repressor gene, metJ, on the Escherichia coli transcriptome. We assayed the effects in the knockout strain of supplying wild-type or mutant MetJ repressors from an expression plasmid, thus establishing a rapid assay for in vivo effects of mutations characterized previously in vitro. Repression is largely restricted to known genes involved in the biosynthesis and uptake of methionine. However, we identified a number of additional genes that are significantly up-regulated in the absence of repressor. Sequence analysis of the 5' promoter regions of these genes identified plausible matches to met-box sequences for three of these, and subsequent electrophoretic mobility-shift assay analysis showed that for two such loci their repressor affinity is higher than or comparable with the known metB operator, suggesting that they are directly regulated. This can be rationalized for one of the loci, folE, by the metabolic role of its encoded enzyme; however, the links to the other regulated loci are unclear, suggesting both an extension to the known met regulon and additional complexity to the role of the repressor. The plasmid gene replacement system has been used to examine the importance of protein-protein co-operativity in operator saturation using the structurally characterized mutant repressor, Q44K. In vivo, there are detectable reductions in the levels of regulation observed, demonstrating the importance of balancing protein-protein and protein-DNA affinity.

Coarse-grained model of entropic allostery

Many signaling functions in molecular biology require proteins to bind to substrates such as DNA in response to environmental signals such as the simultaneous binding to a small molecule. Examples are repressor proteins which may transmit information via a conformational change in response to the ligand binding. An alternative entropic mechanism of "allostery" suggests that the inducer ligand changes the intramolecular vibrational entropy, not just the mean static structure. We present a quantitative, coarse-grained model of entropic allostery, which suggests design rules for internal cohesive potentials in proteins employing this effect. It also addresses the issue of how the signal information to bind or unbind is transmitted through the protein. The model may be applicable to a wide range of repressors and also to signaling in trans-membrane proteins.

E. coli trp repressor forms a domain-swapped array in aqueous alcohol

The E. coli trp repressor (trpR) homodimer recognizes its palindromic DNA binding site through a pair of flexible helix-turn-helix (HTH) motifs displayed on an intertwined helical core. Flexible N-terminal arms mediate association between dimers bound to tandem DNA sites. The 2.5 A X-ray structure of trpR crystallized in 30% (v/v) isopropanol reveals a substantial conformational rearrangement of HTH motifs and N-terminal arms, with the protein appearing in the unusual form of an ordered 3D domain-swapped supramolecular array. Small angle X-ray scattering measurements show that the self-association properties of trpR in solution are fundamentally altered by isopropanol.

Evidence for assembly-dependent folding of protein and RNA in an icosahedral virus

Ordered nucleic acid in an icosahedral virus was first visualized in the X-ray structure of the Picorna-like plant virus, Bean pod mottle virus (BPMV). Virus particles containing the 3500 nucleotide segment of the BPMV bipartite RNA genome (middle component) had nearly 20% of the genome ordered. Here we report the refined structures of the middle component, bottom component (particles containing the 5800 nucleotide segment of the genome), and top component (empty particles of BPMV capsid protein). The bottom component particles contain ordered RNA in the same location as middle component. Although the ordered RNA density in both nucleoprotein particles is the average of the contents of 60 icosahedral asymmetric units, both nucleoprotein components show that the base density for the first two nucleotides is predominantly purine, while the next five appear to be predominantly pyrimidine. The empty capsid demonstrates that RNA dictates the order of the N-terminal 19 residues of the large subunit because these residues are invisible in the top component.

Amino acid discrimination by a class I aminoacyl-tRNA synthetase specified by negative determinants

The 2.5 A crystal structure of Escherichia coli glutaminyl-tRNA synthetase in a quaternary complex with tRNA(Gln), an ATP analog and glutamate reveals that the non-cognate amino acid adopts a distinct binding mode within the active site cleft. In contrast to the binding of cognate glutamine, one oxygen of the charged glutamate carboxylate group makes a direct ion-pair interaction with the strictly conserved Arg30 residue located in the first half of the dinucleotide fold domain. The nucleophilic alpha-carboxylate moiety of glutamate is mispositioned with respect to both the ATP alpha-phosphate and terminal tRNA ribose groups, suggesting that a component of amino acid discrimination resides at the catalytic step of the reaction. Further, the other side-chain carboxylate oxygen of glutamate is found in a position identical to that previously proposed to be occupied by the NH(2) group of the cognate glutamine substrate. At this position, the glutamate oxygen accepts hydrogen bonds from the hydroxyl moiety of Tyr211 and a water molecule. These findings demonstrate that amino acid specificity by GlnRS cannot arise from hydrogen bonds donated by the cognate glutamine amide to these same moieties, as previously suggested. Instead, Arg30 functions as a negative determinant to drive binding of non-cognate glutamate into a non-productive orientation. The poorly differentiated cognate amino acid-binding site in GlnRS may be a consequence of the late emergence of this enzyme from the eukaryotic lineage of glutamyl-tRNA synthetases.

Correlation between normal modes in the 20-200 cm-1 frequency range and localized torsion motions related to certain collective motions in proteins

In certain biologically relevant collective motions, such as protein domain motions and sub-domain motions, large amplitude movements are localized in one or a few flexible regions consisting of a small number of residues. This paper explores the possible use of normal mode analysis in probing localized vibrational torsion motions in these flexible regions that may be related to certain collective motions. The normal modes of 10 structures of five proteins in different conformation (TRP repressor, calmodulin, calbindin D(9k), HIV-1 protease and troponin C), known to have shear or hinge domain or sub-domain motion, respectively, are analyzed. Our study identifies, for each structure, unique normal modes in the 20-200 cm-1 frequency range, whose corresponding motions are primarily concentrated in the region where large amplitude torsion movements of a known domain or sub-domain motion occur. This suggests possible correlation between normal modes at 20-200 cm-1 frequency range and initial fluctuational motions leading to localized collective motions in proteins, and thus the potential application of normal mode analysis in facilitating the study of biologically important localized motions in biomolecules.

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.

Experimental and theoretical high pressure strategies for investigating protein-nucleic acid assemblies

A method was developed to investigate the stability of protein-nucleic acid complexes using hydrostatic pressure during electrophoretic gel mobility shift analysis. The initial system probed by this technique was the well-characterized cognate BamHI-DNA complex. Band shift analysis at several elevated pressures found the equilibrium dissociation (K(d)) constant to be dependent on pressure, which allowed the volume change of dissociation (deltaV) to be calculated. In order to describe the effects of pressure on the specific BamHI-DNA complex at the molecular level, molecular dynamics simulations at both ambient and elevated pressure was performed. Comparison of the simulation trajectories identified several individual BamHI-DNA contacts that are disrupted due to pressure. The disruption of these contacts can be attributed to an observed pressure-induced increase in hydration at the protein-DNA interface during the elevated pressure simulation.

Advancing our knowledge in biochemistry, genetics, and microbiology through studies on tryptophan metabolism

I was fortunate to practice science during the last half of the previous century, when many basic biological and biochemical concepts could be experimentally addressed for the first time. My introduction to research involved isolating and identifying intermediates in the niacin biosynthetic pathway. These studies were followed by investigations focused on determining the properties of genes and enzymes essential to metabolism and examining how they were alterable by mutation. The most challenging problem I initially attacked was establishing the colinear relationship between gene and protein. Subsequent research emphasized identification and characterization of regulatory mechanisms that microorganisms use to control gene expression. An elaborate regulatory strategy, transcription attenuation, was discovered that is often based on selection between alternative RNA structures. Throughout my career I enjoyed the excitement of solving basic scientific problems. Most rewarding, however, was the feeling that I was helping young scientists experience the pleasure of performing creative research.

Protein and drug interactions in the minor groove of DNA

Interactions between proteins, drugs, water and B-DNA minor groove have been analyzed in crystal structures of 60 protein-DNA and 14 drug-DNA complexes. It was found that only purine N3, pyrimidine O2, guanine N2 and deoxyribose O4' are involved in the interactions, and that contacts to N3 and O2 are most frequent and more polar than contacts to O4'. Many protein contacts are mediated by water, possibly to increase the DNA effective surface. Fewer water-mediated contacts are observed in drug complexes. The distributions of ligands around N3 are significantly more compact than around O2, and distributions of water molecules are the most compact. Distributions around O4' are more diffuse than for the base atoms but most distributions still have just one binding site. Ligands bind to N3 and O2 atoms in analogous positions, and simultaneous binding to N3 and N2 in guanines is extremely rare. Contacts with two consecutive nucleotides are much more frequent than base-sugar contacts within one nucleotide. The probable reason for this is the large energy of deformation of hydrogen bonds for the one nucleotide motif. Contacts of Arg, the most frequent amino acid ligand, are stereochemically indistinguishable from the binding of the remaining amino acids except asparagine (Asn) and phenylalanine (Phe). Asn and Phe bind in distinct ways, mostly to a deformed DNA, as in the complexes of TATA-box binding proteins. DNA deformation concentrates on dinucleotide regions with a distinct deformation of the delta and epsilon backbone torsion angles for the Asn and delta, epsilon, zeta and chi for the Phe-contacted regions.

Structural basis for cooperative DNA binding by two dimers of the multidrug-binding protein QacR

The Staphylococcus aureus multidrug-binding protein QacR represses transcription of the qacA multidrug transporter gene and is induced by multiple structurally dissimilar drugs. QacR is a member of the TetR/CamR family of transcriptional regulators, which share highly homologous N-terminal DNA-binding domains connected to seemingly non-homologous ligand-binding domains. Unlike other TetR members, which bind approximately 15 bp operators, QacR recognizes an unusually long 28 bp operator, IR1, which it appears to bind cooperatively. To elucidate the DNA-binding mechanism of QacR, we determined the 2.90 A resolution crystal structure of a QacR-IR1 complex. Strikingly, our data reveal that the DNA recognition mode of QacR is distinct from TetR and involves the binding of a pair of QacR dimers. In this unique binding mode, recognition at each IR1 half-site is mediated by a complement of DNA contacts made by two helix-turn-helix motifs. The inferred cooperativity does not arise from cross-dimer protein-protein contacts, but from the global undertwisting and major groove widening elicited by the binding of two QacR dimers.

Testing water-mediated DNA recognition by the Hin recombinase

The Hin recombinase specifically recognizes its DNA-binding site by means of both major and minor groove interactions. A previous X-ray structure, together with new structures of the Hin DNA-binding domain bound to a recombination half-site that were solved as part of the present study, have revealed that two ordered water molecules are present within the major groove interface. In this report, we test the importance of these waters directly by X-ray crystal structure analysis of complexes with four mutant DNA sequences. These structures, combined with their Hin-binding properties, provide strong support for the critical importance of one of the intermediate waters. A lesser but demonstrable role is ascribed to the second water molecule. The mutant structures also illustrate the prominent roles of thymine methyls both in stabilizing intermediate waters and in interfering with water or amino acid side chain interactions with DNA.

Dissecting the molecular origins of specific protein-nucleic acid recognition: hydrostatic pressure and molecular dynamics

The fundamental processes by which proteins recognize and bind to nucleic acids are critical to understanding cellular function. To explore the factors involved in protein-DNA recognition, we used hydrostatic pressure to perturb the binding of the BamHI endonuclease to cognate DNA, both in experiment and in molecular dynamic simulations. A new technique of high-pressure gel mobility shift analysis was used to test the effects of elevated hydrostatic pressure on the binding of BamHI to its cognate recognition sequence. Upon application of a pressure of 500 bar, the equilibrium dissociation constant of BamHI binding to the cognate site was found to increase nearly 10-fold. A challenge has been to link this type of pure thermodynamic measurement to functional events occurring at the molecular level. Thus, we used molecular dynamic simulations at both ambient and elevated pressures to reveal details of the direct and water-mediated interactions between BamHI and cognate DNA, which allow explanation of the effects of pressure on site-specific protein-DNA binding and complex stability.

Recognition of specific DNA sequences

Proteins that recognize specific DNA sequences play a central role in the regulation of transcription. The tremendous increase in structural information on protein-DNA complexes has uncovered a remarkable structural diversity in DNA binding folds, while at the same time revealing common themes in binding to target sites in the genome.

Mass spectrometric study of the Escherichia coli repressor proteins, Ic1R and Gc1R, and their complexes with DNA

In Escherichia coli, the IclR protein regulates both the aceBAK operon and its own synthesis. Database homology searches have identified many IclR-like proteins, now known as the IclR family, which can be identified by a conserved C-terminal region. We have cloned and purified one of these proteins, which we have named GclR (glyoxylate carboligase repressor). Although purification is straightforward, both the IclR and GclR proteins are difficult to manipulate, requiring high salt (up to 0.6 M KCl) for solubility. With the advent of nanospray ionization, we could transfer the proteins into much higher concentrations of volatile buffer than had been practical with ordinary electrospray. In 0.5 M ammonium bicarbonate buffer, both proteins were stable as tetramers, with a small amount of dimer. In a separate experiment, we found that IclR protein selected from a random pool a sequence which matched exactly that of the presumed binding region of the GclR protein, although IclR does not regulate the gcl gene. We designed a 29 bp synthetic DNA to which IclR and GclR bind, and with which we were able to form noncovalent DNA-protein complexes for further mass spectrometry analysis. These complexes were far more stable than the proteins alone, and we have evidence of a stoichiometry which has not been described previously with (protein monomer : dsDNA) = (4 : 1).

Amino acid-base interactions: a three-dimensional analysis of protein-DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid-base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein-DNA complex structures. DNA-backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base-amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein-DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein-DNA interface. In conclusion, the majority of amino acid-base interactions observed follow general principles that apply across all protein-DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are 'universal' and 'context-dependent'. An interactive Web-based atlas of side chain-base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.

Helix bending as a factor in protein/DNA recognition

Normal vectors perpendicular to individual base pairs are a powerful tool for studying the bending behavior of B-DNA, both in the form of normal vector plots and in matrices that list angles between vectors for all possible base pair combinations. A new analysis program, FREEHELIX, has been written for this purpose, and applied to 86 examples of sequence-specific protein/DNA complexes whose coordinates are on deposit in the Nucleic Acid Data Base. Bends in this sample of 86 structures almost invariably follow from roll angles between adjacent base pairs; tilt makes no net contribution. Roll in a direction compressing the broad major groove is much more common than that which compresses the minor groove. Three distinct types of B-DNA bending are observed, each with a different molecular origin: (1) Localized kinking is produced by large roll at single steps or at two steps separated by one turn of helix. (2) Smooth, planar curvature is produced by positive and negative roll angles spaced a half-turn apart, with random side-to-side zigzag roll at intermediate points, rather than a tilt contribution that might have been expected theoretically. (3) Three-dimensional writhe results from significant roll angles at a continuous series of steps. Writhe need not change the overall direction of helix axis, if it is continued indefinitely or for an integral number of helical turns. A-DNA itself can be formally considered as possessing uniform, continuous writhe that yields no net helix bending. Smooth curvature is the most intricate deformation of the three, and is least common. Writhe is the simplest deformation and is most common; indeed, a low level of continuous writhe is the normal condition of an otherwise unbent B-DNA helix of general sequence. With one exception, every example of major kinking in this sample of 86 structures involves a pyrimidine-purine step: C-A/T-G, T-A, or C-G. Purine-purine steps, especially A-A, show the least tendency toward roll deformations.

Conformational flexibility of B-DNA at 0.74 A resolution: d(CCAGTACTGG)(2)

The affinity and specificity of a ligand for its DNA site is a function of the conformational changes between the isolated and complexed states. Although the structures of a hydroxypyrrole-imidazole-pyrrole polyamide dimer with 5'-CCAGTACTGG-3' and the trp repressor recognizing the sequence 5'-GTACT-3' are known, the baseline conformation of the DNA site would contribute to our understanding of DNA recognition by these ligands. The 0.74 A resolution structure of a B-DNA double helix, 5'-CCAGTACTGG-3', has been determined by X-ray crystallography. Six of the nine phosphates, two of four bound calcium ions and networks of water molecules hydrating the oligonucleotide have alternate conformations. By contrast, nine of the ten bases have a single, unique conformation with hydrogen atoms visible in most cases. The polyamide molecules alter the geometry of the phosphodiester backbone, and the water molecules mediating contacts in the trp repressor/operator complex are conserved in the unliganded DNA. Furthermore, the multiple conformational states, ions and hydration revealed by this ultrahigh resolution structure of a B-form oligonucleotide are potentially general considerations for understanding DNA-binding affinity and specificity by ligands.

An overview of the structures of protein-DNA complexes

On the basis of a structural analysis of 240 protein-DNA complexes contained in the Protein Data Bank (PDB), we have classified the DNA-binding proteins involved into eight different structural/functional groups, which are further classified into 54 structural families. Here we present this classification and review the functions, structures and binding interactions of these protein-DNA complexes.

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.

Oligomerization of the EK18 mutant of the trp repressor of Escherichia coli as observed by NMR spectroscopy

The regulation of the trp repressor system of Escherichia coli is frequently modeled by a single equilibrium, that between the aporepressor (TR) and the corepressor, l-tryptophan (Trp), at their intracellular concentrations. The actual mechanism, which is much more complex and more finely tuned, involves multiple equilibria: TR and Trp association, TR oligomerization, specific and nonspecific binding of various states of TR to DNA, and interactions between these various species and ions. TR in isolation exists primarily as a homodimer, but the state of oligomerization increases as the TR concentration goes up and/or the salt concentration goes down, leading to species with lower affinity for DNA. We have used multinuclear, multidimensional NMR spectroscopy to investigate structural changes that accompany the oligomerization of TR. For these investigations, the superrepressor mutant EK18 (TR with Glu 18 replaced by Lys) was chosen because it exhibits less severe oligomerization at higher protein concentration than other known variants; this made it possible to study the dimer to tetramer oligomerization step by NMR. The NMR results suggest that the interaction between TR dimers is structurally linked to folding of the DNA binding domain and that it likely involves direct contacts between the C-terminal residues of the C-helix of one dimer with the next dimer. This implies that oligomerization can compete with DNA binding and thus serves as a factor in the fine-tuning of gene expression.

The hydration of nucleic acid duplexes as assessed by a combination of volumetric and structural techniques

Using high precision densimetric and ultrasonic measurements, we have determined, at 25 degrees C, the apparent molar volumes PhiV and the apparent molar compressibilities PhiK(S) of four nucleic acid duplexes-namely, the DNA duplex, poly(dIdC)poly(dIdC); the RNA duplex, poly(rA)poly(rU); and the two DNA/RNA hybrid duplexes, poly(rA)poly(dT) and poly(dA)poly(rU). Using available fiber diffraction data on these duplexes, we have calculated the molecular volumes as well as the solvent-accessible surface areas of the constituent charged, polar, and nonpolar atomic groups. We found that the hydration properties of these nucleic acid duplexes do not correlate with the extent and the chemical nature of the solvent-exposed surfaces, thereby suggesting a more specific set of duplex-water interactions beyond general solvation effects. A comparative analysis of our volumetric data on the four duplexes, in conjunction with available structural information, suggests the following features of duplex hydration: (a) The four duplexes exhibit different degrees of hydration, in the order poly(dIdC)poly(dIdC) > poly(dGdC)poly(dGdC) > poly(dAdT)poly(dAdT) approximately poly(dA)poly(dT). (b) Repetitive AT and IC sequences within a duplex are solvated beyond general effects by a spine of hydration in the minor groove, with this sequence-specific water network involving about 8 additional water molecules from the second and, perhaps, even the third hydration layers. (c) Repetitive GC and IC sequences within a duplex are solvated beyond general effects by a "patch of hydration" in the major groove, with this water network involving about 13 additional water molecules from the second and, perhaps, even the third hydration layers. (d) Random sequence, polymeric DNA duplexes, which statistically lack extended regions of repetitive AT, GC, or IC sequences, do not experience such specific enhancements of hydration. Consequently, consistent with our previous observations (T. V. Chalikian, A. P. Sarvazyan, G. E. Plum, and K. J. Breslauer, Biochemistry, 1994, Vol. 33, pp. 2394-2401), duplexes with approximately 50% AT content exhibit the weakest hydration, while an increase or decrease from this AT content causes enhancement of hydration, either due to stronger hydration of the minor groove (an increase in AT content) or due to stronger hydration of the major groove (an increase in GC content). (e) In dilute aqueous solutions, a B-DNA duplex is more hydrated than an A-DNA duplex, a volumetric-based conclusion that is in agreement with previous results obtained on crystals, fibers, and DNA solutions in organic solvent-water mixtures. (f) the A-like, RNA duplex poly(rA)poly(rU) and the structurally similar A-like, hybrid duplex poly(rA)poly(dT), exhibit similar hydration properties, while the structurally distinct A-like, hybrid duplex poly(rA)poly(dT) and non-A-like, hybrid duplex poly(dA)poly(rU) exhibit differential hydration properties, consistent with structural features dictating hydration characteristics. We discuss how volumetric characterizations, in conjunction with structural studies, can be used to describe, define, and resolve the general and sequence/conformation-specific hydration properties of nucleic acid duplexes.