Complement assembly of two fragments of the streptococcal protein G B1 domain in aqueous solution

GB Tags: Small Covalent Peptide Tags Based on Protein Fragment Reconstitution

Developing peptide tags that can bind target proteins covalently under mild conditions is of great importance for a myriad of applications, ranging from chemical biology to biotechnology. Here we report the development of a small covalent peptide tag system, termed as GB tags, that can covalently label the target protein with high specificity and high yield under oxidizing conditions. The GB tags consist of a pair of short peptides, GN and GC (GN contains 45 residues and GC contains 19 residues). GN and GC, which are split from a parent protein GB1, can undergo protein fragment reconstitution to reconstitute the folded structure of the parent protein spontaneously. The engineered cysteines in GN and GC can readily form a disulfide bond oxidized by air oxygen after protein reconstitution. Using thermally stable variants of GB1, we identified two pairs of GB tags that display improved thermodynamic stability and binding affinity. They can serve as efficient covalent peptide tags for various applications, including specific labeling of mammalian cell surface receptors. We anticipate that these new GB tags will find applications in biochemical labeling as well as biomaterials, such as protein hydrogels.

Life in Phases: Intra- and Inter- Molecular Phase Transitions in Protein Solutions

Proteins, these evolutionarily-edited biological polymers, are able to undergo intramolecular and intermolecular phase transitions. Spontaneous intramolecular phase transitions define the folding of globular proteins, whereas binding-induced, intra- and inter- molecular phase transitions play a crucial role in the functionality of many intrinsically-disordered proteins. On the other hand, intermolecular phase transitions are the behind-the-scenes players in a diverse set of macrosystemic phenomena taking place in protein solutions, such as new phase nucleation in bulk, on the interface, and on the impurities, protein crystallization, protein aggregation, the formation of amyloid fibrils, and intermolecular liquid-liquid or liquid-gel phase transitions associated with the biogenesis of membraneless organelles in the cells. This review is dedicated to the systematic analysis of the phase behavior of protein molecules and their ensembles, and provides a description of the major physical principles governing intramolecular and intermolecular phase transitions in protein solutions.

Engineering protein polymers of ultrahigh molecular weight via supramolecular polymerization: towards mimicking the giant muscle protein titin

Utilizing protein fragment reconstitution of a small protein GB1, we developed an efficient, supramolecular polymerization strategy to engineer protein polymers with ultrahigh molecular weight that mimic the giant muscle protein titin.

The giant muscle protein titin is the largest protein in cells and responsible for the passive elasticity of muscles. Titin, made of hundreds of individually folded globular domains, is a protein polymer with folded globular domains as its macromonomers. Due to titin's ultrahigh molecular weight, it has been challenging to engineer high molecular weight artificial protein polymers that mimic titin. Taking advantage of protein fragment reconstitution (PFR) of a small protein GB1, which can be reconstituted from its two split fragments G_N and G_C, here we report the development of an efficient, PFR-based supramolecular polymerization strategy to engineer protein polymers with ultrahigh molecular weight. We found that the engineered bifunctional protein macromonomers (G_C-macromonomer-G_N) can undergo supramolecular polymerization, in a way similar to condensation polymerization, via the reconstitution of G_N and G_C to produce protein polymers with ultrahigh molecular weight (with an average molecular weight of 0.5 MDa). Such high molecular weight linear protein polymers closely mimic titin and provide protein polymer building blocks for the construction of biomaterials with improved physical and mechanical properties.

Autonomous Sequences in Myoglobin Emerging from X-ray Structure of Holomyoglobin

The proposed continuous folding structure units are fundamental to analyze protein structure. Here, we could elucidate for the first time two types of hydrophobic core networks in apomyoglobin using continuous folding structure units. In myoglobin, two autonomous sequences emerged clearly. We could thus characterize the autonomous sequences using well-defined hydrophobic core networks within respective semifolds. A hydrophobic core is defined as a pair of topology-local hydrophobic amino acids in different folding structures. Hydrophobic core formation is indispensable to stabilize the different folding structures via an efficient hydrophobic interaction. Autonomous sequences in myoglobin encode tertiary structure information for semifolds. These sequences fold autonomously into small sets of continuous folding structure units to grow separate semifolds on each separate framework. The autonomous sequence can be defined as the local sequence assigned to the small set of continuous folding structure units. They create the discrete hydrophobic region in a semifold by assembly of their hydrophobic regions. Semifolds were characterized by discrete hydrophobic regions stabilized by respective type I hydrophobic core networks, which were present within each semifold. The discrete hydrophobic region of a semifold propagated itself with that of a different semifold by hydrophobic interactions in type II hydrophobic core network, which was present between different semifolds, as observed by the X-ray structures of semifolds. The most significant feature of semifolds in apomyoglobin was that they could be verified by the X-ray structure of holomyoglobin regardless of the instability of folds characteristic to autonomous sequence fragments. This work presents the first description of autonomous sequences.

Hydrogels Constructed from Engineered Proteins

Due to their various potential biomedical applications, hydrogels based on engineered proteins have attracted considerable interest. Benefitting from significant progress in recombinant DNA technology and protein engineering/design techniques, the field of protein hydrogels has made amazing progress. The latest progress of hydrogels constructed from engineered recombinant proteins are presented, mainly focused on biorecognition-driven physical hydrogels as well as chemically crosslinked hydrogels. The various bio-recognition based physical crosslinking strategies are discussed, as well as chemical crosslinking chemistries used to engineer protein hydrogels, and protein hydrogels' various biomedical applications. The future perspectives of this fast evolving field of biomaterials are also discussed.

Amino-acid-dependent main-chain torsion-energy terms for protein systems

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters for main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for ψ (N-C(α)-C-N) and ζ (C(β)-C(α)-C-N) angles for each amino acid by using our optimization method, which is one of the knowledge-based approach. In order to test the validity of the new force-field parameters, we then performed folding simulations of α-helical and β-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff03 force field.

PROTEIN FORCE-FIELD PARAMETERS OPTIMIZED WITH THE PROTEIN DATA BANK II: COMPARISONS OF FORCE FIELDS BY FOLDING SIMULATIONS OF SHORT PEPTIDES

In Paper I of this series, the formulations of the optimization method of existing force-field parameters for protein systems have been presented. We then applied it to five sets of force-field parameters, namely, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM version 22, and OPLS-AA. In order to test the validity of these force fields, the folding simulations of α-helical and β-hairpin peptides have been performed with each of the original and optimized force-field parameters. We found that all five modified force-field parameters gave both α-helical and β-hairpin structures more consistent with the experimental implications than the original force fields.

pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK(a) values. However, pK(a) values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK(a) values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK(a) values were determined for carboxyl groups by monitoring their pH-dependent (13)C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK(a) values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK(a) values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK(a) values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.

A novel force field parameter optimization method based on LSSVR for ECEPP

In this paper, we propose a novel force field parameter optimization method based on LSSVR and optimize the torsion energy parameters of ECEPP force field. In this method force field parameter optimization problem is turned into a support vector regression problem. Protein samples for regression model training are chosen from Protein Data Bank. The experiments show that the optimized force-field parameters make both α-helix and β-hairpin structures more consistent with the experimental implications than the original parameters.

In vivo protein stabilization based on fragment complementation and a split GFP system

Protein stabilization was achieved through in vivo screening based on the thermodynamic linkage between protein folding and fragment complementation. The split GFP system was found suitable to derive protein variants with enhanced stability due to the correlation between effects of mutations on the stability of the intact chain and the effects of the same mutations on the affinity between fragments of the chain. PGB1 mutants with higher affinity between fragments 1 to 40 and 41 to 56 were obtained by in vivo screening of a library of the 1 to 40 fragments against wild-type 41 to 56 fragments. Colonies were ranked based on the intensity of green fluorescence emerging from assembly and folding of the fused GFP fragments. The DNA from the brightest fluorescent colonies was sequenced, and intact mutant PGB1s corresponding to the top three sequences were expressed, purified, and analyzed for stability toward thermal denaturation. The protein sequence derived from the top fluorescent colony was found to yield a 12 °C increase in the thermal denaturation midpoint and a free energy of stabilization of -8.7 kJ/mol at 25 °C. The stability rank order of the three mutant proteins follows the fluorescence rank order in the split GFP system. The variants are stabilized through increased hydrophobic effect, which raises the free energy of the unfolded more than the folded state; as well as substitutions, which lower the free energy of the folded more than the unfolded state; optimized van der Waals interactions; helix stabilization; improved hydrogen bonding network; and reduced electrostatic repulsion in the folded state.

Protein GB1 folding and assembly from structural elements

Folding of the Protein G B1 domain (PGB1) shifts with increasing salt concentration from a cooperative assembly of inherently unstructured subdomains to an assembly of partly pre-folded structures. The salt-dependence of pre-folding contributes to the stability minimum observed at physiological salt conditions. Our conclusions are based on a study in which the reconstitution of PGB1 from two fragments was studied as a function of salt concentrations and temperature using circular dichroism spectroscopy. Salt was found to induce an increase in β-hairpin structure for the C-terminal fragment (residues 41 – 56), whereas no major salt effect on structure was observed for the isolated N-terminal fragment (residues 1 – 41). In line with the increasing evidence on the interrelation between fragment complementation and stability of the corresponding intact protein, we also find that salt effects on reconstitution can be predicted from salt dependence of the stability of the intact protein. Our data show that our variant (which has the mutations T2Q, N8D, N37D and reconstitutes in a manner similar to the wild type) displays the lowest equilibrium association constant around physiological salt concentration, with higher affinity observed both at lower and higher salt concentration. This corroborates the salt effects on the stability towards denaturation of the intact protein, for which the stability at physiological salt is lower compared to both lower and higher salt concentrations. Hence we conclude that reconstitution reports on molecular factors that govern the native states of proteins.

Conformational studies of the alpha-helical 28-43 fragment of the B3 domain of the immunoglobulin binding protein G from Streptococcus

To determine whether the alpha-helix in the B3 immunoglobulin binding domain of protein G from group G Streptococcus has conformational stability as an isolated fragment, we carried out a CD and NMR study of the 16-residue peptide in solution corresponding to this alpha-helix. Based on two-dimensional H-NMR spectra recorded at three different temperatures (283, 305, and 313 K), it was found that this peptide is mostly unstructured in water at these temperatures. Weak signals corresponding to i,i+3 or i,i+4 interactions, which are characteristic of formation of turn-like structures, were observed in the ROE spectra at all temperatures. The absence of a stable three-dimensional structure of the investigated peptide supports an earlier study (Blanco and Serrano, Eur J Biochem 1995, 230, 634-649) of a possible mechanism for folding of other (B1 and B2) immunoglobulin binding domains of Protein G.

Protein reconstitution and three-dimensional domain swapping: benefits and constraints of covalency

The phenomena of protein reconstitution and three-dimensional domain swapping reveal that highly similar structures can be obtained whether a protein is comprised of one or more polypeptide chains. In this review, we use protein reconstitution as a lens through which to examine the range of protein tolerance to chain interruptions and the roles of the primary structure in related features of protein structure and folding, including circular permutation, natively unfolded proteins, allostery, and amyloid fibril formation. The results imply that noncovalent interactions in a protein are sufficient to specify its structure under the constraints imposed by the covalent backbone.

Minimization and optimization of designed beta-hairpin folds

Minimized beta hairpins have provided additional data on the geometric preferences of Trp interactions in TW-loop-WT motifs. This motif imparts significant fold stability to peptides as short as 8 residues. High-resolution NMR structures of a 16- (KKWTWNPATGKWTWQE, DeltaG(U)(298) >or= +7 kJ/mol) and 12-residue (KTWNPATGKWTE, DeltaG(U)(298) = +5.05 kJ/mol) hairpin reveal a common turn geometry and edge-to-face (EtF) packing motif and a cation-pi interaction between Lys(1) and the Trp residue nearest the C-terminus. The magnitude of a CD exciton couplet (due to the two Trp residues) and the chemical shifts of a Trp Hepsilon3 site (shifted upfield by 2.4 ppm due to the EtF stacking geometry) provided near-identical measures of folding. CD melts of representative peptides with the -TW-loop-WT- motif provided the thermodynamic parameters for folding, which reflect enthalpically driven folding at laboratory temperatures with a small DeltaC(p) for unfolding (+420 J K(-)(1)/mol). In the case of Asx-Pro-Xaa-Thr-Gly-Xaa loops, mutations established that the two most important residues in this class of direction-reversing loops are Asx and Gly: mutation to alanine is destabilizing by about 6 and 2 kJ/mol, respectively. All indicators of structuring are retained in a minimized 8-residue construct (Ac-WNPATGKW-NH(2)) with the fold stability reduced to DeltaG(U)(278) = -0.7 kJ/mol. NMR and CD comparisons indicate that -TWXNGKWT- (X = S, I) sequences also form the same hairpin-stabilizing W/W interaction.

Cooperative folding mechanism of a β-hairpin peptide studied by a multicanonical replica-exchange molecular dynamics simulation

G-peptide is a 16-residue peptide of the C-terminal end of streptococcal protein G B1 domain, which is known to fold into a specific beta-hairpin within 6 micros. Here, we study molecular mechanism on the stability and folding of G-peptide by performing a multicanonical replica-exchange (MUCAREM) molecular dynamics simulation with explicit solvent. Unlike the preceding simulations of the same peptide, the simulation was started from an unfolded conformation without any experimental information on the native conformation. In the 278-ns trajectory, we observed three independent folding events. Thus MUCAREM can be estimated to accelerate the folding reaction more than 60 times than the conventional molecular dynamics simulations. The free-energy landscape of the peptide at room temperature shows that there are three essential subevents in the folding pathway to construct the native-like beta-hairpin conformation: (i) a hydrophobic collapse of the peptide occurs with the side-chain contacts between Tyr45 and Phe52, (ii) then, the native-like turn is formed accompanying with the hydrogen-bonded network around the turn region, and (iii) finally, the rest of the backbone hydrogen bonds are formed. A number of stable native hydrogen bonds are formed cooperatively during the second stage, suggesting the importance of the formation of the specific turn structure. This is also supported by the accumulation of the nonnative conformations only with the hydrophobic cluster around Tyr45 and Phe52. These simulation results are consistent with high phi-values of the turn region observed by experiment.

Intra- versus intermolecular interactions in monellin: contribution of surface charges to protein assembly

The relative significance of weak non-covalent interactions in biological context has been much debated. Here, we have addressed the contribution of Coulombic interactions to protein stability and assembly experimentally. The sweet protein monellin, a non-covalently linked heterodimeric protein, was chosen for this study because of its ability to spontaneously reconstitute from separated fragments. The reconstitution of monellin mutants containing large surface charge perturbations was compared to the thermostability of structurally equivalent single-chain monellin containing the same sets of mutations under varying salt concentrations. The affinity between monellin fragments is found to correlate with the thermostability of single chain monellin, indicating the involvement of the same underlying Coulombic interactions. This confirms that there are no principal differences in the interactions involved in folding and binding. Based on comparison with a previous mutational study involving hydrophobic core residues, the relative contribution of Coulombic interactions to stability and affinity is modest. However, the Coulombic perturbations only affect the association rates of reconstitution in contrast to perturbations involving hydrophobic residues, which affect primarily the dissociation rates. These results indicate that Coulombic interactions are likely to be of main importance for the association of protein assembly, relevant for functions of proteins.

Reconstitution of calmodulin from domains and subdomains: Influence of target peptide

Reconstitution studies of a protein from domain fragments can furnish important insights into the distinctive role of particular domain interactions and how they affect biophysical properties important for function. Using isothermal titration calorimetry (ITC) and a number of spectroscopic and chromatographic tools, including CD, fluorescence and NMR spectroscopy, size-exclusion chromatography and non-denaturing agarose gel electrophoresis, we have investigated the reconstitution of the ubiquitous Ca2+-sensor protein calmodulin (CaM) and its globular domains from fragments comprising one or two EF-hands. The studies were carried out with and without the target peptide from smooth muscle myosin light chain kinase (smMLCKp). The CaM-target complex can be reconstituted from the three components consisting of the target peptide and the globular domains TR1C and TR2C. In the absence of peptide, there is no evidence for association of the globular domains. The globular domains can further be reconstituted from their corresponding native subdomains. The dissociation constant, K(D), in 2 mM Tris-HCl (pH 7.5), for the subdomain complexes, EF1:EF2 and EF3:EF4, was determined with ITC to 9.3 x 10(-7) M and 5.9 x 10(-8) M, respectively. Thus, the affinity between the two C-terminal subdomains, located within TR2C, is stronger by a factor of 16 than that between the corresponding subdomains within TR1C. These observations are corroborated by the spectroscopic and chromatographic investigations.

Multi-method global analysis of thermodynamics and kinetics in reconstitution of monellin

Accurate and precise determinations of thermodynamic parameters of binding are important steps toward understanding many biological mechanisms. Here, a multi-method approach to binding analysis is applied and a detailed error analysis is introduced. Using this approach, the binding thermodynamics and kinetics of the reconstitution of the protein monellin have been quantitatively determined in detail by simultaneous analysis of data collected with fluorescence spectroscopy, surface plasmon resonance and isothermal titration calorimetry at 25 degrees C, pH 7.0 and 150 mM NaCl. Monellin is an intensely sweet protein composed of two peptide chains that form a single globular domain. The kinetics of the reconstitution reaction are slow, with an association rate constant, k(on) of 8.8 x 10(3) M(-1) s(-1) and a dissociation rate constant, k(off) of 3.1 x 10(-4) s(-1). The equilibrium constant K(A) is 2.8 x 10(7) M(-1) corresponding to a standard free energy of association, DeltaG degrees , of -42.5 kJ/mol. The enthalpic component, DeltaH degrees , is -18.7 kJ/mol and the entropic contribution, DeltaS degrees , is 79.8 J mol(-1) K(-1) (-TDeltaS degrees = -23.8 kJ/mol). The association of monellin is therefore a bimolecular intra-protein association whose energetics are slightly dominated by entropic factors.

Hairpin folding rates reflect mutations within and remote from the turn region

Hairpins play a central role in numerous protein folding and misfolding scenarios. Prior studies of hairpin folding, many conducted with analogs of a sequence from the B1 domain of protein G, suggest that faster folding can be achieved only by optimizing the turn propensity of the reversing loop. Based on studies using dynamic NMR, the native GB1 sequence is a slow folding hairpin (k(F)(278)=1.5 x 10(4)/s). GB1 hairpin analogs spanning a wide range of thermodynamic stabilities (DeltaG(U)(298)=-3.09 to+3.25 kJ/mol) were examined. Fold-stabilizing changes in the reversing loop can act either by accelerating folding or retarding unfolding; we present examples of both types. The introduction of an attractive side-chain/side-chain Coulombic interaction at the chain termini further stabilizes this hairpin. The 1.9-fold increase in folding rate constant observed for this change at the chain termini implies that this Coulombic interaction contributes before or at the transition state. This observation is difficult to rationalize by "zipper" folding pathways that require native turn formation as the sole nucleating event; it also suggests that Coulombic interactions should be considered in the design of systems intended to probe the protein folding speed limit.

High-affinity fragment complementation of a fibronectin type III domain and its application to stability enhancement

The tenth fibronectin type III (FN3) domain of human fibronectin (FNfn10), a prototype of the ubiquitous FN3 domain, is a small, monomeric beta-sandwich protein. In this study, we have bisected FNfn10 in each loop to generate a total of six fragment pairs. We found that fragment pairs bisected at multiple loops of FNfn10 show complementation in vivo as tested with a yeast two-hybrid system. The dissociation constant of these fragment pairs determined in vitro were as low as 3 nM, resulting in one of the tightest fragment complementation systems reported so far. Furthermore, we show that the affinity of fragment complementation is correlated with the stability of the uncut parent protein. Exploring this correlation, we screened a yeast two-hybrid library of one fragment and identified mutations that suppress the effect of a destabilizing mutation in the other fragment. One of the identified mutations significantly increased the stability of the uncut wild-type protein, proving that fragment complementation can be used as a novel strategy for the selection of proteins with enhanced stability.

Probing Fragment Complementation by Rigid-Body Docking: in Silico Reconstitution of Calbindin D9k

Fragment complementation is gaining an increasing impact as a nonperturbing method to probe noncovalent interactions within protein supersecondary structures. In this study, the fast Fourier transform rigid-body docking algorithm ZDOCK has been employed for in silico reconstitution of the calcium binding protein calbindin D9k, from its two EF-hands subdomains, namely, EF1 (residues 1-43) and EF2 (residues 44-75). The EF1 fragment has been used both in its wild type and in nine mutant forms, in line with in vitro experiments. Consistent with in vitro data, ZDOCK reconstituted the proper fold of wild-type and mutated calbindin, locating the nativelike structures (i.e., holding a root-mean-square deviation < 1 A with respect to the X-ray structure) among the first 10 top-scored solutions out of 4000. Moreover, the three independent in silico reconstitutions of wild-type calbindin ranked a nativelike structure at the top of the output list, that is, the best scored one. The algorithm has been also successfully challenged in reconstituting the EF2 homodimer from two identical copies of the monomer. Furthermore, quantitative models consisting of linear correlations between thermodynamic data and ZDOCK scores were built, providing a tested tool for very fast in silico predictions of the free energy of association of protein-protein complexes solved at the atomic level and known to not undergo significant conformational changes upon binding.

10 Residue Folded Peptide Designed by Segment Statistics

We have designed a peptide termed chignolin, consisting of only 10 amino acid residues (GYDPETGTWG), on the basis of statistics derived from more than 10,000 protein segments. The peptide folds into a unique structure in water and shows a cooperative thermal transition, both of which may be hallmarks of a protein. Also, the experimentally determined beta-hairpin structure was very close to what we had targeted. The performance of the short peptide not only implies that the methodology employed here can contribute toward development of novel techniques for protein design, but it also yields insights into the raison d'etre of an autonomous element involved in a natural protein. This is of interest for the pursuit of folding mechanisms and evolutionary processes of proteins.

A Captured Folding Intermediate Involved in Dimerization and Domain-swapping of GB1

Immunoglobulin binding domain B1 of streptococcal protein G (GB1), a small (56 residues), stable, single domain protein, is one of the most extensively used model systems in the area of protein folding and design. The recently determined NMR structure of a quadruple mutant (HS#124F26A, L5V/F30V/Y33F/A34F) revealed a domain-swapped dimer that dissociated into a partially folded, monomeric species at low micromolar protein concentrations. Here, we have characterized this monomeric, partially folded species by NMR and show that extensive conformational heterogeneity for a substantial portion of the polypeptide chain exists. Exchange between the conformers within the monomer ensemble on the microsecond to millisecond timescale renders the majority of backbone amide resonances broadened beyond detection. Despite these extensive temporal and spatial fluctuations, the overall architecture of the monomeric mutant protein resembles that of wild-type GB1 and not the monomer unit of the domain-swapped dimer.

Heat capacities and a snapshot of the energy landscape in protein GB1 from the pre-denaturation temperature dependence of backbone NH nanosecond fluctuations

Protein stability is usually characterized calorimetrically by a melting temperature and related thermodynamic parameters. Despite its importance, the microscopic origin of the melting transition and the relationship between thermodynamic stability and dynamics remains a mystery. Here, NMR relaxation parameters were acquired for backbone 15NH groups of the 56 residue immunoglobulin-binding domain of streptococcal protein G over a pre-denaturation temperature range of 5-50 degrees C. Relaxation data were analyzed using three methods: the standard three-Lorentzian model free approach; the F(omega)=2omegaJ(omega) spectral density approach that yields motional correlation time distributions, and a new approach that determines frequency-dependent order parameters. Regardless of the method of analysis, the temperature dependence of internal motional correlation times and order parameters is essentially the same. Nanosecond time-scale internal motions are found for all NHs in the protein, and their temperature dependence yields activation energies ranging up to about 33kJ/mol residue. NH motional barrier heights are structurally correlated, with the largest energy barriers being found for residues in the most "rigid" segments of the fold: beta-strands 1 and 4 and the alpha-helix. Trends in this landscape also parallel the free energy of folding-unfolding derived from hydrogen-deuterium (H-D) exchange measurements, indicating that the energetics for internal motions occurring on the nanosecond time-scale mirror those occurring on the much slower time-scale of H-D exchange. Residual heat capacities, derived from the temperature dependence of order parameters, range from near zero to near 100J/mol K residue and correlate with this energy landscape. These results provide a unique picture of this protein's energy landscape and a relationship between thermodynamic stability and dynamics that suggests thermosensitive regions in the fold that could initiate the melting process.

Reconstitution of a native-like SH2 domain from disordered peptide fragments examined by multidimensional heteronuclear NMR

The N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3' kinase is cleaved specifically into 9- and 5-kD fragments by limited proteolytic digestion with trypsin. The noncovalent SH2 domain complex and its constituent tryptic peptides have been investigated using high-resolution heteronuclear magnetic resonance (NMR). These studies have established the viability of the SH2 domain as a fragment complementation system. The individual peptide fragments are predominantly unstructured in solution. In contrast, the noncovalent 9-kD + 5-kD complex shows a native-like (1)H-(15)N HSQC spectrum, demonstrating that the two fragments fold into a native-like structure on binding. Chemical shift analysis of the noncovalent complex compared to the native SH2 domain reveals that the highest degree of perturbation in the structure occurs at the cleavage site within a flexible loop and along the hydrophobic interface between the two peptide fragments. Mapping of these chemical shift changes on the structure of the domain reveals changes consistent with the reduction in affinity for the target peptide ligand observed in the noncovalent complex relative to the intact protein. The 5-kD fragment of the homologous Src protein is incapable of structurally complementing the p85 9-kD fragment, either in complex formation or in the context of the full-length protein. These high-resolution structural studies of the SH2 domain fragment complementation features establish the suitability of the system for further protein-folding and design studies.

Probing the structure and stability of a hybrid protein: the human-E. coli thioredoxin chimera

The structure and stability of a hybrid protein composed of N-terminal human and C-terminal E. coli thioredoxin domains were investigated by NMR, fluorescence, and circular dichroism spectroscopy. We demonstrate that the chimeric protein is correctly folded and exhibits the common thioredoxin architecture. However, the stability of the hybrid protein toward thermal and chemical denaturation is clearly reduced when compared with both parent proteins. Altogether, our data indicate that the interface between the two folding units of thioredoxin is tolerant toward changes in exact interdigitation of side chains, allowing for the formation of the unique overall thioredoxin fold. Further, the gene encoding the human-E. coli chimera was tested in vivo whether it supports the assembly of filamentous phages. No complementation of a thioredoxin-deficient E. coli mutant for the replication of the phages M13 or fd was observed, suggesting that parts of the overall protein structure in the N-terminal domain are crucial for this activity.

Cooperative assembly of a nativelike ubiquitin structure through peptide fragment complexation: energetics of peptide association and folding

Peptide fragments corresponding to the N- and C-terminal portions of bovine ubiquitin, U(1-35) and U(36-76), are shown by NMR to associate in solution to form a complex of modest stability (Kassn approximately 1.4 x 10(5) M(-1) at pH 7.0), with NMR features characteristic of a nativelike structure. The complex undergoes cold denaturation, with temperature-dependent estimates of stability from NMR indicating a DeltaC(p) degrees for fragment complexation in good agreement with that determined for native ubiquitin, suggesting that fragment association results in the burial of a similar hydrophobic surface area. The stability of the complex shows appreciable pH dependence, suggesting that ionic interactions on the surface of the protein contribute significantly. However, denaturation studies of native ubiquitin in the presence of guanidine hydrochloride (Gdn.HCl) show little pH dependence, suggesting that ionic interactions may be "screened" by the denaturant, as recently suggested. Examination of the conformation of the isolated peptide fragments has shown evidence for a low population of nativelike structure in the N-terminal beta-hairpin (residues 1-17) and weak nascent helical propensity in the helical fragment (residues 21-35). In contrast, the C-terminal peptide (36-76) shows evidence in aqueous solution, from some Halpha chemical shifts, for nonnative phi and psi angles; nonnative alpha-helical structure is readily induced in the presence of organic cosolvents, indicating that tertiary interactions in both native ubiquitin and the folded fragment complex strongly dictate its structural preference. The data suggest that the N-terminal fragment (1-35), where interaction between the helix and hairpin requires the minimum loss of conformational entropy, may provide the nucleation site for fragment complexation.

Nanosecond dynamics of the single tryptophan reveals multi-state equilibrium unfolding of protein GB1

Unfolding of the immunoglobulin binding domain B1 of streptococcal protein G (GB1) was induced by guanidine hydrochloride (GdnHCl) and studied by circular dichroism, steady-state, and time-resolved fluorescence spectroscopy. The fluorescence methods employed the single tryptophan residue of GB1 as an intrinsic reporter. While the transitions monitored by circular dichroism and steady-state fluorescence coincided with each other, the transitions followed by dynamic fluorescence were markedly different. Specifically, fluorescence anisotropy data showed that a relaxation spectrum of tryptophan contained a slow motion with relaxation times of 9 ns in the native state and 4 ns in the unfolded state in 6 M GdnHCl. At intermediate GdnHCl concentrations of 3.8-4.2 M, however, the slow relaxation time increased to 18 ns. The fast nanosecond motion had an average time of 0.8 ns and showed no dependence on the formation of native structure. Overall, dynamic fluorescence revealed two preliminary stages in GB1 folding, which are equated with the formation of local structure in the beta(3)-strand hairpin and the initial collapse. Both stages exist without alpha-helix formation, i. e., before the appearance of any ordered secondary structure detectable by circular dichroism. Another stage in GB1 folding might exist at very low ( approximately 1 M) GdnHCl concentrations.

The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (eta(xy)), and steady state [1H]-15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 degrees C for 89-100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/eta(xy) was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to approximately 0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be approximately 9 degrees C lower (78 degrees C rather than 87 degrees C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.

Role of side-chains in the cooperative beta-hairpin folding of the short C-terminal fragment derived from streptococcal protein G

A short C-terminal fragment of immunoglobulin-binding domain of streptococcal protein G is known to form nativelike beta-hairpin at physiological conditions. To understand the cooperative folding of the short peptide, eight Ala-substituted mutants of the fragment were investigated with respect to their structural stabilities by analyzing temperature dependence of NMR signals. On comparison of the obtained thermodynamic parameters, we found that the nonpolar residues Tyr45 and Phe52 and the polar residues Asp46 and Thr49 are crucial for the beta-hairpin folding. The results suggest a strong interaction between the nonpolar side chains that participates in a putative hydrophobic cluster and that the polar side chains form a fairly rigid conformation around the loop (46-51). We also investigated the complex formation of the mutants with N-terminal fragment at the variety of temperature to get their thermal unfolding profiles and found that the mutations on the residues Asp46 and Thr49 largely destabilized the complexes, while substitution of Asp47 slightly stabilized the complex. From these results, we deduced that both the hydrophobic cluster formation and the rigidity of the loop (46-51) cooperatively stabilize the beta-hairpin structure of the fragment. These interactions which form a stable beta-hairpin may be the initial structural scaffold which is important in the early folding events of the whole domain.

Thermodynamics of a beta-hairpin structure: evidence for cooperative formation of folding nucleus

To elucidate early nucleation stages in protein folding, multi-probed thermodynamic characterization was applied to the beta-hairpin structural formation of G-peptide, which is a C-terminal fragment of the B1 domain of streptococcal protein G. The segment corresponding to the sequence of G-peptide is believed to act as a nucleus during the folding process of the B1 domain. In spite of the broad thermal transition of G-peptide, nuclear magnetic resonance (NMR) melting measurements combined with our original analytical theory enabled us to obtain the thermodynamic properties of the beta-hairpin formation with considerable accuracy. Additionally, all the thermodynamic properties determined by every NMR probe on both the main-chain and the side-chains were quite similar, and also comparable to the values that were independently determined by calorimetric analysis of G-peptide. These results demonstrate that G-peptide folds cooperatively throughout the molecule. In other words, the formation of the beta-hairpin is interpreted as the fashion of a first-order phase transition between two states without any distinguishable intermediates. This cooperative formation of the short linear peptide consisting of only 16 residues provides insight into not only the first folding events of the B1 domain, but also the general principles of proteins in terms of structural hierarchy, stability and folding mechanism.

Unusual stability of a multiply nicked form of Plasmodium falciparum triosephosphate isomerase

BACKGROUND

The limited proteolytic cleavage of proteins can result in distinct polypeptides that remain noncovalently associated so that the structural and biochemical properties of the 'nicked' protein are virtually indistinguishable from those of the native protein. The remarkable observation that rabbit muscle triosephosphate isomerase (TIM) can be multiply nicked by subtilisin and efficiently religated in the presence of an organic solvent formed the stimulus for our study on a homologous system, Plasmodium falciparum triosephosphate isomerase (PfTIM).

RESULTS

The subtilisin nicked form of PfTIM was prepared by limited proteolysis using subtilisin and the major fragments identified using electrospray ionization mass spectrometry. The order of susceptibility of the peptide bonds in the protein was also determined. The structure of the nicked form of TIM was investigated using circular dichroism, fluorescence and gel filtration. The nicked enzyme exhibited remarkable stability to denaturants, although significant differences were observed with the wild-type enzyme. Efficient religation could be achieved by addition of an organic cosolvent, such as acetonitrile, in the presence of subtilisin. Religation was also demonstrated after dissociation of the proteolytic fragments in guanidinium chloride, followed by reassembly after removal of the denaturant.

CONCLUSIONS

The eight-stranded beta8/alpha8 barrel is a robust, widely used protein structural motif. This study demonstrates that the TIM barrel can withstand several nicks in the polypeptide backbone with a limited effect on its structure and stability.

Fragment reconstitution of a small protein: disulfide mutant of a short C-terminal fragment derived from streptococcal protein G

Hierarchical studies on the folding of protein G B1 domain have shown that the C-terminal fragment (C16) has a considerable amount of beta-hairpin structure that exchanges between the folded and unfolded states at room temperature, and that the C16 fragment binds noncovalently to an N-terminal fragment (N40) under physiological conditions. Those studies have led us to the hypothesis that the amphipathic beta-hairpin structure of C16 initiates folding of the domain. To obtain a more detailed understanding of the folding mechanism of the domain, we designed a mutant of C16 (SS16ox) with a disulfide bond between residues 41 and 56, and then examined the interaction of the mutant with N40 by surface plasmon resonance (SPR) and by thermal denaturation studies using circular dichroism. SS16ox strongly interacted with N40, with an equilibrium constant, KD, that was 7-fold higher than wild-type. The association rate constant, kon, of SS16ox was 8.7-fold higher than that of wild-type. This strong interaction can be explained by the entropic effect of the disulfide bond. The introduction of the disulfide bond into C16 stabilizes the beta-hairpin structure of C16, accelerates the association rate with N40, and then stabilizes the whole complex. These results support a hypothetical folding mechanism of protein G where the amphipathic beta-hairpin structure of C16 acts as a nucleus and accelerates folding of the whole molecule.

Fragment reconstitution of a small protein: folding energetics of the reconstituted immunoglobulin binding domain B1 of streptococcal protein G

To elucidate early stages in protein folding, we have adopted a fragment reconstitution method for small proteins. This approach is expected to provide nuclei for protein folding and to allow us to investigate folding mechanisms. In previous work [Kobayashi, N., et al. (1995) FEBS Lett. 366, 99-103.] we demonstrated the association of two complementary fragments, derived from the immunoglobulin G-binding domain B1 of streptococcal Protein G, and showed the structural similarity between the reconstituted domain and the uncleaved wild-type domain. In this work we have further characterized the reconstituted domain as well as the uncleaved domain thermodynamically by means of differential scanning calorimetry (DSC) and circular dichroism (CD) measurements. Although composed of short peptide fragments not linked by covalent bonds, the reconstituted domain showed a typical folding/unfolding curve in both DSC and CD melting measurements and behaved like a globular protein. The domain was not very stable, and the small value of the Gibbs free energy corresponded to the class of the weakest protein-protein binding systems. The denaturation temperature of 0. 78 mM solution was 313 K at pH 5.9 as measured by DSC, which was more than 40 degrees lower than the uncleaved domain. This apparent instability was primarily caused by entropic disadvantage attributed to a bimolecular reaction. The temperature dependence of the enthalpy change from the folded to the unfolded state was almost identical for the reconstituted domain and the uncleaved one. This indicates that most of the noncovalent intramolecular interactions stabilizing the native structure, such as hydrogen bonding and hydrophobic interactions, are regenerated in the reconstituted domain. By comparing the equilibrium constants of the reconstituted and uncleaved domains, we determined the effective concentration to be approximately 6 M at 298 K. Structure-based estimation of the thermodynamic properties from the values of accessible surface areas showed that approximately 35% of the total heat capacity change and approximately 25% of the total enthalpy change can be attributed to the interchain interaction at 298 K. Furthermore, the folding/unfolding equilibrium of beta-hairpin structure of the fragment 41-56 alone was also characterized. These analyses allow us to envision the microdomain folding mechanism of the Protein G B1 domain, in which segment 41-56 first forms a stable beta-hairpin structure and then collides with segment 1-40, followed by spontaneous folding of the whole molecule.

Conformational analysis of peptide fragments derived from the peripheral subunit-binding domain from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus: evidence for nonrandom structure in the unfolded state

There is currently a great deal of interest in the early events in protein folding. Two issues that have generated particular interest are the nature of the unfolded state under native conditions and the role of local interactions in folding. Here, we report the results of a study of a set of peptides derived from a small two-helix protein, the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex. Five peptides of overlapping sequence were prepared, including sequences corresponding to each of the helices and to the region connecting them. The peptides were characterized by CD and, where possible, nmr. A peptide corresponding to the second helix is between 12 and 17% helical at neutral pH. CD also indicates a lower percentage of helical structure in the peptide corresponding to the first alpha-helix, although the values of the alpha-proton chemical shifts suggest some preference for nonrandom structure. Peptides corresponding to the interhelical loop, which in the full domain contains two overlapping beta-turns and a 5-residue 3(10)-helix, are less structured. There is no significant change in the helicity of any of these peptides with pH. To test for fragment complementation, CD spectra of the two peptides derived from each helix and the long connecting peptide were compared to the spectra of each possible pair, as well as to a mixture containing all three. No increase in structure was observed. We complement our peptide studies by characterizing a point mutant, D34V, which disrupts a critical hydrogen bonding network. This mutant is unable to fold and provides a useful model of the denatured state. The mutant is between 9 and 16% helical as judged by CD. The modest amount of helical structure formed in some of the peptide fragments and in the point mutant suggests that the denatured state of the peripheral subunit binding domain is not completely unstructured. This may contribute to the very rapid folding observed for the intact protein.