Crowding-induced Elongated Conformation of Urea-unfolded Apoazurin: Investigating the Role of Crowder Shape In Silico
Fabio C. Zegarra, Dirar Homouz, Andrei G. Gasic, Lucas Babel, Michael Kovermann, Pernilla Wittung-Stafshede, Margaret S. Cheung
1 Crowding-induced Elongated Conformation of Urea-unfolded Apoazurin: Investigating the Role of Crowder Shape In Silico
Fabio C. Zegarra &1 , Dirar Homouz &1, 2, 3 , Andrei G. Gasic , Lucas Babel , Michael Kovermann , Pernilla Wittung-Stafshede , Margaret S. Cheung & Equal contribution 1. Department of Physics, University of Houston, USA 2. Department of Physics, Khalifa University of Science and Technology, Abu Dhabi, UAE 3. Center for Theoretical Biological Physics, Rice University, USA 4. Department of Chemistry, Universität Konstanz, Germany 5. Research School Chemical Biology (KoRS-CB), Universität Konstanz, Germany 6. Biology and Biological Engineering Department, Chalmers University, Sweden CORRESPONDENCE *Dr. Margaret Cheung [email protected] *Dr. Pernilla Wittung-Stafshede [email protected] *Dr. Michael Kovermann [email protected] Abstract
Here, we show by solution nuclear magnetic resonance measurements that the urea-unfolded protein apoazurin becomes elongated when the synthetic crowding agent dextran 20 is present, in contrast to the prediction from the macromolecular crowding effect based on the argument of volume exclusion. To explore the complex interactions beyond volume exclusion, we employed coarse-grained molecular dynamics simulations to explore the conformational ensemble of apoazurin in a box of monodisperse crowders under strong chemically denaturing conditions. The elongated conformation of unfolded apoazurin appears to result from the interplay of the effective attraction between the protein and crowders and the shape of the crowders. With a volume-conserving crowder model, we show that the crowder shape provides an anisotropic direction of the depletion force, in which a bundle of surrounding rod-like crowders stabilize an elongated conformation of unfolded apoazurin in the presence of effective attraction between the protein and crowders. I. INTRODUCTION
The cellular environment is highly crowded with different types of macromolecules that can occupy volumes up to 40 % of the total volume inside cells . To properly understand how proteins fold and function in vivo , it is necessary to consider the volume exclusion from surrounding macromolecules, known as the macromolecular crowding effect , which influences the structures and dynamics of a biopolymer in this restricted space. To test the crowding effect in vitro , sugar-based, branched polymers such as Ficoll and dextran are typically utilized as synthetic crowder agents to investigate the macromolecular crowding effect because they are considered to be chemically inactive agents in solution. Under the assumption that the dominant interaction between a crowding agent and a biopolymer is steric repulsion, the macromolecular crowding is shown to stabilize the conformation of compact biopolymers as shown in proteins , nucleic acids , and chromatins . However, here we made the puzzling discovery by solution nuclear magnetic resonance (NMR) spectroscopy that our model protein, apoazurin, unfolded by a high urea concentration adopts more elongated conformations in the presence of dextran 20 than without a crowding agent (bulk). This is a deviation from a previous study from the Wittung-Stafshede group where they have shown that apoazurin increases in stability with synthetic crowders . The thermodynamic stability of the two-state folding protein apoazurin increases with the increased volume fraction of crowders, and this effect was independent of the crowder size when various dextran sizes were analyzed, suggesting that a cylindrical representation from Minton’s model is appropriate to model dextran that interacts with apoazurin through steric repulsions . To find an explanation to this counterintuitive result, we considered in silico coarse-grained molecular dynamics simulations, allowing more complex modeling of crowders than a simple hard sphere model for volume exclusion. The question central to all computational modeling is how to choose an appropriate level of the coarse graining of crowders that captures essential features without a prohibitive computational cost. All-atom molecular dynamics simulations have shown an unprecedented resolution of protein and protein crowders that allows investigation of subtle chemical interactions affecting protein stability and kinetics . These simulations demand massive computational resources and still require advanced sampling techniques to explore a wide range of phase space in a system for the computation of thermodynamic properties . Coarse-grained models, although losing atomistic details, allow for greater sampling than in all-atom models and provide key polymeric insight into the distribution of protein conformations and the folding energy landscape under the macromolecular crowding effect . To investigate the structural ensemble of apoazurin in 6 M urea in a periodic box of dextran crowders with coarse-grained molecular simulations, we must consider the appropriate dextran model. It is challenging to determine the geometry and property of dextran in a coarse-grained model because of its branched topology and polydisperse spatial distribution in solution . Early experiments based on filtration have characterized that the shapes of dextran crowders are asymmetrical by showing that dextrans pass through a pore much smaller than their presumed diameters . Sasahara et al have further modeled the shape of Dextrans as rods to fit the data from experimental measurement of the protein stability in dextrans with circular dichroism . Our previous work has adopted a spherocylindral model for dextrans to explain the features of folding from various protein mutants . Nevertheless, these studies have assumed dextrans to interact with a protein through hard-core interactions. Interestingly, studies by Jiao et al . and, recently, by Ebbinghaus et al . showed that dextran stabilizes a protein through enthalpic interactions similar to its building block glucose, instead of stabilizing through volume exclusion. A suitable choice of dextran models again has become elusive from the past literature. A mixture of dextran with destabilizing urea complicates our analysis of apoazuring folding. The net opposing impact from dextran and urea on protein stability can depend on the composition and become nonadditive . Given only the hydrodynamic radius of dextran in 6.8 M urea from the NMR experiments, our modeling approach was to vary the shape of a volume-conserving crowder model for dextran while changing the effective attraction between urea-unfolded apoazurin and its surrounding dextran crowders. We first showed that a naïve crowder model of hard-core repulsion between dextrans and apoazurin fails to capture the experimental observation. The addition of an effective attraction between dextrans and apoazurin is necessary to shift the protein ensemble towards elongated conformations. These effective attractions are unlikely the ones such as electrostatics and hydrogen bonding interactions found in protein crowders . We are attracted to another possible explanation that compressible (physically soft) crowders or chemically soft repulsion produce effective attractions between a protein and a crowder because of the increased depletion between soft repulsive depletants. It is noteworthy that the soft repulsion from Zaccarelli’s work is due to deformation of polymers, which is an entropic effect. In Harries’s work , however, the soft repulsion from cosolute-solvent mixtures, resulting from changes in the hydrogen bonding in solvent near the first hydration shell, is represented by a step-like function in the interacting potential – an enthalpic effect. Our effective attractive potential can rise from either or both factors. The geometry of an elongated rod model for dextrans further shifts the ensemble of apoazurin towards elongated conformations (favoring urea-unfolded conformations). The attraction between urea-unfolded apoazurin and dextrans facilitates dextrans to form aligned bundled conformations. The shape of the dextran breaks the symmetry of the depletion force, thus favoring an elongated geometry for the unfolded apoazurin. With coarse-grained modeling, we provide a plausible explanation to the puzzling discovery from the NMR experiments. II. MATERIALS AND METHODS II.1 Apoazurin
We here used the apo-form
Pseudomonas aeruginosa azurin, termed apoazurin, to study the excluded-volume effects of dextran 20 on the urea-denatured protein in vitro and in silico . It is a bacterial protein that binds a copper ion that can undergo oxidation-reduction as part of electron transfer chains in bacteria. Apoazurin has 128 residues, and the folded state adopts a β -sheet Greek-key fold with a disulfide bond (Cys3-Cys26) (see Figure 1). In vitro studies on apoazurin have advanced our knowledge on the folding mechanism of sandwich-like proteins . Apoazurin has provided an excellent template to compare experimental measurement with the predictions from protein folding theories and simulations . Figure 1 . (a) Representations of the folded structure of apoazurin (PDB ID: 1E65). The folded structure is illustrated in a cartoon representation and a side-chain C α representation. These protein illustrations were produced using VMD . The algorithm DSSP was used to assign the secondary structures. The residues that form the disulfide bond are shown in yellow. (b) Illustrations of two models for dextran 20 with the same volume: a spherical model (DX20-S), and a spherocylinder model with an aspect ratio of 4 (DX20-AR4). As reference, the folded structure of apoazurin is shown. The reduced unit of length is s . II.2 NMR measurements
All NMR experiments were conducted at a Bruker Avance III HD 850 MHz spectrometer equipped with a z-gradient cryogenic probe using
20 mM NaP, pH 7.0 and 5 % D O (v/v) as sample buffer. All one-dimensional data sets were processed and analysed by Topspin 3.2 whereas two-dimensional data sets were processed with NMRPipe and analysed in NMRView . Water suppression was achieved by presaturation and WATERGATE. II.2.1 Diffusion experiments . NMR diffusion spectra were recorded using pulsed field bipolar gradient stimulated echo experiments at T = 298 K. For each diffusion profile, 21 different gradient strengths G were used for 6 ms along the z axis followed by a 100 ms recovery delay. The diffusion of apoazurin was allowed to proceed for 100 ms ( c d20 = 0 g/l, c d20 = 100 g/l, c d20 = 150 g/l, c d20 = 180 g/l, c d20 = 200 g/l and c urea = 0 M; c d20 = 0 g/l, c d20 = 100 g/l and c urea = 6.8 M), 150 ms ( c d20 = 150 g/l and c urea = 6.8 M) or 200 ms ( c d20 = 180 g/l and c urea = 6.8 M). The calibration of G was performed by a standard protocol . For error estimation, four different gradient strengths were repeated (relative gradient strengths of 1, 10, 40 and 70 %). The measured H NMR spectra were integrated within the aliphatic signal region I ( 0.5... 1.5 ppm) and fitted to Equation (1): I ( G )= I (0)exp(- G g d D ( D - d /3)) (1) where g is the gyromagnetic ratio, d is the gradient length, D is the diffusion time and D is the calculated diffusion coefficient . II.2.2 Determination of viscosity.
The viscosity h of the solvent has been determined by measurement of the diffusion coefficient D of the internal reference molecules dioxane present at different conditions ( c d20 = 0 g/l and c urea = 0 M; c d20 = 0 g/l and c urea = 6.8 M) or dextran 20 ( c d20 = 100 g/l, c d20 = 150 g/l, c d20 = 180 g/l, c d20 = 200 g/l and c urea = 0 M; c d20 = 100 g/l, c d20 = 150 g/l, c d20 = 180 g/l and c urea = 6.8 M), known hydrodynamic radii of dioxane, r Hdioxane = 1.87 Å and r Hd20 = 32.4 Å and by applying Equation (2) D = k B T/ (6 ph r H ) (2) with Boltzmann’s constant k B and temperature T . II.3 Coarse-grained models
We have performed coarse-grained molecular simulations of apoazurin under several crowding conditions and under 0 M and 6 M of urea denaturant. 6 M urea is a concentration that guarantees the denaturation of the protein at room temperature.
II.3.1 Coarse-grained apoazurin model . Our model protein apoazurin (PDB ID: 1E65; Figure 1) is studied employing a coarse-grained model. We coarse-grained each residue into a side-chain and C α bead with exception to glycine, which does not contain a side-chain . The C α bead is located at the position of the α -carbon, and the side-chain bead is located at the center of mass of the side-chain atoms. We have employed a structure-based Hamiltonian for our protein model, i.e., the energy landscape of the protein follows the principle of minimal frustration , by following previous work
44, 48 . The dimensions of the protein are expressed in units of s ( s = 3.8 Å), which is the average distance between two consecutive C α beads in a protein. This reduced unit of length is also used for the dimensions of the systems with crowders. II.3.2 Coarse-grained Dextran 20 model . We employed twelve volume-conserving models on dextran 20, a glucose-based polydisperse polymer of 20 kDa molar mass, to investigate the effects of crowder shapes and effective attraction with a protein (Figure 2). The radius of r Hd20 = 32.4 Å was taken from the current experimental measurement. The rod-like shapes used by the models are spherocylinders with various aspect ratios ( AR ). A spherocylindral model is created by connecting several overlapping spheres where the center-to-center distance between adjacent spheres is the radius of a sphere. The AR of a spherocylinder is defined as the ratio L/D , where L is the distance between the centers of the two spheres located at the ends of the spherocylinder, and D is the diameter of a sphere. For hard-core crowder models, the AR varies from 0 (DX20-S), to 1 (DX20-AR1), 2 (DX20-AR2), 4 (DX20-AR4), 8 (DX20-AR8), and 16 (DX20-AR16), while keeping the volume the same as DX20-S with the radius of r Hd20 = 32.4 Å which is equivalent to 8.3 s . This value corresponds to the hydrodynamics radius of dextran 20 from our in vitro experiment, and it is approximately twice the radius of gyration R g of the folded state ( R gN ≈ σ ). The crowders with DX20-series interact with a protein through repulsion, and the crowders with DX20A-series interact with the protein through attractive interactions. Table 1 shows the number of beads and radii of those beads needed to create the various AR s. Name Aspect ratio, AR Number of beads Radii of beads ( σ ) DX20-S and DX20A-S 0 1 8.3 DX20-AR1 and DX20A-AR1 1 3 6.1 DX20-AR2 and DX20A-AR2 2 5 5.3 DX20-AR4 and DX20A-AR4 4 9 4.4 DX20-AR8 and DX20A-AR8 8 17 3.6 DX20-AR16 and DX20A-AR16 16 33 2.9
Table 1
Physical characteristics of the coarse-grained model of dextran 20. The spherocylinders are characterized by aspect ratio AR , number of beads needed to compose the spherocylinder, and the radii of those beads. DX20-series are crowders with repulsive interactions with a protein. DX20A-series are crowders with attractive interactions with a protein. An apoazurin protein is placed in the center of the box of crowders of the same type. The number of the crowders in a box depends on the volume fraction ϕ c under study. The volume fractions studied are ϕ c = 0 % (or bulk), 5 %, 10 %, and 20 %. II.3.3. The energy function of apoazurin under urea interference in a periodic cubic box of crowders.
The Hamiltonian of the system ( E s ) is the sum of the potential energy of the protein ( E p ), the interaction between protein and crowders ( E pc ), and the interaction between crowders ( E cc ), Figure 2 . PFG-NMR based diffusion profiles obtained for native apoazurin in absence (colored in black, closed circles) and in presence of c d20 = 100 g/l (colored in red, open circles), c d20 = 150 g/l (colored in orange, closed rectangles) and c d20 = 200 g/l (colored in blue, open rectangles). The straight lines are fits of Equation (1) to the experimental data. All data were recorded at T = 298 K. Particular gradient strengths have been repeated (1, 10, 40, and 70 % related to maximum gradient strength possible). i n t eg r a l H = (- . ... . ) pp m -1 π i.e., E s = E p +E cc +E pc . E p adopts a variant of structure-based potentials where only the contacts found in the native states are attractive according to the amino acid types, while the non-native contacts are repulsive. We have incorporated the chemical interference by adopting the “Urea-SC-DIH” model from a previous study where the urea solvent-mediated interactions were included in the interactions between side-chain beads and the dihedral-potential that corresponds to four successive C α beads. The expression of E p can be found in the supplementary information. For E cc the interactions between crowders are repulsive and inversely related to the distance between the two crowders to the 12 th power [Equation (S3)]. We set the interactions between protein and crowders, E pc, to be repulsive or attractive. The expression of repulsion between hard spheres of crowders labeled as DX20-series in Table 1 follows the same expression as E cc . The expression of non-specific attraction between a crowder labeled as DX20A-series in Table 1 and a bead from the apoazurin model follows a Lennard-Jones potential with the depth well of ε s . The relative strength of the attractive interaction is denoted by λ = ε s / ε , where 𝜀 = 0.6 kcal/mol is the solvent-mediated interaction between Thr and Thr residues at 0 M. We used λ = 0.83 and 2.50 for this study. II.4. Molecular simulations
We used an in-house version of the Amber 10 molecular dynamics package to integrate the Langevin equation of motion in the low friction limit . The integration time step was set to 10 -4 𝜏 , where 𝜏 = %(𝑚𝜎 ) 𝜀⁄ ) , 𝑚 is the mass of a C α bead, 𝜀 is the solvent-mediated interaction, and 𝜎 is the van der Waals radius of a C α bead. The replica exchange method (REM) was used in order to efficiently sample the configurational space of samples when required. The exchanges between the neighboring replicas i and j are attempted every 40 𝜏 . The acceptance or rejection of each exchange follows the Metropolis criterion, 𝑚𝑖𝑛 /1, 𝑒𝑥𝑝 56𝛽 − 𝛽 : ; ∙ =𝑈(𝑟 ) − 𝑈6𝑟 : ;@AB , where 𝛽 = 1 𝑘 D 𝑇⁄ , 𝑘 D is the Boltzmann constant, T is the temperature, and 𝑈 ( r ) is the potential energy of the system. Twenty replicas were employed where the temperature range used in a REM simulation varies for each condition. For each replica, more than 100 000 statistically significant conformations were collected. For each condition of the simulations with crowders, there are ten sets of randomly assigned crowder configurations in a periodic boundary cubic box (PBC) as initial conditions. The size of PBC is 120 × × σ for the crowders with AR from 0 to 8, and it is 185 × × σ for the crowders with AR = 16. The canonical simulations were performed at a temperature that corresponds to the in silico room temperature T u . T u is obtained by transferring the difference between the melting temperature ( T m ) and the room temperature (298 K) for an in vitro experiment to the in silico computational study as suggested in a prior study . The in silico room temperature is given by 𝑇 F (𝑖𝑛 𝑠𝑖𝑙𝑖𝑐𝑜) = 𝑇 K (𝑖𝑛 𝑠𝑖𝑙𝑖𝑐𝑜) − [𝑇 M (𝑖𝑛 𝑣𝑖𝑡𝑟𝑜) − 298 K]. (3) Where T f is the in silico folding temperature. In our case the in silico room temperature at 0 M urea is k B T u / ε = 1.14. II.5 Data Analysis
II.4.1 Geometrical properties of a protein . We measured the geometry of apoazurin by using the radius of gyration ( R g ) and asphericity parameter ( Δ ) [see Equation (S7)]. R g measures the size of the protein, and Δ measures the degree of anisotropy of the structure of the protein. Δ has the following range: 0 ≤ Δ ≤
1, where a perfect sphere has Δ = 0. When Δ approaches 1, it is an indication that the protein is fully extended. II.4.2 The geometrical properties of a crowder void . The crowder void is defined as the space between the crowders in which the protein lies. To find this void, we partition the entire simulation box through Voronoi tessellation . From the tessellation, a 3D polygon is created around each bead (both crowders and protein residues). We select the crowder beads whose polygon is in contact with the polygon of the protein beads, and the positions of these selected crowder beads are used to define the vertices of the crowder void. We characterize the crowder void by computing two geometrical descriptors: the radius of gyration R g and the asphericity parameter Δ . III. RESULTS
III.1. Probing apoazurin in vitro
III.1.1. Structural changes
It has been established that high-resolution NMR spectroscopy enables the profound molecular characterization of the structural and dynamic properties of a protein ensemble as it is present in solution. We have made use of that by applying two-dimensional heteronuclear H- N NMR experiments on apoazurin present under native and denatured conditions in a macromolecular crowding environment. Moreover, we have used Pulse-field gradient (PFG)-NMR techniques to decode the potential impact of macromolecular crowding on the hydrodynamic dimension of apoazurin present under both folded and denatured conditions. The addition of dextran 20 to apoazurin up to a concentration of c d20 = 200 g/l neither changes the chemical shifts nor the signal intensities of apoazurin as probed by one-dimensional proton (Figure S1) as well as by two-dimensional H- N NMR HSQC spectroscopy (Figure S2). This indicates that the addition of dextran 20 does not modify the structural properties of folded apoazurin at secondary and tertiary levels. The maximum change in H, N weighted chemical shifts, Δω max , is 0.025 ppm with an average Δω mean of (0.012±0.004) ppm. Note that chemical shifts for aliphatic protons of apoazurin are highly conserved as well by comparing diluted and c d20 = 200 g/l conditions for the folded state. Next, dextran 20 of up to c d20 = 180 g/l was additionally supplemented to apoazurin that has been chemically unfolded using c = 6.8 M of urea (Figure S3 and S4). Again, changes in H, N weighted chemical shifts are small, with Δω max = 0.025 ppm and Δω mean = (0.009±0.006) ppm. Chemical shifts for aliphatic protons are conserved in the presence of c urea = 6.8 M by comparing diluted and c d20 = 180 g/l conditions as observed for folded apoazurin. Taken together, 200 g/l of dextran 20 does not change the structural integrity of the folded state and 180 g/l of dextran 20 does not perturb the chemical shift pattern seen for chemically unfolded apoazurin. Consequently, we conclude that there is no significant solute-protein interaction between dextran 20 and apoazurin. III.1.2. Molecular dimension
NMR diffusion methodology has been successfully employed for the profound characterization of the molecular dimensions of both folded and highly denatured proteins in solution . The initial idea with these experiments was to probe the expected (based on excluded volume theory) compaction of unfolded apoazurin when macromolecular crowding agents have been added. First, we probed the diffusion properties of apoazurin in the absence of urea. The diffusion coefficient, D , obtained for folded apoazurin ranges from D = (11.4±0.3) × -11 m s -1 for diluted conditions to D = (2.44±0.07) × -11 m s -1 present at c d20 = 200 g/l (Figure 2). As the solvent viscosity h increases from h = 0.89±0.01 mPa · s to h = 3.99±0.06 mPa · s (Figure 3), the hydrodynamic radius r H for folded apoazurin remains constant at r H = 21.7±0.8 Å. Notably, r H determined for folded apoazurin is in excellent agreement at all concentrations of dextran 20 with the molecular dimension expected for a globular protein of azurin’s size —i.e., r H ( N aa =128) = (20±5 Å) (Figure 4). Figure 4 . The hydrodynamic radius r H of apoazurin depends on the macromolecular crowding conditions used. Whereas native apoazurin does not show a dependence on the concentration of dextran 20 (closed circles, colored in red), denatured apoazurin increases in the hydrodynamic dimension by increasing the concentration of dextran 20 (open circles, colored in blue). The horizontal lines (dotted mode) show the hydrodynamic radii that are theoretically expected for a globular folded protein (colored in red) and a denatured polypeptide chain of that size (colored in blue) according to . All the data were recorded at T = 298 K. Error bars have been included (and seen). The calculation has been done in the same way as for we did Figure 3. Next, we probed the diffusion coefficient for apoazurin under denaturing conditions of c urea r H apoa z u r i n / e - m Figure 3 . The viscosity h of the solvent depends on the concentrations of dextran 20 and urea (absence of urea: closed circles, colored in red; presence of c = 6.8 M urea: open circles, colored in blue) and has been calculated by determining the diffusion of dioxane and dextran 20 and Equation (2) at T = 298 K. The straight lines have been computed by assuming a polynomial function to guide the eye. On page 11 in the caption of Figure 3:… Error bars have been included. They are too small to be seen. They have been calculated based on the standard deviation of the diffusion coefficient obtained by fitting Eq. 1 to the diffusion data. This error of the diffusion coefficient has been used to determine the error of rH by using Eq. 2. -3 η s o l v en t / P a s = 6.8 M. Again, D decreases with an increasing concentration of dextran 20 from D = (5.22±0.08)e-11 m s -1 ( c d20 = 0 g/l) to D = (0.62±0.06)e-11 m s -1 ( c d20 = 180 g/l). However, the decrease in D was not compensated by the increase in h as found under native conditions. Instead, surprisingly, the hydrodynamic radius for denatured apoazurin increases from r H = 35.8±0.6 Å ( c d20 = 0 g/l) to r H = 46.6±2.5 Å ( c d20 = 180 g/l) (Figure 4). Notably, r H , as determined here for denatured apoazurin under dilute conditions, agrees with the molecular dimensions expected for a denatured polypeptide chain of apoazurin’s size , r H ( N aa =128)=(35±18 Å). Because excluded volume theory predicts that unfolded proteins become compacted in the presence of macromolecular crowding, more factors must be at play. To find a molecular-mechanistic explanation to this in vitro observation, we turned to in silico computations. III.2. In silico analysis of apoazurin unfolding To compare our simulation results to the experimental measurements of apoazurin’s hydrodynamic radii r H under different conditions, we used two order parameters: radius of gyration R g and asphericity ∆ . R g is not an exact comparison to r H because the shape of the object also matters in determining r H . For a perfect sphere, r H can be approximated by R g ; however, as the protein becomes less spherical (or an increase in ∆ ), r H and R g diverge (see Figure S5). r H / R g can exceed 3 when ∆ approaches 1 (a rod-like shape). This is due to the diffusion of spherical molecules being different than that of other elongated molecules. Therefore, r H is in fact measuring a combination of the size and shape of a molecule. III.2.1. Structural changes for the repulsive crowder model and chemical denaturation
We have investigated the ensemble conformations of apoazurin in various crowding contents at a chemically denatured state at 6 M urea and at an in silico room temperature of k B T u / ε = 1.14 (see Methods) with coarse-grained molecular simulations. First, we have assumed that crowders and apoazurin interact through a repulsive potential that accounts for the excluded volume effects at ϕ c = 5 %, 10 %, and 20 %. While conserving the volume of the crowders, we tested the shape effect by varying the aspect ratio of the crowders on apoazurin unfolding. Figure 5(a) shows the average asphericity ( Δ ) normalized by the asphericity ( Δ B ) under the bulk condition for DX20-S, DX20-AR4 and DX20-AR8 against ϕ c . Δ / Δ B are all less than 1 over a wide range of ϕ c , implying that all dextran 20 models with repulsive interactions with apoazurin do not contribute to the elongation of apoazurin. For most of the conditions, the macromolecular crowding effect leads to compaction of apoazurin with R g / R g,B less than 1 [Figure 5(b)]. The crowders with a rod-like geometry show greater compaction than the spherical crowder over all ϕ c because rod-like crowders render a larger covolume than spherical crowders. Because the effects found upon including only volume exclusion cannot explain the experimental results of dextran 20-induced elongation of unfolded apoazurin, we next investigated a dextran 20 model that also provides effective attractive interactions. Figure 5 . (a) Average asphericity ( Δ ) of the protein apoazurin in varying shapes of crowders normalized by the average asphericity ( Δ B ) of apoazurin under the bulk condition against the volume fractions of crowders ( ϕ c ). (b) Average radius of gyration ( R g ) of apoazurin normalized by the average radius of gyration under the bulk condition against ϕ c . The interactions between the protein apoazurin and crowders are repulsive. The concentration of urea is c urea = 6 M. The temperature is k B T u / ε = 1.14, corresponding to room temperature (see Method), for all conditions. Error bars are included. III.2.2. Structural changes of apoazurin with which crowders interact attractively and under chemical denaturation
We hypothesize the existence of non-specific attractive interactions between the unfolded protein and crowders. For this, we tested two types of non-specific attractive interactions with λ = 0.83 and λ = 2.50 (see methods). λ measures the attraction between a protein and crowder relative to the attraction between two Thr residues at 0 M urea. λ = 0.83 represents a weak attraction, and λ = 2.50 represents a strong attraction. We noticed that a box of spherical crowders with a weak attraction (DX20A-S at λ = 0.83) cannot produce elongated ensembles of apoazurin in silico at room temperature with k B T u / ε = 1.14 in Figure 6. The distribution of D for apoazurin, centered at D ≈ ϕ c . We note that, with a strong attraction at λ = 2.5, there was a slight increase in Δ / Δ B or R g / R g , B along ϕ c (Figure S6). However, R g / R g , B is still less than 1 when the aspect ratio of the crowder is less than 4. This result leads us to a working hypothesis that, in addition to weak attractions, the crowders must be rod-like to stabilize an extended unfolded conformation of apoazurin. We showed that, in the system of weakly attractive crowders with a rod-like shape (DX20A-AR4), the distribution of D for the protein apoazurin under the chemically denatured condition becomes broad and comparable to that under the bulk condition in Figure 6. Interestingly, the distribution of D for apoazurin shifts toward D = 1 and forms a large peak at D ≈ ϕ c = 5 % for DX20A-AR8 (Figure 7). This finding implies that more elongated conformations of apoazurin in chemically unfolded states are sampled in a system with the crowder model DX20A-AR8 than with the crowder model DX20A-AR4. However, the height and position of the peak shifts toward D = 0 with increasing ϕ c . At a high ϕ c of 20 %, apoazurin adopts a compact configuration again (Figure 8). This feature was not observed for spherical crowders (DX20A-S) with a weakly attractive interaction with a protein at λ = 0.83 (Figure S6). Figure 6 . Distribution of the asphericity ( D ) at c urea = 6 M in bulk, with DX20A-AR0, and with DX20A-AR4 with a strength of attractive interactions λ = 0.83 at a wide range of volume fraction of crowders ( ϕ c ). The temperature is k B T u / ε = 1.14. Figure 7 . Distribution of the asphericity ( D ) for chemically unfolded apoazurin at 6 M urea in bulk and in a system of weakly attractive crowders (DX20A-AR8 at λ = 0.83) over several volume fractions ( ϕ c ). The temperature for all conditions is k B T u / ε = 1.14. Two representative structures are shown: at the left bottom with D = 0.14 ( ϕ c = 20 %) and the right bottom with D = 0.83 ( ϕ c = 5 %). Figure 8 . (a) Ratio of the average of the asphericity ( Δ ) under a specific condition and average of the asphericity ( Δ B ) under the bulk condition. (b) Ratio of the average of the radius of gyration ( R g ) under a specific condition and average of the radius of gyration ( R g,B ) under the bulk condition. These ratios are shown as a function of the aspect ratio of crowders, AR . The concentration of urea is 6 M, and the temperature is k B T u / ε = 1.14 for all conditions. Error bars are included. III.2.3. Both the rod-like shape of a crowder and attractive interaction with a protein contribute to the structural changes of apoazurin under chemically denatured conditions.
We evaluated Δ / Δ B as well as R g / R g,B of chemically denatured apoazurin at varying aspect ratios ( AR ) of crowders for a wide range of ϕ c in Figure 8. Overall, there is an increase in Δ / Δ B when AR increases [Figure 8(a)]. Δ / Δ B exceeds 1 at AR greater than 8. Interestingly, R g / R g , B also exceeds 1 at the most extreme case of AR = 16 [Figure 8(b)] at ϕ c = 5 %. However, we noted that the increasing trends in Δ / Δ B or R g / R g , B are not uniform with ϕ c s. When AR is less than
4, the protein has the largest Δ and R g values at ϕ c = 20 % compared with the other ϕ c conditions. When AR is greater than 4, the largest Δ and R g values occur at ϕ c = 5 %. Such a crossover behavior at AR = 4 is contributed by the shape of the protein-crowder assembly (Figure 9). The unfolded ensemble of apoazurin conformations adopts a shape according to the available void formed by surrounding crowders [Figure 9(a)]. A rod-like crowder with non-specific attractive interactions provides an elongated surface where several regions of the protein can stick [Figure 9(b)]. This particular geometry of rod-like crowder favors the formation of extended unfolded conformations of apoazurin. As several crowders stick to a protein, they tend to align themselves and restrict the available space for the protein. It appears that, at high ϕ c , there is insufficient space to accommodate the volume of a large assembly of protein crowders. Therefore, the uptick of Δ / Δ B or R g / R g , B at high ϕ c is weaker than that at low ϕ c . Figure 9 . Snapshots of an exemplarily unfolded structure of apoazurin and the void between the protein and closest crowders for (a) spherical crowders (DX20A-S) (b) and rod-like crowders (DX20A-AR8) at ϕ c = 5 % and at c urea = 6 M with weakly attractive interactions ( λ = 0.83). III.2.3. The geometry of crowder-induced voids is modulated by the protein-crowder interactions
Next, we investigated the geometry of voids, in terms of shape and size, formed by crowders of the various models modulated by the protein-crowder interactions. The voids are the depletion zones due to non-overlapping volume exclusion between a protein and crowders. The distribution of the asphericity of the void D void is virtually unchanged between spherical and rod-like crowders without attractive interactions [Figure 10(a)]. Additionally, as ϕ c increases, the peaks of the distributions are unperturbed. These results indicate that, without attractive interactions, the variation in the shape effect from rod-like crowders has been averaged in the ensemble. However, when attractive interactions are added to the crowders, there is a clear shift rightward (towards higher asphericity) of the peak. The spread of the distributions increases with AR when the shape of a crowder becomes elongated [Figure 10(b)]. These effects are greatest for ϕ c = 5 %. We show that the distribution of the size of a void ( R gvoid ) depends on the volume fraction of crowders, the shape of a crowder, and interaction between the protein and crowder. To make a fair comparison between spherical and rod-like crowders, we plotted the distribution of ( R gvoid – R )/ σ in Figure 10(c) and Figure 10(d), where R is the radius of a sphere in a crowder model. The overall size of ( R gvoid – R )/ σ is smaller for rod-like crowders (DX20-AR8) than spherical crowders (DX20-S) by roughly a third at each ϕ c when a protein and crowders interact through hard-core interaction [Figure 10(c)]. However, with a weak attraction between a protein and crowders ( λ = 0.83), ( R gvoid – R )/ σ is reduced by half across all crowder types [Figure 10(d)] compared to those in Figure 10(c). The peaks in the distribution of ( R gvoid – R )/ σ for rod-like crowders (DX20A-AR8) at all ϕ c nearly overlap at the same position around 10. This feature is consistent with the finding of a protein-crowder assembly in Figure 9. Figure 10 . Distribution of the asphericity of the void ( Δ void ) between the protein and crowders with (a) steric repulsive interactions and (b) non-specific attractive interactions between the protein and crowders. The distribution of the radius of gyration of the void ( R gvoid - R )/ σ between the protein and crowders with (c) steric repulsive interactions and (d) weakly attractive interactions between the protein and crowders. IV. DISCUSSION
This in silico investigation was initiated by the anomalous observation from NMR experiments that the addition of 180 g/l of dextran 20 alters the ensemble of urea-unfolded apoazurin in the direction of an extended conformation. Initially, this observation contrasts the physical intuition of the excluded volume theory that proteins generally favor a collapsed, compact conformation in the presence of macromolecular crowders . This intuition is built on a hard-sphere crowder (or depletant) that exerts entropic depletion interactions, or the macromolecular crowding effect , on the protein of interest, that effectively makes the unfolded expanded states of a protein entropically unfavorable, thus favoring a collapsed, compact conformation over an unfolded, extended one . The strength of the macromolecular crowding effect grows with the volume fraction of crowders and is often strongest when the radius of the crowders is smaller than the R g of the test protein . However, recent studies have shown that the sequence-dependent properties of a protein from protein crowders can offset the volume exclusion of protein crowders
24, 59 . The sequence dependent interactions between a protein and protein crowder include chemically soft interactions such as hydrogen bonding, electrostatic interactions, or dispersion. In some cases, these chemical attractions that possibly stick a protein to the surface of crowders further form quinary structures and destabilize the native folded state . Chemically inert crowders such as Ficoll 70, PEG, dextrans do not share the same chemical surfaces as protein crowders. However, a spherocylinder model for dextrans that lacks atomistic details and interacts with a protein through hard-core potentials fails to explain the anomaly from experiments, unless we consider the effective attractions between a protein and dextrans. A recent study by Zhou et al . also suggests that synthetic crowders may attain strong attractions with intrinsically disordered proteins (IDPs) such that IDP polypeptides populate expanded configurations under crowded conditions . The underlying reasons may include soft repulsion between the protein and synthetic crowders, deepening the depletion in the potential of mean force . This physically soft aspect of an interaction (in contrast to chemically soft as mentioned above) is due to the compressibility of soft-matter materials that does not exist in non-overlapping hard spheres. It is possible that, in the mixture of urea and dextran, dextran becomes physically soft. This speculation is supported by a recent essay from Ferreira et al . that, in the mixture of osmolyte and a macromolecular crowder, the solvent properties may not always be additive, depending on the composition of the mixture. We speculate that, at 6 M urea, not only the solvent-mediated interactions between amino acids in the protein are weakened, leading to unfolding, but dextrans also become softened to trap unfolded apoazurin conformations. This scenario can be imagined as spaghetti trapped by a fork. This “fork-spaghetti” mechanism may be facilitated by the non-uniform length of the branches of dextran . The geometric feature of branched dextrans is at odds with the dimensions of a linear dextran model , especially when the molecular weight increases . Noteworthily, the behavior of branched polymers (dextran) is distinct from those of linear polymers (PEG) with a similar molecular weight shown by Soranno et al . who investigates the macromolecular crowding effect on the compaction of intrinsically disorder peptides (IDPs). When large PEG polymers reach a semi dilute region, they unwind and overlap one another at a moderate volume fraction of 10%. The monomers of overlapping PEG polymers still exert volume exclusion effect on the compaction of intrinsically disorder peptides (IDPs). However, the strength of volume exclusion from PEG monomers are not as strong as non-overlapping PEG polymers as predicted from the scaled particle theory . Notably, we showed that the effective attraction alone is not sufficient to bias the ensemble of unfolded apoazurin towards an elongated one as the shape of the crowder must be cylindrical with a sufficiently large aspect ratio. With spherical crowders containing either weak or strong effective attraction, the enhanced depletion interactions still favor a compact unfolded apoazurin configuration [Figure 8(b)]. This contrasts the findings from Zhou’s work in which peptide model samples with both extended and compact conformations interacted strongly with spherical crowders . We speculate that the biphasic distribution might change if the length of the simulation is increased or if a more efficient sampling scheme is introduced. Because both the aspherical shape of a crowder and its effective attraction with a protein increase the correlation time in the simulations, we ensured that our conclusions are independent of the sampling scheme in our simulations. Our results at a volume fraction of 5 % for the DX20A-AR8 model at λ = 0.83 do not depend on the initial unfolded conformation of apoazurin. The temporal evolution of D is shown in Figure 11 at the early stage of two canonical simulations: one starting from a compact unfolded conformation and the other from an extended unfolded conformation. In the first case, D fluctuates and rearranges initially, and then it increases as time increases. In the second case, the conformation becomes compact from a stretched one, and then it rearranges to a have a less extended shape. Thus, the initial conformation is not the reason why the simulations favor the formation of extended conformations with crowders and non-specific attractive interactions. Figure 11 . Initial stage of two canonical simulations for the model DX20A-AR8 at k B T u / ε = 1.14: one starting from (a) a compact configuration and the other from (b) an extended configuration. The snapshots of the protein depict the temporal structural changes. To explain our in silico results, the elongated shape of a crowder must break the symmetry of the isotropic direction of the depletion force, creating an anisotropic force that sandwiches the backbone of the protein. We found that the elongated conformation of apoazurin is stabilized by a bundle of rod-like crowders in the presence of effective attraction. The rod-like crowders are still free to rotate, and there is no evidence of an isotropic-nematic transition in our system. In the absence of effective attraction, the repulsive cylinders do not self-assemble into nematic states. To reach a nematic state for rod-like crowders with an aspect ratio of 5, ϕ c typically must exceed 40 % . In fact, hard spherocylinders contribute to greater compactness of apoazurin because of an increased covolume compared with a volume-conserving hard sphere . The difference here is the effective attraction of apoazurin and several rod-like crowders that collectively form a bundle. The bundles formed around elongated apoazurin maximize the volume accessible to the other crowders (not in a bundle), thus maximizing the total entropy of the system. We showed that elongated apoazurin appears at low ϕ c = 5 %, but not at high ϕ c = 20 %; the latter was observed in in vitro experiments. This difference may be due to the increase in excluded volume interactions between hard spherocylinders instead of the spherocylinder-apoazurin interactions, as the volume fraction of spherocylinders increases. To induce the formation of significant bundles at high ϕ c , the concentration of apoazurin in a binary mixture of spherocylinder and apoazurin should be increased. V. CONCLUSION
New understanding in science generally begins when theoretical expectations and new experimental findings do not agree. Here, the intriguing NMR discovery showing elongated urea-unfolded states in the presence of dextran goes against intuition from conventional macromolecular crowding theory and previous experimental findings. With very little clues from experimental data, we have employed coarse-grained modeling to investigate the underlying mechanism. Even though we cannot definitively exclude other possible explanations for the trends found in the experiments, our work denotes the importance of both crowder shapes and solvent-mediated interactions in manipulating the geometry of a protein. These two properties are pivotal to understand the full complexity of proteins in the cytoplasm. Furthermore, this newfound understanding can potentially be developed into a design principle for probes to accurately measure crowding volume fraction, shapes, and interactions locally in a subcellular environment.
Supporting Information
Hamiltonian of the coarse-grained protein model in the presence of urea, NMR measurements (Figure S1-S4), definition of the asphericity ∆ , ratio of r H and R g as a function of ∆ (Figure S5), changes in the geometry of the protein as a function of volume fraction of crowders (Figure S6), and protein stability at different concentrations of urea (Figure S7) are presented in the supporting information (SI). Acknowledgements
We thank Alexander Christiansen for providing apoazurin and Birgit Köhn for technical assistance. We acknowledge the funding support from the National Science Foundation (MCB-1412532 and ACI-1531814) and the computational resources from the Center for Advanced Computing and Data Science (CACDS) and Research Computing Center (RCC) at the University of Houston. MSC thanks Yossi Eliaz for his intellectual discussion. AGG is supported by a fellowship from the Houston Area Biophysics Program (T32 GM008280). PWS acknowledges funding from Knut and Alice Wallenberg Foundation and the Swedish Research Council. MK acknowledges financial support from the Young Scholar Fund of the Universität Konstanz. Author contribution statement
M.K., P.W.S., and M.S.C. designed research, F.Z., D.H., A.G., L.B., and M.K. performed research work and did data analysis. All authors contributed to discussions and data interpretation. M.S.C., P.W.S., M.K. D.H., F.Z., and A.G. wrote the manuscript.
Conflict of interest statement
There is no conflict of interest.
References
TOC Graphic (1) Ellis, R. J.; Minton, A. P., Join the crowd.
Nature , 27-28. (2) Link, A. J.; Robison, K.; Church, G. M., Comparing the predicted and observed properties of proteins encoded in the genome of Escherichia coli K-12.
Electrophoresis , 1259-1313. (3) Luby-Phelps, K., Cytoarchitecture and physical properties of cytoplasm: Volume, viscosity, diffusion, intracellular surface area. Int. Rev. Cytol. , 189-221. (4) Minton, A. P., Excluded volume as a determinant of macromolecular structure and reactivity.
Biopolymers , 2093-2120. (5) Venturoli, D.; Rippe, B., Ficoll and dextran vs. globular proteins as probes for testing glomerular permselectivity: Effects of molecular size, shape, charge, and deformability. Am. J. Physiol. Renal Physiol. , F605-F613. (6) Homouz, D.; Perham, M.; Samiotakis, A.; Cheung, M. S.; Wittung-Stafshede, P., Crowded, cell-like environment induces shape changes in aspherical protein.
Proc. Natl. Acad. Sci. U. S. A. , 11754-11759. (7) Christiansen, A.; Wang, Q.; Samiotakis, A.; Cheung, M. S.; Wittung-Stafshede, P., Factors defining effects of macromolecular crowding on protein stability: An in vitro/in silico case study using cytochrome c.
Biochemistry , 6519-6530. (8) Dhar, A.; Samiotakis, A.; Ebbinghaus, S.; Nienhaus, L.; Homouz, D.; Gruebele, M.; Cheung, M. S., Structure, function and folding of phosphoglycerate kinase are strongly perturbed by macromolecular crowding. Proc. Natl. Acad. Sci. U. S. A. , 17586-17591. (9) Gupta, A. N.; van der Maarel, J. R. C., Compaction of plasmid DNA by macromolecular crowding. Macromolecules , 1666-1671. (10) Zinchenko, A.; Berezhnoy, N. V.; Chen, Q.; Nordenskiold, L., Compaction of single-molecule megabase-long chromatin under the influence of macromolecular crowding. Biophys. J. , 2326-2335. (11) Christiansen, A.; Wittung-Stafshede, P., Quantification of excluded volume effects on the folding landscape of Pseudomonas aeruginosa Apoazurin in vitro.
Biophys. J. , 1689-1699. (12) Sasahara, K.; McPhie, P.; Minton, A. P., Effect of dextran on protein stability and conformation attributed to macromolecular crowding.
J. Mol. Biol. , 1227-1237. (13) Christiansen, A.; Wittung-Stafshede, P., Synthetic crowding agent dextran causes excluded volume interactions exclusively to tracer protein apoazurin.
FEBS Lett. , 811-814. (14) Yu, I.; Mori, T.; Ando, T.; Harada, R.; Jung, J.; Sugita, Y.; Feig, M., Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm. eLife , e19274. (15) McGuffee, S. R.; Elcock, A. H., Diffusion, crowding and protein stability in dynamic molecular model of bacterial cytoplasm. PLoS Comput. Biol. , e1000694. (16) Dong, H.; Qin, S.; Zhou, H.-X., Effects of macromolecular crowding on protein conformational changes. PLoS Comput. Biol. , e1000833. (17) Qin, S. B.; Minh, D. D. L.; McCammon, J. A.; Zhou, H.-X., Method to predict crowding effects by postprocessing molecular dynamics trajectories: Application to the flap dynamics of HIV-1 protease. J. Phys. Chem. Lett. , 107-110. (18) Cheung, M. S.; Klimov, D.; Thirumalai, D., Molecular crowding enhances native state stability and refolding rates. Proc. Natl. Acad. Sci. U. S. A. , 4753-4758. (19) Mittal, J.; Best, R. B., Dependence of protein folding stability and dynamics on the density and composition of macromolecular crowders.
Biophys. J. , 315-320. (20) Wang, Q.; Cheung, M. S., A physics-based approach of coarse-graining the cytoplasm of E. coli (CGCYTO). Biophys. J. , 2353-2361. (21) Kudlay, A.; Cheung, M. S.; Thirumalai, D., Influence of the shape of crowding particles on the structural transitions in a polymer.
J. Phys. Chem. B , 8513-8522. (22) Granath, K. A., Solution properties of branched dextrans.
J. Colloid Sci. , 308-328. (23) Homouz, D.; Stagg, L.; Wittung-Stafshede, P.; Cheung, M. S., Macromolecular crowding modulates folding mechanism of α / β protein apoflavodoxin. Biophys. J. , 671-680. (24) Jiao, M.; Li, H.-T.; Chen, J.; Minton, A. P.; Liang, Y., Attractive protein-polymer interactions markedly alter the effect of macromolecular crowding on protein association equilibria. Biophys. J. , 914-923. (25) Senske, M.; Tork, L.; Born, B.; Havenith, M.; Herrmann, C.; Ebbinghaus, S., Protein stability by macromolecular crowding through enthalpy rather than entropy. J. Am. Chem. Soc. , 9036-9041. (26) Ferreira, L. A.; Uversky, V. N.; Zaslavsky, B. Y., Role of solvent properties of water in crowding effects induced by macromolecular crowding agents and osmolytes.
Mol. BioSyst. , 2551-2563. (27) Gorensek-Benitez, A. H.; Smith, A. E.; Stadmiller, S. S.; Perez Goncalves, G. M.; Pielak, G. J., Cosolutes, crowding, and protein folding kinetics. J. Phys. Chem. B , 6527-6537. (28) Mu, X.; Choi, S.; Lang, L.; Mowray, D.; Dokholyan, N. V.; Danielsson, J.; Oliveberg, M., Physicochemical code for quinary protein interactions in Escherichia coli.
Proc. Natl. Acad. Sci. U. S. A. , E4556-E4563. (29) Rovigatti, L.; Gnan, N.; Parola, A.; Zaccerelli, E., How soft repulsion enhances the depletion mechanism.
Soft Matter , 692-700. (30) Sapir, L.; Harries, D., Origin of enthalpic depletion forces. J. Phys. Chem. Lett. , 1061-1065. (31) Sapir, L.; Harries, D., Is the depletion force entropic? Molecular crowding beyond steric interactions. Curr. Opin. Colloid Interface Sci. , 3-10. (32) Adman, E. T., Copper protein structures. Adv. Protein Chem. , 145-197. (33) Wilson, C. J.; Apiyo, D.; Wittung-Stafshede, P., Solvation of the folding-transition state in Pseudomonas aeruqinosa azurin is modulated by metal. Protein Sci. , 843-852. (34) Zong, C. H.; Wilson, C. J.; Shen, T. Y.; Wolynes, P.; Wittung-Stafshede, P., φ -value analysis of apo-azurin folding: Comparison between experiment and theory. Biochemistry , 6458-6466. (35) Wang, Q.; Christiansen, A.; Samiotakis, A.; Wittung-Stafshede, P.; Cheung, M. S., Comparison of chemical and thermal protein denaturation by combination of computational and experimental approaches. II J. Chem. Phys. , 175102. (36) Humphrey, W.; Dalke, A.; Schulten, K., VMD: Visual molecular dynamics.
J. Mol. Graphics Modell. , 33-38. (37) Kabsch, W.; Sander, C., Dictionary of protein secondary structure: Pattern recognition of hydrogen - bonded and geometrical features. Biopolymers , 2577-2637. (38) Delaglio, F.; Grzesiek, S.; Vuister, G. W.; Zhu, G.; Pfeifer, J.; Bax, A., NMRPipe: A multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR , 277-293. (39) Johnson, B. A.; Blevins, R. A., NMR View: A computer program for the visualization and analysis of NMR data. J. Biomol. NMR , 603-614. (40) Jones, J. A.; Wilkins, D. K.; Smith, L. J.; Dobson, C. M., Characterisation of protein unfolding by NMR diffusion measurements. J. Biomol. NMR , 199-203. (41) Berger, S.; Braun, S.,
200 and more NMR experiments: A practical course . Wiley-Vch Weinheim: 2004. (42) Wilkins, D. K.; Grimshaw, S. B.; Receveur, V.; Dobson, C. M.; Jones, J. A.; Smith, L. J., Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques.
Biochemistry , 16424-16431. (43) Armstrong, J. K.; Wenby, R. B.; Meiselman, H. J.; Fisher, T. C., The hydrodynamic radii of macromolecules and their effect on red blood cell aggregation. Biophys. J. , 4259-4270. (44) Cheung, M. S.; Finke, J. M.; Callahan, B.; Onuchic, J. N., Exploring the interplay between topology and secondary structural formation in the protein folding problem. J. Phys. Chem. B , 11193-11200. (45) Bryngelson, J. D.; Wolynes, P. G., Spin glasses and the statistical mechanics of protein folding.
Proc. Natl. Acad. Sci. U. S. A. , 7524-7528. (46) Onuchic, J. N.; Luthey-Schulten, Z.; Wolynes, P. G., Theory of protein folding: The energy landscape perspective. Annu. Rev. Phys. Chem. , 545-600. (47) Sheinerman, F. B.; Brooks III, C. L., Calculations on folding of segment B1 of streptococcal protein G. J. Mol. Biol. , 439-456. (48) Samiotakis, A.; Cheung, M. S., Folding dynamics of Trp-cage in the presence of chemical interference and macromolecular crowding. I.
J. Chem. Phys. , 175101 (49) Case, D. A.; Darden, T. A.; Cheatham, T. E.; Simmerling, C. L.; Wang, J.; Duke, R. E.; Luo, R.; Crowley, M.; Walker, R. C.; Zhang, W., et al.
AMBER 10 , University of California, San Francisco, 2008. (50) Veitshans, T.; Klimov, D.; Thirumalai, D., Protein folding kinetics: Timescales, pathways, and energy landscapes in terms of sequence-dependent properties.
Folding Des. , 1-22. (51) Sanbonmatsu, K. Y.; Garcia, A. E., Structure of Met-Enkephalin in explicit aqueous solution using replica exchange molecular dynamics. Proteins: Struct., Funct., Bioinf. , 225-234. (52) Sugita, Y.; Okamoto, Y., Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. , 141-151. (53) Dima, R. I.; Thirumalai, D., Asymmetry in the shapes of folded and denatured states of proteins.
J. Phys. Chem. B , 6564-6570. (54) Barber, C. B.; Dobkin, D. P.; Huhdanpaa, H., The quickhull algorithm for convex hulls.
ACM Trans. Math. Softw. , 469-483. (55) Dyson, H. J.; Wright, P. E., How does your protein fold? Elucidating the apomyoglobin folding pathway. Acc. Chem. Res. , 105-111. (56) Kovermann, M.; Rogne, P.; Wolf-Watz, M., Protein dynamics and function from solution state NMR spectroscopy. Q. Rev. Biophys. , E6. (57) Minton, A. P., Models for excluded volume interaction between an unfolded protein and rigid macromolecular cosolutes: Macromolecular crowding and protein stability revisited Biophys. J. , 971-985. (58) Kang, H.; Pincus, P. A.; Hyeon, C.; Thirumalai, D., Effects of macromolecular crowding on the collapse of biopolymers. Phys. Rev. Lett. , 068303. (59) Miklos, A. C.; Sarkar, M.; Wang, Y. Q.; Pieak, G. J., Protein crowding tunes protein stability.
J. Am. Chem. Soc. , 7116-7120. (60) Gnutt, D.; Timr, S.; Ahlers, J.; Konig, B.; Manderfeld, E.; Heyden, M.; Sterpone, F.; Ebbinghaus, S., Stability effect of quinary interactions reverse by single point mutations.
J. Am. Chem. Soc. , 4660–4669. (61) Cohen, R. D.; Pielak, G. J., Electrostatic contributions to protein quinary structure.
J. Am. Chem. Soc. , 13139–13142. (62) Banks, A.; Qin, S.; Weiss, K. L.; Stanley, C. B.; Zhou, H.-X., Intrinsically disorder protein exhibits both compaction and expansion under macromolecular crowding.
Biophys. J. , 1067-1079. (63) Ebert, K. H.; Brosche, M., Origin of branches in native dextrans.
Biopolymers , 423-430. (64) Soranno, A.; Koenig, I.; Borgia, M. B.; Hoffman, H.; Zosel, F.; Nettles, D.; Schuler, B., Single-molecule spectroscopy reveals polymer effects of disordered proteins in crowded environments. Proc. Natl. Acad. Sci. U. S. A. , 4874-4879. (65) Lago, S.; Cuetos, A.; Martinez-Haya, B.; Rull, L. F., Crowding effects in binary mixtures of rod-like and spherical particles. J. Mol. Recognit. , 417-425. Supplementary Information
Crowding-induced Elongated Conformation of Urea-unfolded Apoazurin: Investigating the Role of Crowder Shape
In Silico
Fabio C. Zegarra &1 , Dirar Homouz &1, 2, 3 , Andrei G. Gasic , Lucas Babel , Michael Kovermann , Pernilla Wittung-Stafshede , Margaret S. Cheung & Equal contribution 1. Department of Physics, University of Houston, USA 2. Department of Physics, Khalifa University of Science and Technology, Abu Dhabi, UAE 3. Center for Theoretical Biological Physics, Rice University, USA 4. Department of Chemistry, Universität Konstanz, Germany 5. Research School Chemical Biology (KoRS-CB), Universität Konstanz, Germany 6. Biology and Biological Engineering Department, Chalmers University, Sweden Corresponding author, email: *Dr. Margaret Cheung [email protected] *Dr. Pernilla Wittung-Stafshede [email protected] *Dr. Michael Kovermann [email protected]
1. The Hamiltonian of a protein at varying concentration of urea
We adopted the “Urea-SC-DIH” model from a previous study where solvent-mediated interactions ε α ij for a particular concentration of urea ( α ) were used in the native side-chain interactions between the side-chain beads i and j and in the dihedral angle formed by four successive C α beads. 𝐸 " = 𝐸 $%&’(%’&)*+’&,) + 𝐸 ./0123+’&,) +𝐸 ./0141 + 𝐸 ./567 , (S1) 𝐸 $%&’(%’&)*+’&,) = 𝐸 + 𝐸 ):<*, + 𝐸 ;.=,;&)*+’&,) + 𝐸 (=.&)* . (S2) The potential energy of the protein ( E p ) is the sum of a structural term ( E structural-urea ) and several non-bonded terms. The list of native pair contacts between the i th and j th side-chain beads was established by the CSU software . The structural potential term ( E structural-urea ) includes a bond-length potential, bond-angle potential, dihedral-angle potential modulated by the concentration of urea, and a chiral potential term that considers the L-form of every side-chain . The equilibrium values, used to estimate the potential terms, were taken from the crystal structure of apoazurin (PDB ID: 1E65) . The disulfide bond between the side-chain beads of Cys3 and Cys26 was modeled by a bond-length potential. The non-bonded potential term considers the native and non-native interactions among residues. The potential term for the native pairs consists of the side-chain-side-chain potential term ( 𝐸 ./0123+’&,) ) and the backbone hydrogen-bonding potential term ( 𝐸 ./0141 ) . The repulsive potential term ( 𝐸 ./567 ) is regarded for interactions between C α beads as well as between non-native side-chain pairs. The repulsive potential term is given by Equation (S3) where 𝜀 = 0.6 kcal/mol and 𝜎 ./ = 0.9C𝜎 . + 𝜎 / D , 𝜎 . and 𝜎 / are the Van der Waals radii of C α or side-chain beads. 𝐸 ./567 = 𝜀 F G HI & HI J KL . (S3) The influence of urea at a certain concentration was introduced in the Hamiltonian by using two terms as indicated by the protocol “Urea-SC-DIH ” from our prior study : the side-chain-side-chain potential term and the dihedral-angle potential term. The side-chain-side-chain potential ( 𝐸 ./0123+’&,) ) for a specific concentration of urea ( α ) is given by a Lennard-Jones potential in Equation (S4) where the solvent-mediated interactions ( 𝜀 ./M ) were derived from all-atom molecular dynamics simulations using a Boltzmann inversion approach . α represents a solution with a concentration of urea of 0M, 6M, or 8M. In the absence of urea (i.e. 0M), the solvent-mediated interactions ( 𝜀 ./M ) for the non-bonded interactions were obtained from the Betancourt–Thirumalai potential , in which case the solvent-mediated interaction between Thr and Thr residues was used to obtain the solvent-mediated interactions of other residue pairs. The distance between the side-chain beads from two different residues is represented by 𝑟 ./ , and the equilibrium position between two beads is given by 𝜎 ./ = 0.9C𝜎 . + 𝜎 / D , where 𝜎 . and 𝜎 / are the radii of the side-chain beads i and j, respectively. 𝐸 ./0123+’&,) = 𝜀 ./M OF G HI & HI J KL − 2 F G HI & HI J R S . (S4) The dihedral-angle potential ( 𝐸 ;.=,;&)*+’&,) ) for a specific concentration of urea ( α ) is given by Equation (S5) where the dihedral force constant depends on the average of the solvent-mediated interactions of four successive C α beads with indexes i , j , k , and l . In the absence the urea, 𝜀̅ = 𝜀 . 𝐸 ;.=,;&)*+’&,) = ∑ 𝜀̅[1 − cos(𝜙 − 𝜙 ] )] + 0.5𝜀̅a1 − cosC3(𝜙 − 𝜙 ] )Dc dH +3 dI +3 de +3 df ;.=,;&)*$ , (S5) where 𝜀̅ is given by 𝜀̅ = C𝜀 ./M + 𝜀 .gM + 𝜀 .*M + 𝜀 /gM + 𝜀 /*M + 𝜀 g*M D 6⁄ . (S6) The interaction between crowder and crowder ( E cc ) of the spherical model of dextran 20 (DX20-S or DX20A-S) is given by a repulsive term as described by Equation (S3). E cc of a spherocylinder crowder (DX20-series or DX20A-series) is composed of the potential energy of crowders and the repulsive interaction between crowders. The potential term for the interaction between crowder and crowder is also described by Equation (S3), and the intra-potential energy of a spherocylinder crowder consists of a bond-length potential, bond-angle potential, and a dihedral-angle potential, which equilibrium values were set to the radius a crowder bead, 180°, and 0, respectively.
2. Results from the NMR measurements
Figure S1.
High-field chemical shifts for one-dimensional H NMR spectra for native apoazurin in absence (colored in black) and in presence of c d20 = 200 g/l (colored in red). Both spectra were recorded at T = 298 K and B = 20 T. i n t en s i t y H / ppm 0 % 20 % Figure S2.
Two-dimensional H- N NMR HSQC spectra for native apoazurin at different concentrations c of dextran 20 ( c d20 = 0 g/l colored in black, c d20 = 100 g/l colored in red and c d20 = 200 g/l colored in blue). All spectra were recorded for at T = 298 K and B = 20 T. The assignment of signals is found in . Figure S3.
Two-dimensional H- N NMR HSQC spectra for denatured apoazurin ( c urea = 6.8 M) at different concentrations c of dextran 20 ( c d20 = 0 g/l colored in black, c d20 = 100 g/l colored in red and c d20 = 180 g/l colored in blue). All spectra were recorded at T = 298 K and B = 20 T. Figure S4.
High-field chemical shifts for one-dimensional H NMR spectra for denatured apoazurin ( c urea = 6.8 M) in absence (colored in black) and in presence of c d20 = 180 g/l (colored in red). Both spectra were recorded at T = 298 K and B = 20 T.
3. The asphericity parameter ∆ The asphericity parameter Δ is defined as ∆= kL l∑ Cm H +mnD opHqr s(%&𝐓) o , (S7) where 𝜆 . is an eigenvalue of the inertia tensor T and 𝜆̅ is the average eigenvalue given by 𝜆̅ = %&𝐓k , (S8) tr is the trace of T . The elements of T are obtained from the following expression, 𝑇 Mx = KL0 o ∑ C𝑟 .M − 𝑟 /M D C𝑟 .x − 𝑟 /x D, (S9) where α and β represent the axes x , y , and z . N is the number of beads of the coarse-grained protein model, and 𝑟 .M is the component of the position of the bead i in the α -direction.
4. Analysis of r H and R g Hydrodynamic radius r H is an effective radius, which assumes the diffusing protein is spherical. However, as the asphericity ∆ of the protein increases, r H and radius of gyration R g diverge. For pedagogical purposes, we use an extreme case to see how both ∆ and R g are needed to describe r H . This example assumes that as the unfolded protein becomes less spherical, it becomes more like a cylindrical rod with increasing aspect ratio AR . We can write an expression 𝑟 z = 𝑅 < 𝑓(Δ) by using the Stokes-Einstein relation 𝐷 $"=,&, = g (cid:127) (cid:128)R(cid:129)(cid:130)& (cid:131) , (S10) where D sphere is the diffusion coefficient of a sphere, 𝑘 (cid:133) is the Boltzmann constant, 𝑇 is temperature, and 𝜂 is the viscosity, and also using the modified relation for a rod , 𝐷 &9; = g (cid:127) (cid:128)R(cid:129)(cid:130)&(cid:135)5 (cid:136)ln(𝐴𝑅) + 𝑎 + + ((cid:135)5 o (cid:141), (S11) i n t en s i t y H / ppm 0 % 18 % where D rod is the diffusion coefficient of a rod, and a, b, c are constants ( a = 0.312, b = 0.565, c = 0.1). Equating D sphere to D rod , we can find the (effective) r H for the rod, which is 𝑟 = &(cid:135)5(cid:136)ln((cid:135)5)(cid:142))(cid:142) (cid:143)(cid:144)(cid:145) (cid:142) (cid:146)(cid:144)(cid:145)o (cid:141) . (S12) By expressing AR as a function of ∆ 𝐴𝑅(Δ) = (cid:147)
L(cid:148)(cid:142)k√(cid:148)(cid:142)KK+(cid:148) , (S13) and letting 𝑅 = 𝑅 < , we can show how 𝑟 z /𝑅 < diverges as ∆ of the protein increases (see Figure S5). Since AR is the ratio of z and x principle components of gyration tensor and since ∆ is a function of the principle components, we can obtain 𝐴𝑅(Δ) by assuming x and y principle components are equal (as in a cylinder). Figure S5.
Hydrodynamic radius r H over radius of gyration R g of a protein with respect to asphericity ∆ . Figure S6. (a) Average asphericity ( Δ ) of protein apoazurin in varying shapes of crowders normalized by the average asphericity ( Δ B ) of protein apoazurin at the bulk condition against the volume fractions of crowders ( ϕ c ). (b) Average radius of gyration ( R g ) of apoazurin normalized by the average radius of gyration under the bulk condition against ϕ c . The interactions between the protein apoazurin and crowders are attractive. The concentration of urea is 6 M. The temperature is k B T u / ε =1.14, corresponding to a room temperature (see Method), for all conditions. Error bars are included. Figure S7. (a) The radius of gyration R g and (b) asphericity Δ of apoazurin as a function of temperature for apoazurin at several concentrations of urea (0M, 2M, 4M, 6M). Error bars are included. References (1) Samiotakis, A.; Cheung, M. S., Folding dynamics of Trp-cage in the presence of chemical interference and macromolecular crowding. I.
J. Chem. Phys. , 175101. (2) Sobolev, V.; Sorokine, A.; Prilusky, J.; Abola, E. E.; Edelman, M., Automated analysis of interatomic contacts in proteins.
Bioinformatics , 327-332. (3) Takada, S.; Luthey-Schulten, Z.; Wolynes, P. G., Folding dynamics with nonadditive forces: A simulation study of a designed helical protein and a random heteropolymer. The Journal of Chemical Physics , 11616-11629. (4) Nar, H.; Messerschmidt, A.; Huber, R.; Vandekamp, M.; Canters, G. W., Crystal-Structure of Pseudomonas-Aeruginosa Apo-Azurin at 1.85 Angstrom Resolution.
FEBS Letters , 119-124. (5) Cheung, M. S.; Finke, J. M.; Callahan, B.; Onuchic, J. N., Exploring the interplay of topology and secondary structural formation in the protein folding problem.
The Journal of Physical Chemistry B , 11193-11200. (6) Samiotakis, A.; Homouz, D.; Cheung, M. S., Multiscale investigation of chemical interference in proteins.
The Journal of Chemical Physics , 175101. (7) Betancourt, M. R.; Thirumalai, D., Pair potentials for protein folding: choice of reference states and sensitivity of predicted native states to variations in the interaction schemes.
Protein Sci. , 361-369. (8) Zaballa, M. E.; Abriata, L. A.; Donaire, A.; Vila, A. J., Flexibility of the metal-binding region in apo-cupredoxins. Proc. Natl. Acad. Sci. U.S.A. , 9254-9259. (9) Dima, R. I.; Thirumalai, D., Asymmetry in the shapes of folded and denatured states of proteins.
J. Phys. Chem. B , 6564-6570. (10) Tirado, M. M.; Martinez, C. L.; Delatorre, J. G., Comparison of Theories for the Translational and Rotational Diffusion-Coefficients of Rod-Like Macromolecules - Application to Short DNA Fragments.
The Journal of Chemical Physics , 2047-2052., 2047-2052.