Structure and stability effects of mutations designed to increase the primary sequence symmetry within the core region of a beta-trefoil

FGF1

Fibroblast Growth Factor 1 (Fgf1), also known as acidic FGF (aFGF), is involved in the regulation of various biological processes, ranging from development to disease pathogenesis. It is a single chain polypeptide and is highly expressed in adult brain and kidney tissues. Its expression has been shown to be directed by multiple tissue-specific promoters, which generate transcripts of varying lengths. During development the Fgf1 gene is widely expressed, including in the neural tube, heart and lung. Mouse mutants for this gene are normal under standard laboratory conditions. However, when Fgf1 mutants are exposed to a high fat diet, an aggressive diabetic phenotype has been reported, along with aberrant adipose tissue expansion. Ongoing research on FGF1 and its signalling pathways holds promise for greater understanding of developmental processes as well as the development of novel therapeutic interventions for diseases including diabetes.

Why we are made of proteins and nucleic acids: Structural biology views on extraterrestrial life

Is it a miracle that life exists on the Earth, or is it a common phenomenon in the universe? If extraterrestrial organisms exist, what are they like? To answer these questions, we must understand what kinds of molecules could evolve into life, or in other words, what properties are generally required to perform biological functions and store genetic information. This review summarizes recent findings on simple ancestral proteins, outlines the basic knowledge in textbooks, and discusses the generally required properties for biological molecules from structural biology viewpoints (e.g., restriction of shapes, and types of intra- and intermolecular interactions), leading to the conclusion that proteins and nucleic acids are at least one of the simplest (and perhaps very common) forms of catalytic and genetic biopolymers in the universe. This review article is an extended version of the Japanese article, On the Origin of Life: Coevolution between RNA and Peptide, published in SEIBUTSU BUTSURI Vol. 61, p. 232-235 (2021).

Engineering the kinetic stability of a β-trefoil protein by tuning its topological complexity

Kinetic stability, defined as the rate of protein unfolding, is central to determining the functional lifetime of proteins, both in nature and in wide-ranging medical and biotechnological applications. Further, high kinetic stability is generally correlated with high resistance against chemical and thermal denaturation, as well as proteolytic degradation. Despite its significance, specific mechanisms governing kinetic stability remain largely unknown, and few studies address the rational design of kinetic stability. Here, we describe a method for designing protein kinetic stability that uses protein long-range order, absolute contact order, and simulated free energy barriers of unfolding to quantitatively analyze and predict unfolding kinetics. We analyze two β-trefoil proteins: hisactophilin, a quasi-three-fold symmetric natural protein with moderate stability, and ThreeFoil, a designed three-fold symmetric protein with extremely high kinetic stability. The quantitative analysis identifies marked differences in long-range interactions across the protein hydrophobic cores that partially account for the differences in kinetic stability. Swapping the core interactions of ThreeFoil into hisactophilin increases kinetic stability with close agreement between predicted and experimentally measured unfolding rates. These results demonstrate the predictive power of readily applied measures of protein topology for altering kinetic stability and recommend core engineering as a tractable target for rationally designing kinetic stability that may be widely applicable.

Functionalization of a symmetric protein scaffold: Redundant folding nuclei and alternative oligomeric folding pathways

Successful de novo protein design ideally targets specific folding kinetics, stability thermodynamics, and biochemical functionality, and the simultaneous achievement of all these criteria in a single step design is challenging. Protein design is potentially simplified by separating the problem into two steps: (a) an initial design of a protein "scaffold" having appropriate folding kinetics and stability thermodynamics, followed by (b) appropriate functional mutation-possibly involving insertion of a peptide functional "cassette." This stepwise approach can also separate the orthogonal effects of the "stability/function" and "foldability/function" tradeoffs commonly observed in protein design. If the scaffold is a protein architecture having an exact rotational symmetry, then there is the potential for redundant folding nuclei and multiple equivalent sites of functionalization; thereby enabling broader functional adaptation. We describe such a "scaffold" and functional "cassette" design strategy applied to a β-trefoil threefold symmetric architecture and a heparin ligand functionality. The results support the availability of redundant folding nuclei within this symmetric architecture, and also identify a minimal peptide cassette conferring heparin affinity. The results also identify an energy barrier of destabilization that switches the protein folding pathway from monomeric to trimeric, thereby identifying another potential advantage of symmetric protein architecture in de novo design.

Cooperative hydrophobic core interactions in the β-trefoil architecture

Symmetric protein architectures have a compelling aesthetic that suggests a plausible evolutionary process (i.e., gene duplication/fusion) yielding complex architecture from a simpler structural motif. Furthermore, symmetry inspires a practical approach to computational protein design that substantially reduces the combinatorial explosion problem, and may provide practical solutions for structure optimization. Despite such broad relevance, the role of structural symmetry in the key area of hydrophobic core-packing cooperativity has not been adequately studied. In the present report, the threefold rotational symmetry intrinsic to the β-trefoil architecture is shown to form a geometric basis for highly-cooperative core-packing interactions that both stabilize the local repeating motif and promote oligomerization/long-range contacts in the folding process. Symmetry in the β-trefoil structure also permits tolerance towards mutational drift that involves a structural quasi-equivalence at several key core positions.

Regenerative responses of rabbit corneal endothelial cells to stimulation by fibroblast growth factor 1 (FGF1) derivatives, TTHX1001 and TTHX1114

Utilising rabbit corneal endothelial cells (CEC) in three different paradigms, two human FGF1 derivatives (TTHX1001 and TTHX1114), engineered to exhibit greater stability, were tested as proliferative agents. Primary CECs and mouse NIH 3T3 cells treated with the two FGF1 derivatives showed equivalent EC₅₀ ranges (3.3-24 vs.1.9-16. ng/mL) and, in organ culture, chemically lesioned corneas regained half of the lost endothelial layer in three days after treatment with the FGF1 derivatives as compared to controls. In vivo, following cryolesioning, the CEC monolayer, as judged by specular microscopy, regenerated 10-11 days faster when treated with TTHX1001. Over two weeks, all treated eyes showed clearing of opacity about twice that of untreated controls. In all three rabbit models, both FGF1 derivatives were effective in inducing CEC proliferation over control conditions, supporting the prediction that these stabilised FGF1 derivatives can potentially regenerate corneal endothelial deficits in humans.

Conserved buried water molecules enable the β-trefoil architecture

Available high-resolution crystal structures for the family of β-trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β-trefoil folds. Three buried waters conserved among all β-trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil-fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B-factor) substantially lower than the average Cα atom, providing essentially ideal H-bonding geometry with three solvent-inaccessible main chain groups, simultaneously serving as a bridging H-bond for three different β-strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core-packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β-trefoil fold is therefore, surprisingly, a buried water molecule.

Oligomerization of a symmetric β-trefoil protein in response to folding nucleus perturbation

Gene duplication and fusion events in protein evolution are postulated to be responsible for the common protein folds exhibiting internal rotational symmetry. Such evolutionary processes can also potentially yield regions of repetitive primary structure. Repetitive primary structure offers the potential for alternative definitions of critical regions, such as the folding nucleus (FN). In principle, more than one instance of the FN potentially enables an alternative folding pathway in the face of a subsequent deleterious mutation. We describe the targeted mutation of the carboxyl-terminal region of the (internally located) FN of the de novo designed purely-symmetric β-trefoil protein Symfoil-4P. This mutation involves wholesale replacement of a repeating trefoil-fold motif with a "blade" motif from a β-propeller protein, and postulated to trap that region of the Symfoil-4P FN in a nonproductive folding intermediate. The resulting protein (termed "Bladefoil") is shown to be cooperatively folding, but as a trimeric oligomer. The results illustrate how symmetric protein architectures have potentially diverse folding alternatives available to them, including oligomerization, when preferred pathways are perturbed.

Ab initio folding of a trefoil-fold motif reveals structural similarity with a β-propeller blade motif

Many protein architectures exhibit evidence of internal rotational symmetry postulated to be the result of gene duplication/fusion events involving a primordial polypeptide motif. A common feature of such structures is a domain-swapped arrangement at the interface of the N- and C-termini motifs and postulated to provide cooperative interactions that promote folding and stability. De novo designed symmetric protein architectures have demonstrated an ability to accommodate circular permutation of the N- and C-termini in the overall architecture; however, the folding requirement of the primordial motif is poorly understood, and tolerance to circular permutation is essentially unknown. The β-trefoil protein fold is a threefold-symmetric architecture where the repeating ~42-mer "trefoil-fold" motif assembles via a domain-swapped arrangement. The trefoil-fold structure in isolation exposes considerable hydrophobic area that is otherwise buried in the intact β-trefoil trimeric assembly. The trefoil-fold sequence is not predicted to adopt the trefoil-fold architecture in ab initio folding studies; rather, the predicted fold is closely related to a compact "blade" motif from the β-propeller architecture. Expression of a trefoil-fold sequence and circular permutants shows that only the wild-type N-terminal motif definition yields an intact β-trefoil trimeric assembly, while permutants yield monomers. The results elucidate the folding requirements of the primordial trefoil-fold motif, and also suggest that this motif may sample a compact conformation that limits hydrophobic residue exposure, contains key trefoil-fold structural features, but is more structurally homologous to a β-propeller blade motif.

Fine-Sampled Photographic Quantitation of Dermal Wound Healing Senescence in Aged BALB/cByJ Mice and Therapeutic Intervention with Fibroblast Growth Factor-1

Objective: To determine quantitative parameters of dermal wound healing senescence in aged BALB/cByJ mice (an important animal model of aging) and to evaluate the potential for therapeutic intervention by fibroblast growth factor-1 (FGF-1). Approach: Utilize a novel noninvasive fine-sampled photographic methodology to quantify wound healing parameters for healing phases from wounding through to wound closure. Results: Parameters associated with key healing phases were quantified and compared between nonaged and aged cohorts of both genders. The results identify a sexual dimorphism in dermal wound healing, with nonaged females exhibiting a greater overall healing efficiency than males. This enhanced healing in females, however, senesces with age such that healing parameters for aged males and females are statistically indistinguishable. Topical application of FGF-1 was identified as an effective therapeutic intervention to treat dermal healing senescence in aged females. Innovation: The FGF intervention is being analyzed using a new recently published model. This approach significantly increases the amount of preclinical animal data obtainable in wound healing studies, minimizes cohort number compared with (lethal) histological studies, and permits a direct statistical comparison between different healing studies. Conclusion: Quantitative parameters of dermal wound healing, obtained from noninvasive fine-sampled photographic data, identify topical FGF-1 as an effective therapeutic to treat the senescence of dermal healing present in aged female BALB/cByJ mice.

Investigating the dynamics and polyanion binding sites of fibroblast growth factor-1 using hydrogen-deuterium exchange mass spectrometry

In this study, we examined the local dynamics of acidic fibroblast growth factor (FGF-1) as well as the binding sites of various polyanions including poly-sulfates (heparin and low MW heparin) and poly-phosphates (phytic acid and ATP) using hydrogen-deuterium exchange mass spectrometry (HX-MS). For local dynamics, results are analyzed at the peptide level as well as in terms of buried amides employing crystallographic B-factors and compared with a residue level heat map generated from HX-MS results. Results show that strand 4 and 5 and the turn between them to be the most flexible regions as was previously seen by NMR. On the other hand, the C-terminal strands 8, 9, and 10 appear to be more rigid which is also consistent with crystallographic B-factors as well as local dynamics studies conducted by NMR. Crystal structures of FGF-1 in complex with heparin have shown that heparin binds to N-terminal Asn18 and to C-terminal Lys105, Tryp107, Lys112, Lys113, Arg119, Pro121, Arg122, Gln127, and Lys128 indicating electrostatic forces as dominant interactions. Heparin binding as determined by HX-MS is consistent with crystallography data. Previous studies have also shown that other polyanions including low MW heparin, phytic acid and ATP dramatically increase the thermal stability of FGF-1. Using HX-MS, we find other poly anions tested bind in a similar manner to heparin, primarily targeting the turns in the lysine rich C-terminal region of FGF-1 along with two distinct N-terminal regions that contains lysines and arginines/histidines. This confirms the interactions between FGF-1 and polyanions are primary directed by electrostatics.

Effect of extension of the heparin binding pocket on the structure, stability, and cell proliferation activity of the human acidic fibroblast growth factor

Acidic human fibroblast growth factor (hFGF1) plays a key role in cell growth and proliferation. Activation of the cell surface FGF receptor is believed to involve the glycosaminoglycan, heparin. However, the exact role of heparin is a subject of considerable debate. In this context, in this study, the correlation between heparin binding affinity and cell proliferation activity of hFGF1 is examined by extending the heparin binding pocket through selective engineering via charge reversal mutations (D82R, D84R and D82R/D84R). Results of biophysical experiments such as intrinsic tryptophan fluorescence and far UV circular dichroism spectroscopy suggest that the gross native structure of hFGF1 is not significantly perturbed by the engineered mutations. However, results of limited trypsin digestion and ANS binding experiments show that the backbone structure of the D82R variant is more flexible than that of the wild type hFGF1. Results of the temperature and urea-induced equilibrium unfolding experiments suggest that the stability of the charge-reversal mutations increases in the presence of heparin. Isothermal titration calorimetry (ITC) data reveal that the heparin binding affinity is significantly increased when the charge on D82 is reversed but not when the negative charge is reversed at both positions D82 and D84 (D82R/D84R). However, despite the increased affinity of D82R for heparin, the cell proliferation activity of the D82R variant is observed to be reduced compared to the wild type hFGF1. The results of this study clearly demonstrate that heparin binding affinity of hFGF1 is not strongly correlated to its cell proliferation activity.

Engineering a Cysteine-Free Form of Human Fibroblast Growth Factor-1 for "Second Generation" Therapeutic Application

Human fibroblast growth factor-1 (FGF-1) has broad therapeutic potential in regenerative medicine but has undesirable biophysical properties of low thermostability and 3 buried cysteine (Cys) residues (at positions 16, 83, and 117) that interact to promote irreversible protein unfolding under oxidizing conditions. Mutational substitution of such Cys residues eliminates reactive buried thiols but cannot be accomplished simultaneously at all 3 positions without also introducing further substantial instability. The mutational introduction of a novel Cys residue (Ala66Cys) that forms a stabilizing disulfide bond (i.e., cystine) with one of the extant Cys residues (Cys83) effectively eliminates one Cys while increasing overall stability. This increase in stability offsets the associated instability of remaining Cys substitution mutations and permits production of a Cys-free form of FGF-1 (Cys16Ser/Ala66Cys/Cys117Ala) with only minor overall instability. The addition of a further stabilizing mutation (Pro134Ala) creates a Cys-free FGF-1 mutant with essentially wild-type biophysical properties. The elimination of buried free thiols in FGF-1 can substantially increase the protein half-life in cell culture. Here, we show that the effective cell survival/mitogenic functional activity of a fully Cys-free form is also substantially increased and is equivalent to wild-type FGF-1 formulated in the presence of heparin sulfate as a stabilizing agent. The results identify this Cys-free FGF-1 mutant as an advantageous "second generation" form of FGF-1 for therapeutic application.

An S116R Phosphorylation Site Mutation in Human Fibroblast Growth Factor-1 Differentially Affects Mitogenic and Glucose-Lowering Activities

Fibroblast growth factor-1 (FGF-1), a potent human mitogen and insulin sensitizer, signals through both tyrosine kinase receptor-mediated autocrine/paracrine pathways as well as a nuclear intracrine pathway. Phosphorylation of FGF-1 at serine 116 (S116) has been proposed to regulate intracrine signaling. Position S116 is located within a ∼17 amino acid C-terminal loop that contains a rich set of functional determinants including heparin∖heparan sulfate affinity, thiol reactivity, nuclear localization, pharmacokinetics, functional half-life, nuclear ligand affinity, stability, and structural dynamics. Mutational targeting of specific functionality in this region without perturbing other functional determinants is a design challenge. S116R is a non-phosphorylatable variant present in bovine FGF-1 and other members of the human FGF family. We show that the S116R mutation in human FGF-1 is accommodated with no perturbation of biophysical or structural properties, and is therefore an attractive mutation with which to elucidate the functional role of phosphorylation. Characterization of S116R shows reduction in NIH 3T3 fibroblast mitogenic stimulation, increase in fibroblast growth factor receptor-1c activation, and prolonged duration of glucose lowering in ob/ob hyperglycemic mice. A novel FGF-1/fibroblast growth factor receptor-1c dimerization interaction combined with non-phosphorylatable intracrine signaling is hypothesized to be responsible for these observed functional effects.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

A single aromatic core mutation converts a designed "primitive" protein from halophile to mesophile folding

The halophile environment has a number of compelling aspects with regard to the origin of structured polypeptides (i.e., proteogenesis) and, instead of a curious niche that living systems adapted into, the halophile environment is emerging as a candidate "cradle" for proteogenesis. In this viewpoint, a subsequent halophile-to-mesophile transition was a key step in early evolution. Several lines of evidence indicate that aromatic amino acids were a late addition to the codon table and not part of the original "prebiotic" set comprising the earliest polypeptides. We test the hypothesis that the availability of aromatic amino acids could facilitate a halophile-to-mesophile transition by hydrophobic core-packing enhancement. The effects of aromatic amino acid substitutions were evaluated in the core of a "primitive" designed protein enriched for the 10 prebiotic amino acids (A,D,E,G,I,L,P,S,T,V)-having an exclusively prebiotic core and requiring halophilic conditions for folding. The results indicate that a single aromatic amino acid substitution is capable of eliminating the requirement of halophile conditions for folding of a "primitive" polypeptide. Thus, the availability of aromatic amino acids could have facilitated a critical halophile-to-mesophile protein folding adaptation-identifying a selective advantage for the incorporation of aromatic amino acids into the codon table.

A novel β-trefoil lectin from the parasol mushroom (Macrolepiota procera) is nematotoxic

Lectins are carbohydrate-binding proteins present in all organisms. Some cytoplasmic lectins from fruiting bodies of dikaryotic fungi are toxic against a variety of parasites and predators. We have isolated, cloned and expressed a novel, single domain lectin from Macrolepiota procera, designated MpL. Determination of the crystal structure revealed that MpL is a ricin B-like lectin with a β-trefoil fold. Biochemical characterization, site-directed mutagenesis, co-crystallization with carbohydrates, isothermal titration calorimetry and glycan microarray analyses show that MpL forms dimers with the carbohydrate-binding site at the α-repeat, with the highest specificity for terminal N-acetyllactosamine and other β-galactosides. A second putative carbohydrate-binding site with a low affinity for galactose is present at the γ-repeat. In addition, a novel hydrophobic binding site was detected in MpL with specificity for molecules other than carbohydrates. The tissue specific distribution of MpL in the stipe and cap tissue of fruiting bodies and its toxicity towards the nematode Caenorhabditis elegans indicate a function of MpL in protecting fruiting bodies against predators and parasites.

DATABASE

Nucleotide sequence data have been deposited in the DDBJ/EMBL/GenBank databases under accession numbers HQ449738 and HQ449739. Structural data have been deposited in the Protein Data Bank under accession codes 4ION, 4IYB, 4IZX and 4J2S.

Symmetric protein architecture in protein design: top-down symmetric deconstruction

Top-down symmetric deconstruction (TDSD) is a joint experimental and computational approach to generate a highly stable, functionally benign protein scaffold for intended application in subsequent functional design studies. By focusing on symmetric protein folds, TDSD can leverage the dramatic reduction in sequence space achieved by applying a primary structure symmetric constraint to the design process. Fundamentally, TDSD is an iterative symmetrization process, in which the goal is to maintain or improve properties of thermodynamic stability and folding cooperativity inherent to a starting sequence (the "proxy"). As such, TDSD does not attempt to solve the inverse protein folding problem directly, which is computationally intractable. The present chapter will take the reader through all of the primary steps of TDSD-selecting a proxy, identifying potential mutations, establishing a stability/folding cooperativity screen-relying heavily on a successful TDSD solution for the common β-trefoil fold.

Creation and structure determination of an artificial protein with three complete sequence repeats

Symfoil-4P is a de novo protein exhibiting the threefold symmetrical β-trefoil fold designed based on the human acidic fibroblast growth factor. First three asparagine-glycine sequences of Symfoil-4P are replaced with glutamine-glycine (Symfoil-QG) or serine-glycine (Symfoil-SG) sequences protecting from deamidation, and His-Symfoil-II was prepared by introducing a protease digestion site into Symfoil-QG so that Symfoil-II has three complete repeats after removal of the N-terminal histidine tag. The Symfoil-QG and SG and His-Symfoil-II proteins were expressed in Eschericha coli as soluble protein, and purified by nickel affinity chromatography. Symfoil-II was further purified by anion-exchange chromatography after removing the HisTag by proteolysis. Both Symfoil-QG and Symfoil-II were crystallized in 0.1 M Tris-HCl buffer (pH 7.0) containing 1.8 M ammonium sulfate as precipitant at 293 K; several crystal forms were observed for Symfoil-QG and II. The maximum diffraction of Symfoil-QG and II crystals were 1.5 and 1.1 Å resolution, respectively. The Symfoil-II without histidine tag diffracted better than Symfoil-QG with N-terminal histidine tag. Although the crystal packing of Symfoil-II is slightly different from Symfoil-QG and other crystals of Symfoil derivatives having the N-terminal histidine tag, the refined crystal structure of Symfoil-II showed pseudo-threefold symmetry as expected from other Symfoils. Since the removal of the unstructured N-terminal histidine tag did not affect the threefold structure of Symfoil, the improvement of diffraction quality of Symfoil-II may be caused by molecular characteristics of Symfoil-II such as molecular stability.

Protein engineering for cardiovascular therapeutics: untapped potential for cardiac repair

A number of new and innovative approaches for repairing damaged myocardium are currently undergoing investigation, with several encouraging results. In addition to the progression of stem cell-based approaches and gene therapy/silencing methods, evidence continues to emerge that protein therapeutics may be used to directly promote cardiac repair and even regeneration. However, proteins are often limited in their therapeutic potential by short local half-lives and insufficient bioavailability and bioactivity, and many academic laboratories studying cardiovascular diseases are more comfortable with molecular and cellular biology than with protein biochemistry. Protein engineering has been used broadly to overcome weaknesses traditionally associated with protein therapeutics and has the potential to specifically enhance the efficacy of molecules for cardiac repair. However, protein engineering as a strategy has not yet been used in the development of cardiovascular therapeutics to the degree that it has been used in other fields. In this review, we discuss the role of engineered proteins in cardiovascular therapies to date. Further, we address the promise of applying emerging protein engineering technologies to cardiovascular medicine and the barriers that must be overcome to enable the ultimate success of this approach.

Emergence of symmetric protein architecture from a simple peptide motif: evolutionary models

Structural symmetry is observed in the majority of fundamental protein folds and gene duplication and fusion evolutionary processes are postulated to be responsible. However, convergent evolution leading to structural symmetry has also been proposed; additionally, there is debate regarding the extent to which exact primary structure symmetry is compatible with efficient protein folding. Issues of symmetry in protein evolution directly impact strategies for de novo protein design as symmetry can substantially simplify the design process. Additionally, when considering gene duplication and fusion in protein evolution, there are two competing models: "emergent architecture" and "conserved architecture". Recent experimental work has shed light on both the evolutionary process leading to symmetric protein folds as well as the ability of symmetric primary structure to efficiently fold. Such studies largely support a "conserved architecture" evolutionary model, suggesting that complex protein architecture was an early evolutionary achievement involving oligomerization of smaller polypeptides.

Experimental support for the foldability-function tradeoff hypothesis: segregation of the folding nucleus and functional regions in fibroblast growth factor-1

The acquisition of function is often associated with destabilizing mutations, giving rise to the stability-function tradeoff hypothesis. To test whether function is also accommodated at the expense of foldability, fibroblast growth factor-1 (FGF-1) was subjected to a comprehensive φ-value analysis at each of the 11 turn regions. FGF-1, a β-trefoil fold, represents an excellent model system with which to evaluate the influence of function on foldability: because of its threefold symmetric structure, analysis of FGF-1 allows for direct comparisons between symmetry-related regions of the protein that are associated with function to those that are not; thus, a structural basis for regions of foldability can potentially be identified. The resulting φ-value distribution of FGF-1 is highly polarized, with the majority of positions described as either folded-like or denatured-like in the folding transition state. Regions important for folding are shown to be asymmetrically distributed within the protein architecture; furthermore, regions associated with function (i.e., heparin-binding affinity and receptor-binding affinity) are localized to regions of the protein that fold after barrier crossing (late in the folding pathway). These results provide experimental support for the foldability-function tradeoff hypothesis in the evolution of FGF-1. Notably, the results identify the potential for folding redundancy in symmetric protein architecture with important implications for protein evolution and design.

Pharmacokinetic properties of 2nd-generation fibroblast growth factor-1 mutants for therapeutic application

Fibroblast growth factor-1 (FGF-1) is an angiogenic factor with therapeutic potential for the treatment of ischemic disease. FGF-1 has low intrinsic thermostability and is characteristically formulated with heparin as a stabilizing agent. Heparin, however, adds a number of undesirable properties that negatively impact safety and cost. Mutations that increase the thermostability of FGF-1 may obviate the need for heparin in formulation and may prove to be useful “2nd-generation” forms for therapeutic use. We report a pharmacokinetic (PK) study in rabbits of human FGF-1 in the presence and absence of heparin, as well as three mutant forms having differential effects upon thermostability, buried reactive thiols, and heparin affinity. The results support the hypothesis that heparan sulfate proteoglycan (HSPG) in the vasculature of liver, kidney and spleen serves as the principle peripheral compartment in the distribution kinetics. The addition of heparin to FGF-1 is shown to increase endocrine-like properties of distribution. Mutant forms of FGF-1 that enhance thermostability or eliminate buried reactive thiols demonstrate a shorter distribution half-life, a longer elimination half-life, and a longer mean residence time (MRT) in comparison to wild-type FGF-1. The results show how such mutations can produce useful 2nd-generation forms with tailored PK profiles for specific therapeutic application.

Designing proteins from simple motifs: opportunities in Top-Down Symmetric Deconstruction

The purpose of this review is to describe the development of 'top-down' approaches to protein design. It will be argued that a diverse number of studies over the past decade, involving many investigators, and focused upon elucidating the role of symmetry in protein evolution and design, are converging into a novel top-down approach to protein design. Top-down design methodologies have successfully produced comparatively simple polypeptide 'building blocks' (typically comprising 40-60 amino acids) useful in generating complex protein architecture, and have produced compelling data in support of macro-evolutionary pathways of protein structure. Furthermore, a distillation of the experimental approaches utilized in such studies suggests the potential for method formalism, one that may accelerate future success in this field.

An empirical phase diagram approach to investigate conformational stability of "second-generation" functional mutants of acidic fibroblast growth factor-1

Acidic fibroblast growth factor-1 (FGF-1) is an angiogenic protein which requires binding to a polyanion such as heparin for its mitogenic activity and physicochemical stability. To evaluate the extent to which this heparin dependence on solution stability could be reduced or eliminated, the structural integrity and conformational stability of 10 selected FGF-1 mutants were examined as a function of solution pH and temperature by a series of spectroscopic methods including circular dichroism, intrinsic and extrinsic fluorescence spectroscopy and static light scattering. The biophysical data were summarized in the form of colored empirical phase diagrams (EPDs). FGF-1 mutants were identified with stability profiles in the absence of heparin comparable to that of wild-type FGF-1 in the presence of heparin while still retaining their biological activity. In addition, a revised version of the EPD methodology was found to provide an information rich, high throughput approach to compare the effects of mutations on the overall conformational stability of proteins in terms of their response to environmental stresses such as pH and temperature.

PROTS: a fragment based protein thermo-stability potential

Designing proteins with enhanced thermo-stability has been a main focus of protein engineering because of its theoretical and practical significance. Despite extensive studies in the past years, a general strategy for stabilizing proteins still remains elusive. Thus effective and robust computational algorithms for designing thermo-stable proteins are in critical demand. Here we report PROTS, a sequential and structural four-residue fragment based protein thermo-stability potential. PROTS is derived from a nonredundant representative collection of thousands of thermophilic and mesophilic protein structures and a large set of point mutations with experimentally determined changes of melting temperatures. To the best of our knowledge, PROTS is the first protein stability predictor based on integrated analysis and mining of these two types of data. Besides conventional cross validation and blind testing, we introduce hypothetical reverse mutations as a means of testing the robustness of protein thermo-stability predictors. In all tests, PROTS demonstrates the ability to reliably predict mutation induced thermo-stability changes as well as classify thermophilic and mesophilic proteins. In addition, this white-box predictor allows easy interpretation of the factors that influence mutation induced protein stability changes at the residue level.

A polypeptide "building block" for the β-trefoil fold identified by "top-down symmetric deconstruction"

Fibroblast growth factor-1, a member of the 3-fold symmetric β-trefoil fold, was subjected to a series of symmetric constraint mutations in a process termed "top-down symmetric deconstruction." The mutations enforced a cumulative exact 3-fold symmetry upon symmetrically equivalent positions within the protein and were combined with a stability screen. This process culminated in a β-trefoil protein with exact 3-fold primary-structure symmetry that exhibited excellent folding and stability properties. Subsequent fragmentation of the repeating primary-structure motif yielded a 42-residue polypeptide capable of spontaneous assembly as a homotrimer, producing a thermostable β-trefoil architecture. The results show that despite pronounced reduction in sequence complexity, pure symmetry in the design of a foldable, thermostable β-trefoil fold is possible. The top-down symmetric deconstruction approach provides a novel alternative means to successfully identify a useful polypeptide "building block" for subsequent "bottom-up" de novo design of target protein architecture.

Experimental support for the evolution of symmetric protein architecture from a simple peptide motif

The majority of protein architectures exhibit elements of structural symmetry, and "gene duplication and fusion" is the evolutionary mechanism generally hypothesized to be responsible for their emergence from simple peptide motifs. Despite the central importance of the gene duplication and fusion hypothesis, experimental support for a plausible evolutionary pathway for a specific protein architecture has yet to be effectively demonstrated. To address this question, a unique "top-down symmetric deconstruction" strategy was utilized to successfully identify a simple peptide motif capable of recapitulating, via gene duplication and fusion processes, a symmetric protein architecture (the threefold symmetric β-trefoil fold). The folding properties of intermediary forms in this deconstruction agree precisely with a previously proposed "conserved architecture" model for symmetric protein evolution. Furthermore, a route through foldable sequence-space between the simple peptide motif and extant protein fold is demonstrated. These results provide compelling experimental support for a plausible evolutionary pathway of symmetric protein architecture via gene duplication and fusion processes.

Symmetric key structural residues in symmetric proteins with beta-trefoil fold

To understand how symmetric structures of many proteins are formed from asymmetric sequences, the proteins with two repeated beta-trefoil domains in Plant Cytotoxin B-chain family and all presently known beta-trefoil proteins are analyzed by structure-based multi-sequence alignments. The results show that all these proteins have similar key structural residues that are distributed symmetrically in their structures. These symmetric key structural residues are further analyzed in terms of inter-residues interaction numbers and B-factors. It is found that they can be distinguished from other residues and have significant propensities for structural framework. This indicates that these key structural residues may conduct the formation of symmetric structures although the sequences are asymmetric.

Structure, evolutionary conservation, and conformational dynamics of Homo sapiens fascin-1, an F-actin crosslinking protein

Eukaryotes have several highly conserved actin-binding proteins that crosslink filamentous actin into compact ordered bundles present in distinct cytoskeletal processes, including microvilli, stereocilia and filopodia. Fascin is an actin-binding protein that is present predominantly in filopodia, which are believed to play a central role in normal and aberrant cell migration. An important outstanding question regards the molecular basis for the unique localization and functional properties of fascin compared with other actin crosslinking proteins. Here, we present the crystal structure of full-length Homo sapiens fascin-1, and examine its packing, conformational flexibility, and evolutionary sequence conservation. The structure reveals a novel arrangement of four tandem beta-trefoil domains that form a bi-lobed structure with approximate pseudo 2-fold symmetry. Each lobe has internal approximate pseudo 2-fold and pseudo 3-fold symmetry axes that are approximately perpendicular, with beta-hairpin triplets located symmetrically on opposite sides of each lobe that mutational data suggest are actin-binding domains. Sequence conservation analysis confirms the importance of hydrophobic core residues that stabilize the beta-trefoil fold, as well as interfacial residues that are likely to stabilize the overall fascin molecule. Sequence conservation also indicates highly conserved surface patches near the putative actin-binding domains of fascin, which conformational dynamics analysis suggests to be coupled via an allosteric mechanism that might have important functional implications for F-actin crosslinking by fascin.

Engineering an improved crystal contact across a solvent-mediated interface of human fibroblast growth factor 1

Large-volume protein crystals are a prerequisite for neutron diffraction studies and their production represents a bottleneck in obtaining neutron structures. Many protein crystals that permit the collection of high-resolution X-ray diffraction data are inappropriate for neutron diffraction owing to a plate-type morphology that limits the crystal volume. Human fibroblast growth factor 1 crystallizes in a plate morphology that yields atomic resolution X-ray diffraction data but has insufficient volume for neutron diffraction. The thin physical dimension has been identified as corresponding to the b cell edge and the X-ray structure identified a solvent-mediated crystal contact adjacent to position Glu81 that was hypothesized to limit efficient crystal growth in this dimension. In this report, a series of mutations at this crystal contact designed to both reduce side-chain entropy and replace the solvent-mediated interface with direct side-chain contacts are reported. The results suggest that improved crystal growth is achieved upon the introduction of direct crystal contacts, while little improvement is observed with side-chain entropy-reducing mutations alone.

X-ray structure and biophysical properties of rabbit fibroblast growth factor 1

The rabbit is an important and de facto animal model in the study of ischemic disease and angiogenic therapy. Additionally, fibroblast growth factor 1 (FGF-1) is emerging as one of the most important growth factors for novel proangiogenic and pro-arteriogenic therapy. However, despite its significance, the fundamental biophysical properties of rabbit FGF-1, including its X-ray structure, have never been reported. Here, the cloning, crystallization, X-ray structure and determination of the biophysical properties of rabbit FGF-1 are described. The X-ray structure shows that the amino-acid differences between human and rabbit FGF-1 are solvent-exposed and therefore potentially immunogenic, while the biophysical studies identify differences in thermostability and receptor-binding affinity that distinguish rabbit FGF-1 from human FGF-1.

Effects of external interactions on protein sequence-structure relations of beta-trefoil fold

Proteins with symmetric structures are ideal models to investigate the sequence-structure relations. We investigate proteins with beta-trefoil fold and find they have different degrees of sequence symmetries although they show similar symmetric structures. To understand this, we calculate the strength of interactions of the beta-trefoil folds with surrounding environments and find the low degrees of sequence symmetries are often correlated with large external interactions. Our results give an additional confirmation of Anfinsen's thermodynamic hypothesis that protein structures are not only determined by their sequences but also by their surrounding environments. We suggest the external interactions should be considered additionally in protein structure prediction through ab initio folding.

Interaction of anti-aggregation agent dimethylethylammonium propane sulfonate with acidic fibroblast growth factor

Prevention of aggregation is critical for analyzing protein structure. Non-detergent sulfobetaines (NDSBs) are known to prevent protein aggregation, but the molecular mechanisms of their anti-aggregation effect are poorly understood. To elucidate the underlying mechanisms, we analyzed the effects of dimethylethylammonium propane sulfonate (NDSB-195) on acidic fibroblast growth factor (aFGF). NDSB-195 (0.5M) increased both aggregation and denaturation temperatures of aFGF by 4 degrees C. Chemical shift perturbation analyses indicated that many affected residues were located at the junction between a beta-strand (or 3(10)-helix) and a loop, irrespective of the chemical properties of the residue. The apparent dissociation constants of the interaction ranged from 0.04 to 3M, indicating weak interactions between NDSB and protein molecules.

Mutagenesis of the crystal contact of acidic fibroblast growth factor

An attempt has been made to improve a crystal contact of human acidic fibroblast growth factor (haFGF; 140 amino acids) to control the crystal growth, because haFGF crystallizes only as a thin-plate form, yielding crystals suitable for X-ray but not neutron diffraction. X-ray crystal analysis of haFGF showed that the Glu81 side chain, located at a crystal contact between haFGF molecules, is in close proximity with an identical residue related by crystallographic symmetry, suggesting that charge repulsion may disrupt suitable crystal-packing interactions. To investigate whether the Glu residue affects the crystal-packing interactions, haFGF mutants in which Glu81 was replaced by Ala, Val, Leu, Ser and Thr were constructed. Although crystals of the Ala and Leu mutants were grown as a thin-plate form by the same precipitant (formate) as the wild type, crystals of the Ser and Thr mutants were grown with increased thickness, yielding a larger overall crystal volume. X-ray structural analysis of the Ser mutant determined at 1.35 A resolution revealed that the hydroxy groups of Ser are linked by hydrogen bonds mediated by the formate used as a precipitant. This approach to engineering crystal contacts may contribute to the development of large protein crystals for neutron crystallography.

Common occurrence of internal repeat symmetry in membrane proteins

Symmetry plays significant roles in protein structure and function. Particularly, symmetric interfaces are known to act as switches for two-state conformational change. Membrane proteins often undergo two-state conformational change during the transport process of ion channels or the active/inactive transitions in receptors. Here, we provide the first comprehensive analyses of internal repeat symmetry in membrane proteins. We examined the known membrane protein structures and found that, remarkably, nearly half of them have internal repeat symmetry. Moreover, we found that the conserved cores of these internal repeats are positioned at the interface of symmetric units when they are mapped on structures. Because of the large sequence divergence that occurs between internal repeats, the inherent symmetry present in protein sequences often has only been detected after structure determination. We therefore developed a sensitive procedure to predict the internal repeat symmetry from sequence information and identified 4653 proteins that are likely to have internal repeat symmetry.

Activation profiles and regulatory cascades of the human kallikrein-related peptidases

The human kallikrein (KLK)-related peptidases are the largest family of serine peptidases, comprising 15 members (KLK1-15) and with the majority (KLK4-15) being identified only within the last decade. Members of this family are associated with important diseased states (including cancer, inflammation, and neurodegeneration) and have been utilized or proposed as clinically important biomarkers or therapeutic targets of interest. All human KLKs are synthesized as prepro-forms that are proteolytically processed to secreted pro-forms via the removal of an amino-terminal secretion signal peptide. The secreted inactive pro-KLKs are then activated extracellularly to mature peptidases by specific proteolytic release of their amino-terminal propeptide. Although a key step in the regulation of KLK function, details regarding the activation of the human pro-KLKs (i.e. the KLK "activome") are unknown, to a significant extent, but have been postulated to involve "activation cascades" with other KLKs and endopeptidases. To characterize more completely the KLK activome, we have expressed from Escherichia coli individual KLK propeptides fused to the amino terminus of a soluble carrier protein. The ability of 12 different mature KLKs to process the 15 different pro-KLK peptide sequences has been determined. Various autolytic and cross-activation relationships identified using this system have subsequently been characterized using recombinant pro-KLK proteins. The results demonstrate the potential for extensive KLK activation cascades and, when combined with available data for the tissue-specific expression of the KLK family, permit the construction of specific regulatory cascades. One such tissue-specific cascade is proposed for the central nervous system.

Spackling the crack: stabilizing human fibroblast growth factor-1 by targeting the N and C terminus beta-strand interactions

The beta-trefoil protein human fibroblast growth factor-1 (FGF-1) is made up of a six-stranded antiparallel beta-barrel closed off on one end by three beta-hairpins, thus exhibiting a 3-fold axis of structural symmetry. The N and C terminus beta-strands hydrogen bond to each other and their interaction is postulated from both NMR and X-ray structure data to be important in folding and stability. Specific mutations within the adjacent N and C terminus beta-strands of FGF-1 are shown to provide a substantial increase in stability. This increase is largely correlated with an increased folding rate constant, and with a smaller but significant decrease in the unfolding rate constant. A series of stabilizing mutations are subsequently combined and result in a doubling of the DeltaG value of unfolding. When taken in the context of previous studies of stabilizing mutations, the results indicate that although FGF-1 is known for generally poor thermal stability, the beta-trefoil architecture appears capable of substantial thermal stability. Targeting stabilizing mutations within the N and C terminus beta-strand interactions of a beta-barrel architecture may be a generally useful approach to increase protein stability. Such stabilized mutations of FGF-1 are shown to exhibit significant increases in effective mitogenic potency, and may prove useful as "second generation" forms of FGF-1 for application in angiogenic therapy.

Structure of rat acidic fibroblast growth factor at 1.4 A resolution

Fibroblast growth factors (FGFs) constitute a family of 22 structurally related heparin-binding polypeptides that are involved in the regulation of cell growth, survival, differentiation and migration. Here, a 1.4 A resolution X-ray structure of rat FGF1 is presented. Two molecules are present in the asymmetric unit of the crystal and they coordinate a total of five sulfate ions. The structures of human, bovine and newt FGF1 have been published previously. Human and rat FGF1 are found to have very similar structures.

Conversion of type I 4:6 to 3:5 beta-turn types in human acidic fibroblast growth factor: effects upon structure, stability, folding, and mitogenic function

Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold, a protein architecture that exhibits a characteristic threefold axis of structural symmetry. FGF-1 contains 11 beta-turns, the majority being type I 3:5; however, a type I 4:6 turn is also found at three symmetry-related locations. The relative uniqueness of the type I 4:6 turn in the FGF-1 structure suggests it may play a key role in the stability, folding, or function of the protein. To test this hypothesis a series of deletion mutations were constructed, the aim of which was to convert existing type I 4:6 turns at two locations into type I 3:5 turns. The results show it is possible to successfully substitute the type I 4:6 turn by a type I 3:5 turn with minimal impact upon protein stability or folding. Thus, these different turn structures, even though they differ in length, exhibit similar energetic properties. Additional sequence swapping mutations within the introduced type I 3:5 turns suggests that the turn sequence primarily affects stability but not turn structure (which appears dictated primarily by the local environment). Although the results suggest that a stable, foldable beta-trefoil protein may be designed utilizing a single turn type (type I 3:5), a type I 4:6 turn at turn 1 of FGF-1 appears essential for efficient mitogenic function.

Highly stable mutants of human fibroblast growth factor-1 exhibit prolonged biological action

Fibroblast growth factor 1 (FGF-1) shows strong angiogenic, osteogenic and tissue-injury repair properties that might be relevant to medical applications. Since FGF-1 is partially unfolded at physiological temperature we decided to increase significantly its conformational stability and test how such an improvement will affect its biological function. Using an homology approach and rational strategy we designed two new single FGF-1 mutations: Q40P and S47I that appeared to be the most strongly stabilizing substitutions among those reported so far, increasing the denaturation temperature by 7.8 deg. C and 9.0 deg. C, respectively. As our goal was to produce highly stable variants of the growth factor, we combined these two mutations with five previously described stabilizing substitutions. The multiple mutants showed denaturation temperatures up to 27 deg. C higher than the wild-type and exhibited full additivity of the mutational effects. All those mutants were biologically competent in several cell culture assays, maintaining typical FGF-1 activities, such as binding to specific cell surface receptors and activation of downstream signaling pathways. Thus, we demonstrate that the low denaturation temperature of wild-type FGF-1 is not related to its fundamental cellular functions, and that FGF-1 action is not affected by its stability. A more detailed analysis of the biological behavior of stable FGF-1 mutants revealed that, compared with the wild-type, their mitogenic properties, as probed by the DNA synthesis assay, were significantly increased in the absence of heparin, and that their half-lives were extensively prolonged. We found that the biological action of the mutants was dictated by their susceptibility to proteases, which strongly correlated with the stability. Mutants which were much more resistant to proteolytic degradation always displayed a significant improvement in the half-life and mitogenesis. Our results show that engineered stable growth factor variants exhibit enhanced and prolonged activity, which can be advantageous in terms of the potential therapeutic applications of FGF-1.

A common sequence-associated physicochemical feature for proteins of beta-trefoil family

Different amino acid sequences can fold into similar tertiary structures but the reasons for it are not very clear. It has been suggested in the literature that these sequences may have some common features associated with them but the exact nature of such shared properties remains largely unknown. We studied a representative sample of proteins from the beta-trefoil family and observed that their amino acid sequences, despite being considerably divergent from each other, can be accounted for by matching to a repetition of three physicochemically similar segments. This observation in turn is consistent with the three-fold pseudo-symmetry in tertiary structures of these proteins.

Redesigning symmetry-related "mini-core" regions of FGF-1 to increase primary structure symmetry: thermodynamic and functional consequences of structural symmetry

Previous reports detailing mutational effects within the hydrophobic core of human acidic fibroblast growth factor (FGF-1) have shown that a symmetric primary structure constraint is compatible with a stably folded protein. In the present report, we investigate symmetrically related pairs of buried hydrophobic residues in FGF-1 (termed "mini-cores") that are not part of the central core. The effect upon the stability and function of FGF-1 mutations designed to increase primary structure symmetry within these "mini-core" regions was evaluated. At symmetry-related positions 22, 64, and 108, the wild-type protein contains either Tyr or Phe side chains. The results show that either residue can be readily accommodated at these positions. At symmetry-related positions 42, 83, and 130, the wild-type protein contains either Cys or Ile side chains. While positions 42 and 130 can readily accommodate either Cys or Ile side chains, position 83 is substantially destabilized by substitution by Ile. Tertiary structure asymmetry in the vicinity of position 83 appears responsible for the inability to accommodate an Ile side chain at this position, and is known to contribute to functional half-life. A mutant form of FGF-1 with enforced primary structure symmetry at positions 22, 64, and 108 (all Tyr) and 42, 83, and 130 (all Cys) is shown to be more stable than the reference FGF-1 protein. The results support the hypothesis that a symmetric primary structure within a symmetric protein superfold represents a solution to achieving a foldable, stable polypeptide, and highlight the role that function may play in the evolution of asymmetry within symmetric superfolds.