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.     

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.     

Circular permutation: a different way to engineer enzyme structure and function

Protein engineers traditionally rely on amino acid substitutions to alter the functional properties of biomacromolecules, yet have largely overlooked the potential benefits of reorganizing the polypeptide chain of a protein by circular permutation (CP). By connecting the native protein termini via a covalent linker and introducing new ends through the cleavage of an existing peptide bond, CP can perturb local tertiary structure and protein dynamics, as well as introduce possible quaternary structure changes. In several recent studies, these effects have successfully been exploited to manipulate protein scaffolds, resulting in improved catalytic activity and altered substrate or ligand binding affinity, as well as enabling the design of novel biocatalysts and biosensors.     

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.     

Construction of a random circular permutation library using an engineered transposon

Circular permutation is an important protein engineering tool used to create sequence diversity of a protein by changing its linear order of amino acid sequence. Circular permutation has proven to be effective in the evolution of proteins for desired properties while maintaining similar three-dimensional structures. Due to the lack of a robust design principle guiding the selection of new termini, construction of a combinatorial library is much preferred for comprehensive evaluation of circular permutation. Unfortunately, the conventional methods used to create random circular permutation libraries cause significant sequence modification at new termini of circular permutants. In addition, these methods impose additional limitations by requiring either relatively inefficient blunt-end ligation during library construction or redesign of transposons for tailored expression of circular permutants. In this study, we present the development of an engineered transposon for facile construction of random circular permutation libraries. We provide evidence that minimal modification at the new termini of the random circular permutants is possible with our engineered transposon. In addition, our method enables the use of sticky-end ligation during library construction and provides external tunability for expression of random circular permutants.     

Circular permutation of ligand‐binding module improves dynamic range of genetically encoded FRET‐based nanosensor

Quantitative measurement of small molecules with high spatiotemporal resolution provides a solid basis for correct understanding and accurate modeling of metabolic regulation. A promising approach toward this goal is the FLIP (fluorescent indicator protein) nanosensor based on bacterial periplasmic binding proteins (PBPs) and fluorescence resonance energy transfer (FRET) between the yellow and cyan variants of green fluorescent protein (GFP). Each FLIP has a PBP module that specifically binds its ligand to induce a conformation change, leading to a change in FRET between the two GFP variant modules attached to the N‐ and C‐termini of the PBP. The larger is the dynamic range the more reliable is the measurement. Thus, we attempted to expand the dynamic range of FLIP by introducing a circular permutation with a hinge loop deletion to the PBP module. All the six circularly permutated PBPs tested, including structurally distinct Type I and Type II PBPs, showed larger dynamic ranges than their respective native forms when used for FLIP. Notably, the circular permutation made three PBPs, which totally failed to show FRET change when used as their native forms, fully capable of functioning as a ligand binding module of FLIP. These FLIPs were successfully used for the determination of amino acid concentration in complex solutions as well as real‐time measurement of amino acid influx in living yeast cells. Thus, the circular permutation strategy would not only improve the performance of each nanosensor but also expand the repertoire of metabolites that can be measured by the FLIP nanosensor technology.     

Solution structures and backbone dynamics of the ribosomal protein S6 and its permutant P54‐55

The ribosomal protein S6 from Thermus thermophilus has served as a model system for the study of protein folding, especially for understanding the effects of circular permutations of secondary structure elements. This study presents the structure of a permutant protein, the 96‐residue P^54‐55, and the structure of its 101‐residue parent protein S6^wt in solution. The data also characterizes the effects of circular permutation on the backbone dynamics of S6. Consistent with crystallographic data on S6^wt, the overall solution structures of both P^54‐55 and S6^wt show a β‐sheet of four antiparallel β‐strands with two α‐helices packed on one side of the sheet. In clear contrast to the crystal data, however, the solution structure of S6^wt reveals a disordered loop in the region between β‐strands 2 and 3 (Leu43‐Phe60) instead of a well‐ordered stretch and associated hydrophobic mini‐core observed in the crystal structure. Moreover, the data for P^54‐55 show that the joined wild‐type N‐ and C‐terminals form a dynamically robust stretch with a hairpin structure that complies with the in silico design. Taken together, the results explain why the loop region of the S6^wt structure is relatively insensitive to mutational perturbations, and why P^54‐55 is more stable than S6^wt: the permutant incision at Lys54‐Asp55 is energetically neutral by being located in an already disordered loop whereas the new hairpin between the wild‐type N‐ and C‐termini is stabilizing.     

Circularly permuted monomeric red fluorescent proteins with new termini in the ��-sheet

Circularly permuted fluorescent proteins (FPs) have a growing number of uses in live cell fluorescence biosensing applications. Most notably, they enable the construction of single fluorescent protein‐based biosensors for Ca²⁺ and other analytes of interest. Circularly permuted FPs are also of great utility in the optimization of fluorescence resonance energy transfer (FRET)‐based biosensors by providing a means for varying the critical dipole–dipole orientation. We have previously reported on our efforts to create circularly permuted variants of a monomeric red FP (RFP) known as mCherry. In our previous work, we had identified six distinct locations within mCherry that tolerated the insertion of a short peptide sequence. Creation of circularly permuted variants with new termini at the locations corresponding to the sites of insertion led to the discovery of three permuted variants that retained no more than 18% of the brightness of mCherry. We now report the extensive directed evolution of the variant with new termini at position 193 of the protein sequence for improved fluorescent brightness. The resulting variant, known as cp193g7, has 61% of the intrinsic brightness of mCherry and was found to be highly tolerant of circular permutation at other locations within the sequence. We have exploited this property to engineer an expanded series of circularly permuted variants with new termini located along the length of the 10th β‐strand of mCherry. These new variants may ultimately prove useful for the creation of single FP‐based Ca²⁺ biosensors.     

Structural Redesign of Lipase B from Candida antarctica by Circular Permutation and Incremental Truncation

Circular permutation of Candida antarctica lipase B yields several enzyme variants with substantially increased catalytic activity. To better understand the structural and functional consequences of protein termini reorganization, we have applied protein engineering and x-ray crystallography to cp283, one of the most active hydrolase variants. Our initial investigation has focused on the role of an extended surface loop, created by linking the native N- and C-termini, on protein integrity. Incremental truncation of the loop partially compensates for observed losses in secondary structure and the permutants' temperature of unfolding. Unexpectedly, the improvements are accompanied by quaternary-structure changes from monomer to dimer. The crystal structures of one truncated variant (cp283Δ7) in the apo-form determined at 1.49 Å resolution and with a bound phosphonate inhibitor at 1.69 Å resolution confirmed the formation of a homodimer by swapping of the enzyme's 35-residue N-terminal region. Separately, the new protein termini at amino acid positions 282/283 convert the narrow access tunnel to the catalytic triad into a broad crevice for accelerated substrate entry and product exit while preserving the native active-site topology for optimal catalytic turnover.     

RosettaRemodel: a generalized framework for flexible backbone protein design

We describe RosettaRemodel, a generalized framework for flexible protein design that provides a versatile and convenient interface to the Rosetta modeling suite. RosettaRemodel employs a unified interface, called a blueprint, which allows detailed control over many aspects of flexible backbone protein design calculations. RosettaRemodel allows the construction and elaboration of customized protocols for a wide range of design problems ranging from loop insertion and deletion, disulfide engineering, domain assembly, loop remodeling, motif grafting, symmetrical units, to de novo structure modeling.     

Structure of a designed tetrahedral protein assembly variant engineered to have improved soluble expression

We recently reported the development of a computational method for the design of coassembling multicomponent protein nanomaterials. While four such materials were validated at high‐resolution by X‐ray crystallography, low yield of soluble protein prevented X‐ray structure determination of a fifth designed material, T33‐09. Here we report the design and crystal structure of T33‐31, a variant of T33‐09 with improved soluble yield resulting from redesign efforts focused on mutating solvent‐exposed side chains to charged amino acids. The structure is found to match the computational design model with atomic‐level accuracy, providing further validation of the design approach and demonstrating a simple and potentially general means of improving the yield of designed protein nanomaterials.     

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.     

Sequence determinants of thermodynamic stability in a WW domain—An all‐β‐sheet protein

The stabilities of 66 sequence variants of the human Pin1 WW domain have been determined by equilibrium thermal denaturation experiments. All 34 residues composing the hPin1 WW three‐stranded β‐sheet structure could be replaced one at a time with at least one different natural or non‐natural amino acid residue without leading to an unfolded protein. Alanine substitutions at only four positions within the hPin1 WW domain lead to a partially or completely unfolded protein—in the absence of a physiological ligand. The side chains of these four residues form a conserved, partially solvent‐inaccessible, continuous hydrophobic minicore comprising the N‐ and C‐termini. Ala mutations at five other residues, three of which constitute the ligand binding patch on the concave side of the β‐sheet, significantly destabilize the hPin1 WW domain without leading to an unfolded protein. The remaining mutations affect protein stability only slightly, suggesting that only a small subset of side chain interactions within the hPin1 WW domain are mandatory for acquiring and maintaining a stable, cooperatively folded β‐sheet structure.     

Modeling large regions in proteins: applications to loops, termini, and folding

Template-based methods for predicting protein structure provide models for a significant portion of the protein but often contain insertions or chain ends (InsEnds) of indeterminate conformation. The local structure prediction "problem" entails modeling the InsEnds onto the rest of the protein. A well-known limit involves predicting loops of ≤12 residues in crystal structures. However, InsEnds may contain as many as ~50 amino acids, and the template-based model of the protein itself may be imperfect. To address these challenges, we present a free modeling method for predicting the local structure of loops and large InsEnds in both crystal structures and template-based models. The approach uses single amino acid torsional angle "pivot" moves of the protein backbone with a C(β) level representation. Nevertheless, our accuracy for loops is comparable to existing methods. We also apply a more stringent test, the blind structure prediction and refinement categories of the CASP9 tournament, where we improve the quality of several homology based models by modeling InsEnds as long as 45 amino acids, sizes generally inaccessible to existing loop prediction methods. Our approach ranks as one of the best in the CASP9 refinement category that involves improving template-based models so that they can function as molecular replacement models to solve the phase problem for crystallographic structure determination.     

Toward a “Structural BLAST”: Using structural relationships to infer function

We outline a set of strategies to infer protein function from structure. The overall approach depends on extensive use of homology modeling, the exploitation of a wide range of global and local geometric relationships between protein structures and the use of machine learning techniques. The combination of modeling with broad searches of protein structure space defines a “structural BLAST” approach to infer function with high genomic coverage. Applications are described to the prediction of protein–protein and protein–ligand interactions. In the context of protein–protein interactions, our structure‐based prediction algorithm, PrePPI, has comparable accuracy to high‐throughput experiments. An essential feature of PrePPI involves the use of Bayesian methods to combine structure‐derived information with non‐structural evidence (e.g. co‐expression) to assign a likelihood for each predicted interaction. This, combined with a structural BLAST approach significantly expands the range of applications of protein structure in the annotation of protein function, including systems level biological applications where it has previously played little role.     

Random circular permutation leading to chain disruption within and near alpha helices in the catalytic chains of aspartate transcarbamoylase: effects on assembly, stability, and function

A collection of circularly permuted catalytic chains of aspartate transcarbamoylase (ATCase) has been generated by random circular permutation of the pyrB gene. From the library of ATCases containing permuted polypeptide chains, we have chosen for further investigation nine ATCase variants whose catalytic chains have termini located within or close to an alpha helix. All of the variants fold and assemble into dodecameric holoenzymes with similar sedimentation coefficients and slightly reduced thermal stabilities. Those variants disrupted within three different helical regions in the wild-type structure show no detectable enzyme activity and no apparent binding of the bisubstrate analog N:-phosphonacetyl-L-aspartate. In contrast, two variants whose termini are just within or adjacent to other alpha helices are catalytically active and allosteric. As expected, helical disruptions are more destabilizing than loop disruptions. Nonetheless, some catalytic chains lacking continuity within helical regions can assemble into stable holoenzymes comprising six catalytic and six regulatory chains. For seven of the variants, continuity within the helices in the catalytic chains is important for enzyme activity but not necessary for proper folding, assembly, and stability of the holoenzyme.     

Effect of circular permutation on the structure and function of type 1 blue copper center in azurin

Type 1 copper (T1Cu) proteins are electron transfer (ET) proteins involved in many important biological processes. While the effects of changing primary and secondary coordination spheres in the T1Cu ET function have been extensively studied, few report has explored the effect of the overall protein structural perturbation on active site configuration or reduction potential of the protein, even though the protein scaffold has been proposed to play a critical role in enforcing the entatic or “rack‐induced” state for ET functions. We herein report circular permutation of azurin by linking the N‐ and C‐termini and creating new termini in the loops between 1^st and 2^nd β strands or between 3^rd and 4^th β strands. Characterization by electronic absorption, electron paramagnetic spectroscopies, as well as crystallography and cyclic voltammetry revealed that, while the overall structure and the primary coordination sphere of the circular permutated azurins remain the same as those of native azurin, their reduction potentials increased by 18 and 124 mV over that of WTAz. Such increases in reduction potentials can be attributed to subtle differences in the hydrogen‐bonding network in secondary coordination sphere around the T1Cu center.     

Effect of genetic circular permutation near the active site on the activity and stability of an enzyme inhibitor

We report here the effect of circular permutation on the structure and function of a model protein tendamistat, a 74 amino acid competitive inhibitor of porcine pancreatic alpha-amylase. The activity and stability of wild type and two permuted tendamistat variants were characterized by measurement of alpha-amylase kinetic and thermodynamic binding parameters and their thermodynamics of unfolding. Our results show that large variations in structure and function can occur upon circularly permuting tendamistat near its active site that are not obvious, a priori, from the structure of the native protein and we propose a structural thermodynamic explanation of the experimental observations.     

Design and structure of two new protein cages illustrate successes and ongoing challenges in protein engineering

In recent years, new protein engineering methods have produced more than a dozen symmetric, self‐assembling protein cages whose structures have been validated to match their design models with near‐atomic accuracy. However, many protein cage designs that are tested in the lab do not form the desired assembly, and improving the success rate of design has been a point of recent emphasis. Here we present two protein structures solved by X‐ray crystallography of designed protein oligomers that form two‐component cages with tetrahedral symmetry. To improve on the past tendency toward poorly soluble protein, we used a computational protocol that favors the formation of hydrogen‐bonding networks over exclusively hydrophobic interactions to stabilize the designed protein–protein interfaces. Preliminary characterization showed highly soluble expression, and solution studies indicated successful cage formation by both designed proteins. For one of the designs, a crystal structure confirmed at high resolution that the intended tetrahedral cage was formed, though several flipped amino acid side chain rotamers resulted in an interface that deviates from the precise hydrogen‐bonding pattern that was intended. A structure of the other designed cage showed that, under the conditions where crystals were obtained, a noncage structure was formed wherein a porous 3D protein network in space group I2₁3 is generated by an off‐target twofold homomeric interface. These results illustrate some of the ongoing challenges of developing computational methods for polar interface design, and add two potentially valuable new entries to the growing list of engineered protein materials for downstream applications.     

Loop prediction for a GPCR homology model: Algorithms and results

We present loop structure prediction results of the intracellular and extracellular loops of four G‐protein‐coupled receptors (GPCRs): bovine rhodopsin (bRh), the turkey β1‐adrenergic (β1Ar), the human β2‐adrenergic (β2Ar) and the human A2a adenosine receptor (A2Ar) in perturbed environments. We used the protein local optimization program, which builds thousands of loop candidates by sampling rotamer states of the loops' constituent amino acids. The candidate loops are discriminated between with our physics‐based, all‐atom energy function, which is based on the OPLS force field with implicit solvent and several correction terms. For relevant cases, explicit membrane molecules are included to simulate the effect of the membrane on loop structure. We also discuss a new sampling algorithm that divides phase space into different regions, allowing more thorough sampling of long loops that greatly improves results. In the first half of the paper, loop prediction is done with the GPCRs' transmembrane domains fixed in their crystallographic positions, while the loops are built one‐by‐one. Side chains near the loops are also in non‐native conformations. The second half describes a full homology model of β2Ar using β1Ar as a template. No information about the crystal structure of β2Ar was used to build this homology model. We are able to capture the architecture of short loops and the very long second extracellular loop, which is key for ligand binding. We believe this the first successful example of an RMSD validated, physics‐based loop prediction in the context of a GPCR homology model. Proteins 2013. © 2012 Wiley Periodicals, Inc.