Solution NMR chemical shift assignment of apo and molybdate-bound ModA at two pHs.

ModA is a soluble periplasmic molybdate-binding protein found in most gram-negative bacteria. It is part of the ABC transporter complex ModABC that moves molybdenum into the cytoplasm, to be used by enzymes that carry out various redox reactions. Since there is no clear analog for ModA in humans, this protein could be a good target for antibacterial drug design. Backbone 1H, 13C and 15N chemical shifts of apo and molybdate-bound ModA from E. coli were assigned at pHs 6.0 and 4.5. In addition, side chain atoms were assigned for apo ModA at pH 6.0. When comparing apo and molybdate-bound ModA at pH 6.0, large chemical shift perturbations are observed, not only in areas near the bound metal, but also in regions that are distant from the metal-binding site. Given the significant conformational change between apo and holo ModA, we might expect the large chemical shift changes to be more widespread; however, since they are limited to specific regions, the residues with large perturbations may reveal allosteric sites that could ultimately be important for the design of antibiotics that target ModA.

Removal of disulfide from acid stress chaperone HdeA does not wholly eliminate structure or function at low pH.

HdeA is an acid-stress chaperone that operates in the periplasm of various strains of pathogenic gram-negative bacteria. Its primary function is to prevent irreversible aggregation of other periplasmic proteins when the bacteria enter the acidic environment of the stomach after contaminated food is ingested; its role is therefore to help the bacteria survive long enough to enter and colonize the intestines. The mechanism of operation of HdeA is unusual in that this helical homodimer is inactive when folded at neutral pH but becomes activated at low pH after the dimer dissociates and partially unfolds. Studies with chemical reducing agents previously suggested that the intramolecular disulfide bond is important for maintaining residual structure in HdeA at low pH and may be responsible for positioning exposed hydrophobic residues together for the purpose of binding unfolded client proteins. In order to explore its role in HdeA structure and chaperone function we performed a conservative cysteine to serine mutation of the disulfide. We found that, although residual structure is greatly diminished at pH 2 without the disulfide, it is not completely lost; conversely, the mutant is almost completely random coil at pH 6. Aggregation assays showed that mutated HdeA, although less successful as a chaperone than wild type, still maintains a surprising level of function. These studies highlight that we still have much to learn about the factors that stabilize residual structure at low pH and the role of disulfide bonds.

Detection of key sites of dimer dissociation and unfolding initiation during activation of acid-stress chaperone HdeA at low pH.

HdeA is a small acid-stress chaperone protein with a unique activity profile. At physiological pH, it forms a folded, but inactive, dimer. Below pH 3.0, HdeA unfolds and dissociates into disordered monomers, utilizing exposed hydrophobic patches to bind other unfolded proteins and prevent their irreversible aggregation. In this way, HdeA has a key role in helping pathogenic bacteria survive our acidic stomach and colonize our intestines, facilitating the spread of dysentery. Despite numerous publications on the topic, there remain questions about the mechanism by which HdeA unfolding and activation are triggered. Previous studies usually assessed HdeA unfolding over pH increments that are too far apart to gain fine detail of the process of unfolding and dimer dissociation, and often employed techniques that prevented thorough evaluation of specific regions of the protein. We used a variety of heteronuclear NMR experiments to investigate changes to backbone and side chain structure and dynamics of HdeA at four pHs between 3.0 and 2.0. We found that the long loop in the dimer interface is an early site of initiation of dimer dissociation, and that a molecular "clasp" near the disulfide bond is broken open at low pH as part, or as a trigger, of unfolding; this process also results in the separation of C-terminal helices and exposure of key hydrophobic client binding sites. Our results highlight important regions of HdeA that may have previously been overlooked because they lie too close to the disulfide bond or are thought to be too dynamic in the folded state to influence unfolding processes.

The complex role of the N-terminus and acidic residues of HdeA as pH-dependent switches in its chaperone function.

HdeA is a small acid-stress chaperone protein found in the periplasm of several pathogenic gram-negative bacteria. In neutral pH environments HdeA is an inactive folded homodimer but when exposed to strong acidic environments it partially unfolds and, once activated, binds to other periplasmic proteins, protecting them from irreversible aggregation. Here we use a combination of hydrogen/deuterium exchange NMR experiments and constant pH molecular dynamics simulations to elucidate the role of HdeA's N-terminus in its activation mechanism. Previous work indicates that the N-terminus is flexible and unprotected at high pH while exhibiting interactions with some HdeA client binding site residues. It, however, becomes partially solvent-protected at pH 2.6 - 2.8 and then loses protection again at pH 2.0. This protection is not due to the appearance of new secondary structure, but rather increased contacts between N-terminal residues and the C-terminus of the other protomer in the dimer, as well as concurrent loosening of its hold on the client binding site residues, priming HdeA for interactions with periplasmic client proteins. This work also uncovers unusual protonation profiles of some titratable residues and suggests their complex role in chaperone function.

Backbone and side chain chemical shift assignments of apolipophorin III from Galleria mellonella.

Apolipophorin III, a 163 residue monomeric protein from the greater wax moth Galleria mellonella (abbreviated as apoLp-IIIGM), has roles in upregulating expression of antimicrobial proteins as well as binding and deforming bacterial membranes. Due to its similarity to vertebrate apolipoproteins there is interest in performing atomic resolution analysis of apoLp-IIIGM as part of an effort to better understand its mechanism of action in innate immunity. In the first step towards structural characterization of apoLp-IIIGM, 99 % of backbone and 88 % of side chain (1)H, (13)C and (15)N chemical shifts were assigned. TALOS+ analysis of the backbone resonances has predicted that the protein is composed of five long helices, which is consistent with the reported structures of apolipophorins from other insect species. The next stage in the characterization of apoLp-III from G. mellonella will be to utilize these resonance assignments in solving the solution structure of this protein.

NMR-monitored titration of acid-stress bacterial chaperone HdeA reveals that Asp and Glu charge neutralization produces a loosened dimer structure in preparation for protein unfolding and chaperone activation.

HdeA is a periplasmic chaperone found in several gram-negative pathogenic bacteria that are linked to millions of cases of dysentery per year worldwide. After the protein becomes activated at low pH, it can bind to other periplasmic proteins, protecting them from aggregation when the bacteria travel through the stomach on their way to colonize the intestines. It has been argued that one of the major driving forces for HdeA activation is the protonation of aspartate and glutamate side chains. The goal for this study, therefore, was to investigate, at the atomic level, the structural impact of this charge neutralization on HdeA during the transition from near-neutral conditions to pH 3.0, in preparation for unfolding and activation of its chaperone capabilities. NMR spectroscopy was used to measure pKa values of Asp and Glu residues and monitor chemical shift changes. Measurements of R2/R1 ratios from relaxation experiments confirm that the protein maintains its dimer structure between pH 6.0 and 3.0. However, calculated correlation times and changes in amide protection from hydrogen/deuterium exchange experiments provide evidence for a loosening of the tertiary and quaternary structures of HdeA; in particular, the data indicate that the dimer structure becomes progressively weakened as the pH decreases. Taken together, these results provide insight into the process by which HdeA is primed to unfold and carry out its chaperone duties below pH 3.0, and it also demonstrates that neutralization of aspartate and glutamate residues is not likely to be the sole trigger for HdeA dissociation and unfolding.

Expression, purification and preliminary NMR characterization of isotopically labeled wild-type human heterotrimeric G protein αi1.

Molecular-level investigation of proteins is increasingly important to researchers trying to understand the mechanisms of signal transmission. Heterotrimeric G proteins control the activation of many critical signal transmission cascades and are also implicated in numerous diseases. As part of a longer-term investigation of intramolecular motions in RGS and Gα proteins in their apo and complexed forms, we have successfully developed a protocol for preparing milligram quantities of highly purified, isotopically labeled wild-type human Gα(i1) (hGα(i1)) subunit for NMR studies. High levels of expression in Escherichia coli can be attributed to the use of the SUMO fusion protein system, a bacterial strain that produces rare codons, supplementation of minimal medium with small quantities of isotopically labeled rich medium and a lowered induction temperature. Purification of hGα(i1) utilized affinity and size exclusion chromatography, and protein activity was confirmed using fluorescence-based GTP-binding studies. Preliminary NMR analysis of hGα(i1) has shown that high-quality spectra can be obtained at near-physiological temperatures, whereas lower temperature spectra display numerous weak and broadened peaks, providing preliminary evidence for widespread µs-ms timescale exchange. In an effort to further optimize the NMR spectra we prepared a truncated form of hGα(i1) (hGα(i1)-Δ31) in which the 31-residue unstructured N-terminus was removed. This resulted in further improvements in spectral quality by eliminating high-intensity peaks that obscured resonances from structured segments of the protein. We plan to use hGα(i1)-Δ31 in future investigations of protein dynamics by NMR spectroscopy to gain insight into the role of these motions in RGS/Gα binding selectivity.

NMR-detected conformational exchange observed in a computationally designed variant of protein Gbeta1.

Detailed biophysical characterization of computationally designed proteins has become increasingly important in order to thoroughly understand the properties of these variants compared with wild-type and to apply this knowledge to future designs. The protein dynamics and structural properties of a computationally designed variant (Delta1.5) of the beta1 domain of streptococcal protein G (Gbeta1) were measured using multinuclear NMR methods. Results from relaxation, diffusion and hydrogen exchange experiments indicate that the variant weakly self-associates at NMR concentrations, with evidence for multiple binding sites. Although comparison of fast (ps-ns) timescale motions shows only small differences in dynamics between Delta1.5 and wild-type, results from the measurement of intermediate (mus-ms) timescale motions are very different. Significant backbone conformational exchange has been observed in the variant at positions all along the sequence, whereas the wild-type Gbeta1 shows little evidence for this type of motion. This increased conformational exchange in Delta1.5 has been attributed to core overpacking resulting from the incorporation of two large hydrophobic side chains and the loss of an aromatic T-stacking interaction. These data highlight, in detail, the potential consequences of incorporating major perturbations in the core of a protein and the need to carry out more detailed analyses of the biophysical properties of designed proteins in order to better understand and predict the effects of mutations.

Full-sequence computational design and solution structure of a thermostable protein variant.

Computational protein design procedures were applied to the redesign of the entire sequence of a 51 amino acid residue protein, Drosophila melanogaster engrailed homeodomain. Various sequence optimization algorithms were compared and two resulting designed sequences were experimentally evaluated. The two sequences differ by 11 mutations and share 22% and 24% sequence identity with the wild-type protein. Both computationally designed proteins were considerably more stable than the naturally occurring protein, with midpoints of thermal denaturation greater than 99 degrees C. The solution structure was determined for one of the two sequences using multidimensional heteronuclear NMR spectroscopy, and the structure was found to closely match the original design template scaffold.

Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure.

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.

Aromatic and methyl NOEs highlight hydrophobic clustering in the unfolded state of an SH3 domain.

The N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) exists in equilibrium between a folded (F(exch)) state and a relatively compact unfolded (U(exch)) state under nondenaturing conditions. Selectively labeled samples of the domain have been analyzed by NOESY NMR experiments to probe residual hydrophobic clustering in the U(exch) state. The labeling strategy included selective protonation of aromatic rings or delta-methyl groups on Ile and Leu residues in a highly deuterated background. Combined with long mixing times, the methods permitted observation of significant numbers of long-range interactions between hydrophobic side chains, providing evidence for multiple conformers involving non-native hydrophobic clusters around the Trp 36 indole. Comparison of these data with previously reported HN-HN NOEs yields structural insight into the diversity of structures within the U(exch) ensemble in the drkN SH3 domain. Many of the HN-HN NOEs are consistent with models containing compact residual nativelike secondary structure and greater exposure of the Trp 36 indole to solvent, similar to kinetic intermediates formed in the hierarchic condensation model of folding. However, the methyl and aromatic NOE data better fit conformations with non-native burial of the Trp indole surrounded by hydrophobic groups and more loosely formed beta-structure; these structural characteristics are more consistent with those of kinetic intermediates formed during the hydrophobic collapse mechanism of folding. This suite of NOE data provides a more complete picture of the structures that span the U(exch) state ensemble, from conformers with non-native structure but long-range contacts to those that are highly nativelike. Together, the results are also consistent with the folding funnel view involving multiple folding pathways for this molecule.

Site-specific contributions to the pH dependence of protein stability.

Understanding protein stability is a significant challenge requiring characterization of interactions within both folded and unfolded states. Of these, electrostatic interactions influence ionization equilibria of acidic and basic groups and diversify their pK(a) values. The pH dependence of the thermodynamic stability (Delta G(FU)) of a protein arises as a consequence of differential pK(a) values between folded and unfolded states. Previous attempts to calculate pH-dependent contributions to stability have been limited by the lack of experimental unfolded state pK(a) values. Using recently developed NMR spectroscopic methods, we have determined residue-specific pK(a) values for a thermodynamically unstable Src homology 3 domain in both states, enabling the calculation of the pH dependence of stability based on simple analytical expressions. The calculated pH stability profile obtained agrees very well with experiment, unlike profiles derived from two current models of electrostatic interactions within unfolded states. Most importantly, per-residue contributions to the pH dependence of Delta G(FU) derived from the data provide insights into specific electrostatic interactions in both the folded and unfolded states and their roles in protein stability. These interactions include a hydrogen bond between the Asp-8 side-chain and the Lys-21 backbone amide group in the folded state, which represents a highly conserved interaction in Src homology 3 domains.

Cooperative interactions and a non-native buried Trp in the unfolded state of an SH3 domain.

The presence of residual structure in the unfolded state of the N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) has been investigated using far- and near-UV circular dichroism (CD), fluorescence, and NMR spectroscopy. The unfolded (U(exch)) state of the drkN SH3 domain is significantly populated and exists in equilibrium with the folded (F(exch)) state under non-denaturing conditions near physiological pH. Denaturation experiments have been performed on the drkN SH3 domain in order to monitor the change in ellipticity, fluorescence intensity, and chemical shift between the U(exch) state and chemically or thermally denatured states. Differences between the unfolded and chemically or thermally denatured states highlight specific areas of residual structure in the unfolded state that are cooperatively disrupted upon denaturation. Results provide evidence for cooperative interactions in the unfolded state involving residues of the central beta-sheet, particularly the beta4 strand. Denaturation as well as hydrogen-exchange experiments demonstrate a non-native burial of the Trp ring within this "cooperative" core of the unfolded state. These findings support the presence of non-native hydrophobic clusters, organised by Trp rings, within disordered states.

Distribution of molecular size within an unfolded state ensemble using small-angle X-ray scattering and pulse field gradient NMR techniques.

The size distribution of molecules within an unfolded state of the N-terminal SH3 domain of drk (drkN SH3) has been studied by small-angle X-ray scattering (SAXS) and pulsed-field-gradient NMR (PFG-NMR) methods. An empirical model to describe this distribution in the unfolded state ensemble has been proposed based on (i) the ensemble-averaged radius of gyration and hydrodynamic radius derived from the SAXS and PFG-NMR data, respectively, and (ii) a histogram of the size distribution of structures obtained from preliminary analyses of structural parameters recorded on the unfolded state. Results show that this unfolded state, U(exch), which exists in equilibrium with the folded state, F(exch), under non-denaturing conditions, is relatively compact, with the average size of conformers within the unfolded state ensemble only 30-40% larger than the folded state structure. In addition, the model predicts a significant overlap in the size range of structures comprising the U(exch) state with those in a denatured state obtained by addition of 2 M guanidinium chloride.