Comparing effects of attractive interactions in crowded systems: nonspecific, hydrophobic, and hydrogen bond interactions
- Published
- Accepted
- Received
- Academic Editor
- Gerrick Lindberg
- Subject Areas
- Theoretical and Computational Chemistry, Biophysical Chemistry, Physical Chemistry (other), Thermodynamics and Statistical Mechanics
- Keywords
- Crowding, Protein folding, Molecular interactions, Hydrophobic effect, Molecular simulation, Monte Carlo
- Copyright
- © 2024 Bazmi and Wallin
- Licence
- This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ Physical Chemistry) and either DOI or URL of the article must be cited.
- Cite this article
- 2024. Comparing effects of attractive interactions in crowded systems: nonspecific, hydrophobic, and hydrogen bond interactions. PeerJ Physical Chemistry 6:e31 https://doi.org/10.7717/peerj-pchem.31
Abstract
The equilibrium stability of a protein is determined by its amino acid sequence and the solution conditions, such as temperature, pH and presence of chemical denaturant. The stability of a single protein in two identical solutions can nonetheless differ if other macromolecules, termed cosolutes or crowders, are present in one of the solutions at concentrations high enough to occupy a substantial fraction of the solution volume. This effect, due to the presence of the crowders, decreases or increases the stability depending on the interactions between the protein and crowders. Hard-core steric repulsions, which are responsible for the reduction in free volume, are expected to entropically stabilize the protein while attractive interactions can be destabilizing. Here we use a coarse-grained protein model to assess the impact of different types of crowder-protein interactions on the stability of a 35-amino acid model sequence folding into a helical bundle. We find that, for the same interaction strength and concentration, spherical crowders with a hydrophobic character are more destabilizing than crowders interacting nonspecifically with the protein. However, the two types of interactions differ in the degree of association between crowders and protein. At an interaction strength for which the attractive interactions roughly counteracts the stabilizing hard-core repulsions, the nonspecific interactions lead to much stronger crowder-protein association than the hydrophobic interactions. Additionally, we study crowders in the form of polypeptide chains, which are capable of hydrogen bonding with the protein. These peptide crowders have a destabilizing effect even at relatively low crowder concentrations, especially if the sequence of the peptide crowders includes hydrophobic amino acids. Our findings emphasize the importance of the interplay between different types of attractive crowder-protein interactions and entropic effects in determining the net effect on protein stability.
Introduction
Most biophysical experiments on biomolecules or biomolecular systems are performed under dilute solution conditions in which the macromolecular concentration rarely exceeds 10 g/L (Theillet et al., 2014). However, the cellular environment is often anything but dilute (Ellis, 2001). For example, the concentration of macromolecules in Escherichia coli can reach up to 300–400 g/L, corresponding to a volume occupancy of around 30–40% (Zimmerman & Trach, 1991). The crowded cellular milieu has been shown to impact a wide range of biophysical processes, including diffusion (Słyk, Skóra & Kondrat, 2022; Wang et al., 2012a), aggregation (Schreck, Bridstrup & Yuan, 2020; Siddiqui & Naeem, 2018), DNA replication (Akabayov et al., 2013), protein fold switching (Zhang et al., 2023; Bazmi, Seifi & Wallin, 2023), and liquid-liquid phase separation (André, Yewdall & Spruijt, 2023; Julius et al., 2019).
One of the major issues in macromolecular crowding, in fact since the inception of the field (Minton, 1980), is the effect of crowding on the equilibrium stability of proteins (Pastore & Temussi, 2022; Speer et al., 2022; Zhou, Rivas & Minton, 2008). In the simplest scenario, the protein populates mainly two states (Sosnick & Barrick, 2011): the structurally coherent native state, N, and the typically more extended and structurally diverse unfolded state, U. Protein stability reflects the relative population of these two states. It is therefore expected that different types of crowder-protein interactions differently impact stability. Hard-core steric repulsions, which arise simply from the fact that two macromolecules cannot simultaneously occupy the same region in space, are expected to stabilize proteins, because conformations of the extended U ensemble become entropically disfavored relative to the compact N under crowded conditions. Experimental studies employing artificial polymer macromolecules as crowder agents typically (Sasahara, McPhie & Minton, 2003; Spencer et al., 2005; Wang et al., 2012b; Hong & Gierasch, 2010; Christiansen & Wittung-Stafshede, 2014), although not always (Nasreen et al., 2018; Malik et al., 2012), stabilize proteins as measured by, e.g., the folding midpoint temperature or the free energy of folding. It was shown, within a coarse-grained model for folding, that volume exclusion by purely repulsive crowders can have a neutral effect on stability, or even lead to a net destabilization, under some conditions (Bazmi & Wallin, 2022).
In addition to steric repulsions, attractive or repulsive interactions between the protein and crowders will in general also impact protein stability. Such interactions, often called soft or “chemical”, arise from various effects, including charge-charge, van der Waals, hydrogen bonding, and hydrophobic interactions (Sarkar, Li & Pielak, 2013a). Soft interactions between protein and crowders are often assumed to be nonspecific (Sarkar, Li & Pielak, 2013a; Guin & Gruebele, 2019; Rivas & Minton, 2022) in contrast to the specific (Jones & Thornton, 1996) and cooperative (Hsu, Yen & Yeang, 2022) interactions that occur between the components of functionally relevant biomolecular complexes. Repulsive soft interactions are expected to enhance the stabilizing effect of hard-core repulsions, i.e., make the protein even more stable. By contrast, attractive soft interactions between protein and crowders should counteract the hard-core stabilization (Sarkar, Smith & Pielak, 2013b). The reason is that, while the conformations of both U and N will make favorable contacts with the crowders, the larger surface area of the U should lead to a net energetic stabilization of this state, and hence provide a destabilizing effect on the protein. Both computational (Feig & Sugita, 2012; Rosen, Kim & Mittal, 2011; Kim & Mittal, 2013; Douglas, Dudowicz & Freed, 2009; Bille, Mohanty & Irbäck, 2016; Tsao & Dokholyan, 2010; Macdonald et al., 2016) and experimental (Phillip, Kiss & Schreiber, 2012; Jiao et al., 2010; Miklos et al., 2010) studies have been performed to study the effects of attractive crowder-protein interactions. In general, the net effect on the stability will be determined by a competition between stabilizing steric repulsions and destabilizing soft attractive interactions. Because soft interactions have an energetic (enthalpic) component, the impact of crowders should be temperature-dependent (Miklos et al., 2010). Indeed, this was observed by Zhou (2004) who studied the stability of the protein CI2 in the presence of protein crowders. CI2 was found to be destabilized at low temperatures and stabilized at high temperatures. This change allowed the authors to define a so-called crossover temperature at which there was no net effect on the stability of CI2 (Timr & Sterpone, 2021).
In this work we aim to delineate the effects of different types of soft attractive interactions on protein stability. We compare the effect of crowders that are geometrically identical but differ in their physical properties. In particular, we consider spherical “nonspecific crowders”, which are capable of energetically favorable interactions with any part of the protein, and spherical hydrophobic crowders, for which these favorable interactions are limited to nonpolar (hydrophobic) amino acids. Hydrophobic crowders are expected to be destabilizing because U, but not N, will be favored by effective attractive interactions with the crowders, assuming a native conformation with a completely buried hydrophobic core. However, it is unclear a priori if hydrophobic crowders should be more or less destabilizing than crowders with nonspecific interactions, because the overall destabilization is determined by the net effect of crowder interactions with U and N, which will depend on the structural details of these states. For comparison, we also consider the results from excluded volume crowders, which are stabilizing due to steric repulsions. In addition to spherical crowders, we consider the crowding effect from short polypeptide chains, which can favorably interact with the protein through hydrophobic interactions and hydrogen bonding via their backbone NH and CO groups.
To this end, we use a coarse-grained model for folding which combines an all-atom backbone geometry with a one-bead sidechain representation, an enlarged Cβ atom. This model relies on a simplified amino acid alphabet with 3 types: polar (p), hydrophobic (h), and turn (t). Folding is driven by backbone-backbone hydrogen bonding and effective hydrophobic interactions (pairwise hh-attractions). Different sequences can be designed using simple principles (Cordes, Davidson & Sauer, 1996) leading the chain to adopt various protein-like folds (Bhattacherjee & Wallin, 2012; Trotter & Wallin, 2020). We focus here on a 35-amino acid sequence that folds into a stable α-helical hairpin fold at low temperatures (Holzgräfe & Wallin, 2015), and is stabilized in the presence of excluded volume crowders (Bazmi & Wallin, 2022). We implement nonspecific interactions by making contacts between crowders and any Cβ atom on the protein (p and h amino acids) and limit the interactions to h amino acids in the case of hydrophobic crowders. For a given strength of the attraction, we find that the hydrophobic crowders are more destabilizing than the crowders interacting nonspecifically with the protein. Moreover, because the hydrophobic crowders have fewer interaction sites on the protein compared to the nonspecific crowders, the overall crowder-protein association is much weaker for the hydrophobic crowders despite their stronger destabilizing effect. Crowders that drive protein destabilization through nonspecific interactions rely on a difference in accessible interaction sites in U and N, and therefore require a rather strong crowder-protein association at the point where they can overcome the stabilizing effect of volume exclusion.
Methods
Coarse-grained model for protein folding
To model protein folding, we use the coarse-grained “Cβ-model” described in Bhattacherjee & Wallin (2012). It is a model with three different types of amino acids: polar (p), hydrophobic (h), and turn (t) amino acids. Geometrically, the protein chain is described using an atomistic backbone (Cα, C′, N, H, Hα1, and O) and simplified sidechains using an enlarged Cβ atom. The t amino acids differ from p and h in that it does not contain a Cβ atom, which is instead replaced by an Hα2 atom. Hence, t is strongly related to glycine and more flexible than p and h. Chain conformations are completely specified by the 2N Ramachandran angles . Hence, bond lengths, bond angles, and dihedral angles (e.g., the peptide plane angle ω = 180°) are held fixed at standard values. The potential energy function of the model, Ep, includes four terms Ep = Eloc + Eexvol + Ehbond + Ehp, described in detail in Bhattacherjee & Wallin (2012). Briefly, the first term, Eexvol, provides atom-atom repulsions with a range σa determined by the sum of the atoms’ van der Waals radii, . The second term, Eloc, represents electrostatic interactions between partial charges of adjacent peptide planes. This term is included to ensure sampled ϕi, ψj angles agree with statistics from real protein structures (Ramachandran plots).
The two last terms, Ehbond and Ehp, represent hydrogen bonding and effective hydrophobic interactions, respectively. These terms are essential for stabilizing and driving the formation of the native state. Hydrogen bonds are modeled using: (1)
where the sum is over pairs of CO, NH groups, rij is the OH distance, σhb = 2.0 Å, and khbond = 3.1. The interaction strength is modified by a sequence-dependent scale factor γij taken to be 1 for hh, hp, and pp pairs, and 0.75 for tt, th, and tp pairs. The reduced hydrogen bonding capacity of t amino acids mimics the secondary structure breaking ability of glycine. The directional dependence is implemented via the factor , where αij and βij are the NHO and HOC′ angles, respectively. Additionally, if either αij < 90°or βij < 90°, the ij contribution is set to zero. The hydrophobic effect is modeled using pairwise additive interactions between hydrophobic amino acids. Specifically, (2) where the sum is over pairs of Cβ atoms of h amino acids, excluding nearest and next-nearest neighbors along the chain, and (3) is a function with a maximum at r = r0. We set σhp = 5 Å and the interaction strength khp = 0.805. In this model, tuning the relative strength of hydrogen bonding and hydrophobic forces is essential to obtain folding into stable and protein-like structures (Irbäck, Sjunnesson & Wallin, 2000; Irbäck, Sjunnesson & Wallin, 2001).
Crowders
Repulsive interactions are modeled using the pair potential (Mittal & Best, 2010) (4) where r is a crowder-crowder or crowder-atom distance. Because V(ρ) = 1 and V → ∞ as r → ρ − σ (for r < ρ − σ, we set V = ∞) the two parameters ρ and σ determine the range and “sharpness” of the repulsion, respectively. We set ρ = 2Rc for crowder-crowder interactions and ρ = ρcp = Rc + σa for crowder-atom interactions, where Rc and σa are the radii of the crowder and atom, respectively. We set σ = 6 Å for crowder-crowder interactions and 3 Å + σa for crowder-atom interactions. A crowder-atom attractive interaction is modeled with the potential V(r) − ϵattg(r; ρcp), where the function g is given by Eq. (3). The form of this potential for different attraction strengths ϵatt is shown in Fig. 1. For a system consisting of a protein with Na atoms and Ncr crowder particles, the total potential energy then becomes E = Ep + Ecc + Ecp, where (5) and (6) are the crowder-crowder and crowder-protein energies, respectively. In Eq. (6), the second sum within square brackets controls which atoms are subject to crowder attraction. For crowders with nonspecific attraction to the protein, the sum is over all Cβ atoms (p and h amino acids). For hydrophobic crowder, the sum is over the Cβ atoms of h amino acids. For excluded volume crowders, ϵatt = 0.
Simulated tempering Monte Carlo
To find the thermodynamic behavior of various protein-crowder systems, as determined by the amino acid sequence, number of crowders, and the energy function E(r), we use simulated tempering Monte Carlo (MC) (Marinari & Parisi, 1992; Lyubartsev et al., 1992; Irbäck & Potthast, 1995). In addition to a random walk in conformational space, as in basic Monte Carlo, simulated tempering also carries out a random walk in temperature while keeping the simulation at equilibrium. This is achieved by defining a set of temperatures, , and simulating the joint probability distribution (7) where βj = 1/kBTj, kB is Boltzmann’s constant, and j has been made a dynamic parameter. The gj’s are M simulation parameters that control the marginal distribution P(j). Jumps between temperatures, j → j′, are accomplished as MC updates, with acceptance probability (8) A common choice for the gj parameters, which we follow here, is to select them such that P(j) is roughly flat, ensuring that sampling of conformations takes place at each temperature Tj.
Simulations and analysis details
Equilibrium behaviors of crowder-peptide systems are determined using simulated tempering Monte Carlo simulations (Marinari & Parisi, 1992). Runs are carried out with either 8 temperatures in the range kBT = 0.48–0.68 or 10 temperatures in the range 0.40–0.70. For each system, 10 independent runs with 5 × 109 elementary MC updates are performed. Initial configurations are obtained by picking random chain conformations and random spherical crowder positions, followed by a relaxation step to remove any hard-core steric clashes. Monte Carlo chain updates are divided equally between the protein chain and the crowder particles. Two different types of chain moves are used: biased Gaussian steps (BGS) (Favrin, Irbäck & Sjunnesson, 2001), which produce approximately local chain deformations, and pivot moves, which produce global changes. In the latter type, a single ϕi or ψi angle is set to a new random value such that the chain rotates around the NCα or CαC′ bonds. The frequencies of chain updates are set so that BGS is most frequent at low Ts and pivot is most frequent at high Ts. Temperature updates are attempted every 100 MC step. Spherical crowder simulations are carried out in a cubic box of side length L = 100 Å while for simulations of polypeptide crowders L = 65 Å. The number of spherical crowding agents are 14, 28, 42, and 56, corresponding to packing fractions ϕc = 0.10, 0.20, 0.30, and 0,40, respectively. The crowder radius is Rc = 12 Å. For simulation carried out with 8 temperatures in the range kBT = 0.48–0.68, the multistate Bennett acceptance ratio (MBAR) technique was applied to determine thermodynamic averages in the range kBT = 0.40–0.70.
Observables
We quantify native state stability using two quantities, ΔF and Tm. The free energy of folding is determined using (9) where FU and FN are the free energies of the unfolded and native states, respectively, and Pnat is the population of the native state. The native state population Pnat is determined as in Ref. Bazmi & Wallin (2022). Conformations are considered part of the native state if Qnat ≥ Qcut where Qnat is the number of formed native contacts and Qcut = 50. The folding midpoint temperature, Tm, is determined by fitting the temperature dependence of Pnat to a two-state model with two fit parameters.
Results
Spherical crowders
Using MC simulations, we determined the thermodynamics behavior of our model protein, α35, in isolation and in the presence of three different types of spherical crowders having the same radius (Rc = 12 Å) but different interactions with the protein chain: (1) excluded volume crowders, capable of only steric repulsions; (2) “nonspecific” crowders, which in addition to steric repulsions, form energetically favorable contacts with polar and nonpolar amino acids; and (3) hydrophobic crowders, which differ from the nonspecific crowders in that they interact favorably with nonpolar amino acids only. We varied the crowder volume fraction, ϕc, and the strength ϵatt of the attractive protein-crowder interactions. In general, we find that the equilibrium stability of the α35 native state, as quantified by the free energy of folding, ΔF = − kBTln[Pnat/(1 − Pnat)], where Pnat is the population of the native state, depends on both ϕc and ϵatt. As an illustration of our results, Fig. 2 shows the temperature dependence of ΔF at ϵatt = 1.5, for different ϕc and crowder types.
In the absence of crowders (solid black curve in Fig. 2), the helical hairpin fold of α35 is highly stable at low T. For example, at the lowest studied temperature (kBT = 0.40), which corresponds to , where is the midpoint temperature of the α35 folding transition for ϕc = 0, ΔF/kBT = − 7.1, which is within the range of stabilities typically found for single-domain proteins (Zeldovich, Chen & Shakhnovich, 2007). The midpoint temperature , i.e., the temperature for which ΔF = 0, is also commonly used as a measure of protein stability. For α35, , as obtained previously (Bazmi & Wallin, 2022) by fitting the of α35 folding curve i.e., Pnat as function of T, to a two-state model with two free parameters. For the values of ϵatt considered here, α35 remains stable in the low T region even when the crowder packing fraction reaches ϕc = 0.40. However, there are clear variations between different crowder types as can be seen in Fig. 2. For example, for ϕc = 0.30 and kBT = 0 .40, the excluded volume crowders give ΔF/kBT = − 8.3, which is a decrease relative to the ϕc = 0 case. For the nonspecific and hydrophobic crowders the corresponding ΔF/kBT values are −6.1 and −3.3, respectively, which, by contrast, are increases relative to ϕc = 0.
At a temperature just below the folding midpoint, , we observe the following trends. Upon the addition of excluded volume crowders, ΔF decreases monotonically with ϕc, as shown in Figs. 3A and 3C, meaning a strict stabilization of the protein. The soft interactions of the nonspecific and hydrophobic crowders are expected to be destabilizing. However, because these two crowder types occupy the same space as the excluded volume crowders, they also provide a stabilizing effect. The net effect will be determined by a competition between the stabilizing steric repulsions and the destabilizing soft interactions. Indeed, at low attraction strength, ϵatt = 1.0, the addition of either nonspecific or hydrophobic crowders leads to an overall stabilization of α35 (see Fig. 3A and 3C). At ϵatt = 1.5, the hydrophobic crowders overcome the stabilizing excluded volume effect leading to a net destabilization, while the nonspecific crowders are still net stabilizing (see Figs. 3B and 3D), although for the nonspecific crowders ΔF ≈ 0 at ϕc = 0.40. The picture obtained is similar if instead Tm is used to quantify native state stability. As shown in Fig. 3C, there is a monotonic increase in Tm for both the excluded volume crowders and for the hydrophobic and nonspecific crowders with ϵatt = 1.0. For ϵatt = 1.5, the nonspecific crowders still exhibit a monotonic increase in Tm with increases ϕc while the hydrophobic crowders instead exhibit a monotonic decrease (see Fig. 3D). Overall, our results show that excluded volume crowders provide a stabilizing effect on α35, which can be counteracted by the presence of soft attractive interactions. The hydrophobic crowders provide a stronger destabilizing effect than the nonspecific crowders.
How can the stronger destabilizing effect of hydrophobic attractions be explained?
To address the question of why hydrophobic attractive interactions are more destabilizing than nonspecific interactions, we consider the interaction energy between crowders and protein, Ecp (see Methods). Ecp is a mix of repulsive and attractive energy contributions. We consider, in particular, and , i.e., the crowder-protein energy determined separately for the U and N states. When the difference there is a net energetic stabilization of U, which should have a destabilizing effect on the protein.
Figures 4A and 4B show the ϕc-dependence of and for two different values of ϵatt. For hydrophobic crowders and ϵatt = 1.0, (open circles) exhibits a slight positive curve and a slight negative curve. These trends are in line with an expanded U with some nonpolar amino acids available for favorable interactions with the crowders, and an N with a hydrophobic core well shielded from the crowders. For ϵatt = 1.5, the trend is similar but now exhibits a slight negative curve, indicating that, at this strength of attractions, N becomes slightly distorted to accommodate favorable contacts with the hydrophobic crowders. However, these crowder-proteins interactions are rather limited. At the highest packing fraction, ϕc = 0.40, , corresponding to ≈ 1–2 fully formed contacts between an nonpolar amino acid and a crowder. In comparison, for the nonspecific case, interactions between crowders and protein are much more prevalent for both U and N. is a sharply decreasing function of ϕc for both ϵatt = 1.0 and 1.5, which will strongly stabilize U at high ϕc. However, the decrease in is nearly fully compensated by a decrease in . Comparing the two types of crowders, we find that the net energetic effect, , which is what drives destabilization, is more negative for the hydrophobic crowders than for the nonspecific crowders, as can be seen in Figs. 4C and 4D.
Hence, an interesting difference between the two crowder types is that, while their net effects on stability are similar, they differ substantially in the degree of association with the protein. While hydrophobic crowders associate weakly with the protein, even in U, the protein-crowder association in the case of nonspecific interactions is greater by approximately an order of magnitude, as quantified by the magnitude of the interaction energies in the native state, (see Figs. 4A and 4B).
Balancing destabilizing soft interactions and stabilizing steric repulsions
The α35 protein is net stabilized by the presence of hydrophobic crowders at ϵatt = 1.0 but net destabilized at ϵatt = 1.5. This suggests a critical attraction strength at which there is no net change in stability. Indeed, for ϵatt = 1.25, the relative change in Tm is at most one percent and the change in folding free energy, ΔΔF, is very close to zero, as shown in Figs. 5A and 5B. Interestingly, this balance between stabilizing and destabilizing effects is largely independent of ϕc. The crowding effect on stability at any given ϵatt is, however, generally dependent on T, as shown in Fig. 5C for ΔΔF. There is a T-dependence also for excluded volume crowders, which are strictly stabilizing (ΔΔF < 0). Such a T-dependence is possible for excluded volume crowders despite their lack of energetically favorable energetic interactions with the protein, because the size of U is not constant but vary with T (Bazmi & Wallin, 2022). For hydrophobic crowders, due to the energetic stabilization of U relative to N, at low enough temperatures, the net effect crosses over from stabilizing to destabilization, i.e., ΔΔF changes from negative to positive. This crossover temperature (Zhou, 2013), Tc, increases with the contact strength ϵatt, as shown in Fig. 5D. Our findings for α35 are similar to previous studies on other proteins, such as ubiquitin, which showed a crossover temperature in the presence of either synthetic polymers or protein crowders (Wang et al., 2012b; Zhou, 2013).
Polypeptide crowders
We turn now to the folding of α35 in the presence of polypeptide chains. We use the same model for the peptides as for α35, with the restriction that energetically favorable interactions (i.e., hydrophobic and hydrogen bond interactions) between crowding peptides are turned off. This avoids the formation of oligomeric peptide structures, which would complicate the analysis. We consider two different 5-amino acid sequences: ppppp (peptide 1) and pphpp (peptide 2). Both peptides can thus interact with the α35 chain through backbone-backbone hydrogen bonding and peptide 2, due to its central h amino acid, can additionally interact with α35 through hydrophobic attractions.
Figures 6A and 6B show the T-dependence of the free energy of folding, ΔF, for different numbers Npep of peptide 1 or peptide 2 chains added to the system. For Npep ≤ 42, peptide 1 induces a weak stabilization at high T and weak destabilization at low T. This behavior is consistent with a competition between entropy-driven stabilization due to volume exclusion by the peptides and energy-driven destabilization due to inter-chain hydrogen bonding. Interestingly, there is a rather broad temperature range around (kBT ≈ 0.45–0.60) with no detectable change in ΔF (see Figs. 6A and 6D). The apparent lack of crowding effects in this temperature range can arise either because (i) the crowding effects are overall weak at the studied peptide concentrations or (ii) the two opposing crowding effects, peptide excluded volume and peptide-protein attractions, are equal in magnitude and therefore cancel out. To determine which scenario holds, we performed additional simulations with peptides that were only allowed to interact with other chains via repulsive interactions. These “excluded volume peptides” significantly stabilize α35 across all Ts (see Figs. 6C and 6D), which means that scenario (ii) above holds. Hence, in a relatively broad range around the midpoint temperature of the protein, the excluded volume effect due to the peptides is almost perfectly counteracted by soft interactions in the form of hydrogen bonds.
Peptide 2 provides a stronger destabilizing effect on α35 than peptide 1 at a given concentration, as seen in Figs. 6B and 6D. We therefore study up to Npep = 28 peptide 2 chains. At T ≈ T−, ΔF monotonically increases with the number of added peptide 2 chains in contrast to the flat ΔF exhibited by peptide 1. At very low T and Npep = 28, peptide 2 even leads α35 to become net unstable (ΔF > 0), although the error bars are larger at the lowest studied Ts. A similar behavior is seen for peptide 1 and Npep = 56. Structurally, low-energy conformations obtained for the peptide 1 and peptide 2 systems exhibit non-native features, including β-sheet structure. An example of a low-energy α35 structure in the precence of peptide 2 crowders is given in Fig. 7A. The destabilization at low T is, at least in part, driven by energetically favorable crowder-peptide interactions that also causes non-native structures. By contrast, low-energy conformations remain native-like in the case of excluded volume peptides (see Fig. 7B). It is instructive to get a rough idea of the volume fraction ϕc occupied by the peptides in our systems. Assuming each atom in our Cβ-model occupies volume according to its van der Waals radius, we obtain ϕc ≈ 0.07 for Npep = 28. Alternatively, if we assume amino acids are spheres with radius 3.8 Å (typical Cα- Cα distances of peptide bonds) we obtain a slightly higher value, ϕc ≈ 0.12. What makes peptide 2 more destabilizing than peptide 1? Because the two peptides are geometrically identical, any difference must derive from soft interactions. Figures 8A and 8B show the peptide concentration dependence of the quantity (Fig. 4), determined separately for peptide (crowder)-protein hydrogen bond () and hydrophobic () energies. For peptide 1, , since this peptide lacks a hydrophobic amino acid. For peptide 2, for all peptide-protein systems. This means that, vis-à-vis hydrophobic interactions, peptide 2 interacts more favorably with U than with N, causing α35 destabilization. The formation of hydrogen bonds between peptides and the protein are also destabilizing because generally . Interestingly, is more negative for peptide 2 than for peptide 1, such that the presence of a hydrophobic amino acid on peptide 2 also enhances the formation of hydrogen bonds with the protein chain. In other words, there is a cooperative effect between hydrophobic and hydrogen bond interactions in peptide 2, which further enhances the destabilizing effect of this peptide relative to peptide 1 which has only polar amino acids.
Discussion
We have used a sequence-based coarse-grained protein model to study crowding-induced changes to the stability of a model protein with a sequence that folds to an α-helical hairpin at low temperatures. In particular, we investigated crowders making different types of soft interactions with the protein, i.e., interactions different from hard-core steric repulsions. We found that crowders with a hydrophobic character, i.e., crowders interacting favorably with only nonpolar amino acids, provide a stronger destabilizing effect compared to crowders that interact nonspecifically with the protein, i.e., make favorable contacts with both polar and nonpolar amino acids. Both types of soft interactions are counteracting the stabilizing effect of excluded volume. For weak attraction strengths, these crowders still increase the stability of the protein relative to the dilute limit where there are no crowding effects. At a critical strength of the attraction, the stabilizing and destabilizing effects cancel leading to a net zero change in protein stability over a wide range of crowder concentrations. Similar results were obtained in a structure-based one-bead-per-amino acid model studied by Mittal et al. (Kim, Bhattacharya & Mittal, 2014).
The destabilizing effect of either hydrophobic or nonspecific crowder-protein interactions arises because interactions with U are generally more energetically favorable than interactions with N. This is manifested by the net crowder-protein interactions, ΔEcp, which is negative for high ϕc for both crowder types. However, perhaps counterintuitively, the overall association of the hydrophobic crowders with the protein is much less pronounced than for the nonspecific crowders, at a given degree of destabilization. This is illustrated schematically in Fig. 9. Quantitatively, we find that the magnitude of the interaction energy between crowders and protein, Ecp, is larger for the nonspecific crowders than for the hydrophobic crowders by roughly an order of magnitude. The reason for the difference is the way the two crowder types achieve destabilization. Nonspecific crowders rely on a difference between two rather large favorable interaction energies for U and N. Hydrophobic crowders, for which favorable interactions with N are almost absent due to hydrophobic amino acids being buried in the native structure, obtain destabilization even for relatively weak favorable interactions with U. This difference between hydrophobic and nonspecific soft interactions might be tested experimentally using protein crowders, and varying the chemical nature of surface-exposed amino acids through mutations.
Our results also confirm that the addition of crowders attracted to a protein can lead to temperature-dependent crowding effects, as demonstrated experimentally (Wang et al., 2012b). Specifically, we found that the our model protein exhibits a crossover temperature, Tcross, below which the crowders are destabilizing and above which the crowders are stabilizing. The presence of such a crossover temperature is well established in the literature and has been observed in the presence of both synthetic polymer crowders and protein crowders (Wang et al., 2012b; Zhou, 2013; Timr & Sterpone, 2021). The origin of this crossover temperature is believed to be due to the fact that energetic effects are becoming more important at low T, while at high T entropic effects dominate (Zhou, 2013; Timr & Sterpone, 2021). It should be noted that the hydrophobic effect is itself partly entropically driven, while in our model it is treated as an effective energetically driven interaction. Destabilization driven entirely by soft hydrophobic interactions are not guaranteed to increase with decreasing T.
We have also investigated how the folding and stability of our model protein are impacted by short polypeptide chains, which differ from our spherical crowder in that they are able to interact with the protein through hydrogen bonding. Results for this system indicate an interplay between hydrogen bonds and hydrophobic interactions. Peptide chains with partial hydrophobic character exhibit much stronger destabilizing effects than peptides with only polar amino acids. Part of this destabilization is due to hydrophobic interactions directly stabilizing U, in the same way as found for the hydrophobic spherical crowders. However, we found that the presence of a hydrophobic amino acid in the sequence of our peptide crowder also tends to promote the formation of additional hydrogen bonds with the protein (see Fig. 8), thus further strengthening the destabilization. It should be pointed out, however, that our polypeptide crowders are relatively short. It would be interesting to test the crowding effect of much longer polypeptide chains in our model, which could mimic macromolecular crowders capable of both substantial excluded volume effects and soft interactions in the form of hydrophobic interactions and hydrogen bonding.
Conclusions
In summary, we have used a sequenced-based model to study the effects of attractive interactions between crowders and protein on the stability of the protein’s native state. Our results highlight the importance of considering both the type and the strength of soft crowder-protein interactions when evaluating the impact of crowding. Our main conclusions are as follows: (1) soft attractive interactions have a generally destabilizing effect on protein stability in the context of both spherical and polypeptide chain crowders; (2) hydrophobic protein-crowder interactions provide a stronger destabilizing effect in comparison with entirely nonspecific interactions; and, (3) at a given strength of the soft interactions, hydrophobic crowders exhibit a smaller degree of association with the protein compared to the crowders that bind nonspecifically to the protein. In light of our results, it would be very interesting to see experiments that probe the stability of a test protein in the presence of protein crowders with a variable number of hydrophobic surface amino acids.
Supplemental Information
Raw data for individual simulation trajectories
The total energy, crowder-protein energy, number of native contacts, and fraction of native conformations, as a function of temperature, obtained for each individual simulation run. Results are given separately for repulsive, nonspecific and hydrophobic crowders, and different crowder-protein interaction strengths.