AlloRep: A Repository of Sequence, Structural and Mutagenesis Data for the LacI/GalR Transcription Regulators

Protein families evolve functional variation by accumulating point mutations at functionally important amino acid positions. Homologs in the LacI/GalR family of transcription regulators have evolved to bind diverse DNA sequences and allosteric regulatory molecules. In addition to playing key roles in bacterial metabolism, these proteins have been widely used as a model family for benchmarking structural and functional prediction algorithms. We have collected manually curated sequence alignments for >ﾠ3000 sequences, in vivo phenotypic and biochemical data for >ﾠ5750 LacI/GalR mutational variants, and noncovalent residue contact networks for 65 LacI/GalR homolog structures. Using this rich data resource, we compared the noncovalent residue contact networks of the LacI/GalR subfamilies to design and experimentally validate an allosteric mutant of a synthetic LacI/GalR repressor for use in biotechnology. The AlloRep database (freely available at www.AlloRep.org) is a key resource for future evolutionary studies of LacI/GalR homologs and for benchmarking computational predictions of functional change

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

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]

Engineering substrate promiscuity in halophilic alcohol dehydrogenase (HvADH2) by in silico design

An alcohol dehydrogenase from the halophilic archaeon Haloferax volcanii (HvADH2) has been engineered by rational design to broaden its substrate scope towards the conversion of a range of aromatic substrates, including flurbiprofenol, that is an intermediate of the non-steroidal anti-inflammatory drug, flurbiprofen. Wild-type HvADH2 showed minimal activity with flurbiprofenol (11.1 mU/mg). A homology model of HvADH2 was built and docking experiments with this substrate revealed that the biphenyl rings of flurbiprofenol formed strong interactions with residues F85 and F108, preventing its optimal binding in the active site. Mutations at position 85 however did not increase activity. Site directed mutagenesis at position F108 allowed the identification of three variants showing a significant (up to 2.3-fold) enhancement of activity towards flurbiprofenol, when compared to wild-type HvADH2. Interestingly, F108G variant did not show the classic inhibition in the presence of (R)-enantiomer when tested with rac-1-phenylethanol, underling its potential in racemic resolution of secondary alcohols

Cryo-EM structure of lysenin pore elucidates membrane insertion by an aerolysin family protein

Lysenin from the coelomic fluid of the earthworm Eisenia fetida belongs to the aerolysin family of small β-pore-forming toxins (β-PFTs), some members of which are pathogenic to humans and animals. Despite efforts, a high-resolution structure of a channel for this family of proteins has been elusive and therefore the mechanism of activation and membrane insertion remains unclear. Here we determine the pore structure of lysenin by single particle cryo-EM, to 3.1 Å resolution. The nonameric assembly reveals a long β-barrel channel spanning the length of the complex that, unexpectedly, includes the two pre-insertion strands flanking the hypothetical membrane-insertion loop. Examination of other members of the aerolysin family reveals high structural preservation in this region, indicating that the membrane-insertion pathway in this family is conserved. For some toxins, proteolytic activation and pro-peptide removal will facilitate unfolding of the pre-insertion strands, allowing them to form the β-barrel of the channel

Lysine120 Interactions with p53 Response Elements can Allosterically Direct p53 Organization

p53 can serve as a paradigm in studies aiming to figure out how allosteric perturbations in transcription factors (TFs) triggered by small changes in DNA response element (RE) sequences, can spell selectivity in co-factor recruitment. p53-REs are 20-base pair (bp) DNA segments specifying diverse functions. They may be located near the transcription start sites or thousands of bps away in the genome. Their number has been estimated to be in the thousands, and they all share a common motif. A key question is then how does the p53 protein recognize a particular p53-RE sequence among all the similar ones? Here, representative p53-REs regulating diverse functions including cell cycle arrest, DNA repair, and apoptosis were simulated in explicit solvent. Among the major interactions between p53 and its REs involving Lys120, Arg280 and Arg248, the bps interacting with Lys120 vary while the interacting partners of other residues are less so. We observe that each p53-RE quarter site sequence has a unique pattern of interactions with p53 Lys120. The allosteric, DNA sequence-induced conformational and dynamic changes of the altered Lys120 interactions are amplified by the perturbation of other p53-DNA interactions. The combined subtle RE sequence-specific allosteric effects propagate in the p53 and in the DNA. The resulting amplified allosteric effects far away are reflected in changes in the overall p53 organization and in the p53 surface topology and residue fluctuations which play key roles in selective co-factor recruitment. As such, these observations suggest how similar p53-RE sequences can spell the preferred co-factor binding, which is the key to the selective gene transactivation and consequently different functional effects

Investigating Homology between Proteins using Energetic Profiles

Accumulated experimental observations demonstrate that protein stability is often preserved upon conservative point mutation. In contrast, less is known about the effects of large sequence or structure changes on the stability of a particular fold. Almost completely unknown is the degree to which stability of different regions of a protein is generally preserved throughout evolution. In this work, these questions are addressed through thermodynamic analysis of a large representative sample of protein fold space based on remote, yet accepted, homology. More than 3,000 proteins were computationally analyzed using the structural-thermodynamic algorithm COREX/BEST. Estimated position-specific stability (i.e., local Gibbs free energy of folding) and its component enthalpy and entropy were quantitatively compared between all proteins in the sample according to all-vs.-all pairwise structural alignment. It was discovered that the local stabilities of homologous pairs were significantly more correlated than those of non-homologous pairs, indicating that local stability was indeed generally conserved throughout evolution. However, the position-specific enthalpy and entropy underlying stability were less correlated, suggesting that the overall regional stability of a protein was more important than the thermodynamic mechanism utilized to achieve that stability. Finally, two different types of statistically exceptional evolutionary structure-thermodynamic relationships were noted. First, many homologous proteins contained regions of similar thermodynamics despite localized structure change, suggesting a thermodynamic mechanism enabling evolutionary fold change. Second, some homologous proteins with extremely similar structures nonetheless exhibited different local stabilities, a phenomenon previously observed experimentally in this laboratory. These two observations, in conjunction with the principal conclusion that homologous proteins generally conserved local stability, may provide guidance for a future thermodynamically informed classification of protein homology