High pressure studies on transthyretin

High hydrostatic pressure (HHP) is a powerful tool to study protein folding and the dynamics and structure of folding intermediates. Aggregates and amyloids, derived from partially folding intermediates at the junction between productive and off-pathway folding, have been studied as well, which promises better understanding of the protein misfolding diseases. Here is summarized the recent data we have collected with transthyretin under pressure.     

The NMR structures of the major intermediates of the two-domain tick carboxypeptidase inhibitor reveal symmetry in its folding and unfolding pathways

There is a lack of experimental structural information about folding intermediates of multidomain proteins. Tick carboxypeptidase inhibitor (TCI) is a small, disulfide-rich protein consisting of two domains that fold and unfold autonomously through the formation of two major intermediates, IIIa and IIIb. Each intermediate contains three native disulfide bonds in one domain and six free cysteines in the other domain. Here we have determined the NMR structures of these two intermediates trapped and isolated at acidic pH in which they are stable and compared their structures with that of the native protein analyzed under the same conditions. Both IIIa and IIIb were found to contain a folded region that corresponds to the N- and C-terminal domains of TCI, respectively, with structures very similar to the corresponding regions of the native protein. The remainder of the polypeptide chains of the intermediates was shown to be unfolded in a random coil conformation. Solvent exchange measurements further indicated that the two protein domains are not completely independent, but affect each other in terms of dynamics and stability, in agreement with reported inhibitory activity data. The derived results provide structural evidence for symmetric TCI folding and unfolding mechanisms that converge in IIIa and IIIb and reveal the structural basis that accounts for the strong and simultaneous accumulation of both intermediates. Altogether, this work has important implications for a better understanding of the folding mechanisms of multidomain, disulfide-rich proteins.     

New insights into the mechanisms of protein misfolding and aggregation in amyloidogenic diseases derived from pressure studies

Hydrostatic pressure is a robust tool for studying the thermodynamics of protein folding and protein interactions, as well as the dynamics and structure of folding intermediates. One of the main innovations obtained from using high pressure is the stabilization of folding intermediates such as molten-globule conformations, thus providing a unique opportunity for characterizing their structure and dynamics. Equally important is the prospect of understanding protein misfolding diseases by using pressure to populate partially folded intermediates at the junction between productive and off-pathway folding, which may give rise to misfolded proteins, aggregates, and amyloids. High hydrostatic pressure (HHP) has also been used to dissociate nonamyloid aggregates and inclusion bodies. In many proteins, the competition between correct folding and misfolding can lead to formation of insoluble aggregates, an important problem for the biotechnology industry and for human pathologies such as amyloidosis, Alzheimer's, Parkinson's, prion's, and tumor diseases. The diversity of diseases that result from protein misfolding has made this theme an important research focus for pharmaceutical and biotechnology companies. The use of high-pressure promises to contribute to the identification of the mechanisms behind these defects and creation of therapies against these diseases.     

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.     

Protein Dynamics Governed by Interfaces of High Polarity and Low Packing Density

The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.     

Folding of a model three-helix bundle protein: a thermodynamic and kinetic analysis.

The kinetics and thermodynamics of an off-lattice model for a three-helix bundle protein are investigated as a function of a bias gap parameter that determines the energy difference between native and non-native contacts. A simple dihedral potential is used to introduce the tendency to form right-handed helices. For each value of the bias parameter, 100 trajectories of up to one microsecond are performed. Such statistically valid sampling of the kinetics is made possible by the use of the discrete molecular dynamics method with square-well interactions. This permits much faster simulations for off-lattice models than do continuous potentials. It is found that major folding pathways can be defined, although ensembles with considerable structural variation are involved. The large gap models generally fold faster than those with a smaller gap. For the large gap models, the kinetic intermediates are non-obligatory, while both obligatory and non-obligatory intermediates are present for small gap models. Certain large gap intermediates have a two-helix microdomain with one helix extended outward (as in domain-swapped dimers); the small gap intermediates have more diverse structures. The importance of studying the kinetic, as well as the thermodynamics, of folding for an understanding of the mechanism is discussed and the relation between kinetic and equilibrium intermediates is examined. It is found that the behavior of this model system has aspects that encompass both the "new" view and the "old" view of protein folding.     

High-pressure NMR spectroscopy for characterizing folding intermediates and denatured states of proteins

Extensive structural studies using high-pressure NMR spectroscopy have recently been carried out on proteins, which potentially contribute to our understanding of the mechanisms of protein folding. Pressure shifts the conformational equilibrium from higher to lower volume conformers. If the pressure is varied, starting from the folded native structure, in many cases we observe intermediate conformers before the onset of total unfolding. This enables the investigation of details of the structure and thermodynamic characteristics of various intermediate conformers of proteins under equilibrium conditions. We can also examine pressure effects on the structure and stability of some typical denatured states such as helical denatured, molten globule, and unfolded states. The high-pressure NMR method can also be used to investigate association/dissociation equilibria of oligomeric or aggregated proteins. Beside direct observation of kinetic intermediates upon pressure jump, NMR structural investigations of equilibrium conformers under pressure provide information about the structures of kinetic intermediates during folding/unfolding reactions.     

Repeat-protein folding: new insights into origins of cooperativity, stability, and topology

Although our understanding of globular protein folding continues to advance, the irregular tertiary structures and high cooperativity of globular proteins complicates energetic dissection. Recently, proteins with regular, repetitive tertiary structures have been identified that sidestep limitations imposed by globular protein architecture. Here we review recent studies of repeat-protein folding. These studies uniquely advance our understanding of both the energetics and kinetics of protein folding. Equilibrium studies provide detailed maps of local stabilities, access to energy landscapes, insights into cooperativity, determination of nearest-neighbor interaction parameters using statistical thermodynamics, relationships between consensus sequences and repeat-protein stability. Kinetic studies provide insight into the influence of short-range topology on folding rates, the degree to which folding proceeds by parallel (versus localized) pathways, and the factors that select among multiple potential pathways. The recent application of force spectroscopy to repeat-protein unfolding is providing a unique route to test and extend many of these findings.     

Prion disease--the propagation of infectious protein topologies

Prion propagation is associated with accumulation of a conformational isomer of host encoded cellular prion protein, PrP(C). Solution structures of several mammalian PrPs have now been reported and they have stimulated a significant advance in our understanding of the folding dynamics of PrP. Studies on recombinant PrP have shown the polypeptide chain is able to adopt different topologies in different solvent conditions. Concomitantly, advances in the analysis of the abnormal isoform, PrP(Sc), have expanded our knowledge on the molecular basis of prion strains and have done much to reinforce the protein-only hypothesis of prion replication.     

Major kinetic traps for the oxidative folding of leech carboxypeptidase inhibitor

The leech carboxypeptidase inhibitor (LCI) is a 66-amino acid protein, containing four disulfides that stabilize its structure. This polypeptide represents an excellent model for the study and understanding of the diversity of folding pathways in small, cysteine-rich proteins. The pathway of oxidative folding of LCI has been elucidated in this work, using structural and kinetic analysis of the folding intermediates trapped by acid quenching. Reduced and denatured LCI refolds through a rapid, sequential flow of one- and two-disulfide intermediates and reaches a rate-limiting step in which a mixture of three major three-disulfide species and a heterogeneous population of non-native four-disulfide (scrambled) isomers coexist. The three three-disulfide intermediates have been identified as major kinetic traps along the folding pathway of LCI, and their disulfide structures have been elucidated in this work. Two of them contain only native disulfide pairings, and one contains one native and two non-native disulfide bonds. The coexistence of three-disulfide kinetic traps adopting native disulfide bonds together with a significant proportion of fully oxidized scrambled isomers shows that the folding pathway of LCI features properties exhibited by both the bovine pancreatic trypsin inhibitor and hirudin, two diverse models with extreme folding characteristics. The results further demonstrate the large diversity of disulfide folding pathways.     

Theory of protein folding

Protein folding should be complex. Proteins organize themselves into specific three-dimensional structures, through a myriad of conformational changes. The classical view of protein folding describes this process as a nearly sequential series of discrete intermediates. In contrast, the energy landscape theory of folding considers folding as the progressive organization of an ensemble of partially folded structures through which the protein passes on its way to the natively folded structure. As a result of evolution, proteins have a rugged funnel-like landscape biased toward the native structure. Connecting theory and simulations of minimalist models with experiments has completely revolutionized our understanding of the underlying mechanisms that control protein folding.     

Single-molecule approaches to prion protein misfolding

The structural conversion of the prion protein PrP into a transmissible, misfolded form is the central element of prion disease, yet there is little consensus as to how it occurs. Key aspects of conversion into the diseased state remain unsettled, from details about the earliest stages of misfolding such as the involvement of partially- or fully-unfolded intermediates to the structure of the infectious state. Part of the difficulty in understanding the structural conversion arises from the complexity of the underlying energy landscapes. Single molecule methods provide a powerful tool for probing complex folding pathways as in prion misfolding, because they allow rare and transient events to be observed directly. We discuss recent work applying single-molecule probes to study misfolding in prion proteins, and what it has revealed about the folding dynamics of PrP that may underlie its unique behavior. We also discuss single-molecule studies probing the interactions that stabilize non-native structures within aggregates, pointing the way to future work that may help identify the microscopic events triggering pathogenic conversion. Although single-molecule approaches to misfolding are relatively young, they have a promising future in prion science.     

Free Energy Landscape and Multiple Folding Pathways of an H-Type RNA Pseudoknot

How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.     

Energetic Frustrations in Protein Folding at Residue Resolution: A Homologous Simulation Study of Im9 Proteins

Energetic frustration is becoming an important topic for understanding the mechanisms of protein folding, which is a long-standing big biological problem usually investigated by the free energy landscape theory. Despite the significant advances in probing the effects of folding frustrations on the overall features of protein folding pathways and folding intermediates, detailed characterizations of folding frustrations at an atomic or residue level are still lacking. In addition, how and to what extent folding frustrations interact with protein topology in determining folding mechanisms remains unclear. In this paper, we tried to understand energetic frustrations in the context of protein topology structures or native-contact networks by comparing the energetic frustrations of five homologous Im9 alpha-helix proteins that share very similar topology structures but have a single hydrophilic-to-hydrophobic mutual mutation. The folding simulations were performed using a coarse-grained Gō-like model, while non-native hydrophobic interactions were introduced as energetic frustrations using a Lennard-Jones potential function. Energetic frustrations were then examined at residue level based on φ-value analyses of the transition state ensemble structures and mapped back to native-contact networks. Our calculations show that energetic frustrations have highly heterogeneous influences on the folding of the four helices of the examined structures depending on the local environment of the frustration centers. Also, the closer the introduced frustration is to the center of the native-contact network, the larger the changes in the protein folding. Our findings add a new dimension to the understanding of protein folding the topology determination in that energetic frustrations works closely with native-contact networks to affect the protein folding.     

Autonomous subdomains in protein folding.

Proteolytic dissection of native trp repressor and horse heart cytochrome c has been used to infer some of the steps in the folding pathways of the intact proteins. For both proteins, small fragments are capable of undergoing spontaneous noncovalent association to form subdomains with native-like secondary and/or tertiary structural features, suggesting that dissection/reassembly may be a general method to gain insight into the structures of folding intermediates. The importance of this approach is its simplicity and potential applicability to studying the folding pathways of a wide range of proteins. The proteases report on the structure and dynamics of the native state, circumventing the need for prior knowledge of the structures of folding intermediates. The observation that small fragments of proteins can associated noncovalently suggests that protein folding can be viewed as an intramolecular "recognition" process. The results imply that substantial information about protein structure and folding is encoded at the level of subdomains, and that chain connectivity has only a minor role in determining the fold.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

Structure-Based Derivation of Protein Folding Intermediates and Energies from Optical Tweezers

Optical tweezers (OTs) measure the force-dependent time-resolved extension of a single macromolecule tethered between two trapped beads. From this measurement, it is possible to determine the folding intermediates, energies, and kinetics of the macromolecule. Previous data analysis generally has used the extension as a reaction coordinate to characterize the observed folding transitions. Despite its convenience, the extension poorly describes folding in the absence of force. Here, we chose the contour length of the unfolded polypeptide as a reaction coordinate and modeled the extensions of protein structures along their predicted folding pathways based on high-resolution structures of the proteins in their native states. We included the extension in our model to calculate the total extensions, energies, and transition rates of the proteins as a function of force. We fit these calculations to the corresponding experimental measurements and obtained the best-fit conformations and energies of proteins in different folding states. We applied our method to analyze single-molecule trajectories of two representative protein complexes responsible for membrane fusion, the HIV-1 glycoprotein 41 and the synaptic SNARE proteins, which involved transitions between two and five states, respectively. Nonlinear fitting of the model to the experimental data revealed the structures of folding intermediates and transition states and their associated energies. Our results demonstrate that the contour length is a useful reaction coordinate to characterize protein folding and that intrinsic extensions of protein structures should be taken into account to properly derive the conformations and energies of protein folding intermediates from single-molecule manipulation experiments.     

Highly fluctuating protein structures revealed by variable-pressure nuclear magnetic resonance

Although our knowledge of basic folded structures of proteins has dramatically improved, the extent of our corresponding knowledge of higher-energy conformers remains extremely slim. The latter information is crucial for advancing our understanding of mechanisms of protein function, folding, and conformational diseases. Direct spectroscopic detection and analysis of structures of higher-energy conformers are limited, particularly under physiological conditions, either because their equilibrium populations are small or because they exist only transiently in the folding process. A new experimental strategy using pressure perturbation in conjunction with multidimensional NMR spectroscopy is being used to overcome this difficulty. A number of rare conformers are detected under pressure for a variety of proteins such as the Ras-binding domain of RalGDS, beta-lactoglobulin, dihydrofolate reductase, ubiquitin, apomyoglobin, p13(MTCP1), and prion, which disclose a rich world of protein structure between basically folded and globally unfolded states. Specific structures suggest that these conformers are designed for function and are closely identical to kinetic intermediates. Detailed structural determination of higher-energy conformers with variable-pressure NMR will extend our knowledge of protein structure and conformational fluctuation over most of the biologically relevant conformational space.     

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering,Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.