Functional impact of polar and acidic substitutions in the lactose repressor hydrophobic monomer.monomer interface with a buried lysine

Rheostat functional outcomes occur when substitutions are introduced at nonconserved positions that diverge with speciation

When amino acids vary during evolution, the outcome can be functionally neutral or biologically-important. We previously found that substituting a subset of nonconserved positions, "rheostat" positions, can have surprising effects on protein function. Since changes at rheostat positions can facilitate functional evolution or cause disease, more examples are needed to understand their unique biophysical characteristics. Here, we explored whether "phylogenetic" patterns of change in multiple sequence alignments (such as positions with subfamily specific conservation) predict the locations of functional rheostat positions. To that end, we experimentally tested eight phylogenetic positions in human liver pyruvate kinase (hLPYK), using 10-15 substitutions per position and biochemical assays that yielded five functional parameters. Five positions were strongly rheostatic and three were non-neutral. To test the corollary that positions with low phylogenetic scores were not rheostat positions, we combined these phylogenetic positions with previously-identified hLPYK rheostat, "toggle" (most substitution abolished function), and "neutral" (all substitutions were like wild-type) positions. Despite representing 428 variants, this set of 33 positions was poorly statistically powered. Thus, we turned to the in vivo phenotypic dataset for E. coli lactose repressor protein (LacI), which comprised 12-13 substitutions at 329 positions and could be used to identify rheostat, toggle, and neutral positions. Combined hLPYK and LacI results show that positions with strong phylogenetic patterns of change are more likely to exhibit rheostat substitution outcomes than neutral or toggle outcomes. Furthermore, phylogenetic patterns were more successful at identifying rheostat positions than were co-evolutionary or eigenvector centrality measures of evolutionary change.

What Mutagenesis Can and Cannot Reveal About Allostery

Allosteric regulation of protein function is recognized to be widespread throughout biology; however, knowledge of allosteric mechanisms, the molecular changes within a protein that couple one binding site to another, is limited. Although mutagenesis is often used to probe allosteric mechanisms, we consider herein what the outcome of a mutagenesis study truly reveals about an allosteric mechanism. Arguably, the best way to evaluate the effects of a mutation on allostery is to monitor the allosteric coupling constant (Qax), a ratio of the substrate binding constants in the absence versus presence of an allosteric effector. A range of substitutions at a given residue position in a protein can reveal when a particular substitution causes gain-of-function, which addresses a key challenge in interpreting mutation-dependent changes in the magnitude of Qax. Thus, whole-protein mutagenesis studies offer an acceptable means of identifying residues that contribute to an allosteric mechanism. With this focus on monitoring Qax, and keeping in mind the equilibrium nature of allostery, we consider alternative possibilities for what an allosteric mechanism might be. We conclude that different possible mechanisms (rotation-of-solid-domains, movement of secondary structure, side-chain repacking, changes in dynamics, etc.) will result in different findings in whole-protein mutagenesis studies.

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].

Optimized expression and purification of biophysical quantities of Lac repressor and Lac repressor regulatory domain

The recombinant production of Lac repressor (LacI) in Escherichia coli is complicated by its ubiquitous use as a regulatory element in commercially-available expression vectors and host strains. While LacI-regulated expression systems are often used to produce recombinant LacI, the product can be heterogeneous and unsuitable for some studies. Alternative approaches include using unregulated vectors which typically suffer from low yield or vectors with promoters induced by metabolically active sugars which can dilute isotope labels necessary for certain biophysical studies. Here, an optimized expression system and isolation protocol for producing various constructs of LacI is introduced which eliminates these complications. The expression vector is an adaptation of the pASK backbone wherein expression of the lacI gene is regulated by an anhydrotetracyline inducible tetA promoter and the host strain lacks the lacI gene. Typical yields in highly deuterated minimal medium are nearly 2-fold greater than those previously reported. Notably, the new expression system is also able to produce the isolated regulatory domain of LacI without co-expression of the full-length protein and without any defects in cell viability, eliminating the inconvenient requirement for strict monitoring of cell densities during pre-culturing. Typical yields in highly deuterated minimal medium are significantly greater than those previously reported. Characterization by solution NMR shows that LacI constructs produced using this expression system are highly homogenous and functionally active.

Sequence dependence of transcription factor-mediated DNA looping

DNA is subject to large deformations in a wide range of biological processes. Two key examples illustrate how such deformations influence the readout of the genetic information: the sequestering of eukaryotic genes by nucleosomes and DNA looping in transcriptional regulation in both prokaryotes and eukaryotes. These kinds of regulatory problems are now becoming amenable to systematic quantitative dissection with a powerful dialogue between theory and experiment. Here, we use a single-molecule experiment in conjunction with a statistical mechanical model to test quantitative predictions for the behavior of DNA looping at short length scales and to determine how DNA sequence affects looping at these lengths. We calculate and measure how such looping depends upon four key biological parameters: the strength of the transcription factor binding sites, the concentration of the transcription factor, and the length and sequence of the DNA loop. Our studies lead to the surprising insight that sequences that are thought to be especially favorable for nucleosome formation because of high flexibility lead to no systematically detectable effect of sequence on looping, and begin to provide a picture of the distinctions between the short length scale mechanics of nucleosome formation and looping.

Structure of myostatin·follistatin-like 3: N-terminal domains of follistatin-type molecules exhibit alternate modes of binding

TGF-β family ligands are involved in a variety of critical physiological processes. For instance, the TGF-β ligand myostatin is a staunch negative regulator of muscle growth and a therapeutic target for muscle-wasting disorders. Therefore, it is important to understand the molecular mechanisms of TGF-β family regulation. One form of regulation is through inhibition by extracellular antagonists such as the follistatin (Fst)-type proteins. Myostatin is tightly controlled by Fst-like 3 (Fstl3), which is the only Fst-type molecule that has been identified in the serum bound to myostatin. Here, we present the crystal structure of myostatin in complex with Fstl3. The structure reveals that the N-terminal domain (ND) of Fstl3 interacts uniquely with myostatin as compared with activin A, because it utilizes different surfaces on the ligand. This results in conformational differences in the ND of Fstl3 that alter its position in the type I receptor-binding site of the ligand. We also show that single point mutations in the ND of Fstl3 are detrimental to ligand binding, whereas corresponding mutations in Fst have little effect. Overall, we have shown that the NDs of Fst-type molecules exhibit distinctive modes of ligand binding, which may affect overall affinity of ligand·Fst-type protein complexes.

Positions 94-98 of the lactose repressor N-subdomain monomer-monomer interface are critical for allosteric communication

The central region of the LacI N-subdomain monomer-monomer interface includes residues K84, V94, V95, V96, S97, and M98. The side chains of these residues line the β-strands at this interface and interact to create a network of hydrophobic, charged, and polar interactions that significantly rearranges in different functional states of LacI. Prior work showed that converting K84 to an apolar residue or converting V96 to an acidic residue impedes the allosteric response to inducer. Thus, we postulated that a disproportionate number of substitutions in this region of the monomer-monomer interface would alter the complex features of the LacI allosteric response. To explore this hypothesis, acidic, basic, polar, and apolar mutations were introduced at positions 94-98. Despite their varied locations along the β-strands that flank the interface, ∼70% of the mutations impact allosteric behavior, with the most significant effects found for charged substitutions. Of note, many of the LacI variants with minor functional impact exhibited altered stability to urea denaturation. The results confirm the critical role of amino acids 94-98 and indicate that this N-subdomain interface forms a primary pathway in LacI allosteric response.

Tetramer opening in LacI-mediated DNA looping

Lactose repressor protein (LacI) controls transcription of the genes involved in lactose metabolism in bacteria. Essential to optimal LacI-mediated regulation is its ability to bind simultaneously to two operators, forming a loop on the intervening DNA. Recently, several lines of evidence (both theoretical and experimental) have suggested various possible loop structures associated with different DNA binding topologies and LacI tetramer structural conformations (adopted by flexing about the C-terminal tetramerization domain). We address, specifically, the role of protein opening in loop formation by employing the single-molecule tethered particle motion method on LacI protein mutants chemically cross-linked at different positions along the cleft between the two dimers. Measurements on the wild-type and uncross-linked LacI mutants led to the observation of two distinct levels of short tether length, associated with two different DNA looping structures. Restricting conformational flexibility of the protein by chemical cross-linking induces pronounced effects. Crosslinking the dimers at the level of the N-terminal DNA binding head (E36C) completely suppresses looping, whereas cross-linking near the C-terminal tetramerization domain (Q231C) results in changes of looping geometry detected by the measured tether length distributions. These observations lead to the conclusion that tetramer opening plays a definite role in at least a subset of LacI/DNA loop conformations.

Allostery in the LacI/GalR family: variations on a theme

The lactose repressor protein (LacI) was among the very first genetic regulatory proteins discovered, and more than 1000 members of the bacterial LacI/GalR family are now identified. LacI has been the prototype for understanding how transcription is controlled using small metabolites to modulate protein association with specific DNA sites. This understanding has been greatly expanded by the study of other LacI/GalR homologues. A general picture emerges in which the conserved fold provides a scaffold for multiple types of interactions - including oligomerization, small molecule binding, and protein-protein binding - that in turn influence target DNA binding and thereby regulate mRNA production. Although many different functions have evolved from this basic scaffold, each homologue retains functional flexibility: For the same protein, different small molecules can have disparate impact on DNA binding and hence transcriptional outcome. In turn, binding to alternative DNA sequences may impact the degree of allosteric response. Thus, this family exhibits a symphony of variations by which transcriptional control is achieved.