Sequence Determinants of the Folding Free-Energy Landscape of beta alpha-Repeat Proteins: A Dissertation

The most common structural platform in biology, the βα-repeat classes of proteins, are represented by the (βα)8TIM barrel topology and the α/β/α sandwich, CheY-like topology. Previous studies on the folding mechanisms of several members of these proteins have suggested that the initial event during refolding involves the formation of a kinetically trapped species that at least partially unfolds before the native conformation can be accessed. The simple topologies of these proteins are thought to permit access to locally folded regions that may coalesce in non-native ways to form stable interactions leading to misfolded intermediates. In a pair of TIM barrel proteins, αTS and sIGPS, it has been shown that the core of the off-pathway folding intermediates is comprised of locally connected clusters of isoleucine, leucine and valine (ILV) residues. These clusters of Branched Aliphatic Side Chains (BASiC) have the unique ability to very effectively prevent the penetration of water to the underlying hydrogen bond networks. This property retards hydrogen exchange with solvent, strengthening main chain hydrogen bonds and linking tertiary and secondary structure in a cooperative network of interactions. This property would also promote the rapid formation of collapsed species during refolding. From this viewpoint, the locally connected topology and the appropriate distribution of ILV residues in the sequence can modulate the energy landscapes of TIM barrel proteins. Another sequence determinant of protein stability that can significantly alter the structure and stability of TIM barrels is the long-range main chain-side chain hydrogen bond. Three of these interactions have been shown to form the molecular underpinnings for the cooperative access to the native state in αTS. Global analysis results presented in Chapter II and Chapter III, suggest that the off-pathway mechanism is common to three proteins of the CheY-like topology, namely CheY, NT-NtrC and Spo0F. These results are corroborated by Gō-simulations that are able to identify the minimal structure of kinetically trapped species during the refolding of CheY and Spo0F. The extent of transient, premature structure appears to correlate with the number of ILV side chains involved in a large sequence-local cluster that is formed between the central β-sheet and helices α2, α3 and α4. The failure of Gō-simulations to detect off-pathway species during the refolding of NT-NtrC may reflect the smaller number of ILV side chains in its corresponding hydrophobic cluster. In Chapter IV, comparison of the location of large ILV clusters with the hydrogen exchange protected regions in 19 proteins, suggest that clusters of BASiC residues are the primarily determinants of the stability cores of globular proteins. Although the location of the ILV clusters is sufficient to determine a majority of the protected amides in a protein structure, the extent of protection is over predicted by the ILV cluster method. The survey of 71 TIM barrel proteins presented in Chapter V, suggests that a specific type of long-range main chain-side chain hydrogen bond, termed “βα hairpin clamp” is a common feature in the βα-repeat proteins. The location and sequence patterns observed demonstrate an evolutionary signature of the βαβ modules that are the building blocks of several βα-repeat protein families. In summary, the work presented in this thesis recognizes the role of sequence in modulating the folding free energy landscapes of proteins. The formation of off-pathway folding intermediates in three CheY-like proteins and the differences in the proposed extent of structure formed in off-pathway intermediates of these three proteins, suggest that both topology and sequence play important and concerted roles in the folding of proteins. Locally connected ILV can clusters lead to off-pathway traps, whereas the formation of the productive folding path requires the development of long-range nativelike topological features to form the native state. The ability of ILV clusters to link secondary and tertiary structure formation enables them to be at the core of this cooperative folding process. Very good correlations between the locations of ILV clusters and both strong protection against exchange and the positions of folding nuclei for a variety of proteins reported in the literature support the generality of the BASiC hypothesis. Finally, the discovery of a novel pattern of H-bond interactions in the TIM barrel architecture, between the amide hydrogen of a core ILV residue with a polar side chain, bracketing βαβ modules, suggests a means for establishing cooperativity between different types of side chain interactions towards formation of the native structure. See Additional Files for copies of the source code for the global analysis program and the cluster analysis program

Computer design of microfluidic mixers for protein/RNA folding studies

Kinetic studies of biological macromolecules increasingly use microfluidic mixers to initiate and monitor reaction progress. A motivation for using microfluidic mixers is to reduce sample consumption and decrease mixing time to microseconds. Some applications, such as small-angle x-ray scattering, also require large ( \u3e 10 micron) sampling areas to ensure high signal-to-noise ratios and to minimize parasitic scattering. Chaotic to marginally turbulent mixers are well suited for these applications because this class of mixers provides a good middle ground between existing laminar and turbulent mixers. In this study, we model various chaotic to marginally turbulent mixing concepts such as flow turning, flow splitting, and vortex generation using computational fluid dynamics for optimization of mixing efficiency and observation volume. Design iterations show flow turning to be the best candidate for chaotic/marginally turbulent mixing. A qualitative experimental test is performed on the finalized design with mixing of 10 M urea and water to validate the flow turning unsteady mixing concept as a viable option for RNA and protein folding studies. A comparison of direct numerical simulations (DNS) and turbulence models suggests that the applicability of turbulence models to these flow regimes may be limited

Sub-millisecond time-resolved SAXS using a continuous-flow mixer and X-ray micro-beam

Small-angle X-ray scattering (SAXS) is a well established technique to probe the nanoscale structure and interactions in soft matter. It allows one to study the structure of native particles in near physiological environments and to analyze structural changes in response to variations in external conditions. The combination of microfluidics and SAXS provides a powerful tool to investigate dynamic processes on a molecular level with sub-millisecond time resolution. Reaction kinetics in the sub-millisecond time range has been achieved using continuous-flow mixers manufactured using micromachining techniques. The time resolution of these devices has previously been limited, in part, by the X-ray beam sizes delivered by typical SAXS beamlines. These limitations can be overcome using optics to focus X-rays to the micrometer size range providing that beam divergence and photon flux suitable for performing SAXS experiments can be maintained. Such micro-SAXS in combination with microfluidic devices would be an attractive probe for time-resolved studies. Here, the development of a high-duty-cycle scanning microsecond-timeresolution SAXS capability, built around the Kirkpatrick–Baez mirror-based microbeam system at the Biophysics Collaborative Access Team (BioCAT) beamline 18ID at the Advanced Photon Source, Argonne National Laboratory, is reported. A detailed description of the microbeam small-angle-scattering instrument, the turbulent flow mixer, as well as the data acquisition and control and analysis software is provided. Results are presented where this apparatus was used to study the folding of cytochrome c. Future prospects for this technique are discussed

Atomistic structural ensemble refinement reveals non-native structure stabilizes a sub-millisecond folding intermediate of CheY

The dynamics of globular proteins can be described in terms of transitions between a folded native state and less-populated intermediates, or excited states, which can play critical roles in both protein folding and function. Excited states are by definition transient species, and therefore are difficult to characterize using current experimental techniques. Here, we report an atomistic model of the excited state ensemble of a stabilized mutant of an extensively studied flavodoxin fold protein CheY. We employed a hybrid simulation and experimental approach in which an aggregate 42 milliseconds of all-atom molecular dynamics were used as an informative prior for the structure of the excited state ensemble. This prior was then refined against small-angle X-ray scattering (SAXS) data employing an established method (EROS). The most striking feature of the resulting excited state ensemble was an unstructured N-terminus stabilized by non-native contacts in a conformation that is topologically simpler than the native state. Using these results, we then predict incisive single molecule FRET experiments as a means of model validation. This study demonstrates the paradigm of uniting simulation and experiment in a statistical model to study the structure of protein excited states and rationally design validating experiments

Modulation of frustration in folding by sequence permutation

Folding of globular proteins can be envisioned as the contraction of a random coil unfolded state toward the native state on an energy surface rough with local minima trapping frustrated species. These substructures impede productive folding and can serve as nucleation sites for aggregation reactions. However, little is known about the relationship between frustration and its underlying sequence determinants. Chemotaxis response regulator Y (CheY), a 129-amino acid bacterial protein, has been shown previously to populate an off-pathway kinetic trap in the microsecond time range. The frustration has been ascribed to premature docking of the N- and C-terminal subdomains or, alternatively, to the formation of an unproductive local-in-sequence cluster of branched aliphatic side chains, isoleucine, leucine, and valine (ILV). The roles of the subdomains and ILV clusters in frustration were tested by altering the sequence connectivity using circular permutations. Surprisingly, the stability and buried surface area of the intermediate could be increased or decreased depending on the location of the termini. Comparison with the results of small-angle X-ray-scattering experiments and simulations points to the accelerated formation of a more compact, on-pathway species for the more stable intermediate. The effect of chain connectivity in modulating the structures and stabilities of the early kinetic traps in CheY is better understood in terms of the ILV cluster model. However, the subdomain model captures the requirement for an intact N-terminal domain to access the native conformation. Chain entropy and aliphatic-rich sequences play crucial roles in biasing the early events leading to frustration in the folding of CheY

βα-Hairpin Clamps Brace βαβ Modules and Can Make Substantive Contributions to the Stability of TIM Barrel Proteins

Non-local hydrogen bonding interactions between main chain amide hydrogen atoms and polar side chain acceptors that bracket consecutive βα or αβ elements of secondary structure in αTS from E. coli, a TIM barrel protein, have previously been found to contribute 4–6 kcal mol−1 to the stability of the native conformation. Experimental analysis of similar βα-hairpin clamps in a homologous pair of TIM barrel proteins of low sequence identity, IGPS from S. solfataricus and E. coli, reveals that this dramatic enhancement of stability is not unique to αTS. A survey of 71 TIM barrel proteins demonstrates a 4-fold symmetry for the placement of βα-hairpin clamps, bracing the fundamental βαβ building block and defining its register in the (βα)8 motif. The preferred sequences and locations of βα-hairpin clamps will enhance structure prediction algorithms and provide a strategy for engineering stability in TIM barrel proteins

A tightly packed hydrophobic cluster directs the formation of an off-pathway sub-millisecond folding intermediate in the alpha subunit of tryptophan synthase, a TIM barrel protein

Protein misfolding is now recognized as playing a crucial role in both normal and pathogenic folding reactions. An interesting example of misfolding at the earliest state of a natural folding reaction is provided by the alpha-subunit of tryptophan synthase, a (beta/alpha)(8) TIM barrel protein. The molecular basis for the formation of this off-pathway misfolded intermediate, I(BP), and a subsequent on-pathway intermediate, I1, was probed by mutational analysis of 20 branched aliphatic side-chains distributed throughout the sequence. The elimination of I(BP) and the substantial destabilization of I1 by replacement of a selective set of the isoleucine, leucine or valine residues (ILV) with alanine in a large ILV cluster external-to-the-barrel and spanning the N and C termini (cluster 2) implies tight-packing at most sites in both intermediates. Differential effects on I(BP) and I1 for replacements in alpha3, beta4 and alpha8 at the boundaries of cluster 2 suggest that their incorporation into I1 but not I(BP) reflects non-native folds at the edges of the crucial (beta/alpha)(1-2)beta(3) core in I(BP). The retention of I(BP) and the smaller and consistent destabilization of both I(BP) and I1 by similar replacements in an internal-to-the-barrel ILV cluster (cluster 1) and a second external-to-the-barrel ILV cluster (cluster 3) imply molten globule-like packing. The tight packing inferred, in part, for I(BP) or for all of I1 in cluster 2, but not in clusters 1 and 3, may reflect the larger size of cluster 2 and/or the enhanced number of isoleucine, leucine and valine self-contacts in and between contiguous elements of secondary structure. Tightly packed ILV-dominated hydrophobic clusters could serve as an important driving force for the earliest events in the folding and misfolding of the TIM barrel and other members of the (beta/alpha)(n) class of proteins

Kinetic Traps in the Folding of βα-Repeat Proteins: CheY Initially Misfolds before Accessing the Native Conformation

The beta alpha-repeat class of proteins, represented by the (beta alpha)(8) barrel and the alpha/beta/alpha sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of beta alpha-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Go-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY

Non-Native Structure Appears in Microseconds during the Folding of E. coli RNase H

The folding pathway of Escherichia coli RNase H is one of the best experimentally characterized for any protein. In spite of this, spectroscopic studies have never captured the earliest events. Using continuous-flow microfluidic mixing, we have now observed the first several milliseconds of folding by monitoring the tryptophan fluorescence lifetime (60 μs dead time). Two folding intermediates are observed, the second of which is the previously characterized I(core) millisecond intermediate. The new earlier intermediate is likely on-pathway and appears to have long-range non-native structure, providing a rare example of such non-native structure formation in a folding pathway. The tryptophan fluorescence lifetimes also suggest a deviation from native packing in the second intermediate, I(core). Similar results from a fragment of RNase H demonstrate that only half of the protein is significantly involved in this early structure formation. These studies give us a view of the formation of tertiary structure on the folding pathway, which complements previous hydrogen-exchange studies that monitored only secondary structure and observed sequential native structure formation. Our results provide detailed folding information on both a timescale and a size-scale accessible to all-atom molecular dynamics simulations of protein folding

Clusters of branched aliphatic side chains serve as cores of stability in the native state of the HisF TIM barrel protein

Imidazole-3-glycerol phosphate synthase is a heterodimeric allosteric enzyme that catalyzes consecutive reactions in imidazole biosynthesis through its HisF and HisH subunits. The unusually slow unfolding reaction of the isolated HisF TIM barrel domain from the thermophilic bacteria, Thermotoga maritima, enabled an NMR-based site-specific analysis of the main-chain hydrogen bonds that stabilize its native conformation. Very strong protection against exchange with solvent deuterium in the native state was found in a subset of buried positions in alpha-helices and pervasively in the underlying beta-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the alpha7(betaalpha)8(betaalpha)1-2 and alpha2(betaalpha)3-6beta7 segments of the (betaalpha)8 barrel. The most densely packed region of the large cluster, alpha3(betaalpha)4-6beta7, correlates closely with the core of stability previously observed in computational, protein engineering and NMR dynamics studies, demonstrating a key role for this cluster in determining the thermodynamic and structural properties of the native state of HisF. When considered with the results of previous studies where ILV clusters were found to stabilize the hydrogen-bonded networks in folding intermediates for other TIM barrel proteins, it appears that clusters of branched aliphatic side chains can serve as cores of stability across the entire folding reaction coordinate of one of the most common motifs in biology