Long-range side-chain-main-chain interactions play crucial roles in stabilizing the (betaalpha)8 barrel motif of the alpha subunit of tryptophan synthase

Hydrophobic Clusters Regulate Surface Hydration Dynamics of Bacillus subtilis Lipase A

The surface hydration diffusivity of Bacillus subtilis Lipase A (BSLA) has been characterized by low-field Overhauser dynamic nuclear polarization (ODNP) relaxometry using a series of spin-labeled constructs. Sites for spin-label incorporation were previously designed via an atomistic computational approach that screened for surface exposure, reflective of the surface hydration comparable to other proteins studied by this method, as well as minimal impact on protein function, dynamics, and structure of BSLA by excluding any surface site that participated in greater than 30% occupancy of a hydrogen bonding network within BSLA. Experimental ODNP relaxometry coupling factor results verify the overall surface hydration behavior for these BSLA spin-labeled sites similar to other globular proteins. Here, by plotting the ODNP parameters of relative diffusive water versus the relative bound water, we introduce an effective "phase-space" analysis, which provides a facile visual comparison of the ODNP parameters of various biomolecular systems studied to date. We find notable differences when comparing BSLA to other systems, as well as when comparing different clusters on the surface of BSLA. Specifically, we find a grouping of sites that correspond to the spin-label surface location within the two main hydrophobic core clusters of the branched aliphatic amino acids isoleucine, leucine, and valine cores observed in the BSLA crystal structure. The results imply that hydrophobic clustering may dictate local surface hydration properties, perhaps through modulation of protein conformations and samplings of the unfolded states, providing insights into how the dynamics of the hydration shell is coupled to protein motion and fluctuations.

A conserved folding nucleus sculpts the free energy landscape of bacterial and archaeal orthologs from a divergent TIM barrel family

Orthologous proteins from the three superkingdoms have conserved their structures and functions over evolutionary time. We ask whether their folding mechanisms and the structures of their partially folded states are similarly conserved, using bacterial and archaeal representatives of the IGPS TIM barrel enzyme. Comparison of circular dichroism and fluorescence spectroscopic studies reveal a highly conserved mechanism, and hydrogen–deuterium exchange mass spectrometry analyses highlight similar cores of stability in regions dominated by clusters of branched aliphatic side chains. A bioinformatics analysis of hundreds of IGPS sequences from each superkingdom shows a very highly conserved sequence, V/ILLI, that nucleates the formation of a misfolded, microsecond intermediate and has existed since the last universal common ancestor of the IGPS family of proteins.

The amino acid sequences of proteins have evolved over billions of years, preserving their structures and functions while responding to evolutionary forces. Are there conserved sequence and structural elements that preserve the protein folding mechanisms? The functionally diverse and ancient (βα)_1–8 TIM barrel motif may answer this question. We mapped the complex six-state folding free energy surface of a ∼3.6 billion y old, bacterial indole-3-glycerol phosphate synthase (IGPS) TIM barrel enzyme by equilibrium and kinetic hydrogen–deuterium exchange mass spectrometry (HDX-MS). HDX-MS on the intact protein reported exchange in the native basin and the presence of two thermodynamically distinct on- and off-pathway intermediates in slow but dynamic equilibrium with each other. Proteolysis revealed protection in a small (α1β2) and a large cluster (β5α5β6α6β7) and that these clusters form cores of stability in I_a and I_bp. The strongest protection in both states resides in β4α4 with the highest density of branched aliphatic side chain contacts in the folded structure. Similar correlations were observed previously for an evolutionarily distinct archaeal IGPS, emphasizing a key role for hydrophobicity in stabilizing common high-energy folding intermediates. A bioinformatics analysis of IGPS sequences from the three superkingdoms revealed an exceedingly high hydrophobicity and surprising α-helix propensity for β4, preceded by a highly conserved βα-hairpin clamp that links β3 and β4. The conservation of the folding mechanisms for archaeal and bacterial IGPS proteins reflects the conservation of key elements of sequence and structure that first appeared in the last universal common ancestor of these ancient proteins.

Conservation of the structure and function of bacterial tryptophan synthases

Tryptophan biosynthesis is one of the most characterized processes in bacteria, in which the enzymes from Salmonella typhimurium and Escherichia coli serve as model systems. Tryptophan synthase (TrpAB) catalyzes the final two steps of tryptophan biosynthesis in plants, fungi and bacteria. This pyridoxal 5'-phosphate (PLP)-dependent enzyme consists of two protein chains, α (TrpA) and β (TrpB), functioning as a linear αββα heterotetrameric complex containing two TrpAB units. The reaction has a complicated, multistep mechanism resulting in the β-replacement of the hydroxyl group of l-serine with an indole moiety. Recent studies have shown that functional TrpAB is required for the survival of pathogenic bacteria in macrophages and for evading host defense. Therefore, TrpAB is a promising target for drug discovery, as its orthologs include enzymes from the important human pathogens Streptococcus pneumoniae, Legionella pneumophila and Francisella tularensis, the causative agents of pneumonia, legionnaires' disease and tularemia, respectively. However, specific biochemical and structural properties of the TrpABs from these organisms have not been investigated. To fill the important phylogenetic gaps in the understanding of TrpABs and to uncover unique features of TrpAB orthologs to spearhead future drug-discovery efforts, the TrpABs from L. pneumophila, F. tularensis and S. pneumoniae have been characterized. In addition to kinetic properties and inhibitor-sensitivity data, structural information gathered using X-ray crystallo-graphy is presented. The enzymes show remarkable structural conservation, but at the same time display local differences in both their catalytic and allosteric sites that may be responsible for the observed differences in catalysis and inhibitor binding. This functional dissimilarity may be exploited in the design of species-specific enzyme inhibitors.

Coordinated Network Changes across the Catalytic Cycle of Alpha Tryptophan Synthase

Networks of noncovalent interactions are important for protein structural dynamics. We used nuclear magnetic resonance chemical shift covariance analyses on an inactive variant of the alpha subunit of tryptophan synthase to map amino acid interaction networks across its catalytic cycle. Although some network connections were common to every enzyme state, many of the network connections strengthened or weakened over the catalytic cycle; these changes were highly coordinated. These results suggest a higher level of network organization. Our analyses identified periodic, second-order networks that show highly coordinated interaction changes across the catalytic cycle. These periodic networks may help synchronize the sequence of structural transitions necessary for enzyme function. Molecular dynamics simulations identified interaction changes across the catalytic cycle, including those involving the catalytic residue Glu49, which may help drive other interaction changes throughout the enzyme structure. Similar periodic networks may direct structural transitions and allosteric interactions in other proteins.

The protonation state of an evolutionarily conserved histidine modulates domainswapping stability of FoxP1

Forkhead box P (FoxP) proteins are members of the versatile Fox transcription factors, which control the timing and expression of multiple genes for eukaryotic cell homeostasis. Compared to other Fox proteins, they can form domain-swapped dimers through their DNA-binding –forkhead– domains, enabling spatial reorganization of distant chromosome elements by tethering two DNA molecules together. Yet, domain swapping stability and DNA binding affinity varies between different FoxP proteins. Experimental evidence suggests that the protonation state of a histidine residue conserved in all Fox proteins is responsible for pH-dependent modulation of these interactions. Here, we explore the consequences of the protonation state of another histidine (H59), only conserved within FoxM/O/P subfamilies, on folding and dimerization of the forkhead domain of human FoxP1. Dimer dissociation kinetics and equilibrium unfolding experiments demonstrate that protonation of H59 leads to destabilization of the domain-swapped dimer due to an increase in free energy difference between the monomeric and transition states. This pH–dependence is abolished when H59 is mutated to alanine. Furthermore, anisotropy measurements and molecular dynamics evidence that H59 has a direct impact in the local stability of helix H3. Altogether, our results highlight the relevance of H59 in domain swapping and folding stability of FoxP1.

βαβ Super-Secondary Motifs: Sequence, Structural Overview, and Pursuit of Potential Autonomously Folding βαβ Sequences from (β/α)8/TIM Barrels

βαβ super-secondary structures constitute the basic building blocks of (β/α)₈ class of proteins. Despite the success in designing super-secondary structures, till date, there is not a single example of a natural βαβ sequence known to fold in isolation. In this chapter, to address the finding the "needles" in the haystack scenario, we have combined the sequence preferences and structural features of independent βαβ motifs, dictated by natural selection, with rationally derived parameters from a designed βαβ motif adopting stable fold in solution. Guided by this approach, a set of potential βαβ sequences from (β/α)₈/TIM barrels are proposed as likely candidates for autonomously folding based on the assessment of their foldability.

Millisecond Timescale Motions Connect Amino Acid Interaction Networks in Alpha Tryptophan Synthase

Tryptophan synthase is a model system for understanding allosteric regulation within enzyme complexes. Amino acid interaction networks were previously delineated in the isolated alpha subunit (αTS) in the absence of the beta subunit (βTS). The amino acid interaction networks were different between the ligand-free enzyme and the enzyme actively catalyzing turnover. Previous X-ray crystallography studies indicated only minor localized changes when ligands bind αTS, and so, structural changes alone could not explain the changes to the amino acid interaction networks. We hypothesized that the network changes could instead be related to changes in conformational dynamics. As such, we conducted nuclear magnetic resonance relaxation studies on different substrate- and products-bound complexes of αTS. Specifically, we collected ¹⁵N R₂ relaxation dispersion data that reports on microsecond-to-millisecond timescale motion of backbone amide groups. These experiments indicated that there are conformational exchange events throughout αTS. Substrate and product binding change specific motional pathways throughout the enzyme, and these pathways connect the previously identified network residues. These pathways reach the αTS/βTS binding interface, suggesting that the identified dynamic networks may also be important for communication with the βTS subunit.

Octarellin VI: using rosetta to design a putative artificial (β/α)8 protein

The computational protein design protocol Rosetta has been applied successfully to a wide variety of protein engineering problems. Here the aim was to test its ability to design de novo a protein adopting the TIM-barrel fold, whose formation requires about twice as many residues as in the largest proteins successfully designed de novo to date. The designed protein, Octarellin VI, contains 216 residues. Its amino acid composition is similar to that of natural TIM-barrel proteins. When produced and purified, it showed a far-UV circular dichroism spectrum characteristic of folded proteins, with α-helical and β-sheet secondary structure. Its stable tertiary structure was confirmed by both tryptophan fluorescence and circular dichroism in the near UV. It proved heat stable up to 70°C. Dynamic light scattering experiments revealed a unique population of particles averaging 4 nm in diameter, in good agreement with our model. Although these data suggest the successful creation of an artificial α/β protein of more than 200 amino acids, Octarellin VI shows an apparent noncooperative chemical unfolding and low solubility.

Long-range interactions in the α subunit of tryptophan synthase help to coordinate ligand binding, catalysis, and substrate channeling

The α-subunit of tryptophan synthase (αTS) catalyzes the conversion of indole-3-glycerol phosphate to d-glyceraldehyde-3-phosphate and indole. We propose that allosteric networks intrinsic to αTS are modulated by the binding of the β-subunit to regulate αTS function. Understanding these long-range amino acid networks in αTS thus gives insight into the coordination of the two active sites within TS. In this study, we have used Ala residues as probes for structural and dynamic changes of αTS throughout its catalytic cycle, in the absence of the β-subunit. Projection analysis of the chemical shift changes by site-specific amino acid substitutions and ligand titrations indicates that αTS has three important conformational states: ligand-free, glyceraldehyde-3-phosphate-bound(like), and the active states. The amino acid networks within these conformations are different, as suggested by chemical shift correlation analysis. In particular, there are long-range connections, only in the active state, between Ala47, which reports on structural and dynamic changes associated with the general acid/base Glu49, and residues within the β2α2 loop, which contains the catalytically important Asp60 residue. These long-range interactions are likely important for coordinating chemical catalysis. In the free state, but not in the active state, there are connections between the β2α2 and β6α6 loops that likely help to coordinate substrate binding. Changes in the allosteric networks are also accompanied by protein dynamic changes. During catalytic turnover, the protein becomes more rigid on the millisecond timescale and the active-site dynamics are driven to a faster nanosecond timescale.

The structure of a truncated phosphoribosylanthranilate isomerase suggests a unified model for evolution of the (βα)8 barrel fold

The (βα)(8) barrel is one of the most common protein folds, and enzymes with this architecture display a remarkable range of catalytic activities. Many of these functions are associated with ancient metabolic pathways, and phylogenetic reconstructions suggest that the (βα)(8) barrel was one of the very first protein folds to emerge. Consequently, there is considerable interest in understanding the evolutionary processes that gave rise to this fold. In particular, much attention has been focused on the plausibility of (βα)(8) barrel evolution from homodimers of half barrels. However, we previously isolated a three-quarter-barrel-sized fragment of a (βα)(8) barrel, termed truncated phosphoribosylanthranilate isomerase (trPRAI), that is soluble and almost as thermostable as full-length N-(5'-phosphoribosyl)anthranilate isomerase (PRAI). Here, we report the NMR-derived structure of trPRAI. The subdomain is monomeric, is well ordered and adopts a native-like structure in solution. Side chains from strands β(1) (Glu3 and Lys5), β(2) (Tyr25) and β(6) (Lys122) of trPRAI repack to shield the hydrophobic core from the solvent. This result demonstrates that three-quarter barrels were viable intermediates in the evolution of the (βα)(8) barrel fold. We propose a unified model for (βα)(8) barrel evolution that combines our data, previously published work and plausible scenarios for the emergence of (initially error-prone) genetic systems. In this model, the earliest proto-cells contained diverse pools of part-barrel subdomains. Combinatorial assembly of these subdomains gave rise to many distinct lineages of (βα)(8) barrel proteins, that is, our model excludes the possibility that there was a single (βα)(8) barrel from which all present examples are descended.

Betaalpha-hairpin clamps brace betaalphabeta 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⁻¹ 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 (βα)₈ 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.

Structural analysis of kinetic folding intermediates for a TIM barrel protein, indole-3-glycerol phosphate synthase, by hydrogen exchange mass spectrometry and Gō model simulation

The structures of partially folded states appearing during the folding of a (betaalpha)(8) TIM barrel protein, the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS), was assessed by hydrogen exchange mass spectrometry (HX-MS) and Gō model simulations. HX-MS analysis of the peptic peptides derived from the pulse-labeled product of the sub-millisecond folding reaction from the urea-denatured state revealed strong protection in the (betaalpha)(4) region, modest protection in the neighboring (betaalpha)(1-3) and (betaalpha)(5)beta(6) segments and no significant protection in the remaining N and C-terminal segments. These results demonstrate that this species is not a collapsed form of the unfolded state under native-favoring conditions nor is it the native state formed via fast-track folding. However, the striking contrast of these results with the strong protection observed in the (betaalpha)(2-5)beta(6) region after 5 s of folding demonstrates that these species represent kinetically distinct folding intermediates that are not identical as previously thought. A re-examination of the kinetic folding mechanism by chevron analysis of fluorescence data confirmed distinct roles for these two species: the burst-phase intermediate is predicted to be a misfolded, off-pathway intermediate, while the subsequent 5 s intermediate corresponds to an on-pathway equilibrium intermediate. Comparison with the predictions using a C(alpha) Gō model simulation of the kinetic folding reaction for sIGPS shows good agreement with the core of the structure offering protection against exchange in the on-pathway intermediate(s). Because the native-centric Gō model simulations do not explicitly include sequence-specific information, the simulation results support the hypothesis that the topology of TIM barrel proteins is a primary determinant of the folding free energy surface for the productive folding reaction. The early misfolding reaction must involve aspects of non-native structure not detected by the Gō model simulation.