Glycine insertion in the hinge region of lactose repressor protein alters DNA binding

Ligand-specific changes in conformational flexibility mediate long-range allostery in the lac repressor

Biological regulation ubiquitously depends on protein allostery, but the regulatory mechanisms are incompletely understood, especially in proteins that undergo ligand-induced allostery with few structural changes. Here we used hydrogen-deuterium exchange with mass spectrometry (HDX/MS) to map allosteric effects in a paradigm ligand-responsive transcription factor, the lac repressor (LacI), in different functional states (apo, or bound to inducer, anti-inducer, and/or DNA). Although X-ray crystal structures of the LacI core domain in these states are nearly indistinguishable, HDX/MS experiments reveal widespread differences in flexibility. We integrate these results with modeling of protein-ligand-solvent interactions to propose a revised model for allostery in LacI, where ligand binding allosterically shifts the conformational ensemble as a result of distinct changes in the rigidity of secondary structures in the different states. Our model provides a mechanistic basis for the altered function of distal mutations. More generally, our approach provides a platform for characterizing and engineering protein allostery.

A small-molecule chemical interface for molecular programs

In vitro molecular circuits, based on DNA-programmable chemistries, can perform an increasing range of high-level functions, such as molecular level computation, image or chemical pattern recognition and pattern generation. Most reported demonstrations, however, can only accept nucleic acids as input signals. Real-world applications of these programmable chemistries critically depend on strategies to interface them with a variety of non-DNA inputs, in particular small biologically relevant chemicals. We introduce here a general strategy to interface DNA-based circuits with non-DNA signals, based on input-translating modules. These translating modules contain a DNA response part and an allosteric protein sensing part, and use a simple design that renders them fully tunable and modular. They can be repurposed to either transmit or invert the response associated with the presence of a given input. By combining these translating-modules with robust and leak-free amplification motifs, we build sensing circuits that provide a fluorescent quantitative time-response to the concentration of their small-molecule input, with good specificity and sensitivity. The programmability of the DNA layer can be leveraged to perform DNA based signal processing operations, which we demonstrate here with logical inversion, signal modulation and a classification task on two inputs. The DNA circuits are also compatible with standard biochemical conditions, and we show the one-pot detection of an enzyme through its native metabolic activity. We anticipate that this sensitive small-molecule-to-DNA conversion strategy will play a critical role in the future applications of molecular-level circuitry.

Condensin-Mediated Chromosome Folding and Internal Telomeres Drive Dicentric Severing by Cytokinesis

In Saccharomyces cerevisiae, dicentric chromosomes stemming from telomere fusions preferentially break at the fusion. This process restores a normal karyotype and protects chromosomes from the detrimental consequences of accidental fusions. Here, we address the molecular basis of this rescue pathway. We observe that tandem arrays tightly bound by the telomere factor Rap1 or a heterologous high-affinity DNA binding factor are sufficient to establish breakage hotspots, mimicking telomere fusions within dicentrics. We also show that condensins generate forces sufficient to rapidly refold dicentrics prior to breakage by cytokinesis and are essential to the preferential breakage at telomere fusions. Thus, the rescue of fused telomeres results from a condensin- and Rap1-driven chromosome folding that favors fusion entrapment where abscission takes place. Because a close spacing between the DNA-bound Rap1 molecules is essential to this process, Rap1 may act by stalling condensins.

Price of disorder in the lac repressor hinge helix

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

The mechanism and high-free-energy transition state of lac repressor-lac operator interaction

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)–lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (k_a) and dissociation (k_d) rate constants. RO-destabilizing solutes (urea, KCl) reduce k_a comparably (urea) or more than (KCl) they increase k_d, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce k_d without affecting k_a significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Data on publications, structural analyses, and queries used to build and utilize the AlloRep database

The AlloRep database (www.AlloRep.org) (Sousa et al., 2016) [1] compiles extensive sequence, mutagenesis, and structural information for the LacI/GalR family of transcription regulators. Sequence alignments are presented for >3000 proteins in 45 paralog subfamilies and as a subsampled alignment of the whole family. Phenotypic and biochemical data on almost 6000 mutants have been compiled from an exhaustive search of the literature; citations for these data are included herein. These data include information about oligomerization state, stability, DNA binding and allosteric regulation. Protein structural data for 65 proteins are presented as easily-accessible, residue-contact networks. Finally, this article includes example queries to enable the use of the AlloRep database. See the related article, “AlloRep: a repository of sequence, structural and mutagenesis data for the LacI/GalR transcription regulators” (Sousa et al., 2016) [1].

DNA supercoiling, a critical signal regulating the basal expression of the lac operon in Escherichia coli

Escherichia coli lac repressor (LacI) is a paradigmatic transcriptional factor that controls the expression of lacZYA in the lac operon. This tetrameric protein specifically binds to the O1, O2 and O3 operators of the lac operon and forms a DNA loop to repress transcription from the adjacent lac promoter. In this article, we demonstrate that upon binding to the O1 and O2 operators at their native positions LacI constrains three (−) supercoils within the 401-bp DNA loop of the lac promoter and forms a topological barrier. The stability of LacI-mediated DNA topological barriers is directly proportional to its DNA binding affinity. However, we find that DNA supercoiling modulates the basal expression from the lac operon in E. coli. Our results are consistent with the hypothesis that LacI functions as a topological barrier to constrain free, unconstrained (−) supercoils within the 401-bp DNA loop of the lac promoter. These constrained (−) supercoils enhance LacI’s DNA-binding affinity and thereby the repression of the promoter. Thus, LacI binding is superhelically modulated to control the expression of lacZYA in the lac operon under varying growth conditions.

Binding of Multiple Rap1 Proteins Stimulates Chromosome Breakage Induction during DNA Replication

Telomeres, the ends of linear eukaryotic chromosomes, have a specialized chromatin structure that provides a stable chromosomal terminus. In budding yeast Rap1 protein binds to telomeric TG repeat and negatively regulates telomere length. Here we show that binding of multiple Rap1 proteins stimulates DNA double-stranded break (DSB) induction at both telomeric and non-telomeric regions. Consistent with the role of DSB induction, Rap1 stimulates nearby recombination events in a dosage-dependent manner. Rap1 recruits Rif1 and Rif2 to telomeres, but neither Rif1 nor Rif2 is required for DSB induction. Rap1-mediated DSB induction involves replication fork progression but inactivation of checkpoint kinase Mec1 does not affect DSB induction. Rap1 tethering shortens artificially elongated telomeres in parallel with telomerase inhibition, and this telomere shortening does not require homologous recombination. These results suggest that Rap1 contributes to telomere homeostasis by promoting chromosome breakage.

Telomere length is maintained primarily through equilibrium between telomerase-mediated lengthening and the loss of telomeric sequence through the end-replication problem. In budding yeast Rap1 protein binds to telomeric TG repeat and negatively regulates telomerase recruitment in a dosage-dependent manner. In this paper we provide evidence suggesting an alternative Rap1-dependent telomere shortening mechanism in which binding of multiple Rap1 proteins mediates DNA break induction during DNA replication. This process does not involve recombination events; therefore, it is distinct from loop-mediated telomere trimming.

The length of glycine-rich linker in DNA-binding domain is critical for optimal functioning of quorum-sensing master regulatory protein HapR

HapR is a quorum-sensing master regulatory protein in Vibrio cholerae. Though many facts are known regarding its structural and functional aspects, much still can be learnt from natural variants of this wild-type protein. While unraveling the underlying cause of functional inertness of a natural variant (HapRV2), the significance of a conserved glycine residue at position 39 in a glycine-rich linker in DNA-binding domain comes into light. This work aims at investigating how the length of glycine-rich linker (R(33)GIGRGG(39)) bridging helices α1 and α2 modulates the functionality of HapR. In pursuit of our interest, glycine residues were inserted after terminal glycine (G39) of the linker in a sequential manner. To evaluate functionality, all the glycine linker variants were subjected to a battery of performance tests under various conditions. Combined in vitro and in vivo results clearly demonstrated a gradual functional impairment of HapR linker variants coupled with increasing length of glycine-rich linker and finally, linker variant harboring four glycine residues resulted in a functionally compromised protein with significant loss of communication with cognate DNAs. Molecular dynamics studies of modeled HapR linker variants in complex with cognate promoter region show that residues namely Ser50, Thr53 and Asn56 are involved in varying degree of interactions with different nucleotides of HapR-DNA complex. The diminished functionality between variants and DNA appears to result from reduced or no interactions between Phe55 and nucleotides of cognate DNA as observed during simulations.

Alginate-dependent gene expression mechanism in Sphingomonas sp. strain A1

Sphingomonas sp. strain A1, a Gram-negative bacterium, directly incorporates alginate polysaccharide into the cytoplasm through a periplasmic alginate-binding protein-dependent ATP-binding cassette transporter. The polysaccharide is degraded to monosaccharides via the formation of oligosaccharides by endo- and exotype alginate lyases. The strain A1 proteins for alginate uptake and degradation are encoded in both strands of a genetic cluster in the bacterial genome and inducibly expressed in the presence of alginate. Here we show the function of the alginate-dependent transcription factor AlgO and its mode of action on the genetic cluster and alginate oligosaccharides. A putative gene within the genetic cluster seems to encode a transcription factor-like protein (AlgO). Mutant strain A1 (ΔAlgO mutant) cells with a disrupted algO gene constitutively produced alginate-related proteins. DNA microarray analysis indicated that wild-type cells inducibly transcribed the genetic cluster only in the presence of alginate, while ΔAlgO mutant cells constitutively expressed the genetic cluster. A gel mobility shift assay showed that AlgO binds to the specific intergenic region between algO and algS (algO-algS). Binding of AlgO to the algO-algS intergenic region diminished with increasing alginate oligosaccharides. These results demonstrated a novel alginate-dependent gene expression mechanism. In the absence of alginate, AlgO binds to the algO-algS intergenic region and represses the expression of both strands of the genetic cluster, while in the presence of alginate, AlgO dissociates from the algO-algS intergenic region via binding to alginate oligosaccharides produced through the lyase reaction and subsequently initiates transcription of the genetic cluster. This is the first report on the mechanism by which alginate regulates the expression of the gene cluster.

Scoring and manipulating gene position and dynamics using FROS in budding yeast

The spatial organization of the genome within the nucleus is now seen as a key contributor to genome function. Studying chromatin dynamics in living cells has been rendered possible by the development of fast microscopy coupled with fluorescent repressor operator systems (FROS). In these systems, arrays of protein-binding sites integrated at specific loci by homologous recombination are monitored through the fluorescence of tagged DNA-binding proteins. In the budding yeast, where homologous recombination is efficient, this technique, combined with targeting assay and genetic analysis, has been extremely powerful for studying the determinants and function of chromatin dynamics in living cells. However, issues have been recurrently raised in different species regarding the use of these systems. Here we discuss the different uses of gene tagging with FROS and their limitations, focusing in budding yeast as a model organism.

Tight protein-DNA interactions favor gene silencing

The heterochromatin-like structure formed by the yeast silent information regulator complex (SIR) represses transcription at the silent mating type loci and telomeres. Here, we report that tight protein-DNA complexes induce ectopic recruitment of the SIR complex, promoting gene silencing and changes in subnuclear localization when cis-acting elements are nearby. Importantly, lack of the replication fork-associated helicase Rrm3 enhances this induced gene repression. Additionally, Sir3 and Sir4 are enriched genome-wide at natural replication pause sites, including tRNA genes. Consistently, inserting a tRNA gene promotes SIR-mediated silencing of a nearby gene. These results reveal that replication stress arising from tight DNA-protein interactions favors heterochromatin formation.

Monitoring DNA binding to Escherichia coli lactose repressor using quartz crystal microbalance with dissipation

Lactose repressor protein (LacI) functions as a negative transcription regulator in Escherichia coli by binding to the operator DNA sequence. Our understanding of the immobilized LacI function and the effect of ligand binding on the conformation of LacI-DNA complexes remains poorly understood. Here, we have examined the difference in functionality of wild-type and mutant LacI binding to the target DNA using quartz crystal microbalance with dissipation (QCM-D). To direct the orientation of LacI binding to the gold surface, residue 334 was substituted with cysteine (T334C) to generate a sulfur-gold linkage. Position 334 is located on the surface opposite the DNA-binding domain and remote from the site for inducer binding. With T334C immobilized on the gold surface, our sensors successfully detect operator binding as well as the release of the operator DNA from the repressor in the presence of inducer isopropyl-β-D-thiogalactoside (IPTG). Besides the natural operator DNA sequence (O(1)), a symmetric high-affinity DNA sequence (O(sym)), and a non-specific DNA (O(ns)) sequence with low affinity were also used. In addition, the impact of anti-inducer o-nitrophenyl-beta-d-fucoside (ONPF), which stabilizes LacI operator binding, was examined. The results from immobilized mutant LacI are in good agreement with known solution parameters for LacI-ligand interactions, demonstrating that QCM-D provides a rapid and efficient measurement of DNA binding and impact of ligands upon binding for this complex oligomeric protein.

Flexibility in the inducer binding region is crucial for allostery in the Escherichia coli lactose repressor

Lactose repressor protein (LacI) utilizes an allosteric mechanism to regulate transcription in Escherichia coli, and the transition between inducer- and operator-bound states has been simulated by targeted molecular dynamics (TMD). The side chains of amino acids 149 and 193 interact and were predicted by TMD simulation to play a critical role in the early stages of the LacI conformational change. D149 contacts IPTG directly, and variations at this site provide the opportunity to dissect its role in inducer binding and signal transduction. Single mutants at D149 or S193 exhibit a minimal change in operator binding, and alterations in inducer binding parallel changes in operator release, indicating normal allosteric response. The observation that the double mutant D149A/S193A exhibits wild-type properties excludes the requirement for inter-residue hydrogen bond formation in the allosteric response. The double mutant D149C/S193C purified from cell extracts shows decreased sensitivity to inducer binding while retaining wild-type binding affinities and kinetic constants for both operator and inducer. By manipulating cysteine oxidation, we show that the more reduced state of D149C/S193C responds to inducer more like the wild-type protein, whereas the more oxidized state displays diminished inducer sensitivity. These features of D149C/S193C indicate that the novel disulfide bond formed in this mutant impedes the allosteric transition, consistent with the role of this region predicted by TMD simulation. Together, these results establish the requirement for flexibility in the spatial relationship between D149 and S193 rather than a specific D149-S193 interaction in the LacI allosteric response to inducer.

lac repressor is an antivirulence factor of Salmonella enterica: its role in the evolution of virulence in Salmonella

The genus Salmonella includes many pathogens of great medical and veterinary importance. Bacteria belonging to this genus are very closely related to those belonging to the genus Escherichia. lacZYA operon and lacI are present in Escherichia coli, but not in Salmonella enterica. It has been proposed that Salmonella has lost lacZYA operon and lacI during evolution. In this study, we have investigated the physiological and evolutionary significance of the absence of lacI in Salmonella enterica. Using murine model of typhoid fever, we show that the expression of LacI causes a remarkable reduction in the virulence of Salmonella enterica. LacI also suppresses the ability of Salmonella enterica to proliferate inside murine macrophages. Microarray analysis revealed that LacI interferes with the expression of virulence genes of Salmonella pathogenicity island 2. This effect was confirmed by RT-PCR and Western blot analysis. Interestingly, we found that SBG0326 of Salmonella bongori is homologous to lacI of Escherichia coli. Salmonella bongori is the only other species of the genus Salmonella and it lacks the virulence genes of Salmonella pathogenicity island 2. Overall, our results demonstrate that LacI is an antivirulence factor of Salmonella enterica and suggest that absence of lacI has facilitated the acquisition of virulence genes of Salmonella pathogenicity island 2 in Salmonella enterica making it a successful systemic pathogen.

Concentration and length dependence of DNA looping in transcriptional regulation

In many cases, transcriptional regulation involves the binding of transcription factors at sites on the DNA that are not immediately adjacent to the promoter of interest. This action at a distance is often mediated by the formation of DNA loops: Binding at two or more sites on the DNA results in the formation of a loop, which can bring the transcription factor into the immediate neighborhood of the relevant promoter. These processes are important in settings ranging from the historic bacterial examples (bacterial metabolism and the lytic-lysogeny decision in bacteriophage), to the modern concept of gene regulation to regulatory processes central to pattern formation during development of multicellular organisms. Though there have been a variety of insights into the combinatorial aspects of transcriptional control, the mechanism of DNA looping as an agent of combinatorial control in both prokaryotes and eukaryotes remains unclear. We use single-molecule techniques to dissect DNA looping in the lac operon. In particular, we measure the propensity for DNA looping by the Lac repressor as a function of the concentration of repressor protein and as a function of the distance between repressor binding sites. As with earlier single-molecule studies, we find (at least) two distinct looped states and demonstrate that the presence of these two states depends both upon the concentration of repressor protein and the distance between the two repressor binding sites. We find that loops form even at interoperator spacings considerably shorter than the DNA persistence length, without the intervention of any other proteins to prebend the DNA. The concentration measurements also permit us to use a simple statistical mechanical model of DNA loop formation to determine the free energy of DNA looping, or equivalently, the for looping.

Comparison of deterministic and stochastic models of the lac operon genetic network

The lac operon has been a paradigm for genetic regulation with positive feedback, and several modeling studies have described its dynamics at various levels of detail. However, it has not yet been analyzed how stochasticity can enrich the system's behavior, creating effects that are not observed in the deterministic case. To address this problem we use a comparative approach. We develop a reaction network for the dynamics of the lac operon genetic switch and derive corresponding deterministic and stochastic models that incorporate biological details. We then analyze the effects of key biomolecular mechanisms, such as promoter strength and binding affinities, on the behavior of the models. No assumptions or approximations are made when building the models other than those utilized in the reaction network. Thus, we are able to carry out a meaningful comparison between the predictions of the two models to demonstrate genuine effects of stochasticity. Such a comparison reveals that in the presence of stochasticity, certain biomolecular mechanisms can profoundly influence the region where the system exhibits bistability, a key characteristic of the lac operon dynamics. For these cases, the temporal asymptotic behavior of the deterministic model remains unchanged, indicating a role of stochasticity in modulating the behavior of the system.

Functional consequences of exchanging domains between LacI and PurR are mediated by the intervening linker sequence

Homologue function can be differentiated by changing residues that affect binding sites or long-range interactions. LacI and PurR are two proteins that represent the LacI/GalR family (>500 members) of bacterial transcription regulators. All members have distinct DNA-binding and regulatory domains linked by approximately 18 amino acids. Each homologue has specificity for different DNA and regulatory effector ligands; LacI and PurR also exhibit differences in allosteric communication between DNA and effector binding sites. A comparative study of LacI and PurR suggested that alterations in the interface between the regulatory domain and linker are important for differentiating their functions. Four residues (equivalent to LacI positions 48, 55, 58, and 61) appear particularly important for creating a unique interface and were predicted to be necessary for allosteric regulation. However, nearby residues in the linker interact with DNA ligand. Thus, differences observed in interactions between linker and regulatory domain may be the cause of altered function or an effect of the two proteins binding different DNA ligands. To separate these possibilities, we created a chimeric protein with the LacI DNA-binding domain/linker and the PurR regulatory domain (LLhP). If the interface requires homologue-specific interactions in order to propagate the signal from effector binding, then LLhP repression should not be allosterically regulated by effector binding. Experiments show that LLhP is capable of repression from lacO1 and, contrary to expectation, allosteric response is intact. Further, restoring the potential for PurR-like interactions via substitutions in the LLhP linker tends to diminish repression. These effects are especially pronounced for residues 58 and 61. Clearly, binding affinity of LLhP for the lacO1 DNA site is sensitive to long-range changes in the linker. This result also raises the possibility that mutations at positions 58 and 61 co-evolved with changes in the DNA-binding site. In addition, repression measured in the absence and presence of effector ligand shows that allosteric response increases for several LLhP variants with substitutions at positions 48 and 55. Thus, while side chain variation at these sites does not generally dictate the presence or absence of allostery, the nature of the amino acid can modulate the response to effector.

DNA Binding to Protein−Gold Nanocrystal Conjugates

The E. coli DNA binding protein lac repressor (LacI) and a derivative with a designed thiol (T334C) were developed as gold nanocrystal conjugates to assess the effects of conjugation on DNA binding function. The designed derivative was engineered with a solvent-accessible thiol to promote oriented conjugation, avoiding obstruction of the DNA-binding domain by the nanocrystal. Analytical ultracentrifugation (AU) and electrophoretic mobility shift assays (EMSA) were used to evaluate the ability of conjugated repressors to bind the natural operator DNA sequence O(1). The results show that LacI does not retain significant DNA binding function when conjugated to gold nanocrystals, presumably because the basic DNA-binding domain is the site for nonspecific conjugation. T334C, with the potential for both directed and nonspecific conjugation, shows enhanced interaction with O(1) when conjugated. Interestingly, the order of component addition is a key factor in producing functional lac repressor conjugates.

Lac repressor hinge flexibility and DNA looping: single molecule kinetics by tethered particle motion

The tethered particle motion (TPM) allows the direct detection of activity of a variety of biomolecules at the single molecule level. First pioneered for RNA polymerase, it has recently been applied also to other enzymes. In this work we employ TPM for a systematic investigation of the kinetics of DNA looping by wild-type Lac repressor (wt-LacI) and by hinge mutants Q60G and Q60 + 1. We implement a novel method for TPM data analysis to reliably measure the kinetics of loop formation and disruption and to quantify the effects of the protein hinge flexibility and of DNA loop strain on such kinetics. We demonstrate that the flexibility of the protein hinge has a profound effect on the lifetime of the looped state. Our measurements also show that the DNA bending energy plays a minor role on loop disruption kinetics, while a strong effect is seen on the kinetics of loop formation. These observations substantiate the growing number of theoretical studies aimed at characterizing the effects of DNA flexibility, tension and torsion on the kinetics of protein binding and dissociation, strengthening the idea that these mechanical factors in vivo may play an important role in the modulation of gene expression regulation.

Extrinsic interactions dominate helical propensity in coupled binding and folding of the lactose repressor protein hinge helix

A significant number of eukaryotic regulatory proteins are predicted to have disordered regions. Many of these proteins bind DNA, which may serve as a template for protein folding. Similar behavior is seen in the prokaryotic LacI/GalR family of proteins that couple hinge-helix folding with DNA binding. These hinge regions form short alpha-helices when bound to DNA but appear to be disordered in other states. An intriguing question is whether and to what degree intrinsic helix propensity contributes to the function of these proteins. In addition to its interaction with operator DNA, the LacI hinge helix interacts with the hinge helix of the homodimer partner as well as to the surface of the inducer-binding domain. To explore the hierarchy of these interactions, we made a series of substitutions in the LacI hinge helix at position 52, the only site in the helix that does not interact with DNA and/or the inducer-binding domain. The substitutions at V52 have significant effects on operator binding affinity and specificity, and several substitutions also impair functional communication with the inducer-binding domain. Results suggest that helical propensity of amino acids in the hinge region alone does not dominate function; helix-helix packing interactions appear to also contribute. Further, the data demonstrate that variation in operator sequence can overcome side chain effects on hinge-helix folding and/or hinge-hinge interactions. Thus, this system provides a direct example whereby an extrinsic interaction (DNA binding) guides internal events that influence folding and functionality.

Ligand interactions with lactose repressor protein and the repressor-operator complex: the effects of ionization and oligomerization on binding

Specific interactions between proteins and ligands that modify their functions are crucial in biology. Here, we examine sugars that bind the lactose repressor protein (LacI) and modify repressor affinity for operator DNA using isothermal titration calorimetry and equilibrium DNA binding experiments. High affinity binding of the commonly-used inducer isopropyl-beta,D-thiogalactoside is strongly driven by enthalpic forces, whereas inducer 2-phenylethyl-beta,D-galactoside has weaker affinity with low enthalpic contributions. Perturbing the dimer interface with either pH or oligomeric state shows that weak inducer binding is sensitive to changes in this distant region. Effects of the neutral compound o-nitrophenyl-beta,D-galactoside are sensitive to oligomerization, and at elevated pH this compound converts to an anti-inducer ligand with slightly enhanced enthalpic contributions to the binding energy. Anti-inducer o-nitrophenyl-beta,D-fucoside exhibits slightly enhanced affinity and increased enthalpic contributions at elevated pH. Collectively, these results both demonstrate the range of energetic consequences that occur with LacI binding to structurally-similar ligands and expand our insight into the link between effector binding and structural changes at the subunit interface.

Integrated insights from simulation, experiment, and mutational analysis yield new details of LacI function

Protein structural change underlies many signal transduction processes. Although end-state structures are known for various allosteric proteins, intermediates are difficult to observe. Recently, targeted molecular dynamics simulation (TMD) was used to examine the conformational transition and predict relevant intermediates for wild-type lactose repressor (LacI). A catalog of involved residues suggests that the transition of this homodimer is asymmetric and that K84 is a prominent participant in the dynamic N-subdomain interface. Previous experiments indicated that hydrophobic substitutions at position 84 engender slowed, biphasic inducer binding kinetics, which might reflect the same phenomena observed in TMD. Here, we report biochemical confirmation that DNA and inducer binding remain allosterically linked in K84A and K84L, albeit with a differential smaller than that found in wild-type LacI. Other features of these mutant proteins are consistent with an allosteric conformational shift that approximates that of the wild type. As a consequence, these repressors can be utilized to explore an unanswered question about LacI function: How many inducers (one or two per dimer) are required to diminish operator affinity? The biphasic natures of the K84L and K84A inducer association rates allow direct correlation between the two distinct inducer binding events and operator release. Indeed, the kinetics of operator release for the K84A and K84L closely parallel those for the second inducer binding event. Together with implications from previous equilibrium results for wild-type and mutant proteins, these kinetic data demonstrate that binding of two inducers per dimeric DNA binding unit is required to release the operator in these variant LacI proteins.

The lac repressor

Few proteins have had such a strong impact on a field as the lac repressor has had in Molecular Biology. Over 40 years ago, Jacob and Monod [Genetic regulatory mechanisms in the synthesis of proteins, J. Mol. Biol. 3 (1961) 318] proposed a model for gene regulation, which survives essentially unchanged in contemporary textbooks. It is a cogent depiction of how a set of 'structural' genes may be coordinately transcribed in response to environmental conditions and regulates metabolic events in the cell. In bacteria, the genes required for lactose utilization are negatively regulated when a repressor molecule binds to an upstream cis activated operator. The repressor and its operator together form a genetic switch, the lac operon. The switch functions when inducer molecules alter the conformation of the repressor in a specific manner. In the presence of a particular metabolite, the repressor undergoes a conformational change that reduces its affinity for the operator. The structures of the lac repressor and its complexes with operator DNA and effector molecules have provided a physical platform for visualizing at the molecular level the different conformations the repressor and the molecular basis for the switch. The structures of lac repressor, bound to its operator and inducer, have also been invaluable for interpreting a plethora of biochemical and genetic data.

Investigating the putative glycine hinge in Shaker potassium channel

The crystal structure of an open potassium channel reveals a kink in the inner helix that lines the pore (Jiang, Y.X., A. Lee, J.Y. Chen, M. Cadene, B.T. Chait, and R. MacKinnon. 2002. Nature 417:523–526). The putative hinge point is a highly conserved glycine residue. We examined the role of the homologous residue (Gly466) in the S6 transmembrane segment of Shaker potassium channels. The nonfunctional alanine mutant G466A will assemble, albeit poorly, with wild-type (WT) subunits, suppressing functional expression. To test if this glycine residue is critical for activation gating, we did a glycine scan along the S6 segment in the background of G466A. Although all of these double mutants lack the higher-level glycosylation that is characteristic of mature Shaker channels, one (G466A/V467G) is able to generate voltage-dependent potassium current. Surface biotinylation shows that functional and nonfunctional constructs containing G466A express at comparable levels in the plasma membrane. Compared with WT channels, the shifted-glycine mutant has impairments in voltage-dependent channel opening, including a right-shifted activation curve and a decreased rate of activation. The double mutant has relatively normal open-channel properties, except for a decreased affinity for intracellular blockers, a consequence of the loss of the side chain of Val467. Control experiments with the double mutants M440A/G466A and G466A/V467A suggest that the flexibility provided by Gly466 is more important for channel function than its small size. Our results support roles for Gly466 both in biogenesis of the channel and as a hinge in activation gating.

Coupling between catalytic site and collective dynamics: a requirement for mechanochemical activity of enzymes

Growing evidence supports the view that enzymatic activity results from a subtle interplay between chemical kinetics and molecular motions. A systematic analysis is performed here to delineate the type and level of coupling between catalysis and conformational mechanics. The dynamics of a set of 98 enzymes representative of different EC classes are analyzed with the Gaussian network model (GNM) and compared with experimental data. In more than 70% of the examined enzymes, the global hinge centers predicted by the GNM are found to be colocalized with the catalytic sites experimentally identified. Low translational mobility (< 7%) is observed for the catalytic residues, consistent with the fine-tuned design of enzymes to achieve precise mechanochemical activities. Ligand binding sites, while closely neighboring catalytic sites, enjoy a moderate flexibility to accommodate the ligand binding. These findings could serve as additional criteria for assessing drug binding residues and could lessen the computational burden of substrate docking searches.

Perturbation from a distance: mutations that alter LacI function through long-range effects

Allosteric modification of ligand binding is central to LacI transcription control. Recently, the conformational change between LacI operator- and inducer-bound states was simulated with targeted molecular dynamics (TMD) [Flynn, T. C., Swint-Kruse, L., Kong, Y., Booth, C., Matthews, K. S., and Ma, J. (2003) Protein Sci., 12, 2523-2541]. Atomic-level analyses of TMD results indicate the structural importance of the core pivot region that connects the N- and C-subdomains flanking the inducer-binding site. Further, a number of LacI mutations in the core pivot have been identified recently by their altered behaviors in phenotypic screens. Biochemical characterization of three of these variants-L148F, S151P, and P320A-provides an opportunity to directly explore the role of the core pivot in repressor function. For L148F, inducer IPTG binding affinity is strengthened, whereas O(1) operator DNA binding is diminished approximately 30-fold. In contrast, O(1) binding is increased for S151P, whereas IPTG binding is decreased. UV-difference spectroscopy and urea denaturation indicate long-range effects in both variants. Interestingly, P320A binds to DNA approximately 4-fold more tightly than wild-type, yet inducer binding is unaffected. To examine linkage between the core pivot and DNA binding domains, the L148F substitution was combined with Q60G, a previously known mutant with enhanced operator affinity. The double mutant exhibits the properties of both parent proteins, resulting in near wild-type DNA binding affinity and enhanced inducer sensitivity. These features may render Q60G/L148F more cost-effective in technological applications than wild-type repressor. As a group, the behaviors of the core pivot mutants are consistent with the allosteric structural role predicted for this region by TMD and reflect the significant long-range impact that single substitutions can elicit on protein function.

Detection and prevention of protein aggregation before, during, and after purification

The use of proteins for in vitro studies or as therapeutic agents is frequently hampered by protein aggregation during expression, purification, storage, or transfer into requisite assay buffers. A large number of potential protein stabilizers are available, but determining which are appropriate can take days or weeks. We developed a solubility assay to determine the best cosolvent for a given protein that requires very little protein and only a few hours to complete. This technique separates native protein from soluble and insoluble aggregates by filtration and detects both forms of protein by SDS-PAGE or Western blotting. Multiple buffers can be simultaneously screened to determine conditions that enhance protein solubility. The behavior of a single protein in mixtures and crude lysates can be analyzed with this technique, allowing testing prior to and throughout protein purification. Aggregated proteins can also be assayed for conditions that will stabilize native protein, which can then be used to improve subsequent purifications. This solubility assay was tested using both prokaryotic and eukaryotic proteins that range in size from 17 to 150 kDa and include monomeric and multimeric proteins. From the results presented, this technique can be applied to a variety of proteins.

Arginine operator binding by heterologous and chimeric ArgR repressors from Escherichia coli and Bacillus stearothermophilus

Bacillus stearothermophilus ArgR binds efficiently to the Escherichia coli carAB operator, whereas the E. coli repressor binds very poorly to the argCo operator of B. stearothermophilus. In order to elucidate this contradictory behavior between ArgRs, we constructed chimeric proteins by swapping N-terminal DNA-binding and C-terminal oligomerization domains or by exchanging the linker peptide. Chimeras carrying the E. coli DNA-binding domain and the B. stearothermophilus oligomerization domain showed sequence-nonspecific rather than sequence-specific interactions with arg operators. Chimeras carrying the B. stearothermophilus DNA-binding domain and E. coli oligomerization domain exhibited a high DNA-binding affinity for the B. stearothermophilus argCo and E. coli carAB operators and repressed the reporter-gene transcription from the B. stearothermophilus PargCo control region in vitro; arginine had no effect on, and indeed even decreased, their DNA-binding affinity. With the protein array method, we showed that the wild-type B. stearothermophilus ArgR and derivatives of it containing only the exchanged linker from E. coli ArgR or carrying the B. stearothermophilus DNA-binding domain along with the linker and the alpha4 regions were able to bind argCo containing the single Arg box. This binding was weaker than binding to the two-box operator but was no longer arginine dependent. Several lines of observations indicate that the alpha4 helix in the oligomerization domain and the linker peptide can contribute to the recognition of single or double Arg boxes and therefore to the operator DNA-binding specificity in similar but not identical ArgR repressors from two distant bacteria.

Fine-tuning function: correlation of hinge domain interactions with functional distinctions between LacI and PurR

LacI and PurR are highly homologous proteins. Their functional units are homodimers, with an N-terminal DNA binding domain that comprises the helix-turn-helix (HTH), N-linker, and hinge regions from both monomers. Hinge structural changes are known to occur upon DNA dissociation but are difficult to monitor experimentally. The initial steps of hinge unfolding were therefore examined using molecular dynamics simulations, utilizing a truncated, chimeric protein comprising the LacI HTH/N-linker and PurR hinge. A terminal Gly-Cys-Gly was added to allow "dimerization" through disulfide bond formation. Simulations indicate that differences in LacI and PurR hinge primary sequence affect the quaternary structure of the hinge x hinge' interface. However, these alternate hinge orientations would be sterically restricted by the core domain. These results prompted detailed comparison of recently available DNA-bound structures for LacI and truncated LacI(1-62) with the PurR structure. Examination revealed that different N-linker and hinge contacts to the core domain of the partner monomer (which binds effector molecule) affect the juxtapositions of the HTH, N-linker, and hinge regions in the DNA binding domain. In addition, the two full-length repressors exhibit significant differences in the interactions between the core and the C-linker connection to the DNA binding domain. Both linkers and the hinge have been implicated in the allosteric response of these repressors. Intriguingly, one functional difference between these two proteins is that they exhibit opposite allosteric response to effector. Simulations and observed structural distinctions are correlated with mutational analysis and sequence information from the LacI/GalR family to formulate a mechanism for fine-tuning individual repressor function.

Engineered disulfide linking the hinge regions within lactose repressor dimer increases operator affinity, decreases sequence selectivity, and alters allostery

The hinge domain encompasses amino acids 51-60 of lactose repressor (LacI) and plays an important role in its regulatory interaction with operator DNA. This segment makes both hinge-DNA and hinge-hinge' contacts that are critical to DNA binding. Furthermore, this small region serves as a central element in communicating the allosteric response to inducer. Introducing a disulfide bond between partner hinges within a dimer via the mutation V52C results in a protein that has increased affinity for O(1) operator DNA compared to wild-type LacI and abolishes allosteric response to inducer [Falcon, C. M., Swint-Kruse, L., and Matthews, K. S. (1997) J. Biol. Chem. 272, 26818]. We have established that this high affinity is maintained for the disulfide-linked protein even when symmetry and half-site spacing within the operator region are altered, whereas binding by the reduced protein, as for wild-type LacI, is severely diminished by these alterations. Interestingly, the allosteric response to inducer for V52C-oxidized remains intact for a small group of operator variants. Temperature studies demonstrate that the presence of the disulfide alters the thermodynamics of the protein-DNA interaction, with a DeltaC(p) of significantly smaller magnitude compared to wild-type LacI. The results presented here establish the hinge region as an important element not only for LacI high-affinity operator binding but also for the essential communication between ligand binding domains. Moreover, the results confirm that DNA sequence/conformation can profoundly influence allostery for this prototypic regulatory protein.

The position of the alpha and beta subunits in a single chain variant of human chorionic gonadotropin affects the heterodimeric interaction of the subunits and receptor-binding epitopes

The glycoprotein hormone family represents a class of heterodimers, which include the placental hormone human chorionic gonadotropin (CG) and the anterior pituitary hormones follitropin, lutropin, and thyrotropin. They are composed of common alpha subunit and a hormone-specific beta subunit. Based on the CG crystal structure, it was suggested that the quaternary subunit interactions are crucial for biological activity. However, recent observations using single chain glycoprotein hormone analogs, where the beta and alpha subunits are linked (NH(2)-CGbeta-alpha; CGbetaalpha orientation), implied that the heterodimeric-like quaternary configuration is not a prerequisite for receptor binding/signal transduction. To study the heterodimeric alignment of the two subunit domains in a single chain and its role in the intracellular behavior and biological action of the hormone, a single chain CG variant was constructed in which the carboxyl terminus of alpha was fused to the CGbeta amino terminus (NH(2)-alpha-CGbeta; alphaCGbeta orientation). The secretion rate of alphaCGbeta from transfected Chinese hamster ovary cells was less than that seen for CGbetaalpha. The alphaCGbeta tether was not recognized by dimer-specific monoclonal antibodies and did not bind to lutropin/CG receptor. To define if one or both subunit domains were modified in alphaCGbeta, it was co-transfected with a monomeric alpha or CGbeta gene. In each case, alphaCGbeta/alpha and alphaCGbeta/CGbeta complexes were formed indicating that CG dimer-specific epitopes were established. The alphaCGbeta/alpha complex bound to receptor indicating that the beta domain in the alphaCGbeta tether was still functional. In contrast, no significant receptor binding of alphaCGbeta/CGbeta was observed indicating a major perturbation in the alpha domain. These results suggest that although dimeric-like determinants are present in both alphaCGbeta/alpha and alphaCGbeta/CGbeta complexes, the receptor binding determinants in the alpha domain of the tether are absent. These results show that generating heterodimeric determinants do not necessarily result in a bioactive molecule. Our data also indicate that the determinants for biological activity are distinct from those associated with intracellular behavior.

Ion concentration and temperature dependence of DNA binding: comparison of PurR and LacI repressor proteins

Purine repressor (PurR) binding to specific DNA is enhanced by complexing with purines, whereas lactose repressor (LacI) binding is diminished by interaction with inducer sugars despite 30% identity in their protein sequences and highly homologous tertiary structures. Nonetheless, in switching from low- to high-affinity DNA binding, these proteins undergo a similar structural change in which the hinge region connecting the DNA and effector binding domains folds into an alpha-helix and contacts the DNA minor groove. The differences in response to effector for these proteins should be manifest in the polyelectrolyte effect which arises from cations displaced from DNA by interaction with positively charged side chains on a protein and is quantitated by measurement of DNA binding affinity as a function of ion concentration. Consistent with structural data for these proteins, high-affinity operator DNA binding by the PurR-purine complex involved approximately 15 ion pairs, a value significantly greater than that for the corresponding state of LacI (approximately 6 ion pairs). For both proteins, however, conversion to the low-affinity state results in a decrease of approximately 2-fold in the number of cations released per dimeric DNA binding site. Heat capacity changes (DeltaC(p)) that accompany DNA binding, derived from buried apolar surface area, coupled folding, and restriction of motional freedom of polar groups in the interface, also reflect the differences between these homologous repressor proteins. DNA binding of the PurR-guanine complex is accompanied by a DeltaC(p) (-2.8 kcal mol(-1) K(-1)) more negative than that observed previously for LacI (-0.9 to -1.5 kcal mol(-1) K(-1)), suggesting that more extensive protein folding and/or enhanced structural rigidity may occur upon DNA binding for PurR compared to DNA binding for LacI. The differences between these proteins illustrate plasticity of function despite high-level sequence and structural homology and undermine efforts to predict protein behavior on the basis of such similarities.

Mutations in the cytoplasmic tail of murine leukemia virus envelope protein suppress fusion inhibition by R peptide

During viral maturation, the cytoplasmic tail of the murine leukemia virus (MuLV) envelope (Env) protein undergoes proteolytic cleavage by the viral protease to release the 16-amino-acid R peptide, and this cleavage event activates the Env protein's fusion activity. We introduced Gly and/or Ser residues at different positions upstream of the R peptide in the cytoplasmic tail of the Friend MuLV Env protein and investigated their effects on fusion activity. Expression in HeLa T4 cells of a mutant Env protein with a single Gly insertion after I619, five amino acids upstream from the R peptide, induced syncytium formation with overlaid XC cells. Env proteins containing single or double Gly-Ser insertions after F614, 10 amino acids upstream from the R peptide, induced syncytium formation, and mutant proteins with multiple Gly insertions induced various levels of syncytium formation between HeLa T4 and XC cells. Immunoprecipitation and surface biotinylation assays showed that most of the mutants had surface expression levels comparable to those of the wild-type or R peptide-truncated Env proteins. Fluorescence dye redistribution assays also showed no hemifusion in the Env proteins which did not induce fusion. Our results indicate that insertion mutations in the cytoplasmic tail of the MuLV Env protein can suppress the inhibitory effect of the R peptide on membrane fusion and that there are differences in the effects of insertions in two regions in the cytoplasmic tail upstream of the R peptide.

The Lac repressor: a second generation of structural and functional studies

In the past year, the crystal structure of a dimeric version of the Escherichia coli Lac repressor bound to operator DNA was determined at 2.6A resolution, providing a closer view of the operator-bound conformation of the repressor. Refined NMR studies of the DNA-binding portion of the repressor complexed to operator DNA have revealed further details of the unique DNA-binding interactions of the repressor. The structural studies have been complemented by continued biochemical studies, with the overall goal of understanding the mechanism of allosteric regulation.

Operator DNA sequence variation enhances high affinity binding by hinge helix mutants of lactose repressor protein

The mechanism by which genetic regulatory proteins discern specific target DNA sequences remains a major area of inquiry. To explore in more detail the interplay between DNA and protein sequence, we have examined binding of variant lac operator DNA sequences to a series of mutant lactose repressor proteins (LacI). These proteins were altered in the C-terminus of the hinge region that links the N-terminal DNA binding and core sugar binding domains. Variant operators differed from the wild-type operator, O(1), in spacing and/or symmetry of the half-sites that contact the LacI N-terminal DNA binding domain. Binding of wild-type and mutant proteins was affected differentially by variations in operator sequence and symmetry. While the mutant series exhibits a 10(4)-fold range in binding affinity for O(1) operator, only a approximately 20-fold difference in affinity is observed for a completely symmetric operator, O(sym), used widely in studies of the LacI protein. Further, DNA sequence influenced allosteric response for these proteins. Binding of this LacI mutant series to other variant operator DNA sequences indicated the importance of symmetry-related bases, spacing, and the central base pair sequence in high affinity complex formation. Conformational flexibility in the DNA and other aspects of the structure influenced by the sequence may establish the binding environment for protein and determine both affinity and potential for allostery.