Small angle x-ray scattering of the hemoglobin from Biomphalaria glabrata

The hemoglobin from Biomphalaria glabrata is an extracellular respiratory protein of high molecular mass composed by subunits of 360 kDa, each one containing two 180 kDa chains linked by disulfide bridges. In this work, small angle x-ray scattering (SAXS) measurements were performed with the hemoglobin at pH 5.0 and 7.5. Radii of gyration of 98.6 +/- 0.5 and 101.8 +/- 0.2 A and maximum diameters of 300 +/- 10 and 305 +/- 10 A, respectively, were obtained from Guinier plot extrapolation and analytical curve fitting. The pair distance distribution functions p(r) corresponded to globular particles with a somewhat anisotropic shape for both preparations. Computer analysis of the low angle part of the scattering curve led to the determination of the low resolution envelope of the protein, revealing a P(222) symmetry. Shape reconstruction from ab initio calculations using the complete scattering curve furnished a compact prolate three-dimensional (3D) bead model for the protein. Hydrodynamic parameters were obtained from experiments and theoretical calculations using the 3D model. The results of the structural and biochemical studies reported herein indicate that the multisubunit structure of this hemoglobin is compatible with a tetrameric arrangement.     

Relationship between chain collapse and secondary structure formation in a partially folded protein

Chain collapse and secondary structure formation are frequently observed during the early stages of protein folding. Is the chain collapse brought about by interactions between secondary structure units or is it due to polymer behavior in a poor solvent (coil‐globule transition)? To answer this question, we measured small‐angle X‐ray scattering for a series of β‐lactoglobulin mutants under conditions in which they assume a partially folded state analogous to the folding intermediates. Mutants that were designed to disrupt the secondary structure units showed the gyration radii similar to that of the wild type protein, indicating that chain collapse is due to coil‐globule transitions. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 651–658, 2014.     

Simulation and experiment conspire to reveal cryptic intermediates and a slide from the nucleation-condensation to framework mechanism of folding

There is a change from three-state to two-state kinetics of folding across the homeodomain superfamily of proteins as the mechanism slides from framework to nucleation-condensation. The tendency for framework folding in this family correlates with inherent helical propensity. The cellular myeloblastis protein (c-Myb) falls in the mechanistic transition region. An earlier, preliminary report of protein engineering experiments and molecular dynamics simulations (MD) showed that the folding mechanism for this protein has aspects of both the nucleation-condensation and framework models. In the more in-depth analysis of the MD trajectories presented here, we find that folding may be attributed to both of these mechanisms in different regions of the protein. The folding of the loop, middle helix, and turn is best described by nucleation-condensation, whereas folding of the N and C-terminal helices may be described by the framework model. Experimentally, c-Myb folds by apparent two-state kinetics, but the MD simulations predict that the kinetics hide a high-energy intermediate. We stabilized this hypothetical folding intermediate by deleting a residue (P174) in the loop between its second and third helices, and the mutant intermediate is long-lived in the simulations. Equilibrium and kinetic experiments demonstrate that folding of the DeltaP174 mutant is indeed three-state. The presence and shape of the intermediate observed in the simulations were confirmed by small angle X-ray scattering experiments.     

Domain interactions determine the conformational ensemble of the periplasmic chaperone SurA

SurA is thought to be the most important periplasmic chaperone for outer membrane protein (OMP) biogenesis. Its structure is composed of a core region and two peptidylprolyl isomerase domains, termed P1 and P2, connected by flexible linkers. As such these three independent folding units are able to adopt a number of distinct spatial positions with respect to each other. The conformational dynamics of these domains are thought to be functionally important yet are largely unresolved. Here we address this question of the conformational ensemble using sedimentation equilibrium, small‐angle neutron scattering, and folding titrations. This combination of orthogonal methods converges on a SurA population that is monomeric at physiological concentrations. The conformation that dominates this population has the P1 and core domains docked to one another, for example, “P1‐closed” and the P2 domain extended in solution. We discovered that the distribution of domain orientations is defined by modest and favorable interactions between the core domain and either the P1 or the P2 domains. These two peptidylprolyl domains compete with each other for core‐binding but are thermodynamically uncoupled. This arrangement implies two novel insights. Firstly, an open conformation must exist to facilitate P1 and P2 exchange on the core, indicating that the open client‐binding conformation is populated at low levels even in the absence of client unfolded OMPs. Secondly, competition between P1 and P2 binding paradoxically occludes the client binding site on the core, which may serve to preserve the reservoir of binding‐competent apo‐SurA in the periplasm.     

Oligomerization of a symmetric β‐trefoil protein in response to folding nucleus perturbation

Gene duplication and fusion events in protein evolution are postulated to be responsible for the common protein folds exhibiting internal rotational symmetry. Such evolutionary processes can also potentially yield regions of repetitive primary structure. Repetitive primary structure offers the potential for alternative definitions of critical regions, such as the folding nucleus (FN). In principle, more than one instance of the FN potentially enables an alternative folding pathway in the face of a subsequent deleterious mutation. We describe the targeted mutation of the carboxyl‐terminal region of the (internally located) FN of the de novo designed purely‐symmetric β‐trefoil protein Symfoil‐4P. This mutation involves wholesale replacement of a repeating trefoil‐fold motif with a “blade” motif from a β‐propeller protein, and postulated to trap that region of the Symfoil‐4P FN in a nonproductive folding intermediate. The resulting protein (termed “Bladefoil”) is shown to be cooperatively folding, but as a trimeric oligomer. The results illustrate how symmetric protein architectures have potentially diverse folding alternatives available to them, including oligomerization, when preferred pathways are perturbed.     

Protein folding observed by time-resolved synchrotron x-ray scattering. A feasibility study.

A solution to the "protein folding" problem, the successful prediction of tertiary and quaternary protein structure from amino acid or gene sequence, would be a major advance in biology and biotechnology. Knowledge of any intermediate structure between fully unwound and folded would aid folding calculations. The use of high intensity synchrotron x-rays from the SUNY X21 beamline at National Synchrotron Light Source has been investigated as a probe of structural changes during protein folding and unfolding in solution. A temperature jump apparatus was used to study thermally-induced folding and unfolding. Scattering of solutions of myoglobin in the angular range 20 = 1-50 mrad. was measured during temperature jumps between 26 and 76 degrees C. There are clear signs of time/temperature-dependent structural changes, in the small angle region, consistent with those from other equilibrium techniques. Analysis indicates that this experimental technique can be extended to the higher angle region where theoretical calculations indicate more detailed structural information, for example when alpha-helix formation, is present.     

Comparison of protein fragments identified by limited proteolysis and by computational cutting of proteins

Here we present a comparison between protein fragments produced by limited proteolysis and those identified by computational cutting based on the building block folding model. The principles upon which the two methods are based are different. Limited proteolysis of natively folded proteins occurs at flexible sites and never at the level of chain segments of regular secondary structure such as alpha-helices. Therefore, the targets for limited proteolysis are locally unfolded regions. In contrast, the computational cutting algorithm considers the compactness of the fragments, their nonpolar buried surface area, and their isolatedness, that is, the surface area which was buried prior to the cutting and becomes exposed subsequently. Despite the different criteria, there is an overall correspondence between sites or regions of limited proteolysis with those identified by computational cutting. The computational cutting method has been applied to several model proteins for which detailed limited proteolysis data are available, namely apomyoglobin, cytochrome c, ribonuclease A, alpha-lactalbumin, and thermolysin. As expected, more cuts are obtained computationally than experimentally and the agreement is better when a number of proteolytic enzymes are used. For example, cytochrome c is cleaved by thermolysin at 56-57, 45-46, and at 80-81, and by proteinase K at 48-49 and 50-51. Incubation of the noncovalent and native-like complex of cytochrome c fragments 1-56 and 57-104 with proteinase K yielded the gapped protein species 1-48/57-104 and finally 1-40/57-104. Computational cutting of cytochrome c reproduced the major experimental observations, with cuts at 47, 64-65 or 65-66 and 80-81 and an unstable 32-47 region not assigned to any building block. The next step, not addressed in this work, is to probe the ability of the generated fragments to fold independently. Since both the computational algorithm and limited proteolysis attempt to dissect the protein folding problem, the general agreement between the two procedures is gratifying. This consistency allows us to propose the use of limited proteolysis to produce protein fragments that can adopt an independent folding and, therefore, to study folding intermediates. The results of the present study appear to validate the building block folding model and are in line with the proposal that protein folding is a hierarchical process, where parts constituting local minima of energy fold first, with their subsequent association and mutual stabilization to finally yield the global fold.     

Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin

The foldon domain constitutes the C-terminal 30 amino acid residues of the trimeric protein fibritin from bacteriophage T4. Its function is to promote folding and trimerization of fibritin. We investigated structure, stability and folding mechanism of the isolated foldon domain. The domain folds into the same trimeric beta-propeller structure as in fibritin and undergoes a two-state unfolding transition from folded trimer to unfolded monomers. The folding kinetics involve several consecutive reactions. Structure formation in the region of the single beta-hairpin of each monomer occurs on the submillisecond timescale. This reaction is followed by two consecutive association steps with rate constants of 1.9(+/-0.5)x10(6)M(-1)s(-1) and 5.4(+/-0.3)x10(6)M(-1)s(-1) at 0.58 M GdmCl, respectively. This is similar to the fastest reported bimolecular association reactions for folding of dimeric proteins. At low concentrations of protein, folding shows apparent third-order kinetics. At high concentrations of protein, the reaction becomes almost independent of protein concentrations with a half-time of about 3 ms, indicating that a first-order folding step from a partially folded trimer to the native protein (k=210 +/- 20 s(-1)) becomes rate-limiting. Our results suggest that all steps on the folding/trimerization pathway of the foldon domain are evolutionarily optimized for rapid and specific initiation of trimer formation during fibritin assembly. The results further show that beta-hairpins allow efficient and rapid protein-protein interactions during folding.     

Residual interactions in unfolded bile acid-binding protein by 19F NMR

The folding initiation mechanism of human bile acid-binding protein (BABP) has been examined by (19) F NMR. Equilibrium unfolding studies of BABP labeled with fluorine at all eight of its phenylalanine residues showed that at least two sites experience changes in solvent exposure at high denaturant concentrations. Peak assignments were made by site-specific 4FPhe incorporation. The resonances for proteins specifically labeled at Phe17, Phe47, and Phe63 showed changes in chemical shift at denaturant concentrations at which the remaining five phenylalanine residues appear to be fully solvent-exposed. Phe17 is a helical residue that was not expected to participate in a folding initiation site. Phe47 and Phe63 form part of a hydrophobic core region that may be conserved as a site for folding initiation in the intracellular lipid-binding protein family.     

Phage P22 tailspike protein: removal of head-binding domain unmasks effects of folding mutations on native-state thermal stability.

A shortened, recombinant protein comprising residues 109-666 of the tailspike endorhamnosidase of Salmonella phage P22 was purified from Escherichia coli and crystallized. Like the full-length tailspike, the protein lacking the amino-terminal head-binding domain is an SDS-resistant, thermostable trimer. Its fluorescence and circular dichroism spectra indicate native structure. Oligosaccharide binding and endoglycosidase activities of both proteins are identical. A number of tailspike folding mutants have been obtained previously in a genetic approach to protein folding. Two temperature-sensitive-folding (tsf) mutations and the four known global second-site suppressor (su) mutations were introduced into the shortened protein and found to reduce or increase folding yields at high temperature. The mutational effects on folding yields and subunit folding kinetics parallel those observed with the full-length protein. They mirror the in vivo phenotypes and are consistent with the substitutions altering the stability of thermolabile folding intermediates. Because full-length and shortened tailspikes aggregate upon thermal denaturation, and their denaturant-induced unfolding displays hysteresis, kinetics of thermal unfolding were measured to assess the stability of the native proteins. Unfolding of the shortened wild-type protein in the presence of 2% SDS at 71 degrees C occurs at a rate of 9.2 x 10(-4) s(-1). It reflects the second kinetic phase of unfolding of the full-length protein. All six mutations were found to affect the thermal stability of the native protein. Both tsf mutations accelerate thermal unfolding about 10-fold. Two of the su mutations retard thermal unfolding up to 5-fold, while the remaining two mutations accelerate unfolding up to 5-fold. The mutational effects can be rationalized on the background of the recently determined crystal structure of the protein.     

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.     

Folding nucleus structure persists in thermally‐aggregated FGF‐1

An efficient protein‐folding pathway leading to target structure, and the avoidance of aggregation, is essential to protein evolution and de novo design; however, design details to achieve efficient folding and avoid aggregation are poorly understood. We report characterization of the thermally‐induced aggregate of fibroblast growth factor‐1 (FGF‐1), a small globular protein, by solid‐state NMR. NMR spectra are consistent with residual structure in the aggregate and provide evidence of a structured region that corresponds to the region of the folding nucleus. NMR data on aggregated FGF‐1 also indicate the presence of unstructured regions that exhibit hydration‐dependent dynamics and suggest that unstructured regions of aggregated FGF‐1 lie outside the folding nucleus. Since it is known that regions outside the folding nucleus fold late in the folding pathway, we postulate that these regions unfold early in the unfolding pathway and that the partially folded state is more prone to intermolecular aggregation. This interpretation is further supported by comparison with a designed protein that shares the same FGF‐1 folding nucleus sequence, but has different 1° structure outside the folding nucleus, and does not thermally aggregate. The results suggest that design of an efficient folding nucleus, and the avoidance of aggregation in the folding pathway, are potentially separable design criteria – the latter of which could principally focus upon the physicochemical properties of 1° structure outside the folding nucleus.     

Experimental support for the foldability–function tradeoff hypothesis: Segregation of the folding nucleus and functional regions in fibroblast growth factor‐1

The acquisition of function is often associated with destabilizing mutations, giving rise to the stability–function tradeoff hypothesis. To test whether function is also accommodated at the expense of foldability, fibroblast growth factor‐1 (FGF‐1) was subjected to a comprehensive φ‐value analysis at each of the 11 turn regions. FGF‐1, a β‐trefoil fold, represents an excellent model system with which to evaluate the influence of function on foldability: because of its threefold symmetric structure, analysis of FGF‐1 allows for direct comparisons between symmetry‐related regions of the protein that are associated with function to those that are not; thus, a structural basis for regions of foldability can potentially be identified. The resulting φ‐value distribution of FGF‐1 is highly polarized, with the majority of positions described as either folded‐like or denatured‐like in the folding transition state. Regions important for folding are shown to be asymmetrically distributed within the protein architecture; furthermore, regions associated with function (i.e., heparin‐binding affinity and receptor‐binding affinity) are localized to regions of the protein that fold after barrier crossing (late in the folding pathway). These results provide experimental support for the foldability–function tradeoff hypothesis in the evolution of FGF‐1. Notably, the results identify the potential for folding redundancy in symmetric protein architecture with important implications for protein evolution and design.     

Effect of deamidation on folding of ribonuclease A

The folding of ribonuclease A (RNase A) has been extensively studied by characterizing the disulfide containing intermediates using different experimental conditions and analytical techniques. So far, some aspects still remain unclear such as the role of the loop 65-72 in the folding pathway. We have studied the oxidative folding of a RNase A derivative containing at position 67 the substitution Asn --> isoAsp where the local structure of the loop 65-72 has been modified keeping intact the C65-C72 disulfide bond. By comparing the folding behavior of this mutant to that of the wild-type protein, we found that the deamidation significantly decreases the folding rate and alters the folding pathway of RNase A. Results presented here shed light on the role of the 65-72 region in the folding process of RNase A and also clarifies the effect of the deamidation on the folding/unfolding processes. On a more general ground, this study represents the first characterization of the intermediates produced along the folding of a deamidated protein.     

Cooperative folding mechanism of a beta-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.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.     

NMR and SAXS characterization of the denatured state of the chemotactic protein CheY: implications for protein folding initiation

The denatured state of a double mutant of the chemotactic protein CheY (F14N/V83T) has been analyzed in the presence of 5 M urea, using small angle X-ray scattering (SAXS) and heteronuclear magnetic resonance. SAXS studies show that the denatured protein follows a wormlike chain model. Its backbone can be described as a chain composed of rigid elements connected by flexible links. A comparison of the contour length obtained for the chain at 5 M urea with the one expected for a fully expanded chain suggests that approximately 25% of the residues are involved in residual structures. Conformational shifts of the alpha-protons, heteronuclear (15)N-[(1)H] NOEs and (15)N relaxation properties have been used to identify some regions in the protein that deviate from a random coil behavior. According to these NMR data, the protein can be divided into two subdomains, which largely coincide with the two folding subunits identified in a previous kinetic study of the folding of the protein. The first of these subdomains, spanning residues 1-70, is shown here to exhibit a restricted mobility as compared to the rest of the protein. Two regions, one in each subdomain, were identified as deviating from the random coil chemical shifts. Peptides corresponding to these sequences were characterized by NMR and their backbone (1)H chemical shifts were compared to those in the intact protein under identical denaturing conditions. For the region located in the first subdomain, this comparison shows that the observed deviation from random coil parameters is caused by interactions with the rest of the molecule. The restricted flexibility of the first subdomain and the transient collapse detected in that subunit are consistent with the conclusions obtained by applying the protein engineering method to the characterization of the folding reaction transition state.     

Time-resolved small-angle X-ray scattering investigation of the folding dynamics of heme oxygenase: implication of the scaling relationship for the submillisecond intermediates of protein folding

Polypeptide collapse is generally observed as the initial folding dynamics of proteins with more than 100 residues, and is suggested to be caused by the coil-globule transition explained by Flory's theory of polymers. To support the suggestion by establishing a scaling behavior between radius of gyration (Rg) and chain length for the initial folding intermediates, the folding dynamics of heme oxygenase (HO) was characterized by time-resolved, small-angle X-ray scattering. HO is a highly helical protein without disulfide bridges, and is the largest protein (263 residues) characterized by the method. The folding process of HO was found to contain a transient oligomerization; however, the conformation within 10 ms was demonstrated to be monomeric and to possess Rg of 26.1(+/-1.1) A. Together with the corresponding data for proteins with different chain lengths, the seven Rg values demonstrated the scaling relationship to chain length with a scaling exponent of 0.35+/-0.11, which is close to the theoretical value of 1/3 predicted for globules in solutions where monomer-monomer interactions are favored over monomer-solvent interactions (poor solvent). The finding indicated that the initial folding dynamics of proteins bears the signature of the coil-globule transition, and offers a clue to explain the folding mechanisms of proteins with different chain lengths.     

小角X-射线散射法解析多结构域蛋白质或复合物的溶液结构

随着同步辐射装置的建设与发展及各种建模方法的产生与完善,小角X-射线散射（small angle X-ray scattering,SAXS）法已经逐渐成为结构生物学中的一种重要的工具。SAXS可以用于研究溶液中生物大分子的结构及构象变化,蛋白质的组装、折叠等动态过程。本文对SAXS的基本原理、常用的研究技术和建模方法及其应用进行了综述。     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.