CPSARST: an efficient circular permutation search tool applied to the detection of novel protein structural relationships

Circular permutation of a protein can be visualized as if the original amino- and carboxyl termini were linked and new ones created elsewhere. It has been well-documented that circular permutants usually retain native structures and biological functions. Here we report CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) to be an efficient database search tool. In this post-genomics era, when the amount of protein structural data is increasing exponentially, it provides a new way to rapidly detect novel relationships among proteins.     

Effect of circular permutations on transient partial unfolding in proteins

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three‐dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N‐termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native‐state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild‐type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild‐type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild‐type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.     

Folding circular permutants of IL-1β: route selection driven by functional frustration

Interleukin-1β (IL-1β) is the cytokine crucial to inflammatory and immune response. Two dominant routes are populated in the folding to native structure. These distinct routes are a result of the competition between early packing of the functional loops versus closure of the β-barrel to achieve efficient folding and have been observed both experimentally and computationally. Kinetic experiments on the WT protein established that the dominant route is characterized by early packing of geometrically frustrated functional loops. However, deletion of one of the functional loops, the β-bulge, switches the dominant route to an alternative, yet, as accessible, route, where the termini necessary for barrel closure form first. Here, we explore the effect of circular permutation of the WT sequence on the observed folding landscape with a combination of kinetic and thermodynamic experiments. Our experiments show that while the rate of formation of permutant protein is always slower than that observed for the WT sequence, the region of initial nucleation for all permutants is similar to that observed for the WT protein and occurs within a similar timescale. That is, even permutants with significant sequence rearrangement in which the functional-nucleus is placed at opposing ends of the polypeptide chain, fold by the dominant WT "functional loop-packing route", despite the entropic cost of having to fold the N- and C- termini early. Taken together, our results indicate that the early packing of the functional loops dominates the folding landscape in active proteins, and, despite the entropic penalty of coalescing the termini early, these proteins will populate an entropically unfavorable route in order to conserve function. More generally, circular permutation can elucidate the influence of local energetic stabilization of functional regions within a protein, where topological complexity creates a mismatch between energetics and topology in active proteins.     

Circular permutation of red fluorescent proteins

Circular permutation of fluorescent proteins provides a substrate for the design of molecular sensors. Here we describe a systematic exploration of permutation sites for mCherry and mKate using a tandem fusion template approach. Circular permutants retaining more than 60% (mCherry) and 90% (mKate) brightness of the parent molecules are reported, as well as a quantitative evaluation of the fluorescence from neighboring mutations. Truncations of circular permutants indicated essential N- and C-terminal segments and substantial flexibility in the use of these molecules. Structural evaluation of two cp-mKate variants indicated no major conformational changes from the previously reported wild-type structure, and cis conformation of the chromophores. Four cp-mKates were identified with over 80% of native fluorescence, providing important new building blocks for sensor and complementation experiments.     

High-resolution structure prediction of a circular permutation loop

Methods for rapid and reliable design and structure prediction of linker loops would facilitate a variety of protein engineering applications. Circular permutation, in which the existing termini of a protein are linked by the polypeptide chain and new termini are created, is one such application that has been employed for decreasing proteolytic susceptibility and other functional purposes. The length and sequence of the linker can impact the expression level, solubility, structure and function of the permuted variants. Hence it is desirable to achieve atomic‐level accuracy in linker design. Here, we describe the use of RosettaRemodel for design and structure prediction of circular permutation linkers on a model protein. A crystal structure of one of the permuted variants confirmed the accuracy of the computational prediction, where the all‐atom rmsd of the linker region was 0.89 Å between the model and the crystal structure. This result suggests that RosettaRemodel may be generally useful for the design and structure prediction of protein loop regions for circular permutations or other structure‐function manipulations.     

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.     

Random circular permutation of DsbA reveals segments that are essential for protein folding and stability

One of the key questions in protein folding is whether polypeptide chains require unique nucleation sites to fold to the native state. In order to identify possible essential polypeptide segments for folding, we have performed a complete circular permutation analysis of a protein in which the natural termini are in close proximity. As a model system, we used the disulfide oxidoreductase DsbA from Escherichia coli, a monomeric protein of 189 amino acid residues. To introduce new termini at all possible positions in its polypeptide chain, we generated a library of randomly circularly permuted dsbA genes and screened for active circularly permuted variants in vivo. A total of 51 different active variants were identified. The new termini were distributed over about 70 % of the polypeptide chain, with the majority of them occurring within regular secondary structures. New termini were not found in approximately 30 % of the DsbA sequence which essentially correspond to four alpha-helices of DsbA. Introduction of new termini into these "forbidden segments" by directed mutagenesis yielded proteins with altered overall folds and strongly reduced catalytic activities. In contrast, all active variants analysed so far show structural and catalytic properties comparable with those of DsbA wild-type. We suggest that random circular permutation allows identification of contiguous structural elements in a protein that are essential for folding and stability.     

Restriction enzyme-free construction of random gene mutagenesis libraries in Escherichia coli

Directed evolution relies on both random and site-directed mutagenesis of individual genes and regulatory elements to create variants with altered activity profiles for engineering applications. Central to these experiments is the construction of large libraries of related variants. However, a number of technical hurdles continue to limit routine construction of random mutagenesis libraries in Escherichia coli, in particular, inefficiencies during digestion and ligation steps. Here, we report a restriction enzyme-free approach to library generation using megaprimers termed MegAnneal. Target DNA is first exponentially amplified using error-prone polymerase chain reaction (PCR) and then linearly amplified with a single 3′ primer to generate long, randomly mutated, single-stranded megaprimers. These are annealed to single-stranded dUTP-containing template plasmid and extended with T7 polymerase to create a complementary strand, and the resulting termini are ligated with T4 DNA ligase. Using this approach, we are able to reliably generate libraries of approximately 10⁷ colony-forming units (cfu)/μg DNA/transformation in a single day. We have created MegAnneal libraries based on three different single-chain antibodies and identified variants with enhanced expression and ligand-binding affinity. The key advantages of this approach include facile amplification, restriction enzyme-free library generation, and a significantly reduced risk of mutations outside the targeted region and wild-type contamination as compared with current methods.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Common motifs and topological effects in the protein folding transition state

Through extensive experiment, simulation, and analysis of protein S6 (1RIS), we find that variations in nucleation and folding pathway between circular permutations are determined principally by the restraints of topology and specific nucleation, and affected by changes in chain entropy. Simulations also relate topological features to experimentally measured stabilities. Despite many sizable changes in phi values and the structure of the transition state ensemble that result from permutation, we observe a common theme: the critical nucleus in each of the mutants share a subset of residues that can be mapped to the critical nucleus residues of the wild-type. Circular permutations create new N and C termini, which are the location of the largest disruption of the folding nucleus, leading to a decrease in both phi values and the role in nucleation. Mutant nuclei are built around the wild-type nucleus but are biased towards different parts of the S6 structure depending on the topological and entropic changes induced by the location of the new N and C termini.     

Use of an engineered ribozyme to produce a circular human exon.

We report the use of an engineered ribozyme to produce a circular human exon in vitro. Specifically, we have designed a derivative of a yeast self-splicing group II intron that is able to catalyze the formation of a circular exon encoding the first kringle domain (K1) of the human tissue plasminogen activator protein. We show that the circular K1 exon is formed with high fidelity in vitro. Furthermore, the system is designed such that the circular exon that is produced consists entirely of human exon sequence. Thus, our results demonstrate that all yeast exon sequences are dispensable for group II intron catalyzed inverse splicing. This is the first demonstration that an engineered ribozyme can be used to create a circular exon containing only human sequences, linked together at a precise desired ligation point. We expect these results to be generalizable, so that similar ribozymes can be designed to precisely create circular derivatives of any nucleotide sequence.     

Autonomous replication of foreign DNA in Histoplasma capsulatum: role of native telomeric sequences.

Genetic transformation of the dimorphic pathogenic fungus Histoplasma capsulatum can result in chromosomal integration of the transforming DNA or the generation of multicopy linear plasmids carrying the transforming DNA. We showed previously that Escherichia coli plasmids do not replicate autonomously in H. capsulatum without significant modifications, one of which is the in vivo addition of Histoplasma telomeres at the termini of linear DNA. To address the requirements for autonomous replication in H. capsulatum, we constructed a circular E. coli plasmid containing adjacent inverted stretches of Histoplasma telomeric repeats separated by a unique restriction site. The linearized plasmid bearing telomeric termini was maintained in H. capsulatum without modification other than the addition of more telomeric sequence. We recovered the original plasmid in E. coli after removal of the telomeric termini by using engineered restriction sites. Thus, no special Histoplasma modification or sequence other than the telomeres was needed for autonomous replication in H. capsulatum. Additionally, this plasmid provides a shuttle vector that replicates autonomously in E. coli (as a circular plasmid) and in H. capsulatum (as a linear plasmid).     

Quantitative in vivo solubility and reconstitution of truncated circular permutants of green fluorescent protein

Several versions of split green fluorescent protein (GFP) fold and reconstitute fluorescence, as do many circular permutants, but little is known about the dependence of reconstitution on circular permutation. Explored here is the capacity of GFP to fold and reconstitute fluorescence from various truncated circular permutants, herein called "leave-one-outs" using a quantitative in vivo solubility assay and in vivo reconstitution of fluorescence. Twelve leave-one-out permutants are discussed, one for each of the 12 secondary structure elements. The results expand the outlook for the use of permuted split GFPs as specific and self-reporting gene encoded affinity reagents.     

Circular permutation in the Ω-loop of TEM-1 β-lactamase results in improved activity and altered substrate specificity

Generating diverse protein libraries that contain improved variants at a sufficiently high frequency is critical for improving the properties of proteins using directed evolution. Many studies have illustrated how random mutagenesis, cassette mutagenesis, DNA shuffling and similar approaches are effective diversity generating methods for directed evolution. Very few studies have explored random circular permutation, the intramolecular relocation of the N- and C-termini of a protein, as a diversity-generating step for directed evolution. We subjected a library of random circular permutations of TEM-1 β-lactamase to selections on increasing concentrations of a variety of β-lactam antibiotics including cefotaxime. We identified two circularly permuted variants that conferred elevated resistance to cefotaxime but decreased resistance to other antibiotics. These variants were circularly permuted in the Ω-loop proximal to the active site. Remarkably, one variant was circularly permuted such that the key catalytic residue Glu166 was located at the N-terminus of the mature protein.     

Analysis of extrachromosomal Ac/Ds transposable elements

The mechanism of transposition of the maize Ac/Ds elements is not well understood. The true transposition intermediates are not known and it has not been possible to distinguish between excision models involving 8-bp staggered cuts or 1-bp staggered cuts followed by hairpin formation. In this work, we have analyzed extrachromosomal excision products to gain insight into the excision mechanism. Plasmid rescue was used to demonstrate that Ds excision is associated with the formation of circular molecules. In addition, we present evidence for the formation of linear extrachromosomal species during Ds excision. Sequences found at the termini of circular and linear elements showed a broad range of nucleotide additions or deletions, suggesting that these species are not true intermediates. Additional nucleotides adjacent to the termini in extrachromosomal elements were compared to the sequence of the original donor site. This analysis showed that: (1) the first nucleotide adjacent to the transposon end was significantly more similar to the first nucleotide flanking the element in the donor site than to a random sequence and (2) the second and farther nucleotides did not resemble the donor site. The implications of these findings for excision models are discussed.     

Acyclic permutants of naturally occurring cyclic proteins. Characterization of cystine knot and beta-sheet formation in the macrocyclic polypeptide kalata B1

Kalata B1 is a prototypic member of the unique cyclotide family of macrocyclic polypeptides in which the major structural features are a circular peptide backbone, a triple-stranded beta-sheet, and a cystine knot arrangement of three disulfide bonds. The cyclotides are the only naturally occurring family of circular proteins and have prompted us to explore the concept of acyclic permutation, i.e. opening the backbone of a cross-linked circular protein in topologically permuted ways. We have synthesized the complete suite of acyclic permutants of kalata B1 and examined the effect of acyclic permutation on structure and activity. Only two of six topologically distinct backbone loops are critical for folding into the native conformation, and these involve disruption of the embedded ring in the cystine knot. Surprisingly, it is possible to disrupt regions of the beta-sheet and still allow folding into native-like structure, provided the cystine knot is intact. Kalata B1 has mild hemolytic activity, but despite the overall structure of the native peptide being retained in all but two cases, none of the acyclic permutants displayed hemolytic activity. This loss of activity is not localized to one particular region and suggests that cyclization is critical for hemolytic activity.     

Chemical synthesis of circular proteins

Circular proteins, once thought to be rare, are now commonly found in plants. Their chemical synthesis, once thought to be difficult, is now readily achievable. The enabling methodology is largely due to the advances in entropic chemical ligation to overcome the entropy barrier in coupling the N- and C-terminal ends of large peptide segments for either intermolecular ligation or intramolecular ligation in end-to-end cyclization. Key elements of an entropic chemical ligation consist of a chemoselective capture step merging the N and C termini as a covalently linked O/S-ester intermediate to permit the subsequent step of an intramolecular O/S-N acyl shift to form an amide. Many ligation methods exploit the supernucleophilicity of a thiol side chain at the N terminus for the capture reaction, which makes cysteine-rich peptides ideal candidates for the entropy-driven macrocyclization. Advances in desulfurization and modification of the thiol-containing amino acids at the ligation sites to other amino acids add extra dimensions to the entropy-driven ligation methods. This minireview describes recent advances of entropy-driven ligation to prepare circular proteins with or without a cysteinyl side chain.     

Split-intein mediated circular ligation used in the synthesis of cyclic peptide libraries in E. coli

Recent advances in chemical biology and the advantages presented by in vivo screening have highlighted the need for a robust and flexible biologically synthesized small-molecule library. Herein we describe a method for the biosynthesis of cyclic peptide libraries of up to 10(8) members in Escherichia coli using split-intein circular ligation of peptides and proteins (SICLOPPS). The method utilizes split-intein chemistry to cyclize randomized peptide sequences. The cyclic peptide library can potentially be of any size and the peptide itself may contain unlimited random residues. However, the library size is limited by the transformation efficiency of E. coli and random residues are generally limited to five, but additional amino acids can be used in the cyclic peptide backbone, varying the structure and ring size of the cyclic peptide. SICLOPPS libraries have been combined with a bacterial reverse two-hybrid system in our labs and used in the identification of inhibitors of several protein-protein interactions. This protocol is expected to take around 3-4 weeks to implement.