Molecular dynamics study of the stabilities of consensus designed ankyrin repeat proteins

Two designed ankyrin repeat (AR) proteins (E3_5 and E3_19) are high homologous (with about 87% sequence identity) and their crystal structures have a Calpha atom-positional root-mean-square difference of about 0.14 nm. However, it was found that E3_5 is considerably more stable than E3_19 in guanidinium hydrochloride and thermal denaturation experiments. With the goal of providing insights into the various factors contributing to the stabilities of the designed AR proteins and suggesting possible mutations to enhance their stabilities, homology modeling and molecular dynamics (MD) simulations with explicit solvent have been performed. Because the crystal structure of E3_19 was solved later than that of E3_5, a homology model of E3_19 based on the crystal structure of E3_5 was also used in the simulations. E3_5 shows a very stable trajectory in both crystal and solution simulations. In contrast, the C-terminal repeat of E3_19 unfolds in the simulations starting from either the modeled structure or the crystal structure, although it has a sequence identical to that of E3_5. A continuum electrostatic model was used to estimate the effect of single mutations on protein stability and to study the interaction between the internal ARs and the C-terminal capping AR. Mutations involving charged residues were found to have large effects on stability. Due to the difference in charge distribution in the internal ARs of E3_19 and E3_5, their interaction with the C-terminal capping AR is less favorable in E3_19. The simulation trajectories suggest that the stability of the designed AR proteins can be increased by optimizing the electrostatic interactions within and between the different repeats.     

Stabilizing ionic interactions in a full-consensus ankyrin repeat protein

Full-consensus designed ankyrin repeat proteins (DARPins), in which randomized positions of the previously described DARPin library have been fixed, are characterized. They show exceptionally high thermodynamic stabilities, even when compared to members of consensus DARPin libraries and even more so when compared to naturally occurring ankyrin repeat proteins. We determined the crystal structure of a full-consensus DARPin, containing an N-capping repeat, three identical internal repeats and a C-capping repeat at 2.05 Å resolution, and compared its structure with that of the related DARPin library members E3_5 and E3_19. This structural comparison suggests that primarily salt bridges on the surface, which arrange in a network with almost crystal-like regularity, increase thermostability in the full-consensus NI₃C DARPin to make it resistant to boiling. In the crystal structure, three sulfate ions complement this network. Thermal denaturation experiments in guanidine hydrochloride directly indicate a contribution of sulfate binding to the stability, providing further evidence for the stabilizing effect of surface-exposed electrostatic interactions and regular charge networks. The charged residues at the place of randomized residues in the DARPin libraries were selected based on sequence statistics and suggested that the charge interaction network is a hidden design feature of this protein family. Ankyrins can therefore use design principles from proteins of thermophilic organisms and reach at least similar stabilities.     

Folding of a designed simple ankyrin repeat protein

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.     

Characterization and further stabilization of designed ankyrin repeat proteins by combining molecular dynamics simulations and experiments

Multiple molecular dynamics simulations with explicit solvent at room temperature and at 400 K were carried out to characterize designed ankyrin repeat (AR) proteins with full-consensus repeats. Using proteins with one to five repeats, the stability of the native structure was found to increase with the number of repeats. The C-terminal capping repeat, originating from the natural guanine-adenine-binding protein, was observed to denature first in almost all high-temperature simulations. Notably, a stable intermediate is found in experimental equilibrium unfolding studies of one of the simulated consensus proteins. On the basis of simulation results, this intermediate is interpreted to represent a conformation with a denatured C-terminal repeat. To validate this interpretation, constructs without C-terminal capping repeat were prepared and did not show this intermediate in equilibrium unfolding experiments. Conversely, the capping repeats were found to be essential for efficient folding in the cell and for avoiding aggregation, presumably because of their highly charged surface. To design a capping repeat conferring similar solubility properties yet even higher stability, eight point mutations adapting the C-cap to the consensus AR and adding a three-residue extension at the C-terminus were predicted in silico and validated experimentally. The in vitro full-consensus proteins were also compared with a previously published designed AR protein, E3_5, whose internal repeats show 80% identity in primary sequence. A detailed analysis of the simulations suggests that networks of salt bridges between β-hairpins, as well as additional interrepeat hydrogen bonds, contribute to the extraordinary stability of the full consensus.     

Designing repeat proteins: well-expressed,soluble and stable proteins from combinatorial libraries of consensus ankyrin repeat proteins

We describe an efficient way to generate combinatorial libraries of stable, soluble and well-expressed ankyrin repeat (AR) proteins. Using a combination of sequence and structure consensus analyses, we designed a 33 amino acid residue AR module with seven randomized positions having a theoretical diversity of 7.2x10(7). Different numbers of this module were cloned between N and C-terminal capping repeats, i.e. ARs designed to shield the hydrophobic core of stacked AR modules. In this manner, combinatorial libraries of designed AR proteins consisting of four to six repeats were generated, thereby potentiating the theoretical diversity. All randomly chosen library members were expressed in soluble form in the cytoplasm of Escherichia coli in amounts up to 200 mg per 1 l of shake-flask culture. Virtually pure proteins were obtained in a single purification step. The designed AR proteins are monomeric and display CD spectra identical with those of natural AR proteins. At the same time, our AR proteins are highly thermostable, with T(m) values ranging from 66 degrees C to well above 85 degrees C. Thus, our combinatorial library members possess the properties required for biotechnological applications. Moreover, the favorable biophysical properties and the modularity of the AR fold may account, partly, for the abundance of natural AR proteins.     

Self‐association of TPR domains: Lessons learned from a designed,consensus‐based TPR oligomer

The tetratricopeptide repeat (TPR) motif is a protein–protein interaction module that acts as an organizing centre for complexes regulating a multitude of biological processes. Despite accumulating evidence for the formation of TPR oligomers as an additional level of regulation there is a lack of structural and solution data explaining TPR self‐association. In the present work we characterize the trimeric TPR‐containing protein YbgF, which is linked to the Tol system in Gram‐negative bacteria. By subtracting previously identified TPR consensus residues required for stability of the fold from residues conserved across YbgF homologs, we identified residues involved in oligomerization of the C‐terminal YbgF TPR domain. Crafting these residues, which are located in loop regions between TPR motifs, onto the monomeric consensus TPR protein CTPR3 induced the formation of oligomers. The crystal structure of this engineered oligomer shows an asymmetric trimer where stacking interactions between the introduced tyrosines and displacement of the C‐terminal hydrophilic capping helix, present in most TPR domains, are key to oligomerization. Asymmetric trimerization of the YbgF TPR domain and CTPR3Y3 leads to the formation of higher order oligomers both in the crystal and in solution. However, such open‐ended self‐association does not occur in full‐length YbgF suggesting that the protein's N‐terminal coiled‐coil domain restricts further oligomerization. This interpretation is borne out in experiments where the coiled‐coil domain of YbgF was engineered onto the N‐terminus of CTPR3Y3 and shown to block self‐association beyond trimerization. Our study lays the foundations for understanding the structural basis for TPR domain self‐association and how such self‐association can be regulated in TPR domain‐containing proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Structural insights into the evolution of a non-biological protein: importance of surface residues in protein fold optimization

Phylogenetic profiling of amino acid substitution patterns in proteins has led many to conclude that most structural information is carried by interior core residues that are solvent inaccessible. This conclusion is based on the observation that buried residues generally tolerate only conserved sequence changes, while surface residues allow more diverse chemical substitutions. This notion is now changing as it has become apparent that both core and surface residues play important roles in protein folding and stability. Unfortunately, the ability to identify specific mutations that will lead to enhanced stability remains a challenging problem. Here we discuss two mutations that emerged from an in vitro selection experiment designed to improve the folding stability of a non-biological ATP binding protein. These mutations alter two solvent accessible residues, and dramatically enhance the expression, solubility, thermal stability, and ligand binding affinity of the protein. The significance of both mutations was investigated individually and together, and the X-ray crystal structures of the parent sequence and double mutant protein were solved to a resolution limit of 2.8 and 1.65 A, respectively. Comparative structural analysis of the evolved protein to proteins found in nature reveals that our non-biological protein evolved certain structural features shared by many thermophilic proteins. This experimental result suggests that protein fold optimization by in vitro selection offers a viable approach to generating stable variants of many naturally occurring proteins whose structures and functions are otherwise difficult to study.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

The 2.1A crystal structure of copGFP, a representative member of the copepod clade within the green fluorescent protein superfamily

The green fluorescent protein (avGFP), its variants, and the closely related GFP-like proteins are characterized structurally by a cyclic tri-peptide chromophore located centrally within a conserved beta-can fold. Traditionally, these GFP family members have been isolated from the Cnidaria although recently, distantly related GFP-like proteins from the Bilateria, a sister group of the Cnidaria have been described, although no representative structure from this phylum has been reported to date. We have determined to 2.1A resolution the crystal structure of copGFP, a representative GFP-like protein from a copepod, a member of the Bilateria. The structure of copGFP revealed that, despite sharing only 19% sequence identity with GFP, the tri-peptide chromophore (Gly57-Tyr58-Gly59) of copGFP adopted a cis coplanar conformation within the conserved beta-can fold. However, the immediate environment surrounding the chromophore of copGFP was markedly atypical when compared to other members of the GFP-superfamily, with a large network of bulky residues observed to surround the chromophore. Arg87 and Glu222 (GFP numbering 96 and 222), the only two residues conserved between copGFP, GFP and GFP-like proteins are involved in autocatalytic genesis of the chromophore. Accordingly, the copGFP structure provides an alternative platform for the development of a new suite of fluorescent protein tools. Moreover, the structure suggests that the autocatalytic genesis of the chromophore is remarkably tolerant to a high degree of sequence and structural variation within the beta-can fold of the GFP superfamily.     

Denatured state is critical in determining the properties of model proteins designed on different folds

The thermodynamics of proteins designed on three common folds (SH3, chymotrypsin inhibitor 2 [CI2], and protein G) is studied with a simplified C(alpha) model and compared with the thermodynamics of proteins designed on random-generated folds. The model allows to design sequences to fold within a dRMSD ranging from 1.2 to 4.2 A from the crystallographic native conformation and to study properties that are hard to be measured experimentally. It is found that the denatured state of all of them is not random but is, to different extents, partially structured. The degree of structure is more abundant for SH3 and protein G, giving rise to a weaker stability but a more efficient folding kinetics than CI2 and, even more, than the random-generated folds. Consequently, the features of the unfolded state seem to be as important in the determination of the thermodynamic properties of these proteins as the features of the native state.     

A fold-recognition approach to loop modeling

A novel approach is proposed for modeling loop regions in proteins. In this approach, a prerequisite sequence-structure alignment is examined for regions where the target sequence is not covered by the structural template. These regions, extended with a number of residues from adjacent stem regions, are submitted to fold recognition. The alignments produced by fold recognition are integrated into the initial alignment to create an alignment between the target sequence and several structures, where gaps in the main structural template are covered by local structural templates. This one-to-many (1:N) alignment is used to create a protein model by existing protein-modeling techniques. Several alternative approaches were evaluated using a set of ten proteins. One approach was selected and evaluated using another set of 31 proteins. The most promising result was for gap regions not located at the C-terminus or N-terminus of a protein, where the method produced an average RMSD 12% lower than the loop modeling provided with the program MODELLER. This improvement is shown to be statistically significant. Figure The method derived from the training set applied to CASP target T0191     

Folding kinetics of the cooperatively folded subdomain of the IκBα ankyrin repeat domain

The ankyrin repeat (AR) domain of IκBα consists of a cooperative folding unit of roughly four ARs (AR1-AR4) and of two weakly folded repeats (AR5 and AR6). The kinetic folding mechanism of the cooperative subdomain, IκBα_67-206, was analyzed using rapid mixing techniques. Despite its apparent architectural simplicity, IκBα_67-206 displays complex folding kinetics, with two sequential on-pathway high-energy intermediates. The effect of mutations to or away from the consensus sequences of ARs on folding behavior was analyzed, particularly the GXTPLHLA motif, which have not been examined in detail previously. Mutations toward the consensus generally resulted in an increase in folding stability, whereas mutations away from the consensus resulted in decreased overall stability. We determined the free energy change upon mutation for three sequential transition state ensembles along the folding route for 16 mutants. We show that folding initiates with the formation of the interface of the outer helices of AR3 and AR4, and then proceeds to consolidate structure in these repeats. Subsequently, AR1 and AR2 fold in a concerted way in a single kinetic step. We show that this mechanism is robust to the presence of AR5 and AR6 as they do not strongly affect the folding kinetics. Overall, the protein appears to fold on a rather smooth energy landscape, where the folding mechanism conforms a one-dimensional approximation. However, we note that the AR does not necessarily act as a single folding element.     

Structural Determinants for Improved Stability of Designed Ankyrin Repeat Proteins with a Redesigned C-Capping Module

Designed ankyrin repeat proteins (DARPins) that specifically bind to almost any target can be obtained by ribosome display or phage display from combinatorial libraries. Although DARPins are already very stable molecules, molecular dynamics simulations, equilibrium denaturation experiments, structural studies, and recent NMR experiments suggested that the unfolding of the original C-terminal capping repeat (C-cap), taken from a natural ankyrin repeat protein, limits the stability of the initial DARPin design. Several point mutations had been introduced to optimize the C-cap and were shown to indeed further increase the stability of DARPins. We now determined crystal structures of DARPins with one or three full-consensus internal repeats (NI₁C or NI₃C) between an N-terminal capping repeat and mutants of the C-cap. An NI₁C mutant, in which the C-cap was only extended by three additional helix-forming residues, showed no structural change but reduced B-factors in the C-cap. An NI₃C C-cap mutant carrying five additional mutations in the interface to the preceding repeat, previously designed by using the consensus sequence as a guide, showed a rigid-body movement of the C-cap towards the internal repeat. This movement results in an increased buried surface area and a superior surface complementarity and explains the improved stability in equilibrium unfolding, compared to the original C-cap. A C-cap mutant with three additional mutations introducing suitably spaced charged residues did not show formation of salt bridges, explaining why its stability was not increased further. These structural studies underline the importance of repeat coupling for stability and help in the further design of this protein family.     

NMR structure of the human doppel protein

The NMR structure of the recombinant human doppel protein, hDpl(24-152), contains a flexibly disordered "tail" comprising residues 24-51, and a globular domain extending from residues 52 to 149 for which a detailed structure was obtained. The globular domain contains four alpha-helices comprising residues 72-80 (alpha1), 101-115 (alpha2(a)), 117-121 (alpha2(b)), and 127-141 (alpha3), and a short two-stranded anti-parallel beta-sheet comprising residues 58-60 (beta1) and 88-90 (beta2). The fold of the hDpl globular domain thus coincides nearly identically with the structure of the murine Dpl protein. There are close similarities with the human prion protein (hPrP) but, similar to the situation with the corresponding murine proteins, hDpl shows marked local differences when compared to hPrP: the beta-sheet is flipped by 180 degrees with respect to the molecular scaffold formed by the four helices, and the beta1-strand is shifted by two residues toward the C terminus. A large solvent-accessible hydrophobic cleft is formed on the protein surface between beta2 and alpha3, which has no counterpart in hPrP. The helix alpha2 of hPrP is replaced by two shorter helices, alpha2(a) and alpha2(b). The helix alpha3 is shortened by more than two turns when compared with alpha3 of hPrP, which is enforced by the positioning of the second disulfide bond in hDpl. The C-terminal peptide segment 144-149 folds back onto the loop connecting beta2 and alpha2. All but four of the 20 conserved residues in the globular domains of hPrP and hDpl appear to have a structural role in maintaining a PrP-type fold. The conservation of R76, E96, N110 and R134 in hDpl, corresponding to R148, E168, N183 and R208 in hPrP suggests that these amino acid residues might have essential roles in the so far unknown functions of PrP and Dpl in healthy organisms.     

Generation of efficient fingerprint for GFP-like fold and computational identification of potential GFP-like homologs

There is a considerable interest in the detection of GFP-like proteins due to their structural stability and functional usefulness. GFP-like proteins share highly conserved beta-barrel fold with 11 beta-strands. However, their low sequence identity hampers efficient identification of their homologous proteins from database. In this study, an attempt was made to generate a fingerprint for efficient detection of GFP-like proteins. Overlapped conserved residues (OCR)-based approach has been used to design a protein fingerprint based on sequentially and structurally conserved residues in secondary structures to detect homologous proteins very efficiently. Therefore, a fingerprint for GFP-like fold was designed using the OCR approach. However, its specificity was too low to be used for the identification of novel proteins. The conserved residues of loop regions were added and optimized to improve its specificity without losing its high sensitivity. The optimized fingerprint was employed to scan NR database. A total of 20 hypothetical proteins were detected, among which nine were validated as potential GFP-like homologs.     

Structural basis for broad neutralization of hepatitis C virus quasispecies

Monoclonal antibodies directed against hepatitis C virus (HCV) E2 protein can neutralize cell-cultured HCV and pseudoparticles expressing envelopes derived from multiple HCV subtypes. For example, based on antibody blocking experiments and alanine scanning mutagenesis, it was proposed that the AR3B monoclonal antibody recognized a discontinuous conformational epitope comprised of amino acid residues 396–424, 436–447, and 523–540 of HCV E2 envelope protein. Intriguingly, one of these segments (436–447) overlapped with hypervariable region 3 (HVR3), a domain that exhibited significant intrahost and interhost genetic diversity. To reconcile these observations, amino-acid sequence variability was examined and homology-based structural modelling of E2 based on tick-borne encephalitis virus (TBEV) E protein was performed based on 413 HCV sequences derived from 18 subjects with chronic hepatitis C. Here we report that despite a high degree of amino-acid sequence variability, the three-dimensional structure of E2 is remarkably conserved, suggesting broad recognition of structural determinants rather than specific residues. Regions 396–424 and 523–540 were largely exposed and in close spatial proximity at the surface of E2. In contrast, region 436–447, which overlaps with HVR3, was >35 Å away, and estimates of buried surface were inconsistent with HVR3 being part of the AR3B binding interface. High-throughput structural analysis of HCV quasispecies could facilitate the development of novel vaccines that target conserved structural features of HCV envelope and elicit neutralizing antibody responses that are less vulnerable to viral escape.     

Solution structure of ribosomal protein S28E from Methanobacterium thermoautotrophicum

The ribosomal protein S28E from the archaeon Methanobacterium thermoautotrophicum is a component of the 30S ribosomal subunit. Sequence homologs of S28E are found only in archaea and eukaryotes. Here we report the three-dimensional solution structure of S28E by NMR spectroscopy. S28E contains a globular region and a long C-terminal tail protruding from the core. The globular region consists of four antiparallel beta-strands that are arranged in a Greek-key topology. Unique features of S28E include an extended loop L2-3 that folds back onto the protein and a 12-residue charged C-terminal tail with no regular secondary structure and greater flexibility relative to the rest of the protein. The structural and surface resemblance to OB-fold family of proteins and the presence of highly conserved basic residues suggest that S28E may bind to RNA. A broad positively charged surface extending over one side of the beta-barrel and into the flexible C terminus may present a putative binding site for RNA.     

Prediction of protein structure from ideal forms

For many years it has been accepted that the sequence of a protein can specify its three-dimensional structure. However, there has been limited progress in explaining how the sequence dictates its fold and no attempt to do this computationally without the use of specific structural data has ever succeeded for any protein larger than 100 residues. We describe a method that can predict complex folds up to almost 200 residues using only basic principles that do not include any elements of sequence homology. The method does not simulate the folding chain but generates many thousands of models based on an idealized representation of structure. Each rough model is scored and the best are refined. On a set of five proteins, the correct fold score well and when tested on a set of larger proteins, the correct fold was ranked highest for some proteins more than 150 residues, with others being close topological variants. All other methods that approach this level of success rely on the use of templates or fragments of known structures. Our method is unique in using a database of ideal models based on general packing rules that, in spirit, is closer to an ab initio approach.