A calcium-sensitive antibody isolates soluble amyloid-ß aggregates and fibrils from Alzheimer's disease brain.

Aqueously soluble oligomers of amyloid-ß peptide may be the principal neurotoxic forms of amyloid-ß in Alzheimer's disease, initiating downstream events that include tau hyperphosphorylation, neuritic/synaptic injury, microgliosis and neuron loss. Synthetic oligomeric amyloid-ß has been studied extensively, but little is known about the biochemistry of natural oligomeric amyloid-ß in human brain, even though it is more potent than simple synthetic peptides and comprises truncated and modified amyloid-ß monomers. We hypothesized that monoclonal antibodies specific to neurotoxic oligomeric amyloid-ß could be used to isolate it for further study. Here we report a unique human monoclonal antibody (B24) raised against synthetic oligomeric amyloid-ß that potently prevents Alzheimer's disease brain oligomeric amyloid-ß-induced impairment of hippocampal long-term potentiation. B24 binds natural and synthetic oligomeric amyloid-ß and a subset of amyloid plaques, but only in the presence of Ca2+. The amyloid-ß N terminus is required for B24 binding. Hydroxyapatite chromatography revealed that natural oligomeric amyloid-ß is highly avid for Ca2+. We took advantage of the reversible Ca2+-dependence of B24 binding to perform non-denaturing immunoaffinity isolation of oligomeric amyloid-ß from Alzheimer's disease brain-soluble extracts. Unexpectedly, the immunopurified material contained amyloid fibrils visualized by electron microscopy and amenable to further structural characterization. B24-purified human oligomeric amyloid-ß inhibited mouse hippocampal long-term potentiation. These findings identify a calcium-dependent method for purifying bioactive brain oligomeric amyloid-ß, at least some of which appears fibrillar.

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.

Sequential Closure of Loop Structures Forms the Folding Nucleus during the Refolding Transition of the Escherichia coli Adenylate Kinase Molecule.

The ensemble of conformers of globular protein molecules immediately following transfer from unfolding to folding conditions is assumed to be collapsed though still disordered, as the first steps of the folding pathway are initiated. In order to test the hypothesis that long loop closure transitions are part of the initiation of the folding pathway, our groups are studying the initiation of the folding transition of a model protein by time-resolved excitation energy transfer (trFRET) detected fast kinetics experiments. Site-specific double labeling is used to study the timing of conformational transitions of individual loop forming chain segments at the microsecond time regime. Previously, it was shown that at least three long loops in the Escherichia coli adenylate kinase (AK) molecule close within the first 5 ms of folding of AK, while the main global folding transition occurs in a time regime of seconds. In order to enhance the time resolution of the kinetics experiments to the microsecond time regime and determine the rate of closure of the two N terminal loops (loop I residues 1-26 and loop II residues 29-72), we applied a continuous flow based double kinetics experiment. These measurements enabled us to obtain a microsecond series of transient time dependent distributions of distances between the ends of the labeled loops. Analysis of the trFRET experiments show that the N terminal loop (loop I) is closed within less than 60 µs after the initiation of refolding. Loop II is also mostly closed within that time step but shows an additional small reduction of the mean end-to-end distance in a second phase at a rate of 0.005 µs(-1). This second phase can either reflect tightening of a loosely closed loop in the ensemble of conformers or may reflect two subpopulations in the ensemble, which differ in the rate of closure of loop II, but not in the rate of closure of loop I. This study shows the very fast closure of long loops in the otherwise disordered backbone and fine details of the very early hidden pretransition state steps that are essential for the fast and efficient folding of the protein molecule.

A Hydrophobic Core Stabilizes the Residual Structure in the RRM2 Intermediate State of the ALS-linked Protein TDP-43.

Folding intermediates mediate both protein folding and the misfolding and aggregation observed in human diseases, including amyotrophic lateral sclerosis (ALS), and are prime targets for therapeutic interventions. In this study, we identified the core nucleus of structure for a folding intermediate in the second RNA recognition motif (RRM2) of the ALS-linked RNA-binding protein, TDP-43 (TAR DNA-binding protein-43), using a combination of experimental and computational approaches. Urea equilibrium unfolding studies revealed that the RRM2 intermediate state consists of collapsed residual secondary structure localized to the N-terminal half of RRM2, while the C-terminus is largely disordered. Steered molecular dynamics simulations and mutagenesis studies yielded key stabilizing hydrophobic contacts that, when mutated to alanine, severely disrupt the overall fold of RRM2. In combination, these findings suggest a role for this RRM intermediate in normal TDP-43 function as well as serving as a template for misfolding and aggregation through the low stability and non-native secondary structure.

A hydrophobic core stabilizes the residual structure in the RRM2 intermediate state of the ALS-linked protein TDP-43.

Folding intermediates mediate both protein folding and the misfolding and aggregation observed in human diseases, including amyotrophic lateral sclerosis (ALS), and are prime targets for therapeutic interventions. In this study, we identified the core nucleus of structure for a folding intermediate in the second RNA recognition motif (RRM2) of the ALS-linked RNA-binding protein, TDP-43, using a combination of experimental and computational approaches. Urea equilibrium unfolding studies revealed that the RRM2 intermediate state consists of collapsed residual secondary structure localized to the N-terminal half of RRM2, while the C-terminus is largely disordered. Steered molecular dynamics simulations and mutagenesis studies yielded key stabilizing hydrophobic contacts that, when mutated to alanine, severely disrupt the overall fold of RRM2. In combination, these findings suggest a role for this RRM intermediate in normal TDP-43 function as well as serving as a template for misfolding and aggregation through the low stability and non-native secondary structure.

Assessment and incorporation of in vitro correlates to pharmacokinetic outcomes in antibody developability workflows.

In vitro assessments for the prediction of pharmacokinetic (PK) behavior of biotherapeutics can help identify corresponding liabilities significantly earlier in the discovery timeline. This can minimize the need for extensive early in vivo PK characterization, thereby reducing animal usage and optimizing resources. In this study, we recommend bolstering classical developability workflows with in vitro measures correlated with PK. In agreement with current literature, in vitro measures assessing nonspecific interactions, self-interaction, and FcRn interaction are demonstrated to have the highest correlations to clearance in hFcRn Tg32 mice. Crucially, the dataset used in this study has broad sequence diversity and a range of physicochemical properties, adding robustness to our recommendations. Finally, we demonstrate a computational approach that combines multiple in vitro measurements with a multivariate regression model to improve the correlation to PK compared to any individual assessment. Our work demonstrates that a judicious choice of high throughput in vitro measurements and computational predictions enables the prioritization of candidate molecules with desired PK properties.

A next-generation Fab library platform directly yielding drug-like antibodies with high affinity, diversity, and developability.

We previously described an in vitro single-chain fragment (scFv) library platform originally designed to generate antibodies with excellent developability properties. The platform design was based on the use of clinical antibodies as scaffolds into which replicated natural complementarity-determining regions purged of sequence liabilities were inserted, and the use of phage and yeast display to carry out antibody selection. In addition to being developable, antibodies generated using our platform were extremely diverse, with most campaigns yielding sub-nanomolar binders. Here, we describe a platform advancement that incorporates Fab phage display followed by single-chain antibody-binding fragment Fab (scFab) yeast display. The scFab single-gene format provides balanced expression of light and heavy chains, with enhanced conversion to IgG, thereby combining the advantages of scFvs and Fabs. A meticulously engineered, quality-controlled Fab phage library was created using design principles similar to those used to create the scFv library. A diverse panel of binding scFabs, with high conversion efficiency to IgG, was isolated against two targets. This study highlights the compatibility of phage and yeast display with a Fab semi-synthetic library design, offering an efficient approach to generate drug-like antibodies directly, facilitating their conversion to potential therapeutic candidates.

Sub-millisecond time-resolved SAXS using a continuous-flow mixer and X-ray microbeam.

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-time-resolution 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.

Advances in turbulent mixing techniques to study microsecond protein folding reactions.

Recent experimental and computational advances in the protein folding arena have shown that the readout of the one-dimensional sequence information into three-dimensional structure begins within the first few microseconds of folding. The initiation of refolding reactions has been achieved by several means, including temperature jumps, flash photolysis, pressure jumps, and rapid mixing methods. One of the most commonly used means of initiating refolding of chemically denatured proteins is by turbulent flow mixing with refolding dilution buffer, where greater than 99% mixing efficiency has been achieved within 10's of microseconds. Successful interfacing of turbulent flow mixers with complementary detection methods, including time-resolved Fluorescence Spectroscopy (trFL), Förster Resonance Energy Transfer, Circular Dichroism, Small-Angle X-ray Scattering, Hydrogen Exchange followed by Mass Spectrometry and Nuclear Magnetic Resonance Spectroscopy, Infrared Spectroscopy (IR), and Fourier Transform IR Spectroscopy, has made this technique very attractive for monitoring various aspects of structure formation during folding. Although continuous-flow (CF) mixing devices interfaced with trFL detection have a dead time of only 30 µs, burst phases have been detected in this time scale during folding of peptides and of large proteins (e.g., CheY and TIM barrels). Furthermore, a major limitation of the CF mixing technique has been the requirement of large quantities of sample. In this brief communication, we will discuss the recent flurry of activity in micromachining and microfluidics, guided by computational simulations, which are likely to lead to dramatic improvements in time resolution and sample consumption for CF mixers over the next few years.

Honing-in antigen-specific cells during antibody discovery: a user-friendly process to mine a deeper repertoire.

Immunization based antibody discovery is plagued by the paucity of antigen-specific B cells. Identifying these cells is akin to finding needle in a haystack. Current and emerging technologies while effective, are limited in terms of capturing the antigen-specific repertoire. We report on the bulk purification of antigen-specific B-cells and the benefits it offers to various antibody discovery platforms. Using five different antigens, we show hit rates of 51-88%, compared to about 5% with conventional methods. We also show that this purification is highly efficient with loss of only about 2% antigen specific cells. Furthermore, we compared clones in which cognate chains are preserved with those from display libraries in which chains either from total B cells (TBC) or antigen-specific B cells (AgSC) underwent combinatorial pairing. We found that cognate chain paired clones and combinatorial clones from AgSC library had higher frequency of functional clones and showed greater diversity in sequence and paratope compared to clones from the TBC library. This antigen-specific B-cell selection technique exemplifies a process improvement with reduced cycle time and cost, by removing undesired clones prior to screening and increasing the chance of capturing desirable and rare functional clones in the repertoire.

The Folding Pathway of 6aJL2.

One-third of the reported cases of light chain amyloidosis are related to the germ line λ6 family; remarkably, healthy individuals express this type of protein in just 2% of the peripheral blood and bone marrow B-cells. The appearance of the disease has been related to the inherent properties of this protein family. A recombinant representative model for λ6 proteins called 6aJL2 containing the amino acid sequence encoded by the 6a and JL2 germ line genes was previously designed and synthesized to study the properties of this family. Previous work on 6aJL2 suggested a simple two-state folding model at 25 °C; no intermediate could be identified either by kinetics or by fluorescence and circular dichroism equilibrium studies, although the presence of an intermediate that is populated at ∼2.4 M urea was suggested by size exclusion chromatography. In this study we employed classic equilibrium and kinetic experiments and analysis to elucidate the detailed folding mechanism of this protein. We identify species that are kinetically accessible and/or are populated at equilibrium. We describe the presence of intermediate and native-like species and propose a five-species folding mechanism at 25 °C at short incubation times, similar to and consistent with those observed in other proteins of this fold. The formation of intermediates in the mechanism of 6aJL2 is faster than that proposed for a Vκ light chain, which could be an important distinction in the amyloidogenic potential of both germ lines.

Bringing the Heavy Chain to Light: Creating a Symmetric, Bivalent IgG-Like Bispecific.

Bispecific molecules are biologically significant, yet their complex structures pose important manufacturing and pharmacokinetic challenges. Nevertheless, owing to similarities with monoclonal antibodies (mAbs), IgG-like bispecifics conceptually align well with conventional expression and manufacturing platforms and often exhibit potentially favorable drug metabolism and pharmacokinetic (DMPK) properties. However, IgG-like bispecifics do not possess target bivalency and current designs often require tedious engineering and purification to ensure appropriate chain pairing. Here, we present a near-native IgG antibody format, the 2xVH, which can create bivalency for each target or epitope and requires no engineering for cognate chain pairing. In this modality, two different variable heavy (VH) domains with distinct binding specificities are grafted onto the first constant heavy (CH1) and constant light (CL) domains, conferring the molecule with dual specificity. To determine the versatility of this format, we characterized the expression, binding, and stability of several previously identified soluble human VH domains. By grafting these domains onto an IgG scaffold, we generated several prototype 2xVH IgG and Fab molecules that display similar properties to mAbs. These molecules avoided the post-expression purification necessary for engineered bispecifics while maintaining a capacity for simultaneous dual binding. Hence, the 2xVH format represents a bivalent, bispecific design that addresses limitations of manufacturing IgG-like bispecifics while promoting biologically-relevant dual target engagement.

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.

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 (>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.

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.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability.

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

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.

Methods for analysis of size-exclusion chromatography-small-angle X-ray scattering and reconstruction of protein scattering.

Size-exclusion chromatography in line with small-angle X-ray scattering (SEC-SAXS) has emerged as an important method for investigation of heterogeneous and self-associating systems, but presents specific challenges for data processing including buffer subtraction and analysis of overlapping peaks. This paper presents novel methods based on singular value decomposition (SVD) and Guinier-optimized linear combination (LC) to facilitate analysis of SEC-SAXS data sets and high-quality reconstruction of protein scattering directly from peak regions. It is shown that Guinier-optimized buffer subtraction can reduce common subtraction artifacts and that Guinier-optimized linear combination of significant SVD basis components improves signal-to-noise and allows reconstruction of protein scattering, even in the absence of matching buffer regions. In test cases with conventional SAXS data sets for cytochrome c and SEC-SAXS data sets for the small GTPase Arf6 and the Arf GTPase exchange factors Grp1 and cytohesin-1, SVD-LC consistently provided higher quality reconstruction of protein scattering than either direct or Guinier-optimized buffer subtraction. These methods have been implemented in the context of a Python-extensible Mac OS X application known as Data Evaluation and Likelihood Analysis (DELA), which provides convenient tools for data-set selection, beam intensity normalization, SVD, and other relevant processing and analytical procedures, as well as automated Python scripts for common SAXS analyses and Guinier-optimized reconstruction of protein scattering.

Microsecond barrier-limited chain collapse observed by time-resolved FRET and SAXS.

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59-heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~27 µs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency>90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.

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 α-helices and pervasively in the underlying ß-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the α7(ßα)8(ßα)1-2 and α2(ßα)3-6ß7 segments of the (ßα)8 barrel. The most densely packed region of the large cluster, α3(ßα)4-6ß7, 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.