The Folding Mechanism of BBL: Plasticity of Transition-State Structure Observed within an Ultrafast Folding Protein Family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Φ-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Φ-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Φ values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Φ values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.     

Downhill versus Barrier-Limited Folding of BBL 1: Energetic and Structural Perturbation Effects upon Protonation of a Histidine of Unusually Low pKa

A dispersion of melting temperatures at pH 5.3 for individual residues of the BBL protein domain has been adduced as evidence for barrier-free downhill folding. Other members of the peripheral subunit domain family fold cooperatively at pH 7. To search for possible causes of anomalies in BBL's denaturation behavior, we measured the pH titration of individual residues by heteronuclear NMR. At 298 K, the pK_a of His142 was close to that of free histidine at 6.47 ± 0.04, while that of the more buried His166 was highly perturbed at 5.39 ± 0.02. Protonation of His166 is thus energetically unfavorable and destabilizes the protein by ∼ 1.5 kcal/mol. Changes in C^α secondary shifts at pH 5.3 showed a decrease in helicity of the C-terminus of helix 2, where His166 is located, which was accompanied by a measured decrease of 1.1 ± 0.2 kcal/mol in stability from pH 7 to 5.3. Protonation of His166 perturbs, therefore, the structure of BBL. Only ∼ 1% of the structurally perturbed state will be present at the biologically relevant pH 7.6. Experiments at pH 5.3 report on a near-equal mixture of the two different native states. Further, at this pH, small changes of pH and pK_a induced by changes in temperature will have near-maximal effects on pH-dependent conformational equilibria and on propagation of experimental error. Accordingly, conventional barrier-limited folding predicts some dispersion of measured thermal unfolding curves of individual residues at pH 5.3.     

Structure, folding and stability of FimA, the main structural subunit of type 1 pili from uropathogenic Escherichia coli strains

Filamentous type 1 pili are responsible for attachment of uropathogenic Escherichia coli strains to host cells. They consist of a linear tip fibrillum and a helical rod formed by up to 3000 copies of the main structural pilus subunit FimA. The subunits in the pilus interact via donor strand complementation, where the incomplete, immunoglobulin-like fold of each subunit is complemented by an N-terminal donor strand of the subsequent subunit. Here, we show that folding of FimA occurs at an extremely slow rate (half-life: 1.6 h) and is catalyzed more than 400-fold by the pilus chaperone FimC. Moreover, FimA is capable of intramolecular self-complementation via its own donor strand, as evidenced by the loss of folding competence upon donor strand deletion. Folded FimA is an assembly-incompetent monomer of low thermodynamic stability (− 10.1 kJ mol^− 1) that can be rescued for pilus assembly at 37 °C because FimC selectively pulls the fraction of unfolded FimA molecules from the FimA folding equilibrium and allows FimA refolding on its surface. Elongation of FimA at the C-terminus by its own donor strand generated a self-complemented variant (FimAa) with alternative folding possibilities that spontaneously adopts the more stable conformation (− 85.0 kJ mol^− 1) in which the C-terminal donor strand is inserted in the opposite orientation relative to that in FimA. The solved NMR structure of FimAa revealed extensive β-sheet hydrogen bonding between the FimA pilin domain and the C-terminal donor strand and provides the basis for reconstruction of an atomic model of the pilus rod.     

Internal bulge and tetraloop of the catalytic domain 5 of a group II intron ribozyme are flexible: implications for catalysis

RNA molecules have an inherent flexibility that enables recognition of other interacting partners through potential disorder-order transitions, yet studies to quantify such motional dynamics remain few. With an increasing database of three-dimensional structures of biologically important RNA molecules, quantifying such motions becomes important to link structural deformations with function. One such system studied intensely is domain 5 (D5) from the self-splicing group II introns, which is at the heart of its catalytic machinery. We report the dynamics of a 36 nucleotide D5 from the Pylaiella littoralis group II intron in the presence and absence of magnesium ions, and at a range of temperatures (298K-318 K). Using high-resolution NMR experiments of heteronuclear nuclear Overhauser enhancement (NOE), spin-lattice (R(1)), and spin-spin (R(2)) (13)C relaxation rates, we determined the rotational diffusion tensor of D5 using the ROTDIF program modified for RNA dynamic analysis (ROTDIF_RNA). The D5 rotational diffusion tensor has an axial symmetric ratio (D(||)/D(perpendicular)) of 1.7+/-0.3, consistent with an estimated overall rotational correlation time of tau(m)=(2D(||)+4D(perpendicular))(-1) of 6.1(+/-0.3) ns at 298 K and 4.1(+/-0.2) ns at 318 K. The measured relaxation data were analyzed with the reduced spectral density mapping formalism using assumed values of the chemical shift anisotropy of the (13)C spins. Both the relaxation data and the values of the spectral density function reveal that the functional groups in D5 implicated in magnesium ion binding and catalysis (catalytic triad, internal bulge, and tetraloop regions) exhibit thermally induced motion on a wide variety of timescales. Because these motions parallel those observed in the intramolecular stem-loop of the U6 element within the spliceosome, we hypothesize that such extensive dynamic disorder likely facilitates D5 engaging both binding and catalytic regions of the ribozyme, and these may be a conserved feature of the catalytic machinery essential for catalysis.     

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.     

Structure and dynamics of pin1 during catalysis by NMR

The link between internal enzyme motions and catalysis is poorly understood. Correlated motions in the microsecond-to-millisecond timescale may be critical for enzyme function. We have characterized the backbone dynamics of the peptidylprolyl isomerase (Pin1) catalytic domain in the free state and during catalysis. Pin1 is a prolyl isomerase of the parvulin family and specifically catalyzes the isomerization of phosphorylated Ser/Thr-Pro peptide bonds. Pin1 has been shown to be essential for cell-cycle progression and to interact with the neuronal tau protein inhibiting its aggregation into fibrillar tangles as found in Alzheimer's disease. (15)N relaxation dispersion measurements performed on Pin1 during catalysis reveal conformational exchange processes in the microsecond timescale. A subset of active site residues undergo kinetically similar exchange processes even in the absence of a substrate, suggesting that this area is already "primed" for catalysis. Furthermore, structural data of the turning-over enzyme were obtained through inter- and intramolecular nuclear Overhauser enhancements. This analysis together with a characterization of the substrate concentration dependence of the conformational exchange allowed the distinguishing of regions of the enzyme active site that are affected primarily by substrate binding versus substrate isomerization. Together these data suggest a model for the reaction trajectory of Pin1 catalysis.     

NMR Structure of a Monomeric Intermediate on the Evolutionarily Optimized Assembly Pathway of a Small Trimerization Domain

Efficient formation of specific intermolecular interactions is essential for self-assembly of biological structures. The foldon domain is an evolutionarily optimized trimerization module required for assembly of the large, trimeric structural protein fibritin from phage T4. Monomers consisting of the 27 amino acids comprising a single foldon domain subunit spontaneously form a natively folded trimer. During assembly of the foldon domain, a monomeric intermediate is formed on the submillisecond time scale, which provides the basis for two consecutive very fast association reactions. Mutation of an intermolecular salt bridge leads to a monomeric protein that resembles the kinetic intermediate in its spectroscopic properties. NMR spectroscopy revealed essentially native topology of the monomeric intermediate with defined hydrogen bonds and side-chain interactions but largely reduced stability compared to the native trimer. This structural preorganization leads to an asymmetric charge distribution on the surface that can direct rapid subunit recognition. The low stability of the intermediate allows a large free-energy gain upon trimerization, which serves as driving force for rapid assembly. These results indicate different free-energy landscapes for folding of small oligomeric proteins compared to monomeric proteins, which typically avoid the transient population of intermediates.     

A Method for Helical RNA Global Structure Determination in Solution Using Small-Angle X-Ray Scattering and NMR Measurements

We report a “top-down” method that uses mainly duplexes' global orientations and overall molecular dimension and shape restraints, which were extracted from experimental NMR and small-angle X-ray scattering data, respectively, to determine global architectures of RNA molecules consisting of mostly A-form-like duplexes. The method is implemented in the G2G (from global measurement to global structure) toolkit of programs. We demonstrate the efficiency and accuracy of the method by determining the global structure of a 71-nt RNA using experimental data. The backbone root-mean-square deviation of the ensemble of the calculated global structures relative to the X-ray crystal structure is 3.0 ± 0.3 Å using the experimental data and is only 2.5 ± 0.2 Å for the three duplexes that were orientation restrained during the calculation. The global structure simplifies interpretation of multidimensional nuclear Overhauser spectra for high-resolution structure determination. The potential general application of the method for RNA structure determination is discussed.     

NMR structure of the Escherichia coli type 1 pilus subunit FimF and its interactions with other pilus subunits

Type 1 pili from uropathogenic Escherichia coli strains mediate bacterial attachment to target receptors on the host tissue. They are composed of up to 3000 copies of the subunit FimA, which form the stiff, helical pilus rod, and the subunits FimF, FimG, and FimH, which form the linear tip fibrillum. All subunits in the pilus interact via donor strand complementation, in which the incomplete immunoglobulin-like fold of each subunit is complemented by insertion of an N-terminal extension from the following subunit. We determined the NMR structure of a monomeric, self-complemented variant of FimF, FimF_F, which has a second FimF donor strand segment fused to its C-terminus that enables intramolecular complementation of the FimF fold. NMR studies on bimolecular complexes between FimF_F and donor strand-depleted variants of FimF and FimG revealed that the relative orientations of neighboring domains in the tip fibrillum cover a wide range. The data provide strong support for the intrinsic flexibility of the tip fibrillum. They lend further support to the hypothesis that this flexibility would significantly increase the probability that the adhesin at the distal end of the fibrillum successfully targets host cell receptors.     

A Small-Molecule Probe Induces a Conformation in HIV TAR RNA Capable of Binding Drug-Like Fragments

The HIV-1 transactivation response (TAR) element-Tat interaction is a potentially valuable target for treating HIV infection, but efforts to develop TAR-binding antiviral drugs have not yet yielded a successful candidate for clinical development. In this work, we describe a novel approach toward screening fragments against RNA that uses a chemical probe to target the Tat-binding region of TAR. This probe fulfills two critical roles in the screen: by locking the RNA into a conformation capable of binding other fragments, it simultaneously allows the identification of proximal binding fragments by ligand-based NMR. Using this approach, we have discovered six novel TAR-binding fragments, three of which were docked relative to the probe-RNA structure using experimental NMR restraints. The consistent orientations of functional groups in our data-driven docked structures and common electrostatic properties across all fragment leads reveal a surprising level of selectivity by our fragment-sized screening hits. These models further suggest linking strategies for the development of higher-affinity lead compounds for the inhibition of the TAR-Tat interaction.     

The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Prokaryotic Ubiquitin-Like Protein Pup Is Intrinsically Disordered

The prokaryotic ubiquitin-like protein Pup targets substrates for degradation by the Mycobacterium tuberculosis proteasome through its interaction with Mpa, an ATPase that is thought to abut the 20S catalytic subunit. Ubiquitin, which is assembled into a polymer to similarly signal for proteasomal degradation in eukaryotes, adopts a stable and compact structural fold that is adapted into other proteins for diverse biological functions. We used NMR spectroscopy to demonstrate that, unlike ubiquitin, the 64-amino-acid protein Pup is intrinsically disordered with small helical propensity in the C-terminal region. We found that the Pup:Mpa interaction involves an extensive contact surface that spans S21-K61 and that the binding is in the “slow exchange” regime on the NMR time scale, thus demonstrating higher affinity than most ubiquitin:ubiquitin receptor pairs. Interestingly, during the titration experiment, intermediate Pup species were observable, suggesting the formation of one or more transient state(s) upon binding. Moreover, Mpa selected one configuration for a region undergoing chemical exchange in the free protein. These findings provide mechanistic insights into Pup's functional role as a degradation signal.     

Energetic Communication between Functional Sites of the Gene-3-Protein during Infection by Phage fd

To initiate infection of Escherichia coli, phage fd uses its gene-3-protein (G3P) to bind first to an F pilus and then to the TolA protein at the cell surface. G3P is normally auto-inhibited because a tight interaction between the two N-terminal domains N1 and N2 buries the TolA binding site. Binding of N2 to the pilus activates G3P by initiating long-range conformational changes that are relayed to the domain interface and to a proline timer. We discovered that the 23–28 loop of the N1 domain is critical for propagating these conformational signals. The analysis of the stability and the folding dynamics of G3P variants with a shortened loop combined with TolA interaction studies and phage infection experiments reveal how the contact between the N2 domain and the 23–28 loop of N1 is energetically linked with the interdomain region and the proline timer and how it affects phage infectivity. Our results illustrate how conformational transitions and prolyl cis/trans isomerization can be coupled energetically and how conformational signals to and from prolines can be propagated over long distances in proteins.     

The Effects of CapZ Peptide (TRTK-12) Binding to S100B-Ca as Examined by NMR and X-ray Crystallography

Structure-based drug design is underway to inhibit the S100B-p53 interaction as a strategy for treating malignant melanoma. X-ray crystallography was used here to characterize an interaction between Ca²⁺-S100B and TRTK-12, a target that binds to the p53-binding site on S100B. The structures of Ca²⁺-S100B (1.5-Å resolution) and S100B-Ca²⁺-TRTK-12 (2.0-Å resolution) determined here indicate that the S100B-Ca²⁺-TRTK-12 complex is dominated by an interaction between Trp7 of TRTK-12 and a hydrophobic binding pocket exposed on Ca²⁺-S100B involving residues in helices 2 and 3 and loop 2. As with an S100B-Ca²⁺-p53 peptide complex, TRTK-12 binding to Ca²⁺-S100B was found to increase the protein's Ca²⁺-binding affinity. One explanation for this effect was that peptide binding introduced a structural change that increased the number of Ca²⁺ ligands and/or improved the Ca²⁺ coordination geometry of S100B. This possibility was ruled out when the structures of S100B-Ca²⁺-TRTK-12 and S100B-Ca²⁺ were compared and calcium ion coordination by the protein was found to be nearly identical in both EF-hand calcium-binding domains (RMSD = 0.19). On the other hand, B-factors for residues in EF2 of Ca²⁺-S100B were found to be significantly lowered with TRTK-12 bound. This result is consistent with NMR ¹⁵N relaxation studies that showed that TRTK-12 binding eliminated dynamic properties observed in Ca²⁺-S100B. Such a loss of protein motion may also provide an explanation for how calcium-ion-binding affinity is increased upon binding a target. Lastly, it follows that any small-molecule inhibitor bound to Ca²⁺-S100B would also have to cause an increase in calcium-ion-binding affinity to be effective therapeutically inside a cell, so these data need to be considered in future drug design studies involving S100B.     

Probing the flexibility of the DsbA oxidoreductase from Vibrio cholerae--a 15N - 1H heteronuclear NMR relaxation analysis of oxidized and reduced forms of DsbA

We have determined the structure of the reduced form of the DsbA oxidoreductase from Vibrio cholerae. The reduced structure shows a high level of similarity to the crystal structure of the oxidized form and is typical of this class of enzyme containing a thioredoxin domain with an inserted alpha-helical domain. Proteolytic and thermal stability measurements show that the reduced form of DsbA is considerably more stable than the oxidized form. NMR relaxation data have been collected and analyzed using a model-free approach to probe the dynamics of the reduced and oxidized states of DsbA. Akaike's information criteria have been applied both in the selection of the model-free models and the diffusion tensors that describe the global motions of each redox form. Analysis of the dynamics reveals that the oxidized protein shows increased disorder on the pico- to nanosecond and micro- to millisecond timescale. Many significant changes in dynamics are located either close to the active site or at the insertion points between the domains. In addition, analysis of the diffusion data shows there is a clear difference in the degree of interdomain movement between oxidized and reduced DsbA with the oxidized form being the more rigid. Principal components analysis has been employed to indicate possible concerted movements in the DsbA structure, which suggests that the modeled interdomain motions affect the catalytic cleft of the enzyme. Taken together, these data provide compelling evidence of a role for dynamics in the catalytic cycle of DsbA.     

Downhill versus Barrier-Limited Folding of BBL 2: Mechanistic Insights from Kinetics of Folding Monitored by Independent Tryptophan Probes

Barrier-free downhill folding has been proposed for the peripheral subunit-binding domain BBL. To date, ultrafast kinetic experiments on BBL, which are crucial for a mechanistic understanding of folding, have been hampered by the lack of good intrinsic spectroscopic probes. Here, we present a detailed kinetic characterization of three single-point tryptophan mutants of BBL that have suitable fluorescence properties for following microsecond and nanosecond folding kinetics using temperature jump fluorescence spectroscopy. Experiments were performed at pH 7, which is optimal for stability and minimizes complications that arise from the presence of an alternative native-state conformation of BBL at lower pH. We examined the dependence of rate and equilibrium constants on concentration of denaturant and found that they follow well-established laws allowing kinetic transients to be related to events in folding and compared with equilibrium data. Logarithms of rate constants versus denaturant concentration yielded plots (chevrons) that are characteristic of barrier-limited folding for all mutants investigated, including a truncated sequence that was previously used in the proposal of downhill folding. The thermodynamic quantities calculated from the rate constants were in excellent agreement with those directly determined from equilibrium denaturation based on empirical two-state equations. We found that sequence truncation of BBL as used in studies proposing downhill folding leads to a large loss in helical content and protein stability, which were exacerbated at the low pH used in those studies. The kinetics and equilibria of folding of BBL fit to conventional barrier-limited kinetics.