Paramita Sahaa,
Samapan Sikdara,
Camelia Mannaa,
Jaydeb Chakrabartiab and
Mahua Ghosh*a
aDepartment of Chemical, Biological and Macromolecular Sciences, S. N. Bose National Centre for Basic Sciences, Block JD, Sector III, Salt lake, Kolkata 700106, India. E-mail: mahuaghosh@bose.res.in; mahua.ghosh@gmail.com
bThe Thematic Unit of Excellence on Computational Materials Science, S. N. Bose National Centre for Basic Sciences, Block JD, Sector III, Salt lake, Kolkata 700106, India
First published on 18th January 2017
STY3178 is a yfdX protein from Salmonella Typhi. yfdX proteins occur ubiquitously in a number of virulent bacteria but their cellular localization is unknown. Our earlier studies have shown that STY3178 is a trimer and can be a periplasmic chaperone protein. In the present study we show the stability of STY3178 in the presence of the bio-mimetic anionic surfactant sodium dodecyl sulphate (SDS). With increasing concentrations of SDS we observe monomeric STY3178 which reversibly forms the trimer upon decreasing the surfactant concentration. Protein tertiary structure is not perturbed in the presence of SDS. We show using molecular dynamics simulation and conformational thermodynamics data that SDS induces stability of the monomer compared to an isolated monomer of STY3178. This supports our experimental observations.
Study of interaction of surfactants with bacterial proteins might provide vital information regarding the host cell surface interactions. It can also provide clue to the intracellular interactions between periplasmic protein and cellular membrane. yfdX proteins are ubiquitously present in virulent bacteria. However, the exact cellular localization for these proteins is still unknown. Importance of these proteins are indicated by the following: (i) its presence in disease causing bacteria, (ii) its up-regulation in E. coli by a multidrug response regulator protein EvgA7 and (iii) its interaction with antibiotics like ciprofloxacin, rifampin and ampicillin for Salmonella Typhi protein STY3178.7 It is also reported that STY3178 show stability in wide range of pH.8 We have shown the thermally stable STY3178 has the chaperone activity capability in absence of ATP like known periplasmic chaperones.9 STY3178 shows reversibility in unfolding in presence of thermal9 and chemical perturbations.8
Here we study the effect of most widely used anionic surfactant sodium dodecyl sulphate (SDS) on the structural stability of STY3178. SDS is known to affect protein structure in a variety of ways. For instance, α-helical and α/β mix proteins appear to be more prone to SDS-induced unfolding than the β-sheeted ones.3 Above the critical micellar concentration (CMC), SDS stabilizes α-helical conformation,10–12 whereas below CMC, it stabilizes β-sheets.12–14 SDS is also reported to alter the tertiary structure keeping the secondary structure unchanged in some proteins such as cytochrome c.12,15–21 It as well affects oligomeric states of various proteins. For instance, EmrE22 from small multidrug resistance (SMR) protein family, BRICHOS23 domain trimer (related to human amyloid disease) and multioligomeric haemoglobin24 are reported to form either monomers or dissociate in presence of SDS. Similarly suppression of aggregation is observed for insulin25 and alcohol dehydrogenase in presence of SDS.
The main objective of the present studies is as follows. The homologues of STY3178 from other virulent bacteria like Salmonella Typhimurium, Escherichia coli, Klebsiella pneumoniae, Salmonella paratyphi yfdX proteins, are all predicted as periplasmic proteins in Cello.26,27 The interactions of yfdX protein with surfactant can possibly mimic the interaction between cell membrane and periplasmic proteins. This leads us to probe the interactions of surfactant with STY3178. Among these yfdX proteins, two are identified as multimers. STY3178 is trimeric in solution as reported in our earlier study.7 The homologue of STY3178 from K. pneumoniae is shown as a tetramer in the deposited structure in PDB. Surfactants at high concentrations are reported to dissociate complex structures of proteins.28–32 Moreover, it is also known that proteins get unfolded in presence of high surfactant concentration.28–32 The key question is: does STY3178 retain its secondary and tertiary structure in SDS?
Structural changes of STY3178 in presence of SDS are monitored using dynamic light scattering (DLS), size exclusion chromatography (SEC), steady state fluorescence, circular dichroism (CD) and isothermal titration calorimetry (ITC) in our present study. DLS measurements show decrease in hydrodynamic size of the protein with increase in SDS concentration, indicating dissociation of STY3178 trimer which is further confirmed by SEC. We find gradual lowering of fluorescence emission intensity without any effect in the peak position in presence of SDS. Near UV-CD spectrum shows loss of broad shoulder in SDS whereas far-UV CD show the α-helical secondary structure remains intact. ITC measurement indicates two regimes during the titration of SDS to the protein, which matches with the DLS observations. To further understand the SDS binding, we perform molecular dynamics simulation. It has been shown that the stabilization of protein–ligand systems can be understood from the conformational free energy and entropy of the system.33,34 Our earlier works show that these quantities can be extracted from the distribution of dihedral angles over the simulated trajectory. Further, the destabilized and disordered residues in a given conformation in reference to another, serve as ligand binding sites in proteins for further activities.35,36 We estimate the conformational thermodynamics changes and observe the SDS induced stability of modeled monomer of the protein compared to that in the absence of SDS. Our data support our experimental observations of SDS interactions with STY3178.
We consider the equilibrated ensemble of the docked complex between 10–40 ns for thermodynamic calculations. The conformational thermodynamic changes of SDS bound STY3178 monomer with respect to isolated monomer is calculated using the histogram based method (HBM),33,35,36,40–42 the detailed description of which can be found in a recent study.8,33 The equilibrated MD trajectory is divided into small time windows (10 ns each) to generate dihedral histograms for every time window. The conformational free energy and entropy have been calculated for each window and averaged over them for each monomer. The error has been estimated from the standard deviations of these values.
Far UV-CD spectrum of SY3178 shows its α-helical nature as reported earlier7 (Fig. 1B). In presence of 69.6 mM SDS (2%), where we obtain a monomer like hydrodynamic size, the CD spectrum is slightly altered. The two minima at 209 and 222 nm show a slight variation in ellipticity. However, the spectrum does not look like an unfolded protein CD signature indicating that the conformation of STY3178 in presence of 69.6 mM SDS (2%) is folded. When we monitor the ellipticity values at 209 and 222 nm for increasing concentration of SDS (0–69.6 mM), we observe ellipticity decreases at 222 nm and increases at 209 nm for 3.48 mM (0.1%) SDS. Beyond this concentration of SDS, no substantial change in ellipticity is observed for these two wavelengths (Fig. 1C). Overall there may be a slight change in helix content compared to the native protein keeping the secondary structure similar. Decrease in helix content in presence of 69.6 mM (2%) SDS compared to native is about 9.3%.
The near UV-CD signature of the protein contains broad shoulder7 in the range 250 nm to 280 nm, which gets altered in presence of 0.696 mM (0.02%) SDS (Fig. 1D). This broad shoulder between 250–280 nm arises due to the presence of phenylalanine and tyrosine in the tertiary structure as reported earlier.7 It has been reported earlier that aromatic residues interact with SDS.43 This might be the reason due to which the broad shoulder of protein is not observed in the near UV-CD spectrum, even at the lowest concentration of SDS. Similarly in fluorescence emission in presence of SDS, we observe decrease in intensity for the three excitation wavelengths 275, 280 and 295 nm (data not shown). There is no change observed in the emission peak position in presence of SDS.
We perform size exclusion chromatography (SEC) using superdex 75 column to confirm the possible dissociation of STY3178 trimer in presence of SDS. SEC is performed in presence of 3.48 mM (0.1%) and 69.6 mM (2%) SDS and elution volumes are compared with that of the native protein. Fig. 1E shows the chromatogram of protein eluted from the column. The elution volume of protein in 3.48 mM SDS is around 120 ml and matches with that of the native STY3178 (Fig. 1E). This confirms that STY3178 remains in trimeric state in presence of 3.48 mM SDS. However, the protein in 69.6 mM SDS elutes at a volume of around 160 ml (Fig. 1E). We compare this elution volume with that of several other standard proteins eluted from the same superdex 75 column. The elution volume of 160 ml corresponds to a protein with molecular weight of about 23 kDa (Fig. 1F), confirming SDS-induced dissociation of trimeric STY3178 into its monomers.
We further perform ITC measurements to determine interaction of STY3178 with SDS. Fig. 2A shows the thermogram of titration of SDS into STY3178 solution at 298 K. The corresponding isotherm is shown in Fig. 2B. There are two regions in the ITC isotherm (Fig. 2B). In the first region (up to 100 molar ratio) up to 14th injection, the final SDS concentration is 33.2 mM (0.95%) (Fig. 2B). In the same SDS concentration range, we find decrease in hydrodynamic diameter from 6.5 nm to 5.6 nm in DLS (Fig. 1A). The second region in ITC isotherm corresponds to SDS concentration of about 63.2 mM (1.8%) at final injection (Fig. 2B), where we observe hydrodynamic diameter (∼4.2 nm) of monomeric protein in DLS (Fig. 1A). The ITC data suggest that SDS binds to the trimeric assembly and the corresponding monomer differently.
Fig. 2 ITC measurement in presence of SDS. (A) Shows the thermogram of titration of 500 mM SDS into STY3178. (B) Shows the isotherm of SDS binding with STY3178. |
To further understand qualitatively the SDS binding and stabilization of the monomer in presence of SDS, we perform molecular dynamics (MD) simulation on model system. We consider a monomeric STY3178 and one SDS molecule in our model studies. A couple of points are noteworthy on the construction of our model system: (1) SDS interaction takes place with both monomeric and trimeric protein. For simplicity of the calculation we restrict only to the monomeric state. (2) We have considered one SDS molecule, even though the experiments indicate that the protein monomers are stable for large SDS concentrations much above the CMC where the SDS molecules form micelles. The SDS molecules are amphiphilic in nature, containing negatively charged head-group and hydrophobic tail. Due to amphiphilic nature, even in the micelles the head groups are solvent exposed which can interact with protein as in a single SDS molecule. Thus the interactions responsible for SDS and the protein interactions have been adequately taken into account in our model calculations.
We have shown earlier a simulated monomer model structure which can explain experimentally observed thermal stability.9 In a recent study we have shown the MD simulated trimeric conformations of STY3178 and compared its stability over the monomer structure.8 We observe that the trimeric conformation is more stable compared to the monomer. However, there are some interfacial residues in the trimer which are destabilized with respect to the monomer.8 The changes in conformational free energy, ΔGconfi and entropy, TΔSconfi of these destabilized interfacial residues between two monomers are shown in Table 1 where the data is averaged over all the three chains. ΔGconfi and TΔSconfi in the table are quoted in units of kJ Mol−1. The estimated error for ΔGconfi is approximately 4% and that for TΔSconfi is 5%. Within this error bar the conformational free energy of the lowest magnitude still remains positive (0.02 ± 0.0008) indicating destabilization, albeit with very small magnitude. The total destabilization conformational free energy of the relevant residues add up to about 1kBT.
Interfacial residues | ΔGconfi | TΔSconfi |
---|---|---|
N114 | 0.1 | −0.5 |
R149 | 0.02 | −2.3 |
Q175 | 0.06 | −0.02 |
S184 | 0.8 | 0.6 |
V185 | 0.7 | 1.2 |
E187 | 0.3 | 0.3 |
The destabilized and disordered polar, basic or hydrophobic residues can bind to amphiphilic SDS molecule. If we consider the destabilized and disordered residues as shown in Table 1, we observe that the majority of the residues are polar, along with basic and hydrophobic ones as well. The destabilized and disordered polar, basic or hydrophobic residues can bind amphiphilic SDS molecule with negatively charged head group. This is because of the amphiphilic nature of SDS containing negatively charged head-group and hydrophobic tail. If we consider the destabilized and disordered residues as shown in Table 1, we observe 66% of them are polar, 17% is basic and 17% is hydrophobic.
We perform biased docking of SDS molecule with a simulated monomer of STY3178 using N114, R149, Q175, S184 and V185 as active residues. The docked monomer–SDS complex is then simulated to obtain ensemble of structures. We estimate conformational thermodynamic changes of SDS bound monomer with respect to SDS-free monomer. The simulated docked complex is shown in Fig. 3A. The SDS molecule after simulation stabilizes near the residues V180, S182, Q183, S184, V185, Q186 and S188, which forms the binding region. The SDS bound monomer shows overall gain in conformational stability, ΔGconf = −100.3 kJ mol−1 and ordering TΔSconf = −431.8 kJ mol−1 with respect to an isolated monomer. A residue wise conformational thermodynamics analysis upon SDS binding show that in the SDS–monomer complex, residues V180, S182, Q183, S184, V185, Q186 and S188 get stabilized. The conformational free energy changes, ΔGconfi, shown in black, in Fig. 3B, are computed for monomer–SDS complex with respect to the free monomer. We compare the ΔGconfi (grey) for these residues in the trimeric assembly, where we compute conformational free energy changes of the trimer with respect to monomer. The residue-wise conformational free energy changes in the trimer is computed by taking an average over the three units (Fig. 3B). ΔGconfi of V180, Q186 and S188 are negative in both the trimeric assembly as well as in the monomer–SDS complex, albeit the degree of stabilization of Q186 is more upon SDS binding. The other residues, S182, Q183, S184 and V185, which are destabilized (ΔGconfi > 0) in the trimeric assembly, undergo conformational stabilization upon SDS binding.
The simulation results indicate that SDS binding stabilizes monomeric STY3178. This is in qualitative agreement to our experimental observations that high concentration of SDS dissociates trimeric STY3178 and stabilizes the monomer. These observations are similar to the interactions between surfactants and multimeric proteins reported earlier.28–32 We also show that monomer–monomer interface of the trimeric assembly is suitable binding surface for SDS. Since the exposed charged head groups are more in number in SDS micelles, the electrostatic interactions with the destabilized polar residues will be enhanced compared to the single SDS molecules. This may lead to difference between the binding of SDS with STY3178 for high and low concentrations of SDS.
This journal is © The Royal Society of Chemistry 2017 |