Symmetric primary and tertiary structure mutations within a symmetric superfold: a solution, not a constraint, to achieve a foldable polypeptide

In previous studies designed to increase the primary structure symmetry within the hydrophobic core of human acidic fibroblast growth factor (FGF-1) a combination of five mutations were accommodated, resulting in structure, stability and folding kinetic properties similar to wild-type (despite the symmetric constraint upon the set of core residues). A sixth mutation in the core, involving a highly conserved Met residue at position 67, appeared intolerant to substitution. Structural analysis suggested that the local packing environment of position 67 involved two regions of apparent insertions that distorted the tertiary structure symmetry inherent in the beta-trefoil architecture. It was postulated that a symmetric constraint upon the primary structure within the core could only be achieved after these insertions had been deleted (concomitantly increasing the tertiary structure symmetry). The deletion of these insertions is now shown to permit mutation of position 67, thereby increasing the primary structure symmetry relationship within the core. Furthermore, despite the imposed symmetric constraint upon both the primary and tertiary structure, the resulting mutant form of FGF-1 is substantially more stable. The apparent inserted regions are shown to be associated with heparin-binding functionality; however, despite a marked reduction in heparin-binding affinity the mutant form of FGF-1 is surprisingly approximately 70 times more potent in 3T3 fibroblast mitogenic assays. The results support the hypothesis that primary structure symmetry within a symmetric protein superfold represents a possible solution, rather than a constraint, to achieving a foldable polypeptide.     

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

Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil hyperfamily and exhibits a characteristic threefold symmetry of the tertiary structure. However, evidence of this symmetry is not readily apparent at the level of the primary sequence. This suggests that while selective pressures may exist to retain (or converge upon) a symmetric tertiary structure, other selective pressures have resulted in divergence of the primary sequence during evolution. Using intra-chain and homologue sequence comparisons for 19 members of this family of proteins, we have designed mutants of FGF-1 that constrain a subset of core-packing residues to threefold symmetry at the level of the primary sequence. The consequences of these mutations regarding structure and stability were evaluated using a combination of X-ray crystallography and differential scanning calorimetry. The mutational effects on structure and stability can be rationalized through the characterization of "microcavities" within the core detected using a 1.0A probe radius. The results show that the symmetric constraint within the primary sequence is compatible with a well-packed core and near wild-type stability. However, despite the general maintenance of overall thermal stability, a noticeable increase in non-two-state denaturation follows the increase in primary sequence symmetry. Therefore, properties of folding, rather than stability, may contribute to the selective pressure for asymmetric primary core sequences within symmetric protein architectures.     

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.     

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.     

An atomic resolution structure for human fibroblast growth factor 1

A 1.10-A atomic resolution X-ray structure of human fibroblast growth factor 1 (FGF-1), a member of the beta-trefoil superfold, has been determined. The beta-trefoil is one of 10 fundamental protein superfolds and is the only superfold to exhibit 3-fold structural symmetry (comprising 3 "trefoil" units). The quality of the diffraction data permits unambiguous assignment of Asn, Gln, and His rotamers, Pro ring pucker, as well as refinement of atomic anisotropic displacement parameters (ADPs). The FGF-1 structure exhibits numerous core-packing defects, detectable using a 1.0-A probe radius. In addition to contributing to the relatively low thermal stability of FGF-1, these defects may also permit domain motions within the structure. The availability of refined ADPs allows a translation/libration/screw (TLS) analysis of putative rigid body domains. The TLS analysis shows that beta-strands 6-12 together form a rigid body, and there is a clear demarcation in TLS motions between the adjacent carboxyl- and amino-termini. Although separate from beta-strands 6-12, the individual beta-strands 1-5 do not exhibit correlated motions; thus, this region appears to be comparatively flexible. The heparin-binding contacts of FGF-1 are located within beta-strands 6-12; conversely, a significant portion of the receptor-binding contacts are located within beta-strands 1-5. Thus, the observed rigid body motion in FGF-1 appears related to the ligand-binding functionalities.     

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.     

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.     

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.     

Extent of structural asymmetry in homodimeric proteins: prevalence and relevance

Most homodimeric proteins have symmetric structure. Although symmetry is known to confer structural and functional advantage, asymmetric organization is also observed. Using a non-redundant dataset of 223 high-resolution crystal structures of biologically relevant homodimers, we address questions on the prevalence and significance of asymmetry. We used two measures to quantify global and interface asymmetry, and assess the correlation of several molecular and structural parameters with asymmetry. We have identified rare cases (11/223) of biologically relevant homodimers with pronounced global asymmetry. Asymmetry serves as a means to bring about 2:1 binding between the homodimer and another molecule; it also enables cellular signalling arising from asymmetric macromolecular ligands such as DNA. Analysis of these cases reveals two possible mechanisms by which possible infinite array formation is prevented. In case of homodimers associating via non-topologically equivalent surfaces in their tertiary structures, ligand-dependent mechanisms are used. For stable dimers binding via large surfaces, ligand-dependent structural change regulates polymerisation/depolymerisation; for unstable dimers binding via smaller surfaces that are not evolutionarily well conserved, dimerisation occurs only in the presence of the ligand. In case of homodimers associating via interaction surfaces with parts of the surfaces topologically equivalent in the tertiary structures, steric hindrance serves as the preventive mechanism of infinite array. We also find that homodimers exhibiting grossly symmetric organization rarely exhibit either perfect local symmetry or high local asymmetry. Binding of small ligands at the interface does not cause any significant variation in interface asymmetry. However, identification of biologically relevant interface asymmetry in grossly symmetric homodimers is confounded by the presence of similar small magnitude changes caused due to artefacts of crystallisation. Our study provides new insights regarding accommodation of asymmetry in homodimers.     

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.     

Assembly of the archaeal box C/D sRNP can occur via alternative pathways and requires temperature-facilitated sRNA remodeling

Archaeal dual-guide box C/D small nucleolar RNA-like RNAs (sRNAs) bind three core proteins in sequential order at both terminal box C/D and internal C'/D' motifs to assemble two ribonuclear protein (RNP) complexes active in guiding nucleotide methylation. Experiments have investigated the process of box C/D sRNP assembly and the resultant changes in sRNA structure or "remodeling" as a consequence of sRNP core protein binding. Hierarchical assembly of the Methanocaldococcus jannaschii sR8 box C/D sRNP is a temperature-dependent process with binding of L7 and Nop56/58 core proteins to the sRNA requiring elevated temperature to facilitate necessary RNA structural dynamics. Circular dichroism (CD) spectroscopy and RNA thermal denaturation revealed an increased order and stability of sRNA folded structure as a result of L7 binding. Subsequent binding of the Nop56/58 and fibrillarin core proteins to the L7-sRNA complex further remodeled sRNA structure. Assessment of sR8 guide region accessibility using complementary RNA oligonucleotide probes revealed significant changes in guide region structure during sRNP assembly. A second dual-guide box C/D sRNA from M. jannaschii, sR6, also exhibited RNA remodeling during temperature-dependent sRNP assembly, although core protein binding was affected by sR6's distinct folded structure. Interestingly, the sR6 sRNP followed an alternative assembly pathway, with both guide regions being continuously exposed during sRNP assembly. Further experiments using sR8 mutants possessing alternative guide regions demonstrated that sRNA folded structure induced by specific guide sequences impacted the sRNP assembly pathway. Nevertheless, assembled sRNPs were active for sRNA-guided methylation independent of the pathway followed. Thus, RNA remodeling appears to be a common and requisite feature of archaeal dual-guide box C/D sRNP assembly and formation of the mature sRNP can follow different assembly pathways in generating catalytically active complexes.     

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.     

Unraveling the symmetry ambiguity in a hexamer: Calculation of the R6 human insulin structure

Crystallographic and NMR studies of insulin have revealed a highly flexible molecule with a range of different aggregation and structural states; the importance of these states for the function of the hormone is still unclear. To address this question, we have studied the solution structure of the insulin R₆ symmetric hexamer using NMR spectroscopy. Structure determination of symmetric oligomers by NMR is complicated due to `symmetry ambiguity' between intra- and intermonomer NOEs, and between different classes of intermonomer NOEs. Hence, to date, only two symmetric tetramers and one symmetric pentamer (VTB, B subunit of verotoxin) have been solved by NMR; there has been no other symmetric hexamer or higher-order oligomer. Recently, we reported a solution structure for R₆ insulin hexamer. However, in that study, a crystal structure was used as a reference to resolve ambiguities caused by the threefold symmetry; the same method was used in solving VTB. Here, we have successfully recalculated R₆ insulin using the symmetry-ADR method, a computational strategy in which ambiguities are resolved using the NMR data alone. Thus the obtained structure is a refinement of the previous R₆ solution structure. Correlated motions in the final structural ensemble were analysed using a recently developed principal component method; this suggests the presence of two major conformational substates. The study demonstrates that the solution structure of higher-order symmetric oligomers can be determined unambiguously from NMR data alone, using the symmetry-ADR method. This success bodes well for future NMR studies of higher-order symmetric oligomers. The correlated motions observed in the structural ensemble suggest a new insight into the mechanism of phenol exchange and the T ₆ R ₆ transition of insulin in solution.     

Topological frustration and the folding of interleukin-1 beta

The cytokine, interleukin-1beta (IL-1beta), adopts a beta-trefoil fold. It is known to be much slower folding than similarly sized proteins, despite having a low contact order. Proteins are sufficiently well designed that their folding is not dominated by local energetic traps. Therefore, protein models that encode only the folded structure and are energetically unfrustrated (Gō-type), can capture the essentials of the folding routes. We investigate the folding thermodynamics of IL-1beta using such a model and molecular dynamics (MD) simulations. We develop an enhanced sampling technique (a modified multicanonical method) to overcome the sampling problem caused by the slow folding. We find that IL-1beta has a broad and high free energy barrier. In addition, the protein fold causes intermediate unfolding and refolding of some native contacts within the protein along the folding trajectory. This "backtracking" occurs around the barrier region. Complex folds like the beta-trefoil fold and functional loops like the beta-bulge of IL-1beta can make some of the configuration space unavailable to the protein and cause topological frustration.     

Thermodynamic characterization of the Saccharomyces cerevisiae telomerase RNA pseudoknot domain in vitro

Recent structural and functional characterization of the pseudoknot in the Saccharomyces cerevisiae telomerase RNA (TLC1) has demonstrated that tertiary structure is present, similar to that previously described for the human and Kluyveromyces lactis telomerase RNAs. In order to biophysically characterize the identified pseudoknot secondary and tertiary structures, UV-monitored thermal denaturation experiments, nuclear magnetic resonance spectroscopy, and native gel electrophoresis were used to investigate various potential conformations in the pseudoknot domain in vitro, in the absence of the telomerase protein. Here, we demonstrate that alternative secondary structures are not mutually exclusive in the S. cerevisiae telomerase RNA, tertiary structure contributes 1.5 kcal mol(-1) to the stability of the pseudoknot (≈ half the stability observed for the human telomerase pseudoknot), and identify additional base pairs in the 3' pseudoknot stem near the helical junction. In addition, sequence conservation in an adjacent overlapping hairpin appears to prevent dimerization and alternative conformations in the context of the entire pseudoknot-containing region. Thus, this work provides a detailed in vitro characterization of the thermodynamic features of the S. cerevisiae TLC1 pseudoknot region for comparison with other telomerase RNA pseudoknots.     

Fold recognition and accurate sequence-structure alignment of sequences directing beta-sheet proteins

The ability to predict structure from sequence is particularly important for toxins, virulence factors, allergens, cytokines, and other proteins of public health importance. Many such functions are represented in the parallel beta-helix and beta-trefoil families. A method using pairwise beta-strand interaction probabilities coupled with evolutionary information represented by sequence profiles is developed to tackle these problems for the beta-helix and beta-trefoil folds. The algorithm BetaWrapPro employs a "wrapping" component that may capture folding processes with an initiation stage followed by processive interaction of the sequence with the already-formed motifs. BetaWrapPro outperforms all previous motif recognition programs for these folds, recognizing the beta-helix with 100% sensitivity and 99.7% specificity and the beta-trefoil with 100% sensitivity and 92.5% specificity, in crossvalidation on a database of all nonredundant known positive and negative examples of these fold classes in the PDB. It additionally aligns 88% of residues for the beta-helices and 86% for the beta-trefoils accurately (within four residues of the exact position) to the structural template, which is then used with the side-chain packing program SCWRL to produce 3D structure predictions. One striking result has been the prediction of an unexpected parallel beta-helix structure for a pollen allergen, and its recent confirmation through solution of its structure. A Web server running BetaWrapPro is available and outputs putative PDB-style coordinates for sequences predicted to form the target folds.     

Conserved and nonconserved features of the folding pathway of hisactophilin, a beta-trefoil protein

Based on previous studies of interleukin-1beta (IL-1beta) and both acidic and basic fibroblast growth factors (FGFs), it has been suggested that the folding of beta-trefoil proteins is intrinsically slow and may occur via the formation of essential intermediates. Using optical and NMR-detected quenched-flow hydrogen/deuterium exchange methods, we have measured the folding kinetics of hisactophilin, another beta-trefoil protein that has < 10% sequence identity and unrelated function to IL-1beta and FGFs. We find that hisactophilin can fold rapidly and with apparently two-state kinetics, except under the most stabilizing conditions investigated where there is evidence for formation of a folding intermediate. The hisactophilin intermediate has significant structural similarities to the IL-1beta intermediate that has been observed experimentally and predicted theoretically using a simple, topology-based folding model; however, it appears to be different from the folding intermediate observed experimentally for acidic FGF. For hisactophilin and acidic FGF, intermediates are much less prominent during folding than for IL-1beta. Considering the structures of the different beta-trefoil proteins, it appears that differences in nonconserved loops and hydrophobic interactions may play an important role in differential stabilization of the intermediates for these proteins.     

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.     

Intrinsic asymmetry of oligomer transitions and biomolecular evolution

Structural transitions in oligomeric proteins due to ligand binding are important in biomolecular regulatory processes. The transitions may occur on the secondary, tertiary or quarternary structure levels. Detailed consideration of the time sequence of ligand binding to the oligomer shows that there is an intrinsic dynamic asymmetry in all oligomer transitions, even if the initial and the final state are completely symmetric. This asymmetry has important bearing on the evolution and the divergence of the primary structure (amino acid sequence) of oligomeric proteins. It may explain (at least in part) the occurrence of oligomeric proteins with similar but not identical protomers. Certain specific groups of oligomers are shown to be under greater evolutionary pressure for protomer structure divergence. The dynamic asymmetry of oligomer transitions also results in higher complexity in reaction kinetics. Some implications on ribosome structural evolution are discussed.