The role of a disulfide bridge in the stability and folding kinetics of Arabidopsis thaliana cytochrome c6A

Cytochrome c_6A is a eukaryotic member of the Class I cytochrome c family possessing a high structural homology with photosynthetic cytochrome c₆ from cyanobacteria, but structurally and functionally distinct through the presence of a disulfide bond and a heme mid-point redox potential of + 71 mV (vs normal hydrogen electrode). The disulfide bond is part of a loop insertion peptide that forms a cap-like structure on top of the core α-helical fold. We have investigated the contribution of the disulfide bond to thermodynamic stability and (un)folding kinetics in cytochrome c_6A from Arabidopsis thaliana by making comparison with a photosynthetic cytochrome c₆ from Phormidium laminosum and through a mutant in which the Cys residues have been replaced with Ser residues (C67/73S). We find that the disulfide bond makes a significant contribution to overall stability in both the ferric and ferrous heme states. Both cytochromes c_6A and c₆ fold rapidly at neutral pH through an on-pathway intermediate. The unfolding rate for the C67/73S variant is significantly increased indicating that the formation of this region occurs late in the folding pathway. We conclude that the disulfide bridge in cytochrome c_6A acts as a conformational restraint in both the folding intermediate and native state of the protein and that it likely serves a structural rather than a previously proposed catalytic role.     

Folding Topology of a Bimolecular DNA Quadruplex Containing a Stable Mini-hairpin Motif within the Diagonal Loop

We describe the NMR structural characterisation of a bimolecular anti-parallel DNA quadruplex d(G₃ACGTAGTG₃)₂ containing an autonomously stable mini-hairpin motif inserted within the diagonal loop. A folding topology is identified that is different from that observed for the analogous d(G₃T₄G₃)₂ dimer with the two structures differing in the relative orientation of the diagonal loops. This appears to reflect specific base stacking interactions at the quadruplex-duplex interface that are not present in the structure with the T₄-loop sequence. A truncated version of the bimolecular quadruplex d(G₂ACGTAGTG₂)₂, with only two core G-tetrads, is less stable and forms a heterogeneous mixture of three 2-fold symmetric quadruplexes with different loop arrangements. We demonstrate that the nature of the loop sequence, its ability to form autonomously stable structure, the relative stabilities of the hairpin loop and core quadruplex, and the ability to form favourable stacking interactions between these two motifs are important factors in controlling DNA G-quadruplex topology.     

Perturbed folding kinetics of circularly permuted RNAs with altered topology

The folding pathway of the Tetrahymena ribozyme correlates inversely with the sequence distance between native interactions, or contact order. The rapidly folding P4-P6 domain has a low contact order, while the slowly folding P3-P7 region has a high contact order. To examine the role of topology and contact order in RNA folding, we screened for circular permutants of the ribozyme that retain catalytic activity. Permutants beginning in the P4-P6 domain fold 5 to 20 times more slowly than the wild-type ribozyme. By contrast, 50% of a permuted RNA that disjoins a non-native interaction in P3 folds tenfold faster than the wild-type ribozyme. Hence, the probability of rapidly folding to the native state depends on the topology of tertiary domains.     

Thermal unfolding studies show the disease causing F508del mutation in CFTR thermodynamically destabilizes nucleotide-binding domain 1

Misfolding and degradation of CFTR is the cause of disease in patients with the most prevalent CFTR mutation, an in-frame deletion of phenylalanine (F508del), located in the first nucleotide-binding domain of human CFTR (hNBD1). Studies of (F508del)CFTR cellular folding suggest that both intra- and inter-domain folding is impaired. (F508del)CFTR is a temperature-sensitive mutant, that is, lowering growth temperature, improves both export, and plasma membrane residence times. Yet, paradoxically, F508del does not alter the fold of isolated hNBD1 nor did it seem to perturb its unfolding transition in previous isothermal chemical denaturation studies. We therefore studied the in vitro thermal unfolding of matched hNBD1 constructs ±F508del to shed light on the defective folding mechanism and the basis for the thermal instability of (F508del)CFTR. Using primarily differential scanning calorimetry (DSC) and circular dichroism, we show for all hNBD1 pairs studied, that F508del lowers the unfolding transition temperature (T_m) by 6–7°C and that unfolding occurs via a kinetically-controlled, irreversible transition in isolated monomers. A thermal unfolding mechanism is derived from nonlinear least squares fitting of comprehensive DSC data sets. All data are consistent with a simple three-state thermal unfolding mechanism for hNBD1 ± F508del: N(±MgATP) ⇄ I_T(±MgATP) → A_T → (A_T)_n. The equilibrium unfolding to intermediate, I_T, is followed by the rate-determining, irreversible formation of a partially folded, aggregation-prone, monomeric state, A_T, for which aggregation to (A_T)_n and further unfolding occur with no detectable heat change. Fitted parameters indicate that F508del thermodynamically destabilizes the native state, N, and accelerates the formation of A_T.     

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 (βα)₈ 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 (βα)₄ region, modest protection in the neighboring (βα)_1-3 and (βα)₅β₆ 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 (βα)_2-5β₆ 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^α 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.     

Common motifs and topological effects in the protein folding transition state

Through extensive experiment, simulation, and analysis of protein S6 (1RIS), we find that variations in nucleation and folding pathway between circular permutations are determined principally by the restraints of topology and specific nucleation, and affected by changes in chain entropy. Simulations also relate topological features to experimentally measured stabilities. Despite many sizable changes in phi values and the structure of the transition state ensemble that result from permutation, we observe a common theme: the critical nucleus in each of the mutants share a subset of residues that can be mapped to the critical nucleus residues of the wild-type. Circular permutations create new N and C termini, which are the location of the largest disruption of the folding nucleus, leading to a decrease in both phi values and the role in nucleation. Mutant nuclei are built around the wild-type nucleus but are biased towards different parts of the S6 structure depending on the topological and entropic changes induced by the location of the new N and C termini.     

Autonomously Folding Protein Fragments Reveal Differences in the Energy Landscapes of Homologous RNases H

An important approach to understanding how a protein sequence encodes its energy landscape is to compare proteins with different sequences that fold to the same general native structure. In this work, we compare E. coli and T. thermophilus homologs of the protein RNase H. Using protein fragments, we create equilibrium mimics of two different potential partially-folded intermediates (I_core and I_core+1) hypothesized to be present on the energy landscapes of these two proteins. We observe that both T. thermophilus RNase H (ttRNH) fragments are folded and have distinct stabilities, indicating that both regions are capable of autonomous folding and that both intermediates are present as local minima on the ttRNH energy landscape. In contrast, the two E. coli RNase H (ecRNH) fragments have very similar stabilities, suggesting that the presence of additional residues in the I_core+1 fragment does not affect the folding or structure as compared to I_core. NMR experiments provide additional evidence that only the I_core intermediate is populated by ecRNH. This is one of the biggest differences that has been observed between the energy landscapes of these two proteins. Additionally, we used a FRET experiment in the background of full-length ttRNH to specifically monitor the formation of the I_core+1 intermediate. We determine that the ttRNH I_core+1 intermediate is likely the intermediate populated prior to the rate-limiting barrier to global folding, in contrast to E. coli RNase H for which I_core is the folding intermediate. This result provides new insight into the nature of the rate-limiting barrier for the folding of RNase H.     

Circular permutants in β-glucosidases (family 3) within a predicted double-domain topology that includes a (β/α)8-barrel

By predicting the general secondary structure for β-glucosidases (family 3), in conjunction with existing knowledge of the circular permutants present in B. fibrisolvens and R. albus, we were able to find the canonical elements of the secondary structure. The way these elements are linked suggests that there is a double-domain topology made up of a (β/α)₈-barrel domain and a “mainly all-β” domain. A number of already known conserved motifs are located within (or near) the C-terminal part of the putative parallel β-strands of the (β/α)₈-barrel, which is consistent with what is known about the location of catalytical sites for enzymes that have this domain topology. Within the circular permutants, two β/α units are located at the N-terminal part of the molecule, whereas the other six β/α units are located at the C-terminal end. In this way, the circular permutants can be seen to have a putative discontinuous double-domain topology. Proteins 31:214–223, 1998. © 1998 Wiley-Liss, Inc.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Acinetobacter baumannii Virulence Is Mediated by the Concerted Action of Three Phospholipases D

Acinetobacter baumannii causes a broad range of opportunistic infections in humans. Its success as an emerging pathogen is due to a combination of increasing antibiotic resistance, environmental persistence and adaptation to the human host. To date very little is known about the molecular basis of the latter. Here we demonstrate that A. baumannii can use phosphatidylcholine, an integral part of human cell membranes, as sole carbon and energy source. We report on the identification of three phospholipases belonging to the PLD superfamily. PLD1 and PLD2 appear restricted to the bacteria and display the general features of bacterial phospholipases D. They possess two PLDc_2 PFAM domains each encompassing the HxKx₄Dx₆GS/GGxN (HKD) motif necessary for forming the catalytic core. The third candidate, PLD3, is found in bacteria as well as in eukaryotes and harbours only one PLDc_2 PFAM domain and one conserved HKD motif, which however do not overlap. Employing a markerless mutagenesis system for A. baumannii ATCC 19606^T, we generated a full set of PLD knock-out mutants. Galleria mellonella infection studies as well as invasion experiments using A549 human lung epithelial cells revealed that the three PLDs act in a concerted manner as virulence factors and are playing an important role in host cell invasion.     

Folding of circular permutants with decreased contact order: general trend balanced by protein stability.

To examine the influence of contact order and stability on the refolding rate constant for two-state proteins, we have analysed the folding kinetics of the small beta-alpha-beta protein S6 and two of its circular permutants with relative contact orders of 0.19, 0.15 and 0.12. Data reveal a small but significant increase of the refolding rate constant (log k(f)) with decreasing contact order. At the same time, the decreased contact order is correlated to losses in global stability and alterations of the folding nucleus. When the differences in stability are accounted for by addition of Na2SO4 or by comparison of the folding kinetics at the transition mid-point, the dependence between log k(f) and contact order becomes stronger and follows the general correlation for two-state proteins. The observation emphasizes the combined action of topology and stability in controlling the rate constant of protein folding.     

Malleability of protein folding pathways: a simple reason for complex behaviour

Although the structures of native proteins are generally unique, the pathways by which they form are often free to vary. Some proteins fold by a multitude of different pathways, whereas others seem restricted to only one choice. An explanation for this variation in folding behaviour has recently emerged from studies of transition state changes: the number of accessible pathways is linked to the number of nucleation motifs contained within the native topology. We refer to these nucleation motifs as 'foldons', as they approach the size of an independent cooperative unit. Thus, with respect to pathway malleability and the composition of the folding funnel, proteins can be seen as modular assemblies of competing foldons. For the split beta-alpha-beta fold, these foldons are two-strand-helix motifs coupled by spatial overlap.     

Residue-Specific Analysis of Frustration in the Folding Landscape of Repeat β/α Protein Apoflavodoxin

Flavodoxin adopts the common repeat β/α topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central β-strand 3 (i.e., at β-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, ?-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the β-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central β-sheet directs initial misfolding, while proper alignment on both sides of β-strand 3 is necessary for productive folding.     

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G_A88 and G_B88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G_A88 and G_B88. We show that, contrary to expectation, G_B88, characterized by a native α+β fold, displays in the denatured state a content of native-like helical structure greater than G_A88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.     

Ab initio folding of a trefoil‐fold motif reveals structural similarity with a β‐propeller blade motif

Many protein architectures exhibit evidence of internal rotational symmetry postulated to be the result of gene duplication/fusion events involving a primordial polypeptide motif. A common feature of such structures is a domain‐swapped arrangement at the interface of the N‐ and C‐termini motifs and postulated to provide cooperative interactions that promote folding and stability. De novo designed symmetric protein architectures have demonstrated an ability to accommodate circular permutation of the N‐ and C‐termini in the overall architecture; however, the folding requirement of the primordial motif is poorly understood, and tolerance to circular permutation is essentially unknown. The β‐trefoil protein fold is a threefold‐symmetric architecture where the repeating ~42‐mer “trefoil‐fold” motif assembles via a domain‐swapped arrangement. The trefoil‐fold structure in isolation exposes considerable hydrophobic area that is otherwise buried in the intact β‐trefoil trimeric assembly. The trefoil‐fold sequence is not predicted to adopt the trefoil‐fold architecture in ab initio folding studies; rather, the predicted fold is closely related to a compact “blade” motif from the β‐propeller architecture. Expression of a trefoil‐fold sequence and circular permutants shows that only the wild‐type N‐terminal motif definition yields an intact β‐trefoil trimeric assembly, while permutants yield monomers. The results elucidate the folding requirements of the primordial trefoil‐fold motif, and also suggest that this motif may sample a compact conformation that limits hydrophobic residue exposure, contains key trefoil‐fold structural features, but is more structurally homologous to a β‐propeller blade motif.     

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3?kcal?mol^?1. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.     

A novel chair-type G-quadruplex formed by a Bombyx mori telomeric sequence

Recently, the human telomeric d[TAGGG(TTAGGG)₃] sequence has been shown to form in K⁺ solution an intramolecular (3+1) G-quadruplex structure, whose G-tetrad core contains three strands oriented in one direction and the fourth in the opposite direction. Here we present a study on the structure of the Bombyx mori telomeric d[TAGG(TTAGG)₃] sequence, which differs from the human counterpart only by one G deletion in each repeat. We found that this sequence adopted multiple G-quadruplex structures in K⁺ solution. We have favored a major G-quadruplex form by a judicious U-for-T substitution in the sequence and determined the folding topology of this form. We showed by NMR that this was a new chair-type intramolecular G-quadruplex which involved a two-layer antiparallel G-tetrad core and three edgewise loops. Our result highlights the effect of G-tract length on the folding topology of G-quadruplexes, but also poses the question of whether a similar chair-type G-quadruplex fold exists in the human telomeric sequences.     

Core requirements for glmS ribozyme self-cleavage reveal a putative pseudoknot structure

The glmS ribozyme is a self-cleaving RNA catalyst that resides in the 5′-untranslated region of glmS mRNA in certain bacteria. The ribozyme is specifically activated by glucosamine-6-phosphate (GlcN6P), the metabolic product of the GlmS protein, and is thus proposed to provide a feedback mechanism of riboswitch regulation. Both phylogenetic and biochemical analyses of the glmS ribozyme have established a highly conserved core sequence and secondary structure required for GlcN6P-dependent self-cleavage. However, the high degree of nucleotide conservation offers few clues regarding the higher-order structural organization of the catalytic core. To further investigate core ribozyme structure, minimal ‘consensus-type’ glmS ribozymes that retain GlcN6P-dependent activity were produced. Mutational analyses of consensus-type glmS ribozymes support a model for core ribozyme folding through a pseudoknot structure formed by the interaction of two highly conserved sequence segments. Moreover, GlcN6P-dependent function is demonstrated for bimolecular constructs in which substrate interaction with the ribozyme is minimally comprised of sequence representing that involved in putative pseudoknot formation. These studies suggest that the glmS ribozyme adopts an intricate multi-strand catalytic core through the formation of a pseudoknot structure, and provide a refined model for further considering GlcN6P interaction and GlcN6P-dependent ribozyme function.