Sedimentation equilibrium (SE) can be used to study protein-protein interactions in a physiological environment. This manuscript describes the use of this technique to determine the effect of pH on the stability of a homo-pentamer formed by the small hydrophobic (SH) protein encoded by the human syncytial respiratory virus (hRSV).
Analytical ultracentrifugation (AUC) can be used to study reversible interactions between macromolecules over a wide range of interaction strengths and under physiological conditions. This makes AUC a method of choice to quantitatively assess stoichiometry and thermodynamics of homo- and hetero-association that are transient and reversible in biochemical processes. In the modality of sedimentation equilibrium (SE), a balance between diffusion and sedimentation provides a profile as a function of radial distance that depends on a specific association model. Herein, a detailed SE protocol is described to determine the size and monomer-monomer association energy of a small membrane protein oligomer using an analytical ultracentrifuge. AUC-ES is label-free, only based on physical principles, and can be used on both water soluble and membrane proteins. An example is shown of the latter, the small hydrophobic (SH) protein in the human respiratory syncytial virus (hRSV), a 65-amino acid polypeptide with a single α-helical transmembrane (TM) domain that forms pentameric ion channels. NMR-based structural data shows that SH protein has two protonatable His residues in its transmembrane domain that are oriented facing the lumen of the channel. SE experiments have been designed to determine how pH affects association constant and the oligomeric size of SH protein. While the pentameric form was preserved in all cases, its association constant was reduced at low pH. These data are in agreement with a similar pH dependency observed for SH channel activity, consistent with a lumenal orientation of the two His residues in SH protein. The latter may experience electrostatic repulsion and reduced oligomer stability at low pH. In summary, this method is applicable whenever quantitative information on subtle protein-protein association changes in physiological conditions have to be measured.
Analytical ultracentrifugation1-5 is one of the most important methods to study interactions of macromolecules under physiological conditions, being accessible to both weak and strong interactions. The method is label-free and uses light absorption or interference, and even fluorescence optical systems can be used to access concentration ranges over several orders of magnitude6.
This method is especially useful since most biochemical processes depend on reversible interactions. The stoichiometry and strength of these interactions have to be quantitatively characterized to understand biological processes, and a number of methods exist for this purpose 7,8. However, transient interactions are difficult to study 9.
The choice of a method to characterize macromolecular interactions depends on its static or dynamic nature. In the first case, sedimentation velocity (SV) is used, where the rate of radial transport is measured and complexes are fractionated on the basis of differences in buoyant mass and shape.
In contrast, dynamic associations that are reversible on the time scale of the experiment cannot be physically separated. In this case, self- or hetero-interactions leading to non-covalent interactions are in an equilibrium that depends on the total protein concentration. These dynamic interactions can be studied by both sedimentation equilibrium (SE) and sedimentation velocity (SV) 10. However, the first method is simpler to perform and is described here. In SE, centrifugation is performed at a sufficiently low speed so that an equilibrium is reached between diffusion and sedimentation. At this point, the equilibrium profile of an optical signal (UV-VIS) as a function of radial distance, can be analyzed using pre-set thermodynamic models for associations11.
In the present paper, a sedimentation equilibrium study is presented of the self-association of a viral membrane protein that forms ion channels. Because of its hydrophobicity, the experiment is run in presence of detergent, and in this case the density of solvent has to be matched to that of the detergent. However, the protocol described would identical in the case of a water soluble protein, except that no solvent density matching would be required.
The protein used is encoded in the human respiratory syncytial virus (hRSV), an enveloped pneumovirus in the paramyxoviridae family that causes lower respiratory tract disease in infants, elderly and immunocompromised populations worldwide 12. Up to 64 million reported cases of hRSV infection and 160,000 deaths occur each year.
The hRSV genome transcribes 11 proteins, including the three membrane proteins F, G, and small hydrophobic (SH). SH protein is involved in the pathogenesis of RSV infection. RSV lacking the SH gene (RSVΔSH) was viable, caused formation of syncytia and grew as well as the wild-type (WT) virus 13-16. However, RSVΔSH virus replicated 10-fold less efficiently than the WT in the upper respiratory tract 15,16. Also, RSVΔSH virus was attenuated in in vivo mouse and chimpanzee models 13,17.
The SH protein is a 64 (RSV subgroup A) or 65 (RSV subgroup B) amino acids long type II integral membrane protein that accumulates mostly at the membranes of the Golgi compartment 18. SH protein has a single predicted a-helical transmembrane (TM) domain 19 which is highly conserved 20,21. The C- and N-terminal extramembrane domains are oriented lumenally/extracellularly and cytoplasmically, respectively.
Both synthetic TM domain (residues 18-43) and full length SH protein have been shown to form homopentamers in a variety of detergents. The homopentameric form is responsible for channel activity in planar lipid bilayers 22,23. The correct orientation of the TM monomers in the lipid bilayer was first determined using site specific infrared dichroism 23, which showed His-22 to be in a lumenal, close to inter-helical, orientation. The same TM domain orientation was confirmed by NMR studies that reconstructed the pentameric a-helical bundle of the full-length protein in dodecylphosphocholine (DPC) micelles 22. In this ‘micelle’ model, a single a- helical TM domain was flanked N-terminally by an a-helix, and C-terminally by an extended b-hairpin. The two protonatable residues of SH protein, His-22 and His-51, are located in the TM domain (lumenally oriented), and at the tip of the extramembrane C-terminal β hairpin (far from the channel pore), respectively. In a bicellar environment, however, the TM α-helix extends up to His-51, and both His residues are accessible to the lumen of the channel24. The channel structure adopts a funnel-like architecture 22, where the narrower region (Ser-29 to Cys-45) 22 is lined with hydrophobic side chains (Ile-32, Ile-36, Ile-40 and Leu-44), and Ile-36 defines the narrowest point in the channel lumen. His-22 is located at the largest opening of this funnel, whereas His-51 is at the tip of the smallest opening.
In the present paper, analytical centrifugation in a sedimentation equilibrium mode has been used to determine if His protonation affects the stability of the SH protein pentamer. In this case, SH protein was solubilized in C14-betaine detergent, which has been used previously to show that SH protein forms pentameric oligomers22.
This protocol is based on the following resources, which are to be referred for more details and special considerations3,25-28.
1. Density matching of detergent micelles with 2H2O
Note: The density of the buffer solution needs to be matched to the density of the detergent micelles. Common density-adjusting agents include 2H2O, H218O, 2H218O, glycerol and sucrose29. H218O has the same density as 2H2O and may be a better choice if deuteration of exchangeable protons in the protein is not desired. In this procedure, the density of 3-(N,N-dimethylmyristylammonio)-propanesulfonate (C14SB) detergent in 50 mM Tris pH 7.3, 100 mM NaCl will be matched with 2H2O. As an initial guess the following concentrations of 2H2O will be used: 10, 30, and 50% v/v.
1.1. Sample preparation
1.2. Assembly of 6-channel AUC cells and sample loading into the cells.
Note: There are two types of AUC cell depending on the sample loading method. Cells without external fill has to be loaded prior to sealing the cell, whereas external-fill cells can be loaded after the cells are sealed. Assembly of an external-fill AUC cell has been described previously3. In this protocol, the assembly of a 6-channel AUC cell without external fill is described. The main difference is that it has screw rings on both sides which need to be tightened separately, and it doesn’t need housing plugs (Fig. 1). The difference in assembly steps are highlighted below.
Figure 1. Exploded view of a 6-channel AUC cell without external fill. This figure has been modified from Beckman Coulter An-50 Ti and An-60 Ti Analytical Rotor, Cells, and Counterbalance user manual.
Figure 2. Exploded view of the window assembly. This figure has been modified from Beckman Coulter An-50 Ti and An-60 Ti Analytical Rotor, Cells, and Counterbalance user manual
1.3. Setting up interference measurement
Figure 3. Result from HeteroAnalysis Match function. The Match function can be used to monitor equilibrium progress by comparing RMSD between successive scans and the last scan. This example shows attainment of equilibrium after 8 hr as indicated by RMSD values asymptotic to X-axis.
1.4. Data analysis
2. Sedimentation equilibrium of SH in C14SB micelles
2.1. Run parameters
2.2. Sample preparations
2.3. Setting up absorbance measurements
Figure 4. Results from HeteroAnalysis Match function. The first and second speeds (top left and right) appear to have reached equilibrium, but it is better to wait a few more hours to be sure. In comparison, the third and fourth speeds (bottom left and right) have clearly reached equilibrium in a shorter time. Please click here to view a larger version of this figure.
2.4. Data analysis in SEDFIT and SEDPHAT
Note: For further details and considerations in data analysis reader is referred to the following website: www.analyticalultracentrifugation.com.
Figure 5. An example on how to fill in Experimental Parameters.
Figure 6. An example on how to fill in Global Parameters for Monomer-n-Mer Self-Association model.
The radial distribution profile of C14SB detergent micelles in 50 mM Tris, 100 mM NaCl pH 7.3 forms a very shallow exponential that could be fitted to a linear model (Figure 7A). The slope of this distribution is inversely correlated to D2O concentration (Figure 7B). The point where the slope is zero, i.e., the matching D2O concentration, was found to be 32.3%.
See Figure 7 Below.
The same experiment was repeated for different buffers: 50 mM sodium phosphate, 100 mM NaCl pH 5.5 and 50 mM phosphate-citrate, 100 mM NaCl pH 3 to obtain matching D2O concentrations of 30.3% and 41.0%, respectively.
Samples of SH wild type (wt) in detergent were exposed to pH 3, 5.5 and 7.3 (total 6 samples), followed by centrifugation at 15,000, 19,000, 23,000, 28,000, 34,000 and 42,000 rpm. Data obtained at lower speed could not be reliably fitted, possibly because equilibrium had not been attained, therefore data obtained from the four higher speeds were used. At pH 7.3, SH WT was found to form pentamers (Figure 8 and Table 1) with apparent log Ka = 21.35 (Table 2). The association constant did not change at pH 5.5, but it was significantly reduced at pH 3. This is consistent with previous reports24, which reported a decrease of conductance at lower pH, with a pK of 4.5.
See Figure 8 Below.
See Table 1 Below.
Figure 7. Density matching of C14SB with D2O. (A) Radial distribution profile formed by C14SB detergent micelles in 50 mM Tris, 100 mM NaCl pH 7.3 with 25, 30, and 35% D2O. The data was separately fitted to a linear model (black). (B) The slopes were plotted against D2O concentration (black squares) and fitted to a linear model (red line). The matching D2O concentration is indicated (red arrow).
Figure 8. Fitting of SH in C14SB to monomer-pentamer model. Radial distribution profile of SH in C14SB at pH 7 (open circles) was found to fit best to the monomer-pentamer self-association model (black solid line). Fitting residual is shown below. Please click here to view a larger version of this figure.
Model (n-mer) | Chi-square (fixed MW) | Chi-square (fitted MW) | Fitted MW |
3 | 15.444 | 1.0492 | 13477 Da |
4 | 3.8094 | 1.0469 | 9889 Da |
5 | 1.0499 | 1.0497 | 7822 Da |
6 | 2.5994 | 1.547 | 6504 Da |
7 | 6.1667 | 1.2112 | 5743 Da |
Table 1. Comparison of global reduced chi-square values of different monomer-n-mer models.
pH | Lower limit (1σ) | Log Ka | Upper limit (1σ) |
3 | 17.432 | 17.576 | 17.737 |
5.5 | 20.064 | 20.419 | 20.839 |
7.3 | 20.687 | 21.052 | 21.492 |
Table 2. Comparison of apparent log Ka values at different pH.
This paper provides an experimental protocol for sample preparation and analysis of oligomerization of a small membrane protein in detergent using equilibrium sedimentation. The protocol described is equally valid –and simpler- for soluble proteins, as the density matching step is not required. Indeed, the system is constituted by a mixture of detergent and protein. To conduct sedimentation studies, the detergent must be invisible to the gravitational field so that it does not contribute to the particle flotation. Thus, the density of the detergent has to be carefully matched by addition of D2O to the buffer, with the limitation that if the detergent used is too dense, e.g., SDS, not even 100 % D2O can match it. The density matching step is not required when working with a water soluble protein, since the sample does not have detergent.
In our application, a small viral protein that forms ion channels has been used. Like all viroporins, SH protein has an α-helical transmembrane domain, and therefore has to be studied in a detergent that will not disrupt its native oligomeric size.
Prior to these studies, therefore, suitable detergents have to be screened. For example, in previous SE studies SH protein formed pentameric oligomers in DPC, C8E5 and C14-betaine, and pentamers were also observed during electrophoresis in the mild detergent PFO22. All these detergents are suitable for SE, since their density can be matched by addition of D2O. Samples should also be tested by sedimentation velocity (SV), which provides the minimum number of species present34.
Also, unsuitable detergents will produce multiple species in SV, and multiple bands in PFO electrophoresis, indicative of multiple non-specific association. Useful information from SE is obtained when the predominant number of species present in the system is not higher than two; in the case of SH protein, the data was fitted to an equilibrium between monomers and pentamers.
It is important to note that the model chosen to fit the SE data, should also be as simple as possible, i.e., a single species of unknown molecular weight should be tried first, followed by reversible equilibria between monomer and oligomers of increasing size. More complex models have associated a higher ambiguity. Also, small populations of other smaller or larger oligomers cannot be detected, and the model will just indicate what are the predominant species present. In the case of viroporins, this is especially evident because oligomeric size changes depending on subtle experimental conditions, e.g., detergent used, pH, protein concentration, centrifugation speed or mutations, e.g, in Hepatitis C Virus p7 35, influenza A M2 36 or SH protein 24.
The results obtained for WT SH protein clearly show that the association constant of the pentamer is reduced at pH 3. These results match channel activity measurements obtained with SH protein in synthetic lipid bilayers, where conductance was markedly reduced at pH 3 (pK ~4.5), whereas it remained constant between pH 7 and 524. Based on these results, it may be possible that SH channel activity is somewhat regulated by lower pH when present in native biological membranes. Although in the Golgi lumen the pH is only one unit below that of the cytoplasm, intravesicular pH drops along the endocytic pathway from pH 6.0–6.5 in early endosomes to pH 4.5–5.5 in late endosomes and lysosomes 37. In the infected cell, the pK for conductance changes or pentameric stability may be higher than those obtained in vitro herein and previously24, therefore pH could indeed play a role in modulating channel activity during the life cycle of the virus. Mutation of protonatable His residues in the context of the infected cell is an interesting avenue for future experiments38.
Finally, the study of SH protein is facilitated by the presence of Trp in its sequence, facilitating UV absorption measurements. However, knowledge of the structure of the oligomer can enable introduction of Trp at insensitive locations, e.g., lipid –or solvent- facing parts of the protein in membrane – or soluble- proteins. Alternatively the protein may be tagged with a Visible-absorbing or fluorescent label.
In summary, this protocol describes the application of SE to determine oligomeric size and association constants of a defined viral channel, when experimental parameters are changed. In this case, pH has been varied to study the effect on stability of His protonation, but many other hypothesis can be tested, such as effects of mutations on the structural integrity of these oligomers.
The authors have nothing to disclose.
This work has been funded by the National Research Foundation grant NRF-CRP4-2008-02 (J.T.) and Tier 1 grant RG 51/13.
Name of Material/ Equipment | Company | Catalog Number | Comments/Description |
3-(N,N-dimethylmyristylammonio)propanesulfonate | Sigma | T0807 | |
Deuterium oxide 99.8% | Cambridge Isotope | DLM-4-99.8 | |
An-50 Ti Rotor, Analytical, 8-Place | Beckman Coulter | 363782 | |
An-60 Ti Rotor, Analytical, 4-Place | Beckman Coulter | 361964 | |
Cell housing | Beckman Coulter | 334784 | |
12 mm six-channel centerpiece, epon charcoal-filled | Beckman Coulter | 331376 | |
Window holder | Beckman Coulter | 305037 | |
Window gasket | Beckman Coulter | 327021 | |
Window liner | Beckman Coulter | 362329 | |
Sapphire window | Beckman Coulter | 307177 | |
Quartz window | Beckman Coulter | 301730 | |
Screw-ring washer | Beckman Coulter | 362328 | |
Screw ring | Beckman Coulter | 301922 | |
Spinkote | Beckman Coulter | 306812 | |
Torque stand assembly | Beckman Coulter | 361318 | |
Counterbalance | Beckman Coulter | 360219 | |
Cell alignment tool | Beckman Coulter | 362340 | |
SEDNTERP | http://bitcwiki.sr.unh.edu/index.php/Main_Page | ||
HeteroAnalysis | http://www.biotech.uconn.edu/auf/?i=aufftp | ||
SEDFIT | http://www.analyticalultracentrifugation .com/sedfit.htm |
||
SEDPHAT | http://www.analyticalultracentrifugation .com/sedphat/default.htm |