Glutamine, alanine or glycine repeats inserted into the loop of a protein have minimal effects on stability and folding rates

The bifunctional 6-phosphofructokinase-2/fructose-2,6-bisphosphatase from the shrimp Litopenaeus vannamei: Molecular characterization and down-regulation of expression in response to severe hypoxia

Hypoxia is a frequent stressor in marine environments with multiple adverse effects on marine species. The white shrimp Litopenaeus vannamei withstands hypoxic conditions by activating anaerobic metabolism with tissue-specific changes in glycolytic and gluconeogenic enzymes. In animal cells, glycolytic/gluconeogenic fluxes are highly controlled by the levels of fructose-2,6-bisphosphate (F-2,6-P₂), a signal metabolite synthesized and degraded by the bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2). PFK-2/FBPase-2 has been studied in vertebrates and some invertebrates, but as far as we know, there are no reports on PFK-2/FBPase-2 from crustaceans. In the present work, we obtained cDNA nucleotide sequences corresponding to two mRNAs for PFK-2/FBPase-2 and named them PFKFBP1 (1644 bp) and PFKFBP2 (1566 bp), from the white shrimp L. vannamei. The deduced PFKFBP1 and PFKFBP2 are 547 and 521 amino acids long, respectively. Both proteins share 99.23% of identity, and only differ in 26 additional amino acids present in the kinase domain of the PFKFBP1. The kinase and phosphatase domains are highly conserved in sequence and structure between both isoforms and other proteins from diverse taxa. Total expression of PFKFBP1-2 is tissue-specific, more abundant in gills than in hepatopancreas and undetectable in muscle. Moreover, severe hypoxia (1 mg/L of DO) decreased expression of PFKFBP1-2 in gills while anaerobic glycolysis was induced, as indicated by accumulation of cellular lactate. These results suggest that negative regulation of PFKFBP1-2 at expression level is necessary to set up anaerobic glycolysis in the cells during the response to hypoxia.

Folding in Place: Design of β-Strap Motifs to Stabilize the Folding of Hairpins with Long Loops

Despite their pivotal role in defining antibody affinity and protein function, β-hairpins harboring long noncanonical loops remain synthetically challenging because of the large entropic penalty associated with their conformational folding. Little is known about the contribution and impact of stabilizing motifs on the folding of β-hairpins with loops of variable length and plasticity. Here, we report a design of minimalist β-straps (strap = strand + cap) that offset the entropic cost of long-loop folding. The judicious positioning of noncovalent interactions (hydrophobic cluster and salt-bridge) within the novel 8-mer β-strap design RW(V/H)W···WVWE stabilizes hairpins with up to 10-residue loops of varying degrees of plasticity (Tm up to 52 °C; 88 ± 1% folded at 18 °C). This "hyper" thermostable β-strap outperforms the previous gold-standard technology of β-strand-β-cap (16-mer) and provides a foundation for producing new classes of long hairpins as a viable and practical alternative to macrocyclic peptides.

Testing the length limit of loop grafting in a helical repeat protein

Alpha-helical repeat proteins such as consensus-designed tetratricopeptide repeats (CTPRs) are exceptionally stable molecules that are able to tolerate destabilizing sequence alterations and are therefore becoming increasingly valued as a modular platform for biotechnology and biotherapeutic applications. A simple approach to functionalize the CTPR scaffold that we are pioneering is the insertion of short linear motifs (SLiMs) into the loops between adjacent repeats. Here, we test the limits of the scaffold by inserting 17 highly diverse amino acid sequences of up to 58 amino acids in length into a two-repeat protein and examine the impact on protein folding, stability and solubility. The sequences include three SLiMs that bind oncoproteins and eleven naturally occurring linker sequences all predicted to be intrinsically disordered but with conformational preferences ranging from compact globules to expanded coils. We show that the loop-grafted proteins retain the native CTPR structure and are thermally stable with melting temperatures above 60 ​°C, despite the longest loop sequence being almost the same size as the CTPR scaffold itself (68 amino acids). Although the main determinant of the effect of stability was found to be loop length and was relatively insensitive to amino acid composition, the relationship between protein solubility and the loop sequences was more complex, with the presence of negatively charged amino acids enhancing the solubility. Our findings will help us to fully realize the potential of the repeat-protein scaffold, allowing a rational design approach to create artificial modular proteins with customized functional capabilities.

Benefits and constrains of covalency: the role of loop length in protein stability and ligand binding

Protein folding is governed by non-covalent interactions under the benefits and constraints of the covalent linkage of the backbone chain. In the current work we investigate the influence of loop length variation on the free energies of folding and ligand binding in a small globular single-domain protein containing two EF-hand subdomains—calbindin D_9k. We introduce a linker extension between the subdomains and vary its length between 1 to 16 glycine residues. We find a close to linear relationship between the linker length and the free energy of folding of the Ca²⁺-free protein. In contrast, the linker length has only a marginal effect on the Ca²⁺ affinity and cooperativity. The variant with a single-glycine extension displays slightly increased Ca²⁺ affinity, suggesting that the slightly extended linker allows optimized packing of the Ca²⁺-bound state. For the extreme case of disconnected subdomains, Ca²⁺ binding becomes coupled to folding and assembly. Still, a high affinity between the EF-hands causes the non-covalent pair to retain a relatively high apparent Ca²⁺ affinity. Our results imply that loop length variation could be an evolutionary option for modulating properties such as protein stability and turnover without compromising the energetics of the specific function of the protein.

Decoupling a tandem-repeat protein: Impact of multiple loop insertions on a modular scaffold

The simple topology and modular architecture of tandem-repeat proteins such as tetratricopeptide repeats (TPRs) and ankyrin repeats makes them straightforward to dissect and redesign. Repeat-protein stability can be manipulated in a predictable way using site-specific mutations. Here we explore a different type of modification - loop insertion - that will enable a simple route to functionalisation of this versatile scaffold. We previously showed that a single loop insertion has a dramatically different effect on stability depending on its location in the repeat array. Here we dissect this effect by a combination of multiple and alternated loop insertions to understand the origins of the context-dependent loss in stability. We find that the scaffold is remarkably robust in that its overall structure is maintained. However, adjacent repeats are now only weakly coupled, and consequently the increase in solvent protection, and thus stability, with increasing repeat number that defines the tandem-repeat protein class is lost. Our results also provide us with a rulebook with which we can apply these principles to the design of artificial repeat proteins with precisely tuned folding landscapes and functional capabilities, thereby paving the way for their exploitation as a versatile and truly modular platform in synthetic biology.

Context-Dependent Energetics of Loop Extensions in a Family of Tandem-Repeat Proteins

Consensus-designed tetratricopeptide repeat proteins are highly stable, modular proteins that are strikingly amenable to rational engineering. They therefore have tremendous potential as building blocks for biomaterials and biomedicine. Here, we explore the possibility of extending the loops between repeats to enable further diversification, and we investigate how this modification affects stability and folding cooperativity. We find that extending a single loop by up to 25 residues does not disrupt the overall protein structure, but, strikingly, the effect on stability is highly context-dependent: in a two-repeat array, destabilization is relatively small and can be accounted for purely in entropic terms, whereas extending a loop in the middle of a large array is much more costly because of weakening of the interaction between the repeats. Our findings provide important and, to our knowledge, new insights that increase our understanding of the structure, folding, and function of natural repeat proteins and the design of artificial repeat proteins in biotechnology.

An NMR Confirmation for Increased Folded State Entropy Following Loop Truncation

Previous studies conducted on flexible loop regions in proteins revealed that the energetic consequences of changing loop length predominantly arise from the entropic cost of ordering a loop during folding. However, in an earlier study of human acylphosphatase (hmAcP) using experimental and computational approaches, we showed that thermodynamic stabilization upon loop truncation can be attributed mainly to the increased entropy of the folded state. Here, using 15N NMR spectroscopy, we studied the effect of loop truncation on hmAcP backbone dynamics on the picosecond-nanosecond timescale with the aim of confirming the effect of folded state entropy on protein stability. NMR-relaxation-derived N-H squared generalized order parameters reveal that loop truncation results in a significant increase in protein conformational flexibility. Comparison of these results with previously acquired all-atom molecular dynamics simulation, analyzed here in terms of squared generalized NMR order parameters, demonstrates general agreement between the two methods. The NMR study not only provides direct evidence for the enhanced conformational entropy of the folded state of hmAcP upon loop truncation but also gives a quantitative measure of the observed effects.

The Cost of Long Catalytic Loops in Folding and Stability of the ALS-Associated Protein SOD1

A conspicuous feature of the amyotrophic lateral sclerosis (ALS)-associated protein SOD1 is that its maturation into a functional enzyme relies on local folding of two disordered loops into a catalytic subdomain. To drive the disorder-to-order transition, the protein employs a single Zn²⁺ ion. The question is then if the entropic penalty of maintaining such disordered loops in the immature apoSOD1 monomer is large enough to explain its unusually low stability, slow folding, and pathological aggregation in ALS. To find out, we determined the effects of systematically altering the SOD1-loop lengths by protein redesign. The results show that the loops destabilize the apoSOD1 monomer by ∼3 kcal/mol, rendering the protein marginally stable and accounting for its aggregation behavior. Yet the effect on the global folding kinetics remains much smaller with a transition-state destabilization of <1 kcal/mol. Notably, this 1/3 transition-state to folded-state stability ratio provides a clear-cut example of the enigmatic disagreement between the Leffler α value from loop-length alterations (typically 1/3) and the "standard" reaction coordinates based on solvent perturbations (typically >2/3). Reconciling the issue, we demonstrate that the disagreement disappears when accounting for the progressive loop shortening that occurs along the folding pathway. The approach assumes a consistent Flory loop entropy scaling factor of c = 1.48 for both equilibrium and kinetic data and has the added benefit of verifying the tertiary interactions of the folding nucleus as determined by phi-value analysis. Thus, SOD1 not only represents a case where evolution of key catalytic function has come with the drawback of a destabilized apo state but also stands out as a well-suited model system for exploring the physicochemical details of protein self-organization.

Exploring Protein Cavities through Rigidity Analysis

The geometry of cavities in the surfaces of proteins facilitates a variety of biochemical functions. To better understand the biochemical nature of protein cavities, the shape, size, chemical properties, and evolutionary nature of functional and nonfunctional surface cavities have been exhaustively surveyed in protein structures. The rigidity of surface cavities, however, is not immediately available as a characteristic of structure data, and is thus more difficult to examine. Using rigidity analysis for assessing and analyzing molecular rigidity, this paper performs the first survey of the relationships between cavity properties, such as size and residue content, and how they correspond to cavity rigidity. Our survey measured a variety of rigidity metrics on 120,323 cavities from 12,785 sequentially non-redundant protein chains. We used VASP-E, a volume-based algorithm for analyzing cavity geometry. Our results suggest that rigidity properties of protein cavities are dependent on cavity surface area.

Control of serine integrase recombination directionality by fusion with the directionality factor

Bacteriophage serine integrases are extensively used in biotechnology and synthetic biology for assembly and rearrangement of DNA sequences. Serine integrases promote recombination between two different DNA sites, attP and attB, to form recombinant attL and attR sites. The 'reverse' reaction requires another phage-encoded protein called the recombination directionality factor (RDF) in addition to integrase; RDF activates attL × attR recombination and inhibits attP × attB recombination. We show here that serine integrases can be fused to their cognate RDFs to create single proteins that catalyse efficient attL × attR recombination in vivo and in vitro, whereas attP × attB recombination efficiency is reduced. We provide evidence that activation of attL × attR recombination involves intra-subunit contacts between the integrase and RDF moieties of the fusion protein. Minor changes in the length and sequence of the integrase-RDF linker peptide did not affect fusion protein recombination activity. The efficiency and single-protein convenience of integrase-RDF fusion proteins make them potentially very advantageous for biotechnology/synthetic biology applications. Here, we demonstrate efficient gene cassette replacement in a synthetic metabolic pathway gene array as a proof of principle.

Mutation Sensitivity Maps: Identifying Residue Substitutions That Impact Protein Structure Via a Rigidity Analysis In Silico Mutation Approach

Understanding how an amino acid substitution affects a protein's structure can aid in the design of pharmaceutical drugs that aim at countering diseases caused by protein mutants. Unfortunately, performing even a few amino acid substitutions in vitro is both time and cost prohibitive, whereas an exhaustive analysis that involves systematically mutating all amino acids in the physical protein is infeasible. Computational methods have been developed to predict the effects of mutations, but even many of them are computationally intensive or are else dependent on homology or experimental data that may not be available for the protein being studied. In this work, we motivate and present a computation pipeline whose only input is a Protein Data Bank file containing the 3D coordinates of the atoms of a biomolecule. Our high-throughput approach uses our ProMuteHT algorithm to exhaustively generate in silico amino acid substitutions at each residue, and it also includes an energy minimization option. This is in contrast to our previous work, where we analyzed the effects of in silico mutations to Alanine, Serine, and Glycine only. We exploit the speed of a fast rigidity analysis approach to analyze our protein variants, and develop a Mutation Sensitivity (MuSe) Map, to permit identifying residues that are most sensitive to mutations. We present a case study to show the degree to which a MuSe Map and whisker plots are able to locate amino acids whose mutations most affect a protein's structure as inferred from a rigidity analysis approach.

Entropic stabilization of a deubiquitinase provides conformational plasticity and slow unfolding kinetics beneficial for functioning on the proteasome

Human ubiquitin C-terminal hydrolyase UCH-L5 is a topologically knotted deubiquitinase that is activated upon binding to the proteasome subunit Rpn13. The length of its intrinsically disordered cross-over loop is essential for substrate recognition. Here, we showed that the catalytic domain of UCH-L5 exhibits higher equilibrium folding stability with an unfolding rate on the scale of 10⁻⁸ s⁻¹, over four orders of magnitudes slower than its paralogs, namely UCH-L1 and -L3, which have shorter cross-over loops. NMR relaxation dynamics analysis confirmed the intrinsic disorder of the cross-over loop. Hydrogen deuterium exchange analysis further revealed a positive correlation between the length of the cross-over loop and the degree of local fluctuations, despite UCH-L5 being thermodynamically and kinetically more stable than the shorter UCHs. Considering the role of UCH-L5 in removing K48-linked ubiquitin to prevent proteasomal degradation of ubiquitinated substrates, our findings offered mechanistic insights into the evolution of UCH-L5. Compared to its paralogs, it is entropically stabilized to withstand mechanical unfolding by the proteasome while maintaining structural plasticity. It can therefore accommodate a broad range of substrate geometries at the cost of unfavourable entropic loss.

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all-atom molecular dynamics simulations to study how gradual shortening a very long or solvent-exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability.

Affinity transfer to the archaeal extremophilic Sac7d protein by insertion of a CDR

Artificially transforming a scaffold protein into binders often consists of introducing diversity into its natural binding region by directed mutagenesis. We have previously developed the archaeal extremophilic Sac7d protein as a scaffold to derive affinity reagents (Affitins) by randomization of only a flat surface, or a flat surface and two short loops with natural lengths. Short loops are believed to contribute to stability of extremophilic proteins, and loop extension has been reported detrimental for the thermal and chemical stabilities of mesophilic proteins. In this work, we wanted to evaluate the possibility of designing target-binding proteins based on Sac7d by using a complementary determining region (CDR). To this aim, we inserted into three different loops a 10 residues CDR from the cAb-Lys3 anti-lysozyme camel antibody. The chimeras obtained were as stable as wild-type (WT) Sac7d at extreme pH and their structural integrity was supported. Chimeras were thermally stable, but with T(m)s from 60.9 to 66.3°C (cf. 91°C for Sac7d) which shows that loop extension is detrimental for thermal stability of Sac7d. The loop 3 enabled anti-lysozyme activity. These results pave the way for the use of CDR(s) from antibodies and/or extended randomized loop(s) to increase the potential of binding of Affitins.

Stabilization of a protein conferred by an increase in folded state entropy

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

A method for the incremental expansion of polyglutamine repeats in recombinant proteins

The polyglutamine diseases are caused by the expansion of CAG repeats. A key step in understanding the disease mechanisms, at the DNA and protein level, is the ability to produce recombinant proteins with specific length glutamine tracts which is a time-consuming first step in setting up in vitro systems to study the effects of polyglutamine expansion. Described here is a PCR-based method for the amplification of CAG repeats, which we used to incrementally extend CAG length by 3-5 repeats per cycle. This method could be translated into various contexts where amplification of repeating elements is necessary.

Protein folding: adding a nucleus to guide helix docking reduces landscape roughness

The elongated three-helix‐bundle spectrin domains R16 and R17 fold and unfold unusually slowly over a rough energy landscape, in contrast to the homologue R15, which folds fast over a much smoother, more typical landscape. R15 folds via a nucleation–condensation mechanism that guides the docking of the A and C-helices. However, in R16 and R17, the secondary structure forms first and the two helices must then dock in the correct register. Here, we use variants of R16 and R17 to demonstrate that substitution of just five key residues is sufficient to alter the folding mechanism and reduce the landscape roughness. We suggest that, by providing access to an alternative, faster, folding route over their landscape, R16 and R17 can circumvent their slow, frustrated wild-type folding mechanism.

Energy landscape of knotted protein folding

Recent experiments have conclusively shown that proteins are able to fold from an unknotted, denatured polypeptide to the knotted, native state without the aid of chaperones. These experiments are consistent with a growing body of theoretical work showing that a funneled, minimally frustrated energy landscape is sufficient to fold small proteins with complex topologies. Here, we present a theoretical investigation of the folding of a knotted protein, 2ouf, engineered in the laboratory by a domain fusion that mimics an evolutionary pathway for knotted proteins. Unlike a previously studied knotted protein of similar length, we see reversible folding/knotting and a surprising lack of deep topological traps with a coarse-grained structure-based model. Our main interest is to investigate how evolution might further select the geometry and stiffness of the threading region of the newly fused protein. We compare the folding of the wild-type protein to several mutants. Similarly to the wild-type protein, all mutants show robust and reversible folding, and knotting coincides with the transition state ensemble. As observed experimentally, our simulations show that the knotted protein folds about ten times slower than an unknotted construct with an identical contact map. Simulated folding kinetics reflect the experimentally observed rollover in the folding limbs of chevron plots. Successful folding of the knotted protein is restricted to a narrow range of temperature as compared to the unknotted protein and fits of the kinetic folding data below folding temperature suggest slow, nondiffusive dynamics for the knotted protein.

Chimeric glutathione S-transferases containing inserts of kininogen peptides: potential novel protein therapeutics

The study of synthetic peptides corresponding to discrete regions of proteins has facilitated the understanding of protein structure-activity relationships. Short peptides can also be used as powerful therapeutic agents. However, in many instances, small peptides are prone to rapid degradation or aggregation and may lack the conformation required to mimic the functional motifs of the protein. For peptides to function as pharmacologically active agents, efficient production or expression, high solubility, and retention of biological activity through purification and storage steps are required. We report here the design, expression, and functional analysis of eight engineered GST proteins (denoted GSHKTs) in which peptides ranging in size from 8 to 16 amino acids and derived from human high molecular weight kininogen (HK) domain 5 were inserted into GST (between Gly-49 and Leu-50). Peptides derived from HK are known to inhibit cell proliferation, angiogenesis, and tumor metastasis, and the biological activity of the HK peptides was dramatically (>50-fold) enhanced following insertion into GST. GSHKTs are soluble and easily purified from Escherichia coli by affinity chromatography. Functionally, these hybrid proteins cause inhibition of endothelial cell proliferation. Crystallographic analysis of GSHKT10 and GSHKT13 (harboring 10- and 13-residue HK peptides, respectively) showed that the overall GST structure was not perturbed. These results suggest that the therapeutic efficacy of short peptides can be enhanced by insertion into larger proteins that are easily expressed and purified and that GST may potentially be used as such a carrier.

Amyloid-like fibril formation by polyQ proteins: a critical balance between the polyQ length and the constraints imposed by the host protein

Nine neurodegenerative disorders, called polyglutamine (polyQ) diseases, are characterized by the formation of intranuclear amyloid-like aggregates by nine proteins containing a polyQ tract above a threshold length. These insoluble aggregates and/or some of their soluble precursors are thought to play a role in the pathogenesis. The mechanism by which polyQ expansions trigger the aggregation of the relevant proteins remains, however, unclear. In this work, polyQ tracts of different lengths were inserted into a solvent-exposed loop of the β-lactamase BlaP and the effects of these insertions on the properties of BlaP were investigated by a range of biophysical techniques. The insertion of up to 79 glutamines does not modify the structure of BlaP; it does, however, significantly destabilize the enzyme. The extent of destabilization is largely independent of the polyQ length, allowing us to study independently the effects intrinsic to the polyQ length and those related to the structural integrity of BlaP on the aggregating properties of the chimeras. Only chimeras with 55Q and 79Q readily form amyloid-like fibrils; therefore, similarly to the proteins associated with diseases, there is a threshold number of glutamines above which the chimeras aggregate into amyloid-like fibrils. Most importantly, the chimera containing 79Q forms amyloid-like fibrils at the same rate whether BlaP is folded or not, whereas the 55Q chimera aggregates into amyloid-like fibrils only if BlaP is unfolded. The threshold value for amyloid-like fibril formation depends, therefore, on the structural integrity of the β-lactamase moiety and thus on the steric and/or conformational constraints applied to the polyQ tract. These constraints have, however, no significant effect on the propensity of the 79Q tract to trigger fibril formation. These results suggest that the influence of the protein context on the aggregating properties of polyQ disease-associated proteins could be negligible when the latter contain particularly long polyQ tracts.

Location trumps length: polyglutamine-mediated changes in folding and aggregation of a host protein

Expanded CAG diseases are progressive neurodegenerative disorders in which specific proteins have an unusually long polyglutamine stretch. Although these proteins share no other sequence or structural homologies, they all aggregate into intracellular inclusions that are believed to be pathological. We sought to determine what impact the position and number of glutamines have on the structure and aggregation of the host protein, apomyoglobin. Variable-length polyQ tracts were inserted either into the loop between the C- and D-helices (Q(n)CD) or at the N-terminus (Q(n)NT). The Q(n)CD mutants lost some α-helix and gained unordered and/or β-sheet in a length-dependent manner. These mutants were partially unfolded and rapidly assembled into soluble chain-like oligomers. In sharp contrast, the Q(n)NT mutants largely retained wild-type tertiary structure but associated into long, fibrillar aggregates. Control proteins with glycine-serine repeats (GS(8)CD and GS(8)NT) were produced. GS(8)CD exhibited similar structural perturbations and aggregation characteristics to an analogously sized Q(16)CD, indicating that the observed effects are independent of amino acid composition. In contrast to Q(16)NT, GS(8)NT did not form fibrillar aggregates. Thus, soluble oligomers are produced through structural perturbation and do not require polyQ, whereas classic fibrils arise from specific polyQ intermolecular interactions in the absence of misfolding.

Introduction of glutamines into the B2-H2 loop promotes prion protein conversion

In prion diseases cellular prion protein (PrP(C)) undergoes conformational transition into the β-sheet-rich form (PrP(Sc)). PrP(C) consists of the disordered N-terminal part and a C-terminal globular domain containing three α-helices (H1, H2, H3) and an antiparallel beta sheet (B1, B2). B2-H2 loop, which has a focal role in the species barrier, contains the highest density of asparagine (N) and glutamine (Q) residues in the whole sequence. Q/N-rich domains are essential for the conversion of yeast prions. We investigated the role of Q/N residues in the B2-H2 loop in PrP conversion. We prepared mouse PrP mutants with increasing number of consecutive Q/N residues in the B2-H2 loop. Stability of the mutants decreased with the increasing number of inserted glutamines. In vitro conversion of mutants yielded fibrils of similar morphology as the wild-type PrP. Q/N mutants accelerated fibrillization in comparison to the wild-type PrP, with mutant containing the most glutamines having the shortest lag phase. The effect of Q/N residues was specific for the B2-H2 loop and was not due to simple increase in flexibility as the introduction of Gly-Ser or Ala residues slowed the conversion despite their decreased stability. Our results thus suggest that Q/N residues in the B2-H2 loop of PrP promote protein conversion and may represent a link to conversion of Q/N-rich prions.

Molecular simulations of mutually exclusive folding in a two-domain protein switch

A major challenge with testing designs of protein conformational switches is the need for experimental probes that can independently monitor their individual protein domains. One way to circumvent this issue is to use a molecular simulation approach in which each domain can be directly observed. Here we report what we believe to be the first molecular simulations of mutually exclusive folding in an engineered two-domain protein switch, providing a direct view of how folding of one protein drives unfolding of the other in a barnase-ubiquitin fusion protein. These simulations successfully capture the experimental effects of interdomain linker length and ligand binding on the extent of unfolding in the less stable domain. In addition, the effect of linker length on the potential for oligomerization, which eliminates switch activity, is in qualitative agreement with analytical ultracentrifugation experiments. We also perform what we believe to be the first study of protein unfolding via progressive localized compression. Finally, we are able to explore the kinetics of mutually exclusive folding by determining the effect of linker length on rates of unfolding and refolding of each protein domain. Our results demonstrate that molecular simulations can provide seemingly novel biological insights on the behavior of individual protein domains, thereby aiding in the rational design of bifunctional switches.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

The construction and use of log-odds substitution scores for multiple sequence alignment

Most pairwise and multiple sequence alignment programs seek alignments with optimal scores. Central to defining such scores is selecting a set of substitution scores for aligned amino acids or nucleotides. For local pairwise alignment, substitution scores are implicitly of log-odds form. We now extend the log-odds formalism to multiple alignments, using Bayesian methods to construct "BILD" ("Bayesian Integral Log-odds") substitution scores from prior distributions describing columns of related letters. This approach has been used previously only to define scores for aligning individual sequences to sequence profiles, but it has much broader applicability. We describe how to calculate BILD scores efficiently, and illustrate their uses in Gibbs sampling optimization procedures, gapped alignment, and the construction of hidden Markov model profiles. BILD scores enable automated selection of optimal motif and domain model widths, and can inform the decision of whether to include a sequence in a multiple alignment, and the selection of insertion and deletion locations. Other applications include the classification of related sequences into subfamilies, and the definition of profile-profile alignment scores. Although a fully realized multiple alignment program must rely upon more than substitution scores, many existing multiple alignment programs can be modified to employ BILD scores. We illustrate how simple BILD score based strategies can enhance the recognition of DNA binding domains, including the Api-AP2 domain in Toxoplasma gondii and Plasmodium falciparum.

Effect of linker length and dockerin position on conversion of a Thermobifida fusca endoglucanase to the cellulosomal mode

We have been developing the cellulases of Thermobifida fusca as a model to explore the conversion from a free cellulase system to the cellulosomal mode. Three of the six T. fusca cellulases (endoglucanase Cel6A and exoglucanases Cel6B and Cel48A) have been converted in previous work by replacing their cellulose-binding modules (CBMs) with a dockerin, and the resultant recombinant "cellulosomized" enzymes were incorporated into chimeric scaffolding proteins that contained cohesin(s) together with a CBM. The activities of the resultant designer cellulosomes were compared with an equivalent mixture of wild-type enzymes. In the present work, a fourth T. fusca cellulase, Cel5A, was equipped with a dockerin and intervening linker segments of different lengths to assess their contribution to the overall activity of simple one- and two-enzyme designer cellulosome complexes. The results demonstrated that cellulose binding played a major role in the degradation of crystalline cellulosic substrates. The combination of the converted Cel5A endoglucanase with the converted Cel48A exoglucanase also exhibited a measurable proximity effect for the most recalcitrant cellulosic substrate (Avicel). The length of the linker between the catalytic module and the dockerin had little, if any, effect on the activity. However, positioning of the dockerin on the opposite (C-terminal) side of the enzyme, consistent with the usual position of dockerins on most cellulosomal enzymes, resulted in an enhanced synergistic response. These results promote the development of more complex multienzyme designer cellulosomes, which may eventually be applied for improved degradation of plant cell wall biomass.

Structural basis for exquisite specificity of affinity clamps, synthetic binding proteins generated through directed domain-interface evolution

We have established a new protein-engineering strategy termed "directed domain-interface evolution" that generates a binding site by linking two protein domains and then optimizing the interface between them. Using this strategy, we have generated synthetic two-domain "affinity clamps" using PDZ and fibronectin type III (FN3) domains as the building blocks. While these affinity clamps all had significantly higher affinity toward a target peptide than the underlying PDZ domain, two distinct types of affinity clamps were found in terms of target specificity. One type conserved the specificity of the parent PDZ domain, and the other increased the specificity dramatically. Here, we characterized their specificity profiles using peptide phage-display libraries and scanning mutagenesis, which suggested a significantly enlarged recognition site of the high-specificity affinity clamps. The crystal structure of a high-specificity affinity clamp showed extensive contacts with a portion of the peptide ligand that is not recognized by the parent PDZ domain, thus rationalizing the improvement of the specificity of the affinity clamp. A comparison with another affinity clamp structure showed that, although both had extensive contacts between PDZ and FN3 domains, they exhibited a large offset in the relative position of the two domains. Our results indicate that linked domains could rapidly fuse and evolve as a single functional module, and that the inherent plasticity of domain interfaces allows for the generation of diverse active-site topography. These attributes of directed domain-interface evolution provide facile means to generate synthetic proteins with a broad range of functions.

Consequences of domain insertion on the stability and folding mechanism of a protein

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the "insert-in-flap" (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDDeltaIF, SlyD* without the chaperone domain (residues 1-69 and 130-165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.

The structural impact of a polyglutamine tract is location-dependent

Polyglutamine (polyQ) expansion leads to protein aggregation and neurodegeneration in Huntington's disease and eight other inherited neurological conditions. Expansion of the polyQ tract beyond a threshold of 37 glutamines leads to the formation of toxic nuclear aggregates. This suggests that polyQ expansion causes a conformational change within the protein, the nature of which is unclear. There is a trend in the disease proteins that the polyQ tract is located external to but not within a structured domain. We have created a model polyQ protein in which the repeat location mimics the flexible environment of the polyQ tract in the disease proteins. Our model protein recapitulates the aggregation features observed with the clinical proteins and allows structural characterization. With the use of NMR spectroscopy and a range of biophysical techniques, we demonstrate that polyQ expansion into the pathological range has no effect on the structure, dynamics, and stability of a domain adjacent to the polyQ tract. To explore the clinical significance of repeat location, we engineered a variant of the model protein with a polyQ tract within the domain, a location that does not mimic physiological context, demonstrating significant destabilization and structural perturbation. These different effects highlight the importance of repeat location. We conclude that protein misfolding within the polyQ tract itself is the driving force behind the key characteristics of polyQ disease, and that structural perturbation of flanking domains is not required.

Biomimetic nanostructures: creating a high-affinity zinc-binding site in a folded nonbiological polymer

One of the long-term goals in developing advanced biomaterials is to generate protein-like nanostructures and functions from a completely nonnatural polymer. Toward that end, we introduced a high-affinity zinc-binding function into a peptoid (N-substituted glycine polymer) two-helix bundle. Borrowing from well-understood zinc-binding motifs in proteins, thiol and imidazole moieties were positioned within the peptoid such that both helices must align in close proximity to form a binding site. We used fluorescence resonance energy transfer (FRET) reporter groups to measure the change of the distance between the two helical segments and to probe the binding of zinc. We systematically varied the position and number of zinc-binding residues, as well as the sequence and size of the loop that connects the two helical segments. We found that certain peptoid two-helix bundles bind zinc with nanomolar affinities and high selectivity compared to other divalent metal ions. Our work is a significant step toward generating biomimetic nanostructures with enzyme-like functions.

Polyglutamine gene function and dysfunction in the ageing brain

The coordinated regulation of gene expression and protein interactions determines how mammalian nervous systems develop and retain function and plasticity over extended periods of time such as a human life span. By studying mutations that occur in a group of genes associated with chronic neurodegeneration, the polyglutamine (polyQ) disorders, it has emerged that CAG/glutamine stretches play important roles in transcriptional regulation and protein-protein interactions. However, it is still unclear what the many structural and functional roles of CAG and other low-complexity sequences in eukaryotic genomes are, despite being the most commonly shared peptide fragments in such proteomes. In this review we examine the function of genes responsible for at least 10 polyglutamine disorders in relation to the nervous system and how expansion mutations lead to neuronal dysfunction, by particularly focusing on Huntington's disease (HD). We argue that the molecular and cellular pathways that turn out to be dysfunctional during such diseases, as a consequence of a CAG expansion, are also involved in the ageing of the central nervous system. These are pathways that control protein degradation systems (including molecular chaperones), axonal transport, redox-homeostasis and bioenergetics. CAG expansion mutations confer novel properties on proteins that lead to a slow-progressing neuronal pathology and cell death similar to that found in other age-related conditions such as Alzheimer's and Parkinson's diseases.

Protein folding transition states probed by loop extension

We propose a new way to characterize protein folding transition states by (1) insertion of one or more residues into an unstructured protein loop, (2) measurement of the effect on protein folding kinetics and thermodynamics, and (3) analysis of the results in terms of a rate-equilibrium free energy relationship, alpha(Loop). alpha(Loop) reports on the fraction of molecules that form the perturbed loop in the transition state. Interpretation of the changes in equilibrium free energy using standard polymer theory can help detect residual structure in the unfolded state. We illustrate our approach with data for the model proteins CI2 and the alpha spectrin SH3 domain.

Loop-closure principles in protein folding

Simple theoretical concepts and models have been helpful to understand the folding rates and routes of single-domain proteins. As reviewed in this article, a physical principle that appears to underly these models is loop closure.

Modeling mutations in protein structures

We describe an automated method for the modeling of point mutations in protein structures. The protein is represented by all non-hydrogen atoms. The scoring function consists of several types of physical potential energy terms and homology-derived restraints. The optimization method implements a combination of conjugate gradient minimization and molecular dynamics with simulated annealing. The testing set consists of 717 pairs of known protein structures differing by a single mutation. Twelve variations of the scoring function were tested in three different environments of the mutated residue. The best-performing protocol optimizes all the atoms of the mutated residue, with respect to a scoring function that includes molecular mechanics energy terms for bond distances, angles, dihedral angles, peptide bond planarity, and non-bonded atomic contacts represented by Lennard-Jones potential, dihedral angle restraints derived from the aligned homologous structure, and a statistical potential for non-bonded atomic interactions extracted from a large set of known protein structures. The current method compares favorably with other tested approaches, especially when predicting long and flexible side-chains. In addition to the thoroughness of the conformational search, sampled degrees of freedom, and the scoring function type, the accuracy of the method was also evaluated as a function of the flexibility of the mutated side-chain, the relative volume change of the mutated residue, and its residue type. The results suggest that further improvement is likely to be achieved by concentrating on the improvement of the scoring function, in addition to or instead of increasing the variety of sampled conformations.

The length dependence of the polyQ-mediated protein aggregation

Polyglutamine (polyQ) repeat disorders are caused by the expansion of CAG tracts in certain genes, resulting in transcription of proteins with abnormally long polyQ inserts. When these inserts expand beyond 35-45 glutamines, affected proteins form toxic aggregates, leading to neuron death. Chymotrypsin inhibitor 2 (CI2) with an inserted glutamine repeat has previously been used to model polyQ-mediated aggregation in vitro. However, polyQ insertion lengths in these studies have been kept below the pathogenic threshold. We perform molecular dynamics simulations to study monomer folding dynamics and dimer formation in CI2-polyQ chimeras with insertion lengths of up to 80 glutamines. Our model recapitulates the experimental results of previous studies of chimeric CI2 proteins, showing high folding cooperativity of monomers as well as protein association via domain swapping. Surprisingly, for chimeras with insertion lengths above the pathogenic threshold, monomer folding cooperativity decreases and the dominant mode for dimer formation becomes interglutamine hydrogen bonding. These results support a mechanism for pathogenic polyQ-mediated aggregation, in which expanded polyQ tracts destabilize affected proteins and promote the formation of partially unfolded intermediates. These unfolded intermediates form aggregates through associations by interglutamine interactions.

Thermodynamic analysis of an antagonistic folding-unfolding equilibrium between two protein domains

A simple model is formulated for analyzing the coupled folding-unfolding equilibrium present in a unique class of molecular switch proteins. We previously fused two single-domain proteins, barnase and ubiquitin, such that the free energy stored in the folded structure of one subunit is used to drive unfolding of the other. Here, we present a thermodynamic test of that mechanism. The antagonistic interaction is represented by a coupling free energy term DeltaGX. DeltaGX is the penalty imposed on folding of one domain by the native structure of the other. If DeltaGX=0, then neither domain senses the other and they fold and unfold independently. If DeltaGX>0, then destabilizing one domain will stabilize the other, and vice versa. In the limit where DeltaGX is greater than the intrinsic stability of either protein, then only one domain can be folded at any given time. We estimate DeltaGX by measuring stability parameters for a series of mutants that destabilize either the barnase or ubiquitin domains. Fitting the data to the model leads to a DeltaGX value of approximately 4 kcal mol(-1). DeltaGX is proposed to depend on both the length of the linker peptides used to join the two proteins, and on the inherent structural plasticity of each domain. We predict that shortening the linkers from their current lengths of two and three amino acid residues will increase structural and thermodynamic coupling.

Stabilization of a binary protein complex by intein-mediated cyclization

The study of protein-protein interactions can be hampered by the instability of one or more of the protein complex components. In this study, we showed that intein-mediated cyclization can be used to engineer an artificial intramolecular cyclic protein complex between two interacting proteins: the largely unstable LIM-only protein 4 (LMO4) and an unstructured domain of LIM domain binding protein 1 (ldb1). The X-ray structure of the cyclic complex is identical to noncyclized versions of the complex. Chemical and thermal denaturation assays using intrinsic tryptophan fluorescence and dynamic light scattering were used to compare the relative stabilities of the cyclized complex, the intermolecular (or free) complex, and two linear versions of the intramolecular complex (in which the interacting domains of LMO4 and ldb1 were fused, via a flexible linker, in either orientation). In terms of resistance to denaturation, the cyclic complex is the most stable variant and the intermolecular complex is the least stable; however, the two linear intramolecular variants show significant differences in stability. These differences appear to be related to the relative contact order (the average distance in sequence between residues that make contacts within a structure) of key binding residues at the interface of the two proteins. Thus, the restriction of the more stable component of a complex may enhance stability to a greater extent than restraining less stable components.

Active TEM-1 beta-lactamase mutants with random peptides inserted in three contiguous surface loops

Engineering of alternative binding sites on the surface of an enzyme while preserving the enzymatic activity would offer new opportunities for controlling the activity by binding of non-natural ligands. Loops and turns are the natural substructures in which binding sites might be engineered with this purpose. We have genetically inserted random peptide sequences into three relatively rigid and contiguous loops of the TEM-1 beta-lactamase and assessed the tolerance to insertion by the percentage of active mutants. Our results indicate that tolerance to insertion could not be correlated to tolerance to mutagenesis. A turn between two beta-strands bordering the active site was observed to be tolerant to random mutagenesis but not to insertions. Two rigid loops comprising rather well-conserved amino acid residues tolerated insertions, although with some constraints. Insertions between the N-terminal helix and the first beta-strand generated active libraries if cysteine residues were included at both ends of the insert, suggesting the requirement for a stabilizing disulfide bridge. Random sequences were relatively well accommodated within the loop connecting the final beta-strand to the C-terminal helix, particularly if the wild-type residue was retained at one of the loops' end. This suggests two strategies for increasing the percentage of active mutants in insertion libraries. The amino acid distribution in the engineered loops was analyzed and found to be less biased against hydrophobic residues than in natural medium-sized loops. The combination of these activity-selected libraries generated a huge library containing active hybrid enzymes with all three loops modified.

Sequential reorganization of beta-sheet topology by insertion of a single strand

Insertions, duplications, and deletions of sequence segments are thought to be major evolutionary mechanisms that increase the structural and functional diversity of proteins. Alternative splicing, for example, is an intracellular editing mechanism that is thought to generate isoforms for 30%-50% of all human genes. Whereas the inserted sequences usually display only minor structural rearrangements at the insertion site, recent observations indicate that they may also cause more dramatic structural displacements of adjacent structures. In the present study we test how artificially inserted sequences change the structure of the beta-sheet region in T4 lysozyme. Copies of two different beta-strands were inserted into two different loops of the beta-sheet, and the structures were determined. Not surprisingly, one insert "loops out" at its insertion site and forms a new small beta-hairpin structure. Unexpectedly, however, the second insertion leads to displacement of adjacent strands and a sequential reorganization of the beta-sheet topology. Even though the insertions were performed at two different sites, looping out occurred at the C-terminal end of the same beta-strand. Reasons as to why a non-native sequence would be recruited to replace that which occurs in the native protein are discussed. Our results illustrate how sequence insertions can facilitate protein evolution through both local and nonlocal changes in structure.

Experimental characterization of the folding kinetics of the notch ankyrin domain

Proteins constructed from linear arrays of tandem repeats provide a simplified architecture for understanding protein folding. Here, we examine the folding kinetics of the ankyrin repeat domain from the Drosophila Notch receptor, which consists of six folded ankyrin modules and a seventh partly disordered N-terminal ankyrin repeat sequence. Both the refolding and unfolding kinetics are best described as a sum of two exponential phases. The slow, minor refolding phase is limited by prolyl isomerization in the denatured state (D). The minor unfolding phase, which appears as a lag during fluorescence-detected unfolding, is consistent with an on-pathway intermediate (I). This intermediate, although not directly detected during refolding, is shown to be populated by interrupted refolding experiments. When plotted against urea, the rate constants for the major unfolding and refolding phases define a single non-linear v-shaped chevron, as does the minor unfolding phase. These two chevrons, along with unfolding amplitudes, are well-fitted by a sequential three-state model, which yields rate constants for the individual steps in folding and unfolding. Based on these fitted parameters, the D to I step is rate-limiting, and closely matches the major observed refolding phase at low denaturant concentrations. I appears to be midway between N and D in folding free energy and denaturant sensitivity, but has Trp fluorescence properties close to N. Although the Notch ankyrin domain has a simple architecture, folding is slow, with the limiting refolding rate constant as much as seven orders of magnitude smaller than expected from topological predictions.

Loop Entropy and Cytochrome c Stability

The loop entropy model proposes that loop closure in a protein becomes entropically more costly as the length of the loop increases. A model protein, cytochrome c, is composed of four loops connecting five helices surrounding a heme-containing core. To test the loop entropy model a series of mutant proteins are constructed with (Gly)n or (Thr)n segments (n = 4-20) inserted between Gly23 and Gly24 of omega loop A of a pseudo wild-type reference protein. Scanning calorimetry shows that protein stability decreases as n increases in the (Gly)n or (Thr)n segment. The dependence of stability on loop length is analyzed with the loop entropy model. Fitting to the model gives a quantitative description of stability differences for the mutant proteins, but with a smaller power-dependence of the probability of loop closure (c-value) than expected from polymer theory. A possible explanation for the discrepancy is that thermodynamically unfavorable loop entropy is partially offset by interactions between the inserted homopolymer and flanking heteropolymer portions of the unfolded protein. The interactions may involve molecular crowding that favors coalescence of the heteropolymer at the insert site and thus closure of the homopolymer loops, possibly as an aspect of the folding code. This may allow use of loop insert mutants to assess the strength of the heteropolymer-encoded folding signals that facilitate loop closure at the insert site.

Site-specific Reflex Response of Ubiquitin to Loop Insertions

Predicting the structural effects of insertions in proteins by homology modeling remains a challenge. To investigate the molecular basis for conformational adaptations to insertions, ten mutants of ubiquitin were generated by introducing five different inserts, varying from five to 11 residues in size, at two different sites. Most insertion sequences were derived from homologous positions in structurally homologous ubiquitin-like proteins; to test sequence specificity, insertions were made into both homologous and non-homologous sites in ubiquitin. Structural inferences from NMR data suggest that each insertion site shows a reflex response to insertions: the sequence of the insertion has much less impact on structural adaptations than does the site of the insertion. Further, each site responds to insertions in a unique but consistent manner. For a given insertion site, different inserted sequences give rise to different stabilities, but the relationship between stability and sequence is not yet clear. However, the change in stability is similar for all insertions in a given site.

Effects of the PT region of EngD and HLD of CbpA on solubility, catalytic activity and purification characteristics of EngD-CBDCbpA fusions from Clostridium cellulovorans

Chimeric proteins combining the catalytic N-terminal region of native EngD with its proline-threonine-threonine (PT) linker region, hydrophilic domain (HLD) and cellulose binding domain (CBD) of cellulose binding protein A (CbpA) from Clostridium cellulovorans were constructed, expressed, and analyzed. The chimeric proteins with CBD(CbpA) all demonstrated strong affinity to Avicel. The chimeric protein with the PT region of EngD and the HLD had the best catalytic activity and the highest estimated percentage of soluble protein amongst the chimeric proteins. Native EngD and two of the chimeric proteins (EngD-PT-HLD-CBD and EngD-CBD) were purified and their characteristics analyzed. Their binding affinities to Avicel as well as their enzymatic activities against various substrates were found to be consistent with the results we saw from protein lysate samples, which was good binding to Avicel but a decrease in solubility and catalytic activities in chimeric proteins without PT and/or HLD. The reasons for these are discussed. These fusion proteins may be important in applications, such as immobilization to solid cellulose substrate for purification of proteins and enrichment/aggregation of protein complexes.

Molecular origins of rapid and continuous morphological evolution

Mutations in cis-regulatory sequences have been implicated as being the predominant source of variation in morphological evolution. We offer a hypothesis that gene-associated tandem repeat expansions and contractions are a major source of phenotypic variation in evolution. Here, we describe a comparative genomic study of repetitive elements in developmental genes of 92 breeds of dogs. We find evidence for selection for divergence at coding repeat loci in the form of both elevated purity and extensive length polymorphism among different breeds. Variations in the number of repeats in the coding regions of the Alx-4 (aristaless-like 4) and Runx-2 (runt-related transcription factor 2) genes were quantitatively associated with significant differences in limb and skull morphology. We identified similar repeat length variation in the coding repeats of Runx-2, Twist, and Dlx-2 in several other species. The high frequency and incremental effects of repeat length mutations provide molecular explanations for swift, yet topologically conservative morphological evolution.

Neurological proteins are not enriched for repetitive sequences

Proteins associated with disease and development of the nervous system are thought to contain repetitive, simple sequences. However, genome-wide surveys for simple sequences within proteins have revealed that repetitive peptide sequences are the most frequent shared peptide segments among eukaryotic proteins, including those of Saccharomyces cerevisiae, which has few to no specialized developmental and neurological proteins. It is therefore of interest to determine if these specialized proteins have an excess of simple sequences when compared to other sets of compositionally similar proteins. We have determined the relative abundance of simple sequences within neurological proteins and find no excess of repetitive simple sequence within this class. In fact, polyglutamine repeats that are associated with many neurodegenerative diseases are no more abundant within neurological specialized proteins than within nonneurological collections of proteins. We also examined the codon composition of serine homopolymers to determine what forces may play a role in the evolution of extended homopolymers. Codon type homogeneity tends to be favored, suggesting replicative slippage instead of selection as the main force responsible for producing these homopolymers.

Polyglutamine Expansion in Ataxin-3 Does Not Affect Protein Stability

Polyglutamine proteins that cause neurodegenerative disease are known to form proteinaceous aggregates, such as nuclear inclusions, in the neurons of affected patients. Although polyglutamine proteins have been shown to form fibrillar aggregates in a variety of contexts, the mechanisms underlying the aberrant conformational changes and aggregation are still not well understood. In this study, we have investigated the hypothesis that polyglutamine expansion in the protein ataxin-3 destabilizes the native protein, leading to the accumulation of a partially unfolded, aggregation-prone intermediate. To examine the relationship between polyglutamine length and native state stability, we produced and analyzed three ataxin-3 variants containing 15, 28, and 50 residues in their respective glutamine tracts. At pH 7.4 and 37 degrees C, Atax3(Q50), which lies within the pathological range, formed fibrils significantly faster than the other proteins. Somewhat surprisingly, we observed no difference in the acid-induced equilibrium and kinetic un/folding transitions of all three proteins, which indicates that the stability of the native conformation was not affected by polyglutamine tract extension. This has led us to reconsider the mechanisms and factors involved in ataxin-3 misfolding, and we have developed a new model for the aggregation process in which the pathways of un/folding and misfolding are distinct and separate. Furthermore, given that native state stability is unaffected by polyglutamine length, we consider the possible role and influence of other factors in the fibrillization of ataxin-3.