for many years , the dominant method for discovering massive clusters of galaxies at @xmath11 , has been via x - ray emission from the hot gas in their dark matter potential wells . over the years
, the _ rosat _ and _ xmm - newton _ telescopes , in particular , have yielded a handful of well - studied examples e.g. , rdcs j0910 + 5422 @xcite , rdcs j1252.9 - 2927 @xcite , rx j0848 + 4452 ( lynx e ; @xcite ) , xmmu j2235 - 25 @xcite , and xmmxcs j2215.9 - 1738 @xcite .
in recent years , it has been realised that , by incorporating deep infrared ( ir ) observations , existing optical imaging techniques can be adapted to successfully detect clusters at redshifts competitive to the existing x - ray surveys .
however , because ir observations covering only modest areas ( 120 arcmin@xmath1 - 7.25 deg@xmath1 ) have been available , the ( cluster or candidate ) systems detected to date have been less massive than those discovered from the x - ray surveys ( @xcite ; although see @xcite ) .
massive clusters of galaxies are rare , and one requires as widefield a survey as possible to detect them . the largest area _
space telescope survey is the 50 square degree seven passband ( 3.6 , 4.5 , 5.8 , 8.0 , 24 , 70 , 160 @xmath12 ) spitzer wide - area infrared extragalactic survey ( swire ) legacy survey @xcite .
we have obtained deep @xmath13band imaging @xcite , and combined this with the pre - existing infrared array camera ( irac ; @xcite ) observations from the swire survey , aiming to select clusters at @xmath14 using a two filter @xmath15 $ ] infrared adaptation of the well - proven optical cluster red sequence ( rs ) method @xcite .
the sparcsgillianw / sparcs ] survey @xcite , is now complete and has an effective area ( defined as the usable overlap with swire after excluding chip gaps , regions near bright stars , etc ) , of 41.9 deg@xmath16 .
the sparcs catalog contains several hundred cluster candidates at @xmath14 , by far the largest , homogeneously selected sample of its kind . in @xcite and @xcite
we presented an overview of the sparcs survey , and reported spectroscopic confirmation of three clusters at @xmath17 and @xmath18 . in this paper , we report spectroscopic confirmation of three additional clusters at @xmath19 and @xmath20 .    the paper is organized as follows : in  [ observations ] , we describe the imaging and spectroscopic observations , and the data reduction ; in  [ analysis ] , we describe the spectroscopic catalog ; in  [ results ] , we present our results for the three individual clusters , and estimate the cluster velocity dispersions and dynamical masses ; and in  [ discussion ] , we present a discussion and the main conclusions based on our results .
we assume a @xmath21-cdm cosmology with @xmath22 and @xmath23 , and h@xmath24 km s@xmath5 mpc@xmath5 .
the sparcs / swire survey is comprised of six fields ( see table 1 of * ? ? ? * ) .
the three clusters described in this paper were identified as high significance cluster candidates in the 7.9 deg@xmath1 elais - n1 field .
a full description of the sparcs data reduction , cluster candidate detection algorithm and catalogs will appear in muzzin et al .
2010 , in preparation .
we provide only a brief overview of the main details here .
sparcs @xmath0 observations of the northern fields were observed by cfht / megacam for a total integration time of 6000s per pointing .
photometry was performed on both the @xmath25 and irac mosaics using the sextractor photometry package @xcite .
the total 3.6 @xmath12 magnitude of all sources in the field was computed using a large aperture , equal to the geometric mean radius of the sextractor isophotal aperture .
further details of the photometric pipeline may be found in @xcite and @xcite .
the depth of the @xmath0 data varies from pointing to pointing depending on the seeing and the sky background , however , the mean 5 @xmath26 depth for extended sources in the elais - n1 field was @xmath27 vega ( 24.2 ab ) .
the @xmath28 color of all sources was computed using an aperture of diameter three irac pixels ( @xmath29 ) .
no aperture corrections were applied to the @xmath25 photometry because the @xmath30 color aperture was much larger than the median seeing . however , the irac point spread function ( psf ) has broad wings compared to a typical ground - based seeing profile , and it was necessary to apply aperture corrections to the measured 3.6 @xmath12 magnitude before computing the color of each galaxy .
galaxy clusters were identified in the sparcs survey , using an algorithm very similar to that described in detail in @xcite , as applied to the spitzer first look survey ( fls ; @xcite ) . for cluster detection in both the fls and sparcs surveys , we used an infrared adaptation of the cluster red - sequence ( crs ) technique @xcite .
this algorithm maps the density of galaxies in a survey within narrow color slices , giving greater weight to brighter galaxies , and flagging the highest overdensities as candidate clusters .
the one important difference here , compared to @xcite , is that the latter paper used an @xmath31 $ ] color to detect clusters at @xmath32 in the fls .
the deeper irac data in swire ( @xmath33s frames ) , compared to the fls ( @xmath34s frames ) , combined with the @xmath15 $ ] color choice , allow sparcs to detect clusters to higher redshift than was possible with the fls dataset .    from analysis of the richnesses and colors , sparcs j161641 + 554513  and sparcs j161037 + 552417
were identified as rich cluster candidates at @xmath11 , and sparcs j161315 + 564930  as an unusually rich cluster candidate at @xmath35 ( table  [ tab_clusters ] ) .
spectroscopic observations of sparcs j161315 + 564930 , sparcs j161641 + 554513  and sparcs j161037 + 552417  were obtained with the low resolution imaging spectrograph ( lris ; @xcite ) on the keck i telescope on the nights of april 3th and april 4th , 2008 .
the seeing on both nights was about @xmath36 .
lris has beam - splitters that separate the light between two sides : red @xcite and blue @xcite . on the red side
, we used the 400 line mm@xmath5 grating which gives a coverage of @xmath37   centered at 8500  , with a dispersion of 1.86  / pixels .
all the masks were designed with slits of 11 width , which gave a resolution of @xmath387 - 9   ( fwhm ) over the wavelength interval @xmath386600 - 9600  . for the redshift range of these particular clusters ,
all of the prominent spectral features e.g. , the 4000  - rest - frame break or [ o@xmath39 emission line , fall in the lris red side wavelength window .
we observed one mask each in the case of clusters sparcs j161315 + 564930  and sparcs j161641 + 554513 , and two in the case of sparcs j161037 + 552417 .
the total exposure time was 6000s ( @xmath40 ) for sparcs j161315 + 564930 , 8400s ( @xmath41 ) for sparcs j161641 + 554513 , and 8400s ( @xmath42 ) and 7000s ( @xmath43 ) for mask @xmath44 and mask @xmath45 of sparcs j161037 + 552417  ( table  1 ) .
the masks contained about 30 slits each ( five of these were alignment stars and the remaining 25 were galaxies , selected according to a sliding scale of priority ) .    to reduce the number of slits placed on obvious foreground galaxies
, we prioritized slits on galaxies with colors near each cluster s red - sequence .
to avoid significant selection bias , we used a very broad cut around the red - sequence , intended to include both star forming and non - star forming systems .
slits were placed on galaxies with priorities from 1 to 4 .
priority 1 was galaxies with colors within 0.2 magnitude of the rs , and with @xmath46 < 17.0 $ ] ( vega ) .
priority 2 was galaxies with colors between 0.2 and 0.6 magnitudes bluer than the rs , and with @xmath47 < 18.7 $ ] .
priority 3 was galaxies with same colors as priority 1 , but with @xmath47 < 18.7 $ ] , and priority 4 was galaxies with same colors as priority 2 , but with @xmath47 < 18.7 $ ] .
priorities 1 to 4 roughly correspond to bright red - sequence , blue cloud , faint red - sequence , and faint blue cloud galaxies respectively .
approximately ten priority 1 galaxies could be accomodated per mask for sparcs j161315 + 564930 , and five per mask for sparcs j161641 + 554513  and sparcs j161037 + 552417 .    to reduce the lris data , we adapted custom software , based on that developed by @xcite to reduce vlt / fors2 data .
overscan and bias corrections were applied to both the calibration ( bias , flats and lamp arcs ) and to the science ( clusters and standard stars ) frames using the iraf _ lrisbias _ and _ lccdproc _ tasks developed by g. d. wirth and c. fassnacht .
after these corrections had been applied , the data were next corrected for geometric distortions across the dispersion axis , and then separated into individual slitlets and reduced using standard long - slit techniques .
each individual slitlet was processed using an algorithm similar to that implemented in _ bogus _ ( developed by d. stern , a. j. bunker and s. a. stanford)stern / homepage / bogus.html ] .
the flatfield slitlets were normalized and applied to the corresponding science slitlets .
the individual 2d spectra were then background - subtracted and fringe - corrected before being stacked to produce a final co - added 2d spectrum .
slits longer than 90 pixels ( @xmath48 ) were further corrected for  long - slit " distortions .
the 1d spectrum of the source ( or sources ) in each slitlet was then extracted using standard iraf tasks .
a set of near lamp exposures were used to wavelength - calibrate the observations .
lris spectra can suffer significant flexure which may introduce wavelength offsets as large as @xmath49   in the wavelength solution , in the most extreme cases .
the offsets in the wavelength calibration due to flexure were corrected by comparing to the wavelengths of known skylines .
we estimate the residual uncertainties in wavelength calibration to be at most 1  , which correspond to a redshift uncertainty of @xmath50 at @xmath51  .
finally , a relative flux calibration was performed , using a sensitivity function spanning the range 5600 - 9400  which was obtained from long - slit observations of the standard stars feige 34 and feige 67 @xcite .
a redshift solution was sought for each object by cross - correlating @xcite each of the lris spectra with observed galaxy templates @xcite using the task _ xcsao _ in the iraf _ rvsao _ package @xcite .
if a redshift could be obtained for a spectrum , it was assigned a quality flag , `` q '' , with one of two possible values : 0 or 1 .
high - confidence redshifts , unambiguously determined by two or more clear features in the continuum or in absorption , or by obvious emission lines such as [ o@xmath39 ( @xmath52 ) were assigned @xmath53 .
lower - confidence redshifts based on ambiguous identification of a few spectral features , often in low s / n spectra , were assigned @xmath54 .
in addition to a quality flag , an emission line flag , `` e '' , was also assigned to each galaxy , based on the presence or absence of emission line features .
spectra displaying one or more emission feature such as h@xmath55 , [ o@xmath56(@xmath574959,@xmath575007 ) or h@xmath58 in the case of low - redshift ( @xmath59 ) galaxies , or [ o@xmath39(@xmath573727 ) in the case of higher redshift galaxies , were assigned @xmath60 .
spectra showing no sign of excess emission above the stellar continuum , were assigned @xmath61 . the total number of galaxies ( cluster members and foreground or background galaxies , excluding stars ) extracted from the fields of sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161641 + 554513  were 22 , 19 and 34 ( see tables [ tab_en1_240 ] , [ tab_en1_349 ] and [ tab_en1_381 ] ) .
occasionally , one or more serendipitous sources fell within a slit .
the breakdown by quality flag was 20 , 13 and 14 @xmath53 galaxies , and 2 , 6 and 20 @xmath54 galaxies , respectively .
example spectra for a subsample of cluster members ( see section  [ cl240 ] for the definition of a `` cluster member '' ) are shown in figure  [ cl_specs ] .
the left column shows examples of sparcs j161315 + 564930  members , the center column shows examples of sparcs j161641 + 554513  members , and the right column shows exampels of sparcs j161037 + 552417  members .
the i d of each object is indicated ( tables [ tab_en1_240 ] , [ tab_en1_349 ] and [ tab_en1_381 ] ) .
fluxes are in relative units , and smoothed by 7 pixels ( 1 pix @xmath62  ) .
prominent spectral features are indicated by vertical lines .
[ results ]      in determining cluster membership , and calculating a velocity dispersion and mass , only high confidence , @xmath53 , galaxies were considered .
definitive cluster membership was determined using the code of @xcite , which is based on the shifting - gap technique of @xcite .
this procedure uses both galaxy angular position and radial velocity information to exclude near - field interlopers .
the squares in the left panel of figure  [ z_histo_240 ] show the 14 @xmath53 galaxies identified as members of cluster sparcs j161315 + 564930  ( see table  [ tab_en1_240 ] ) by the shifting - gap technique . for sparcs j161315 + 564930
, these @xmath53 cluster members fall in the range @xmath63 , indicated by the vertical dashed lines shown in the right panel of figure  [ z_histo_240 ] , which shows the redshift histogram for the sparcs j161315 + 564930  field .
the properties of the cluster members and non - members are summarized in table  [ tab_en1_240 ] . the total 3.6 @xmath12 magnitude ( column 4 in table  [ tab_en1_240 ] ) , and @xmath28 color ( column 5 in table  [ tab_en1_240 ] ) were calculated as described in section  [ photometry ] . in determining membership
, we did not require that any @xmath54 galaxies satisfy the shifting - gap criteria ; if their redshifts fell in the range @xmath63 we included them as  cluster members " in table  [ tab_en1_240 ] , but re - emphasize that we do not utilize them in estimating the cluster mean redshift or velocity dispersion .
the total number of cluster members is 16 ( 14 with @xmath53 and 2 with @xmath54 ) , and six foreground / background galaxies .
all confirmed members of cluster sparcs j161315 + 564930  were passive ( @xmath61 ) galaxies .    in some cases ,
a redshift was obtained which did not correspond to any galaxy in our photometric catalog . in the case of faint galaxies , this was because a spectroscopic redshift was obtained for a strong emission line galaxy whose continuum fell below the detection threshold of the @xmath25 catalog . in the case of bright galaxies ,
this was because of blending issues with the irac @xmath64 @xmath65 psf .
these galaxies were assigned both a @xmath46$]-band magnitude and a @xmath66 $ ] color of 99 in table  [ tab_en1_240 ] .
figure  [ proj_dist_240 ] shows @xmath67 $ ] color composites of cluster sparcs j161315 + 564930 .
the @xmath68 data were obtained from wfc on the isaac newton telescope , and are available with the swire public data release @xcite .
the white squares ( green circles ) overlaid on the right panel show the 16 cluster members ( and foreground / background galaxies in the @xmath69 fov ) with spectroscopically - confirmed redshifts from keck / lris ( see table  [ tab_en1_240 ] ) .
once cluster membership was established , both the redshift and rest - frame velocity dispersion , @xmath70 , of sparcs j003550 - 431224 were then calculated iteratively using the `` robust estimator '' , @xmath71 @xcite .
the robust estimator has been shown to be less sensitive than the standard deviation to outliers which may persist even after rejecting interlopers using the shifting - gap technique .
the actual estimator used depends on the number of cluster members and is either the biweight estimator for datasets with at least 15 members , or , as here in the case of sparcs j161315 + 564930  with 14 @xmath53 members , the gapper estimator .
the gapper estimator is discussed more fully in @xcite , @xcite and @xcite .
the line - of - sight rest - frame velocity dispersion , @xmath72 , was calculated directly from the vector of spectroscopic redshift measurements , , as    @xmath73    where @xmath74 is the speed of light , and @xmath75 is the estimated dispersion of the measured redshifts with respect to the center of the distribution , @xmath76 .    a mean redshift of @xmath77 and a velocity dispersion of @xmath78 km s@xmath5 were calculated for sparcs j161315 + 564930  ( table  [ tab_clusters ] ) .
the uncertainty on the latter was determined using jackknife resampling of the data . for comparison , a gaussian with an rms of 1230
km s@xmath5 has been overlaid on the redshift histogram in the right panel of figure  [ z_histo_240 ] .
the line - of - sight rest - frame velocity dispersion can be used to calculate a dynamical estimate of @xmath79 , the radius at which the mean interior density is 200 times the critical density , @xmath80 , and @xmath81 , the mass contained within @xmath79 . in the spherical collapse model , at redshift @xmath82 , @xmath83 can be calculated from    @xmath84    where @xmath85 , is the hubble parameter at redshift @xmath82 , and    @xmath86    based on its velocity dispersion of @xmath78 km s@xmath5 , we estimate an @xmath87 mpc and a dynamical mass of @xmath88 m@xmath89 for sparcs j161315 + 564930  ( table  [ tab_clusters ] ) .
although this mass is preliminary ( see section  [ richness ] ) , it seems likely that sparcs j161315 + 564930  is an unusually massive cluster , perhaps the most massive cluster in the entire sparcs survey , and comparable in mass to that of cluster ms 1054 - 03 at @xmath90 @xcite .
the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath83 ( @xmath91 mpc ) of the center of sparcs j161315 + 564930  is shown in figure  [ colmag_240 ] .
this radius is approximately equal to the virial radius .
spectroscopically - confirmed galaxies are shown by colored symbols . in this figure ( and in figures  [ colmag_349 ] and [ colmag_381 ] ) , cluster members ( @xmath53 or @xmath54 ) are shown by red circles and foreground / background galaxies by blue squares .
note that there are several cluster members or foreground / background galaxies shown in table  [ tab_en1_240 ] , which fall within a projected radius of @xmath83 of the cluster center , but for which a color could not be determined .
these galaxies do not appear in figure  [ colmag_240 ] ( or in figures  [ colmag_349 ] or [ colmag_381 ] , in the cases of tables  [ tab_en1_349 ] and tables  [ tab_en1_381 ] ) .
seven @xmath53 galaxies were determined to be cluster members of sparcs j161641 + 554513  by the shifting - gap technique .
these galaxies are shown by squares in the left panel of figure  [ z_histo_349 ] .
these galaxies lie in the redshift range @xmath92 , indicated by the vertical dashed lines in the right panel of figure  [ z_histo_349 ] .
table  [ tab_en1_349 ] summarizes the  cluster members " , the ten @xmath53 and @xmath54 galaxies with redshifts in this range , and the nine foreground / background galaxies .
the right panel of figure  [ z_histo_349 ] shows the redshift histogram for the field of sparcs j161641 + 554513 . a gaussian with an rms of 950 km s@xmath5 has been overlaid .
figure  [ proj_dist_349 ] shows @xmath67 $ ] color composites of cluster sparcs j161641 + 554513 .
the white squares ( green circles ) overlaid on the right panel show the cluster members ( foreground / background galaxies ) with spectroscopically - confirmed redshifts ( see table  [ tab_en1_349 ] ) which fall within the @xmath93 fov of the image .    of the seven @xmath53 cluster members , five were classified as passive ( @xmath61 ) , one was classified as emission line ( @xmath60 ) , and one was classifed as an agn ( @xmath94 ) .
the upper left white square in figure  [ proj_dist_349 ] corresponds to the agn .
there were also three passive @xmath54 cluster members ( table  [ tab_en1_349 ] ) .
the spectrum of the confirmed agn ( object i d @xmath95 854892 in table  [ tab_en1_349 ] ) is shown in the lowest panel of the center column in figure  [ cl_specs ] .
this source shows nev(@xmath573346,@xmath573426 ) and neiii(@xmath573869,@xmath573968 ) in emission , which are common in agn , as well as prominent [ oii](@xmath573727 ) emission and high - order balmer lines ( h@xmath96 , h@xmath97 and h6 ) also in emission .
the cluster mean redshift and velocity dispersion , were estimated iteratively from the seven @xmath53 members , using the gapper estimator .
the mean redshift of sparcs j161641 + 554513  was calculated to be @xmath98 .
the velocity dispersion was calculated to be @xmath99 km s@xmath5 , which corresponds to @xmath100 mpc , and a dynamical mass of @xmath101 m@xmath89(table  [ tab_clusters ] ) .
the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath83 ( @xmath102 mpc ) of the center of sparcs j161641 + 554513  is shown in figure  [ colmag_349 ] .
the red circles and blue squares indicate those cluster members and foreground / background galaxies which lie within a projected distance of @xmath83 of the cluster center , and for which a color could be determined .
seven @xmath53 galaxies were also determined to be cluster members of sparcs j161037 + 552417  by the shifting - gap technique .
these galaxies are indicated by squares in the left panel of figure  [ z_histo_381 ] .
two galaxies , indicated by crosses in figure  [ z_histo_381 ] ( i d @xmath95
s 727869 and 734082 in table  [ tab_en1_381 ] ) , were identified as near - field interlopers .
the seven @xmath53 members lie in the redshift range @xmath103 , shown by the vertical dashed lines in the right panel of figure  [ z_histo_381 ] .
table [ tab_en1_381 ] summarizes the  cluster members " , the ten @xmath53 or @xmath54 galaxies with redshifts in this range , and the 24 foreground / background galaxies .
the right panel of figure  [ z_histo_381 ] shows the redshift histogram for the field of sparcs j161037 + 552417 . a gaussian with an rms of 410 km s@xmath5 has been overlaid .
figure  [ proj_dist_381 ] shows @xmath67 $ ] color composites of cluster sparcs j161037 + 552417 .
the white squares ( green circles ) overlaid on the right panel show the cluster members ( foreground / background galaxies ) with spectroscopically - confirmed redshifts ( see table  [ tab_en1_381 ] ) which fall within the @xmath104 fov of the image . of the seven @xmath53 cluster members , all were classified as emission line ( @xmath60 ) .
of the three @xmath54 cluster members , one was classified as passive ( @xmath61 ) , and two as emission line ( @xmath60 ) .
the cluster mean redshift and velocity dispersion , were calculated iteratively from the seven @xmath53 members , using the gapper estimator .
the mean redshift of sparcs j161037 + 552417  was estimated to be @xmath105 and the velocity dispersion , @xmath106 km s@xmath5 , which corresponds to @xmath107 mpc and a dynamical mass of @xmath108 m@xmath89 ( table  [ tab_clusters ] ) .    sparcs j161037 + 552417  is significantly less massive than sparcs j161315 + 564930  or sparcs j161641 + 554513 .
the much smaller @xmath83 estimated for sparcs j161037 + 552417 , combined with the geometric limitations of lris with respect to the redshift number density yield measurable from a single mask , resulted in a yield of only four spectroscopically confirmed members within @xmath83 from two masks .
moreover , because of blending issues with irac s @xmath64 @xmath65 psf , a reliable @xmath66 $ ] color could not be determined for one of these galaxies ( i d @xmath95 4000013 in table  [ tab_en1_381 ] ) .
figure  [ colmag_381 ] shows the @xmath66 $ ] vs. @xmath0 color - magnitude diagram for all galaxies ( gray circles ) within a radius of @xmath109 ( @xmath110 kpc ) of the center of cluster sparcs j161037 + 552417 .
the three red circles denote the cluster members with i d @xmath95 s 723814 , 722784 , and 722712 in table  [ tab_en1_381 ] .
color - magnitude diagrams for sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161037 + 552417  were presented in figures  [ colmag_240 ] , [ colmag_349 ] , and [ colmag_381 ] .
column 2 of table  [ tab_properties ] shows the @xmath15 $ ] color of the red sequence for sparcs j161315 + 564930 , sparcs j161037 + 552417  and sparcs j161037 + 552417 .
the color is calculated from the mean color of the spectroscopically confirmed red - sequence cluster members . also shown in column 2 of table  [ tab_properties ] is the color of the rs for three additional clusters , previously reported in @xcite and @xcite .
these clusters are sparcs j163435 + 402151  at @xmath111 , sparcs j163852 + 403843  at @xmath112 , and sparcs j003550 - 431224  at @xmath113 ( column 6 ) . for all six clusters
a systematic uncertainty in the rs color of 0.15 magnitude has been assumed ( this uncertainty reflects the fact that , at present , we are using the zeropoints provided by elixir for the @xmath0  observations . we expect , in the future , to be able to reduce these photometric uncertainties , using our own internal calibration ) .    columns 3 , 4 and 5 of table  [ tab_properties ] show the redshift that would be estimated for each cluster based on the measured rs color ( column 2 ) , assuming a solar metallicity single burst bc03 model and a formation redshift of either @xmath114 , 4 or 10 . as can be seen from table  [ tab_properties ] ( and the left panel of figure  [ zsandmasses ] ) , the photometric redshift inferred from the measured @xmath66 $ ] color has a slight dependence on one s choice of formation redshift , although differences in color between the models at @xmath115 are fairly small ( @xmath116 between @xmath117 and @xmath118 ) . utilizing the colors from the @xmath119 model ,
the three new clusters presented here were assigned preliminarily redshift estimates of @xmath120 , @xmath121 , and @xmath122 ( table  [ tab_properties ] ) .
these photometrically estimated redshifts are very similar to , albeit slightly lower than , the spectroscopically determined values of @xmath123 , @xmath124 , and @xmath125 ( column 6 in table  [ tab_properties ] ) .
the @xmath126 $ ] color vs. spectroscopic redshift for all six sparcs clusters in table  [ tab_properties ] is plotted in the left panel of figure  [ zsandmasses ] .
the solid , dotted and and dashed lines show the bc03 model colors as a function of redshift for formation redshifts of @xmath114 , 4 and 10 .
it is clear from figure  [ zsandmasses ] that the agreement between the model colors and the observations is very good .    with a larger sample of clusters
, there may turn out to be small but real discrepancies between the models and the measured rs colors . in the left panel of figure  [ zsandmasses ] ,
the @xmath119 model can be seen to be slightly redder than the observations in the case of five clusters ( or equivalantly , the inferred photometric redshift can be seen to be slightly lower than the spectroscopic redshift ) , but slightly bluer in the case of one cluster ( sparcs j003550 - 431224 ) .
these small offsets in color between the models and the observations , can also be seen directly from figures  [ colmag_240 ] , [ colmag_349 ] , and [ colmag_381 ] .
the solid lines in these three figures show the rs color predicted by the bc03 @xmath119 model _ at the spectroscopic redshift of the cluster _ , and can be seen to be slightly redder than the observed color .
a more detailed comparison between the model predictions and the observations will be made in a future paper employing a larger sample of sparcs clusters .    despite the aforementioned caveats , and the issue of degeneracies between the photometric redshift of the clusters and the formation redshifts of their galaxies ,
our overall conclusion is that the inferred one - color photometric redshifts and the spectroscopic redshifts are in excellent ( @xmath127 ) general agreement .
in addition to estimating the masses for sparcs j161315 + 564930 , sparcs j161641 + 554513 , and sparcs j161037 + 552417  from the galaxy line - of - sight velocity dispersion , we also estimated the masses from the richness of the clusters , using the @xmath128 richness parameter .
@xcite introduced @xmath128 , an adaptation of the @xmath129 richness parameter , intended to utilize two - band photometry to increase the contrast of the cluster with the background , and therefore provide a measurement of the richness that is less sensitive to foreground / background large scale structures .
@xmath129 is the amplitude of the three - dimensional , cluster center - galaxy spatial correlation function , @xmath130 @xcite .    instead of counting galaxies in a single passband
, @xmath128 is obtained by counting galaxies in a color slice centered on the location of each cluster s red - sequence in the @xmath66 $ ] color - magnitude diagram . in computing @xmath128
, we used a slice bounded in color by @xmath66 = \pm0.3 $ ] of the best - fit rs color returned by the cluster finding algorithm ( @xcite ) , and bounded in magnitude by @xmath131 , where @xmath132 is the bc03 @xmath119 model prediction of the characteristic magnitude of a galaxy _ at the photometric redshift _ corresponding to that rs color ( table  [ tab_properties ] ) .
the background galaxy counts were determined from the color distribution in the entire 7.9 deg@xmath1 elais - n1 field , minus the regions known to contain galaxy clusters .
the @xmath128 richnesses of the three clusters were computed to be @xmath133 mpc@xmath134 ( sparcs j161315 + 564930 ) , @xmath135 mpc@xmath134 ( sparcs j161641 + 554513 ) , and @xmath136 mpc@xmath134 ( sparcs j161037 + 552417 ) . based on the empirical calibration of @xmath137 vs. @xmath138 determined by @xcite in the k - band for 15 cnoc1 clusters at @xmath139
, these richnesses correspond to @xmath140 for sparcs j161315 + 564930 , @xmath141 for sparcs j161641 + 554513 , and @xmath142 for sparcs j161037 + 552417 .    for comparison ,
columns 7 and 8 of table  [ tab_properties ] show the dynamical mass , @xmath143 , and richness mass , @xmath144 , estimates for all six sparcs clusters .
although the uncertainties associated with both of the mass estimators are large , they are consistent with each other at the 1@xmath145 level for five out of the six clusters and at the 2@xmath145 level for the sixth .
the agreement between the two mass estimators can be seen in the right panel of figure  [ zsandmasses ] .
based on all six sparcs clusters spectroscopically confirmed to date , our conclusion is that , in addition to there being excellent agreement between the photometric and spectroscopic redshifts , there is also reasonable agreement between the cluster dynamical and richness mass estimates .
the dynamical masses estimates should be considered preliminary at this stage .
uncertainties will reduce as more data becomes available from a large spectroscopic follow - up program of sparcs clusters currently being carried out at the gemini telescopes .
collectively , the six sparcs clusters confirmed to date ( the three clusters presented in @xcite and @xcite , plus the three clusters presented here ) , demonstrate that , given the availability of infrared observations , the rs technique is an efficient and effective method of detecting _ bona fide _ massive galaxy clusters at @xmath146 .
moreover , our studies of these six clusters are showing that it is possible to infer fundamental parameters such as cluster redshift and mass _ from the survey data itself _
( see also @xcite ) .    at @xmath147 , both
the optical red - sequence cluster surveys , rcs-1 @xcite and rcs-2 @xcite , and the sloan digital sky survey ( sdss ; @xcite ) have shown that it is feasible to measure cosmological parameters from the evolution of the cluster mass function @xcite . in order to do this efficiently ,
the survey data themselves are used to detect clusters , and also to estimate the redshift and the mass of those clusters .
the redshift is estimated from the red sequence color @xcite , and the mass is estimated from the optical richness @xcite .
the fact that sparcs is also now demonstrating the practicality of estimating redshifts and masses at @xmath146 from the survey data alone is heartening for the current generation of surveys aiming to utilize optical - infrared high redshift cluster observations to constrain cosmological parameters e.g. , sparcs , the ukirt infrared deep - sky survey deep extragalactic survey ( ukidss dxs ; @xcite),and the irac shallow cluster survey ( iscs ; @xcite )
. these optical - ir surveys will provide complementary samples to those selected using the sunyaev - zeldovich effect , e.g. , the south pole telescope survey ( spt ; @xcite ) , the atacama cosmology telescope ( act ; @xcite ) , and the atacama pathfinder experiment ( apex ; @xcite ) .
the complete sparcs catalog contains several hundred cluster candidates at @xmath14 .
with new large , homogeneous , reliable @xmath11 catalogs becoming available from sparcs and other surveys in the very near future , the prospects look bright for high redshift cluster and cluster galaxy evolution studies in the coming years .
+    we thank the referee for useful comments .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .
this work is based in part on archival data obtained with the spitzer space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa .
support for this work was provided by an award issued by jpl / caltech .
gw also gratefully acknowledges support from nsf grant ast-0909198 , and from the college of natural and agricultural sciences at ucr .
andreon , s. , maughan , b. , trinchieri , g. , & kurk , j.  2008 , , 507 , 147 becker , m.  r. , et al .
2007 , , 669 , 905 beers , t.  c. , flynn , k. , & gebhardt , k.  1990 , , 100 , 32 bertin , e. , & arnouts , s.  1996 , , 117 , 393 blindert , k.  2006 , ph.d .
thesis , university of toronto , canada bremer , m.  n. , et al .
2006 , , 371 , 1427 brodwin , m. , et al .
2006 , , 651 , 791 bruzual , g. , & charlot , s.  2003 , , 344 , 1000 carlstrom , j.  e. , et al .
2009 , arxiv:0907.4445 demarco , r. , et al .
2005 , , 432 , 381 demarco , r. , et al .
2007 , , 663 , 164 dobbs , m. , et al .
2006 , new astronomy review , 50 , 960 eisenhardt , p.  r.  m. , et al .
2008 , , 684 , 905 fadda , d. , girardi , m. , giuricin , g. , mardirossian , f. , & mezzetti , m.  1996 , , 473 , 670 fazio , g.  g. , et al .
2004 , , 154 , 10 gilbank , d.  g. , yee , h.  k.  c. , ellingson , e. , gladders , m.  d. , barrientos , l.  f. , & blindert , k.  2007 , , 134 , 282 gioia , i.  m. , braito , v. , branchesi , m. , della ceca , r. , maccacaro , t. , & tran , k .- v .
2004 , , 419 , 517 girardi , m. , biviano , a. , giuricin , g. , mardirossian , f. , & mezzetti , m.  1993 , , 404 , 38 girardi , m. , demarco , r. , rosati , p. , & borgani , s.  2005 , , 442 , 29 gladders , m.  d. , yee , h.  k.  c. , majumdar , s. , barrientos , l.  f. , hoekstra , h. , hall , p.  b. , & infante , l.  2007 , , 655 , 128 gladders , m.  d. , & yee , h.  k.  c.  2000 , , 120 , 2148 gladders , m.  d. , & yee , h.  k.  c.  2005 , , 157 , 1 goto , t. , et al .
2008 , , 60 , 531 hoekstra , h. , franx , m. , & kuijken , k.  2000 , , 532 , 88 jee , m.  j. , white , r.  l. , ford , h.  c. , blakeslee , j.  p. , illingworth , g.  d. , coe , d.  a. , & tran , k .- v .  h.  2005 , , 634 , 813 kinney , a.  l. , calzetti , d. , bohlin , r.  c. , mcquade , k. , storchi - bergmann , t. , & schmitt , h.  r.  1996 , , 467 , 38 koester , b.  p. , et al .
2007 , , 660 , 239 kosowsky , a.  2003 , new astronomy review , 47 , 939 koyama , y. , kodama , t. , tanaka , m. , shimasaku , k. , & okamura , s.  2007 , , 382 , 1719 krick , j.  e. , surace , j.  a. , thompson , d. , ashby , m.  l.  n. , hora , j.  l. , gorjian , v. , & yan , l.  2008 , , 686 , 918 kurk , j. , et al .
2008 , astronomical society of the pacific conference series , 399 , 332 kurtz , m.  j. , mink , d.  j. , wyatt , w.  f. , fabricant , d.  g. , torres , g. , kriss , g.  a. , & tonry , j.  l.  1992 , astronomical data analysis software and systems i , 25 , 432 lacy , m. , et al .
2005 , , 161 , 41 lawrence , a. , et al .
2007 , , 379 , 1599 lonsdale , c.  j. , et al .
2003 , , 115 , 897 mccarthy , j.  k. , et al .
1998 , , 3355 , 81 mccarthy , p.  j. , et al .
2007 , , 664 , l17 mullis , c.  r. , rosati , p. , lamer , g. , bhringer , h. , schwope , a. , schuecker , p. , & fassbender , r.  2005 , , 623 , l85 muzzin , a. , yee , h.  k.  c. , hall , p.  b. , & lin , h.  2007 , , 663 , 150 muzzin , a. , wilson , g. , lacy , m. , yee , h.  k.  c. , & stanford , s.  a.  2008 , , 686 , 966 muzzin , a. , et al .
2009 , , 698 , 1943 oke , j.  b.  1990 , , 99 , 1621 oke , j.  b. , et al .
1995 , , 107 , 375 rosati , p. , stanford , s.  a. , eisenhardt , p.  r. , elston , r. , spinrad , h. , stern , d. , & dey , a.  1999 , , 118 , 76 rosati , p. , et al .  2004 , , 127 , 230 rozo , e. , et al .
2009 , , 703 , 601 rozo , e. , et al .
2009 , arxiv:0902.3702 ruhl , j. , et al .
2004 , , 5498 , 11 shapley , a.  e. , steidel , c.  c. , pettini , m. , & adelberger , k.  l.  2003 , , 588 , 65 shupe , d.  l. , et al .
2008 , , 135 , 1050 stanford , s.  a. , elston , r. , eisenhardt , p.  r. , spinrad , h. , stern , d. , & dey , a.  1997 , , 114 , 2232 stanford , s.  a. , holden , b. , rosati , p. , eisenhardt , p.  r. , stern , d. , squires , g. , & spinrad , h.  2002 , , 123 , 619 stanford , s.  a. , et al .
2005 , , 634 , l129 stanford , s.  a. , et al .
2006 , , 646 , l13 , j. , shupe , d.  l. , fang , f. , lonsdale , c.  j. , & gonzalez - solares , e. 2005 , ssc website release tran , k .- v .
h. , franx , m. , illingworth , g.  d. , van dokkum , p. , kelson , d.  d. , blakeslee , j.  p. , & postman , m.  2007 , , 661 , 750 tonry , j. , & davis , m.  1979 , , 84 , 1511 van breukelen , c. , et al .
2007 , , 382 , 971 wilson , g. , et al .
2009 , , 698 , 1934 yee , h.  k.  c. , & ellingson , e.  2003 , , 585 , 215 yee , h.  k.  c. , & lpez - cruz , o.  1999 , , 117 , 1985 yee , h.  k.  c. , gladders , m.  d. , gilbank , d.  g. , majumdar , s. , hoekstra , h. , ellingson , e. , & the rcs-2 collaboration 2007 , arxiv : astro - ph/0701839 york , d.  g. , et al .
2000 , , 120 , 1579 zatloukal , m. , rser , h .- j .
, wolf , c. , hippelein , h. , & falter , s.  2007 , , 474 , l5      sparcs j161315 + 564930 & 16:13:14.6 & + 56:49:29.9 & 1 & 16 & 14 & @xmath148 & @xmath149 & @xmath150 & @xmath151 + sparcs j161641 + 554513 & 16:16:41.3 & + 55:45:12.5 & 1 & 10 & 7 & @xmath152 & @xmath153 & @xmath154 & @xmath155 + sparcs j161037 + 552417 & 16:10:36.5 & + 55:24:16.6 & 2 & 10 & 7 & @xmath156 & @xmath157 & @xmath158 & @xmath159 + [ tab_clusters ]      1070485 & 243.32851 & 56.821579 & @xmath160 & @xmath161 & @xmath162 & 0 & 0 + 1071085 & 243.31200 & 56.828140 & @xmath163 & @xmath164 & @xmath165 & 0 & 0 + 1071707 & 243.34280 & 56.835220 & @xmath166 & @xmath167 & @xmath168 & 0 & 0 + 1072162 & 243.36000 & 56.840450 & @xmath169 & @xmath170 & @xmath171 & 0 & 0 + 1068630 & 243.33870 & 56.802010 & @xmath172 & @xmath173 & @xmath174 & 0 & 0 + 1072641 & 243.30150 & 56.845970 & @xmath175 & @xmath176 & @xmath177 & 0 & 0 + 2000012 & 243.31107 & 56.826047 & @xmath178 & @xmath178 & @xmath179 & 0 & 0 + 1070775 & 243.31090 & 56.824970 & @xmath180 & @xmath181 & @xmath182 & 0 & 0 + 1066451 & 243.39220 & 56.778770 & @xmath183 & @xmath184 & @xmath185 & 0 & 1 + 2000007 & 243.34483 & 56.772531 & @xmath178 & @xmath178 & @xmath186 & 0 & 0 + 1066059 & 243.35040 & 56.775040 & @xmath187 & @xmath188 & @xmath189 & 0 & 0 + 1069805 & 243.34230 & 56.814430 & @xmath190 & @xmath191 & @xmath192 & 0 & 0 + 1074041 & 243.28720 & 56.862221 & @xmath193 & @xmath194 & @xmath195 & 0 & 1 + 1072364 & 243.29550 & 56.842720 & @xmath196 & @xmath197 & @xmath198 & 0 & 0 + 1069395 & 243.32159 & 56.810188 & @xmath199 & @xmath200 & @xmath201 & 0 & 0 + 1073458 & 243.30260 & 56.855709 & @xmath202 & @xmath173 & @xmath203 & 0 & 0 + 1068289 & 243.33870 & 56.797932 & @xmath204 & @xmath205 & @xmath206 & 1 & 0 + 1071481 & 243.32381 & 56.832401 & @xmath207 & @xmath167 & @xmath208 & 0 & 0 + 1070036 & 243.27380 & 56.817139 & @xmath209 & @xmath210 & @xmath211 & 0 & 0 + 1067689 & 243.27271 & 56.791660 & @xmath212 & @xmath213 & @xmath214 & 1 & 0 + 1065709 & 243.34480 & 56.771900 & @xmath215 & @xmath216 & @xmath217 & 1 & 0 + 2000006 & 243.34517 & 56.771469 & @xmath178 & @xmath178 & @xmath218 & 1 & 0 + [ tab_en1_240 ]      848097 & 244.17630 & 55.748390 & @xmath219 & @xmath220 & @xmath221 & 0 & 1 + 853440 & 244.18610 & 55.764381 & @xmath166 & @xmath222 & @xmath223 & 0 & 0 + 849838 & 244.17340 & 55.753590 & @xmath224 & @xmath225 & @xmath226 & 0 & 0 + 849110 & 244.19260 & 55.751389 & @xmath227 & @xmath228 & @xmath229 & 0 & 1 + 851083 & 244.17180 & 55.757141 & @xmath230 & @xmath231 & @xmath232 & 0 & 0 + 837239 & 244.26920 & 55.715519 & @xmath233 & @xmath234 & @xmath235 & 1 & 0 + 846319 & 244.25529 & 55.743149 & @xmath236 & @xmath237 & @xmath238 & 0 & 1 + 3000009 & 244.26935 & 55.715914 & @xmath178 & @xmath178 & @xmath239 & 0 & 0 + 839226 & 244.27161 & 55.721340 & @xmath240 & @xmath241 & @xmath242 & 0 & 0 + 854892 & 244.21670 & 55.768921 & @xmath243 & @xmath244 & @xmath245 & 2 & 0 + 833647 & 244.16330 & 55.704411 & @xmath246 & @xmath247 & @xmath248 & 1 & 0 + 828383 & 244.22971 & 55.688580 & @xmath249 & @xmath250 & @xmath251 & 1 & 0 + 840274 & 244.26221 & 55.724548 & @xmath252 & @xmath253 & @xmath254 & 1 & 0 + 856147 & 244.18530 & 55.772751 & @xmath255 & @xmath256 & @xmath257 & 1 & 0 + 844301 & 244.21809 & 55.736950 & @xmath258 & @xmath259 & @xmath260 & 1 & 0 + 842789 & 244.15030 & 55.732430 & @xmath261 & @xmath262 & @xmath263 & 0 & 1 + 857305 & 244.18820 & 55.776299 & @xmath264 & @xmath265 & @xmath266 & 0 & 0 + 851794 & 244.17329 & 55.759270 & @xmath267 & @xmath234 & @xmath268 & 0 & 1 + 860365 & 244.18150 & 55.786140 & @xmath209 & @xmath269 & @xmath270 & 0 & 1 + [ tab_en1_349 ]      4000013 & 242.65226 & 55.405353 & @xmath178 & @xmath178 & @xmath271 & 1 & 0 + 723814 & 242.65221 & 55.404598 & @xmath272 & @xmath273 & @xmath274 & 1 & 0 + 722784 & 242.63251 & 55.401890 & @xmath275 & @xmath276 & @xmath277 & 0 & 1 + 718577 & 242.61940 & 55.391060 & @xmath278 & @xmath279 & @xmath280 & 1 & 0 + 722712 & 242.66769 & 55.401711 & @xmath281 & @xmath282 & @xmath283 & 1 & 1 + 4000026 & 242.60676 & 55.475206 & @xmath178 & @xmath178 & @xmath284 & 1 & 1 + 710318 & 242.60490 & 55.370651 & @xmath285 & @xmath286 & @xmath287 & 1 & 0 + 4000010 & 242.67026 & 55.393053 & @xmath178 & @xmath178 & @xmath288 & 1 & 0 + 728525 & 242.69780 & 55.416550 & @xmath289 & @xmath290 & @xmath291 & 1 & 0 + 4000032 & 242.60455 & 55.371836 & @xmath178 & @xmath178 & @xmath292 & 1 & 0 + 744358 & 242.57840 & 55.457119 & @xmath293 & @xmath294 & @xmath295 & 1 & 0 + 4000012 & 242.65180 & 55.404111 & @xmath178 & @xmath178 & @xmath296 & 1 & 0 + 723736 & 242.64661 & 55.404339 & @xmath297 & @xmath298 & @xmath299 & 0 & 1 + 709346 & 242.61549 & 55.368011 & @xmath300 & @xmath301 & @xmath302 & 0 & 1 + 714998 & 242.59080 & 55.381821 & @xmath303 & @xmath304 & @xmath305 & 1 & 1 + 710721 & 242.60280 & 55.371490 & @xmath306 & @xmath307 & @xmath308 & 0 & 0 + 751852 & 242.60730 & 55.476471 & @xmath309 & @xmath310 & @xmath311 & 0 & 1 + 734249 & 242.66650 & 55.430779 & @xmath312 & @xmath313 & @xmath314 & 0 & 1 + 747661 & 242.62151 & 55.465469 & @xmath315 & @xmath316 & @xmath317 & 1 & 1 + 740510 & 242.69310 & 55.447041 & @xmath255 & @xmath318 & @xmath319 & 0 & 1 + 708779 & 242.61410 & 55.366550 & @xmath320 & @xmath321 & @xmath322 & 1 & 1 + 725337 & 242.69600 & 55.408329 & @xmath323 & @xmath324 & @xmath325 & 0 & 1 + 4000015 & 242.60952 & 55.453100 & @xmath178 & @xmath178 & @xmath326 & 0 & 1 + 730800 & 242.66510 & 55.422260 & @xmath327 & @xmath328 & @xmath329 & 0 & 0 + 720317 & 242.64140 & 55.395630 & @xmath330 & @xmath331 & @xmath332 & 0 & 1 + 745152 & 242.66029 & 55.459110 & @xmath333 & @xmath334 & @xmath335 & 0 & 1 + 753137 & 242.63040 & 55.479912 & @xmath336 & @xmath337 & @xmath338 & 0 & 1 + 738292 & 242.63560 & 55.441189 & @xmath339 & @xmath340 & @xmath341 & 0 & 0 + 727869 & 242.68291 & 55.414661 & @xmath333 & @xmath342 & @xmath343 & 1 & 0 + 734082 & 242.71609 & 55.430450 & @xmath344 & @xmath345 & @xmath343 & 1 & 0 + 712763 & 242.63080 & 55.376591 & @xmath346 & @xmath347 & @xmath348 & 1 & 1 + 742510 & 242.67841 & 55.452202 & @xmath349 & @xmath350 & @xmath351 & 1 & 1 + 737166 & 242.70931 & 55.438370 & @xmath352 & @xmath353 & @xmath354 & 1 & 1 + 735789 & 242.69299 & 55.434792 & @xmath355 & @xmath356 & @xmath357 & 1 & 1 + [ tab_en1_381 ]      sparcs j161315 + 564930 & @xmath358 & @xmath359 & @xmath360 & @xmath361 & @xmath148 & @xmath151 & @xmath362 + sparcs j161641 + 554513 & @xmath363 & @xmath364 & @xmath365 & @xmath366 & @xmath152 & @xmath155 & @xmath367 + sparcs j161037 + 552417 & @xmath368 & @xmath369 & @xmath370 & @xmath371 & @xmath156 & @xmath159 & @xmath372 + sparcs j163435 + 402151 & @xmath373 & @xmath374 & @xmath375 & @xmath376 & @xmath377 & @xmath378 & @xmath379 + sparcs j163852 + 403843 & @xmath380 & @xmath381 & @xmath382 & @xmath383 & @xmath384 & @xmath385 & @xmath386 + sparcs j003550 - 431224 & @xmath387 & @xmath388 & @xmath389 & @xmath390 & @xmath391 & @xmath392 & @xmath379 + [ tab_properties ]


cells exist in a crowded environment consisting of macromolecules ( lipids , mrna , ribosome , sugars , etc . ) , making it critical to investigate protein folding in the presence of crowding agents @xcite .
if the interactions between the crowding agents and the protein of interest are short - ranged and non - specific ( as is often the case ) , then the volume excluded by the crowding agents prevents the polypeptide from sampling extended conformations . as a consequence ,
the entropy of the denatured state ensemble ( dse ) decreases relative to the case when the crowding agents are absent .
these arguments suggest that excluded volume of crowding agents should enhance the stability of the folded state provided that the crowding - induced changes in the native state are negligible @xcite .
the entropic stabilization mechanism , described above , has been used in several theoretical models to quantitatively describe the extent of folded protein as a function of the volume fraction , @xmath13 , of the crowding agents @xcite .
more recently , a theory whose origins can be traced to the concept of intra - protein attraction due to depletion of crowding agents near the protein @xcite , predicts that the enhancement in stability , @xmath14 , where @xmath15 is the folding temperature at @xmath13 and @xmath16 is related to the flory exponent that characterizes the size of the protein in the dse @xcite . from this prediction
it follows that crowding affects the dse to a greater extent than the folded state .
although the precise theoretical predictions of the power law change in @xmath17 as @xmath13 changes have not been verified , several experiments using a number of proteins have confirmed that indeed @xmath18 increases with @xmath13 @xcite .
it can not be emphasized enough that the theory described here applies only to cases when the crowding interactions between crowding agents and proteins and between crowding particles themselves are purely repulsive .    while much less is known about the effects of crowding on the folding kinetics , cheung et al .
@xcite predicted that the entropic stabilization also suggests that the folding rates should increase at moderate values of @xmath13 .
they suggest that crowding can enhance folding rates by a factor @xmath19 where @xmath20 @xmath21 is the decrease in the entropy of the dse relative to its value in the bulk . from the arguments of cheung et al .
@xcite it follows that the equilibrium changes in the entropy ( @xmath20 ) of the dse , with respect to the bulk , should also determine rate enhancement provided that neither the barriers to folding @xcite nor the native state is perturbed significantly by crowding particles .
single molecule force spectroscopy , such as atomic force microscopy ( afm ) and laser optical tweezers , have been used to monitor the behavior of biopolymers under tension are ideally suited to probe the enhancement in crowding - induced stability by a direct measurement of @xmath22 .
indeed , ping et al .
@xcite have recently investigated the effect of dextran molecules on the mechanical stability of @xmath23 .
the 8 ub ( fig .
[ fig:1]*a * ) modules were n - c linked ( i.e. , modules i and i+1 were chemically linked together in a head - to - tail manner ) .
they found that the average force required to unfold a module , @xmath24 , increased by 21@xmath25 as the dextran concentration , @xmath26 , was increased from 0 to 300
g / l at @xmath27 pn / s .
similar results have been obtained recently by yuan et al .   at @xmath28
pn / s @xcite .
motivated in part by experiments @xcite , we used simulations to investigate the effects of crowding agents on the mechanical stability of a protein subject to external tension .
we focused on ubiquitin ( ub ) , a 76-residue protein composed of 5 @xmath29-strands and 2 @xmath30-helices ( fig .
[ fig:1]*a * ) , and confined our investigation to non - equilibrium ` force - ramp ' experiments @xcite .
the primary data recorded during such an experiment is a trace of the force exerted on the tip as a function of the extension of the molecule ; a force - extension curve ( fec ) .
when the force exceeds some critical value , the fec displays a sudden increase in length and is often accompanied by a concomitant sharp decrease in force .
presumably , the sharp change corresponds to the unfolding of the protein .
typical afm experiments use tandem arrays of proteins which are chemically linked together ( often through genetic engineering ) .
we use the term module to denote a protein of the array .
the fec resulting from such an experiment reveals several equally spaced peaks punctuated by sharp increases in the extension of the molecule corresponding to the unfolding of individual modules .
the height of these force peaks and their shape depend on the loading rate , @xmath31 , where @xmath32 is the cantilever
s spring constant and @xmath33 is the ( constant ) speed at which the stage is retracted away from the cantilever @xcite .    in order to compare to experiments our simulations
are performed using coarse - grained models for which simulations can be done at @xmath34 that are comparable to those in afm experiments . our work has led to a number of testable results : ( 1 ) at @xmath35 the average unfolding force for ub increases by at most only 7% compared to @xmath36 .
we find that @xmath24 in small crowding agents is greater than in larger particles .
( 2 ) in the presence of crowding agents , secondary structural elements reform multiple times even after initial rupture . (
3 ) although large crowding particles are predicted to have a smaller effect on @xmath24 ( for a given @xmath13 ) , they can profoundly affect the unfolding of poly ub .
we predict that @xmath24 for a given subunit depends on the number of already unfolded portions of the poly protein .
this result is important because many naturally occurring proteins that are subject to tensile stresses exist as tandem arrays of modules .
it further suggests that the existence of such redundancy can more properly be understood in the context of a crowded cellular milieu .
we used a coarse - grained model for proteins to investigate crowding effects on the mechanical stability of ub and @xmath37 at loading rates that are comparable to those used in afm experiments @xcite .
we assumed that ub could be described using the self - organized polymer ( sop ) model ; a model that has been successfully used to make a number of predictions regarding the unfolding of proteins and rna @xcite , allosteric transitions in enzymes @xcite , and movement of molecular motors on polar tracks @xcite .
previous studies @xcite have used more standard go - models @xcite to probe various aspects of forced unfolding of ub .
the sop energy function ( @xmath38 ) for a protein with n amino - acids , specified in terms of the c@xmath39 coordinates @xmath40 @xmath41 , is @xmath42+\sum _ { i=1}^{n-3 } \sum _ { j = i+3}^n \varepsilon _ h \left[\left(\frac{r_{ij}^0}{r_{ij}}\right)^{12}-2\left(\frac{r_{ij}^0}{r_{ij}}\right)^6\right]\delta _ { ij}\\ & + \sum _ { i=1}^{n-2 } \sum _ { j = i+2}^n \varepsilon _ l\left(\frac{\sigma } { r_{ij}}\right)^6\left(1-\delta _ { ij}\right ) , \end{split}\ ] ] where @xmath43 , @xmath44 is the value of @xmath45 in the native structure , @xmath46 kcal/(mol@xmath47@xmath48 ) , @xmath49 kcal / mol , @xmath50 kcal / mol , and @xmath51 nm . note that @xmath52 kcal / mol @xmath53 pn@xmath54 nm . in eq .
( [ eqn:1 ] ) @xmath55 if @xmath56 nm , and @xmath57 otherwise .
native coordinates corresponded to those of the c@xmath39 atoms of the 1.8  resolution protein data bank crystal structure 1ubq @xcite . for ub
@xmath58 and @xmath59 for @xmath37 . the first term in eq . ( [ eqn:1 ] ) is the fene potential @xcite that accounted for chain connectivity . the second ( lennard - jones ) term accounted for the non - bonded interactions that stabilize the native state , and the final ( soft - sphere ) term accounted for excluded - volume interactions ( including those of an angular nature ) .
the sop model is different from the go - model because there are no angular terms in sop , and the connectivity is enforced differently as well .
the sop representation of the polypeptide chain is in the same spirit as other coarse - grained models used in polymers @xcite .
we assumed that the crowding particles are spherical with diameter @xmath60 .
( @xmath61 nm in some simulations , while @xmath62 nm in others . )
crowders interacted amongst themselves and with the protein , respectively , via the following lj potentials : @xmath63    @xmath64    where @xmath65 , @xmath66 nm , @xmath67 , and @xmath68 .
the heaviside functions truncate the potentials at their minima and thereby ensured that only the repulsive portions of eqs .
( [ eqn:2 ] ) and ( [ eqn:3 ] ) contributed to interactions involving crowding agent .
_ mimics of crowding using asakura - oosawa theory : _ even using a coarse - grained sop representation of proteins , it is difficult to carry out converged simulations in the presence of crowding agents .
the reason is that the number of crowding agents can be large .
moreover , the separation in the spatial and temporal scales of the protein and the crowding particles has to be carefully considered to obtain reliable results . in light of these difficulties
, it is of interest to consider the effective attraction between the sites on the protein using the implicit pairwise potential computed by asakura and oosawa .
the intramolecular attraction arises due to the depletion of crowding particles near the protein . to probe the efficacy of these models we employed in some simulations the asakura - oosawa model @xcite of crowding effects . for these simulations , we added the following term to the bare sop hamiltonian ( eqn .
( [ eqn:1 ] ) ) : @xmath69 where @xmath51 nm , @xmath61 nm , @xmath35 , @xmath70 pn@xmath54 nm , and @xmath45 is the distance separating protein beads @xmath71 and @xmath72 .      for simulations involving @xmath37 , residues in different modules interacted via :    @xmath74    where r is the distance separating the two beads .
note that this potential is short - ranged and purely repulsive and that it is the same potential used for non - native intra - protein interactions @xcite .      _
@xmath36 _ : hundreds of simulations of @xmath75 steps ( @xmath76s ) at @xmath77 k were used to generate initial structures for use in the pulling simulations .
the protein was completely free in solution ( i.e. , no forces were applied to either terminus ) , and no crowders were present during the equilibrations .
the n - terminus of the protein was subsequently translated to the origin , and the protein was rotated such that its end - to - end vector , @xmath78 , ( i.e. , the vector pointing from the n - terminal bead to the c - terminal bead ) coincided with the pulling ( + z ) direction .
an unfolding trajectory was initiated by selecting a random initial structure from amongst the set of thermally equilibrated structures , and tethering a harmonic spring to the c - terminal bead .
the n - terminal bead was held fixed throughout the simulations .
tension was applied to the protein by displacing the spring along the + z axis and resulted in application of the following force to the c - terminal bead :    @xmath79-\left[z_s(t)-z_s(0)\right]\right),\ ] ]    where @xmath32 is the spring constant , @xmath80 , and @xmath81 corresponds to the displacement of the end of the spring .
note that @xmath32 was also used to constrain the simulation to the z - axis ; @xmath82 $ ] and @xmath83 $ ] .
the displacement of the spring was updated at every timestep .
we simulated forced - unfolding of monomeric ub at four different @xmath34 ( @xmath84 pn / s , @xmath85 pn / s , @xmath86 pn / s , and @xmath87 pn / s ) , while simulations on n - c - linked @xmath37 were performed at @xmath88 pn / s .
all overdamped force - ramp simulations were performed at the same speed @xmath89 nm / s , and spring constants were varied over a range from 0.3879 pn / nm - 31.032 pn / nm to achieve the aforementioned @xmath34 ( via the relation @xmath90 ) .
our simulations were realistic because they maintained loading rates consistent with experiment and because @xmath34 is the prime determinant of unfolding pathway @xcite .    _
@xmath91 0.0 _ : simulations involving explicit crowders were carried out at a fixed volume fraction @xmath35 and with a fixed number @xmath92 of crowding spheres .
( @xmath93 was fixed to render the problem computationally tractable ) . using the relation @xmath94
, we adjusted the length of a side of the cubic simulation box ( @xmath95 ) to maintain @xmath35 .
thus , @xmath96 nm when @xmath97 nm and @xmath98 nm when @xmath62 nm .
explicit crowders were added to the simulation after loading an equilibrated structure but before the application of tension .
initial crowder positions were chosen randomly and in a serial manner from a uniform distribution .
if the distance between an initial crowder position and that of another crowder or protein bead did not exceed the sum of their radii , then the prospective position was rejected and another random position chosen to avoid highly unfavorable steric overlaps .
periodic boundary conditions ( pbc ) and the minimum image convention @xcite were employed in the simulations .
two sets of coordinates were stored for protein beads at every timestep ; pbc were applied to one set and the other was propagated without pbc .
distances between protein beads were calculated from the uncorrected set of coordinates without minimum imaging , while protein - crowder distances were calculated from the pbc coordinates with minimum imaging .    to improve simulation efficiency ,
a cell list @xcite was used to calculate crowder - crowder and protein - crowder interactions .
the entire simulation volume was partitioned into 64 subvolumes , and it was only necessary to calculate interactions within a subvolume and between beads of the subvolume and those of 13 of its 26 neighbors .
the cell - list was updated at every timestep to ensure the accuracy of the simulations .
the equations of motion in our overdamped simulations ( used in all force - ramp simulations ) were integrated with a timestep @xmath99 ( @xmath100 ps ) using the method of ermak and mccammon @xcite .
the friction coefficient of the crowders , @xmath101 , was determined via the relation @xmath102 , where @xmath60 is the crowder diameter , @xmath51 nm is the diameter of a protein bead , and @xmath103 g / s is the friction coefficient associated with a protein bead of mass @xmath104 g. simulated - times were translated into real - times using @xmath105 @xcite . at @xmath77 k , @xmath106 ps , and since @xmath107 the real - time per step is 5.4306 ps .
at @xmath108 , forced unfolding of ub was simulated at four different @xmath34 ( @xmath87 pn / s , @xmath86 pn / s , @xmath85 pn / s , @xmath84 pn / s ) , where the lowest value corresponds approximately to the value used in the pulling experiments of ping et al .
@xcite , and all @xmath34 are experimentally accessible .    _ force profiles _ : fig .
[ fig:2](*a * and * b * ) provides examples of fec s collected at the highest and lowest @xmath34 .
we used a nominal contour length of ( @xmath109 ) @xmath110 nm and unfolding forces , @xmath111 , to determine contour - length increments , @xmath112 , for each trajectory at @xmath113 pn / s ( fig .
[ fig:2]*a * ) .
we identified @xmath111 with the peak of the fec before the stick - slip transition @xcite .
the average extension @xmath114 nm is in excellent agreement with the experimental result of 24 @xmath115 5 nm found by carrion - vazquez et al .
the projection ( @xmath116 ) of the end - to - end vector at @xmath111 in the z - direction varied between 4.1 nm and 4.7 nm , depending on @xmath34 . since the native end - to - end distance , @xmath117 nm , @xmath118 ranges from 0.4 - 1.0 nm .
the lower end of this range is slightly larger than the 0.25 nm transition - state distance for the mechanical unfolding of the structurally similar titin immunoglobulin domains @xcite .
indeed , we expect @xmath119 nm , because of the non - equilibrium nature of the simulations .
larger @xmath34 typically lead to larger @xmath120 ( @xmath121 ) .
average unfolding forces , @xmath122 , depended approximately logarithmically on @xmath34@xcite ( fig .  [ fig:2]*c * ) , and @xmath123 pn is in fair agreement with the experimental value of @xmath124 pn observed by ping et al .
@xcite at @xmath27 pn / s .
the @xmath125 also showed a logarithmic dependence on @xmath34 , but the difference between the value calculated at @xmath113 pn / s and that calculated at @xmath126 pn / s is small ( @xmath127 2  ) .
although the underlying free - energy landscape is time - dependent in a non - equilibrium force - ramp pulling experiment , these results suggest that the distance from the native state to the transition state is small .
it is likely that at loading rates that are achieved in laser optical tweezer experiments ( @xmath128 10 pn / s ) , the location to the transition state would increase because the response of biopolymers to loading rate changes from being plastic ( low @xmath34 ) to brittle ( high @xmath34 ) @xcite .
_ unfolding pathways _ : in the dominant pathway , unfolding proceeded in a fairly markovian fashion with the primary order of events following the sequence @xmath1291/@xmath1295
@xmath130 @xmath1293/@xmath1295 @xmath130 @xmath1293/@xmath1294 @xmath130 @xmath1291/@xmath1292 ( figs .
[ fig:1]*b * and [ fig:12 ] ) .
this is precisely the same sequence seen in the simulations by li et al.@xcite .
alternative pathways , reminiscent of kinetic partitioning @xcite observed in forced - unfolding of gfp @xcite and lysozyme @xcite , were also infrequently sampled . for example , at @xmath126 pn / s , @xmath131 of the trajectories unfolded as follows : @xmath1291/@xmath1295 @xmath130 @xmath1291/@xmath1292 @xmath130 @xmath1293/@xmath1295 @xmath130 @xmath1293/@xmath1294 ( fig .
[ fig:1]*b * ) , while the remaining @xmath132 followed the dominant pathway . at the highest loading rate ,
one frequently observed the following sequence of events
@xmath1291/@xmath295 @xmath130 @xmath1293/@xmath1295 @xmath130 @xmath1291/@xmath1292 @xmath130 @xmath1291/@xmath1292 ( reform ) @xmath130 @xmath1293/@xmath294 @xmath130 @xmath1291/@xmath1292 , where @xmath1291/@xmath1292 ruptured but then reformed prior to the rupture of @xmath1293/@xmath1294 . as illustrated in fig .
[ fig:2](*a * and * b * ) , unfolding events at smaller @xmath34 tended to result in larger molecular extensions .    the non - equilibrium character of a pulling experiment decreased with decreasing @xmath34 , and smaller @xmath34 resulted in smaller force - drops after the unfolding force , @xmath111 , is reached .
since the force applied by the spring to the end of the protein did not fall off as sharply at lower @xmath34 , more of the protein was extended during an unfolding event .
the smaller extensions following unfolding events at higher @xmath34 ( relative to those observed at lower @xmath34 ) were responsible for the @xmath1291/@xmath1292 unfolding / refolding events mentioned above because the applied tension was very low after the initial rupture event ( fig .
[ fig:2](*a * and * b * ) ) . at lower @xmath34 ,
this situation no longer held because the initial unfolding event resulted in a chain extension that was a significant fraction of the chain s contour length ( fig .
[ fig:2 ] * b * ) .
depletion forces stabilize proteins and shift the folding equilibrium towards more compact states @xcite .
these forces result from an increase in the entropy of the crowding agents that more than compensates for an increase in the free - energy of a protein molecule upon compaction .
simulations of forced - unfolding of ub in the presence of explicit crowders of diameters @xmath61 nm and @xmath62 nm were used to assess the contribution of the depletion forces to mechanical stability .
sixteen trajectories were collected for each @xmath34 investigated .
three @xmath34 ( @xmath86 pn / s , @xmath85 pn / s , @xmath84 pn / s ) , were explored for the @xmath61 nm sized depletants .
only the two highest @xmath34 ( @xmath85 pn / s , @xmath84 pn / s ) , were explored for the @xmath62 nm sized crowders .
_ small crowding particles increase unfolding forces _ :
example fec s collected in the presence of crowders of diameter @xmath61 nm and @xmath62 nm are illustrated in fig .
[ fig:9](*a * and * b * ) .
the @xmath61 nm curves in fig .
[ fig:9 ] look qualitatively very similar to those seen at @xmath108 , while the fec s collected at @xmath62 nm look qualitatively different .
for example , larger @xmath133 than those observed at either @xmath108 or in the presence of the @xmath61 nm crowders are apparent in the fec s collected at @xmath62 nm ( fig .
[ fig:9](*a * and * b * ) ) . indeed , fig .
[ fig:9]*c * reveals that this observation is quantitatively accurate . although the @xmath133 in the presence of the @xmath61 nm crowders were statistically indistinguishable from the @xmath133 at @xmath36 ( compare fig .
[ fig:9]*c * with fig .
[ fig:2]*c * ) , the average unfolding forces in the presence of the @xmath134 nm crowders were statistically greater than those measured at @xmath108 . at @xmath113 pn / s , the @xmath135 in the presence of the @xmath62
nm crowders exceeded that at @xmath108 by 3@xmath25 , while at @xmath136 pn / s the increase was 4@xmath25 ( compare fig .
[ fig:9]*c * with fig .
[ fig:2]*c * ) . in one respect
these results are not surprising ; it follows from the ao theory and eq .
( [ eqn:6 ] ) that smaller crowders stabilize ub more than larger ones . on the other hand ,
the extent of stabilization ( as measured by increases in @xmath133 ) was small .
we should emphasize that although the increase in the unfolding force is small , the stability change upon crowding is significant .
the enhancement in stability is @xmath137 @xmath138 using a 3% increase in the unfolding force , and @xmath139 the location of the unfolded basing @xmath140 nm .    _ crowding leads to transient local refolding _ : the unfolding pathways in the presence of crowding agents of both sizes were very similar to those seen at @xmath108 .
nevertheless , at @xmath35 there tended to be more unfolding / refolding events as the molecule extended past @xmath116 ( fig .
[ fig:12 ] ) than at @xmath108 .
as illustrated , strand - pairing between @xmath1291 and @xmath1292 and between @xmath1293 and @xmath1294 persisted to a greater extent after the initial unfolding event at @xmath35 than at @xmath108 .
as ub passed through the point ( @xmath116,@xmath111 ) , its termini were extended to distances greater than the diameter of even the larger crowders .
depletion forces resulting from the presence of 6.4 nm crowders act on these larger length scales . assuming that the tension applied to the c - terminus is small enough ( e.g. , after rupture events at higher loading rates ) ,
then such depletion forces can promote reformation of contacts between secondary structural elements several times during the course of a trajectory .
if this is indeed the case , then it suggests that unfolding polyub may be different from the unfolding of monomeric ub , because depletion effects should increase with number of modules in the tandem ( see below ) .
_ ao model for forced - unfolding in the presence of crowders _
: we used the asakura - oosawa ao model @xcite ( eq .  ( [ eqn:8 ] ) ) of the depletion interaction to model the effects of a crowded environment on ub .
the ao theory has been successfully used to model the effects of a crowded environment on polymers and colloids @xcite , and we found that it gives qualitatively accurate results for ub .
we used eq .
( [ eqn:8 ] ) to approximate the effective interaction between spherical protein beads immersed in a crowded solution of volume fraction @xmath35 . from the form of the ao potential between protein beads
, it follows that the range of the potential is proportional to @xmath60 , but the strength is greater for smaller crowders .
indeed , as revealed in the previous section , simulations in the presence of explicit crowders showed that 1.0 nm crowders resulted in larger average unfolding forces than 6.0 nm crowders ( fig .
[ fig:9]*c * ) .
although the ao - potential yields qualitatively accurate results , use of the ao - potential of eq .
( [ eqn:8 ] ) to implicitly model non - bonded interactions in ub did not yield quantitatively accurate results . at @xmath113 pn / s @xmath141 pn , while @xmath142 pn .
thus , simulations with the ao potential led to a mean unfolding force that is roughly 50@xmath25 greater than in its absence .
this disagrees sharply with the @xmath143 pn resulting from our own simulations in the presence of explicit crowding agent at @xmath35 at @xmath113 pn / s ( fig .
[ fig:9]*c * ) .
indeed , the result also stands in marked contrast to the experimental results of ping et al .
@xcite on octameric ub , which saw a maximum increase in @xmath133 of 21@xmath25 ( at @xmath144 , @xmath145 pn / s , and with @xmath146 nm ) over the @xmath147 pn at @xmath36 @xcite .
there are a couple of origins to the discrepancy .
first , the ao - potential was derived to understand the equilibrium of colloidal spheres and plates in the presence of smaller - sized spherical crowding agents .
our experiments were of a non - equilibrium nature , so it is somewhat unreasonable to expect such simulations to yield quantitatively accurate unfolding forces or dynamics .
second , as pointed out by shaw and thirumalai @xcite , three - body terms are required to properly model depletion effects even in good solvents let alone in the poor - solvent conditions of our simulations . to elaborate ,
let us consider the volume excluded to crowders by a ub molecule , v@xmath148(ub ) , to be the volume enclosed by a union of spheres of radii @xmath149 = @xmath150 . with an ao potential , the volume excluded to the spherical crowding agents
is :    @xmath151    where @xmath152 is the number of residues in monomeric ub ( 76 ) , @xmath153 is the volume of @xmath149 and @xmath154 is the volume associated with the overlap of @xmath149 and @xmath155 .
( [ eqn:9 ] ) neglects the overlap of three or more spheres .
the importance of such overlaps , for soft - spheres , increases as the crowders become much larger than the protein beads , and the neglect of such overlaps is the reason for the quantitative inaccuracy . as the size of the crowders increases ( i.e. , as the thickness of the depletion layer surrounding the protein beads increases ) , the surface enclosing ub becomes more spherical , loses detail , and undoubtedly changes less in response to changes in the conformation of the molecule . since depletion forces are proportional to the change in @xmath156(ub ) with respect to changes in ub conformation , eq .  ( [ eqn:9 ] )
overestimates the size of depletion forces in the presence of large crowders . despite these limitations the ao model , which is simple , can be used to provide qualitative predictions .
finally , we note that experiments typically use polyproteins to study force induced unfolding ( e.g. , ping et al .
@xcite used ( ub)@xmath157 in their experiments ) . for polyproteins
the quantitative accuracy of the ao theory for unfolding in the presence of crowders of diameter @xmath158 nm is likely to increase , because the volume excluded to the crowders will change significantly with changes in the conformation of the polyprotein .      from arguments based on volume exclusion that lead to crowding - induced entropic stabilization of the folded structures , it follows that crowding effects should be more dramatic on poly ub than the monomer . in order to illustrate the effect of crowding on stretching of ( ub)@xmath0 we chose @xmath61 nm , which had negligible effect on @xmath2 for the monomer .
however , we found significant influence of the large crowding particles when ( ub)@xmath0 was forced to unfold in their presence .
the fec ( fig .  [ fig:14 ] ) shows three peaks that corresponded to unfolding of the three domains , when simulations were performed at crowder volume fractions @xmath108 and @xmath35 at @xmath88 pn / s .
[ fig:13 ] presents average unfolding forces as a function of the unfolding event number .
although the first two unfolding events were statistically indistinguishable at @xmath108 and @xmath35 , the final event occurred at much larger @xmath133 in the presence of crowders than in their absence .
_ order of unfolding was stochastic : _
@xmath37 has 3 chemically identical modules . in the pulling simulations
the n - terminus of module @xmath159 was held fixed , while force was applied to the c - terminus of module @xmath160 .
figure [ fig:20 ] illustrates the time - dependence of contacts between secondary structure elements of modules @xmath159 , @xmath161 , @xmath160 .
it is clear from fig .
[ fig:20]*a * that module @xmath160 unfolded first , followed by module @xmath161 , and finally by module @xmath159 . on the other hand , the order of events in fig .
[ fig:20]*b * was @xmath162 .
the frequency with which the @xmath163 possible permutations of these orders at @xmath36 and at @xmath35 were observed is presented in table [ tbl : moduleorderfrequencies ] .
it is clear that at both @xmath36 and at @xmath35 , the most probably order of events was @xmath164 . at @xmath36
this order was only marginally more probable than the order @xmath165 , while at @xmath35 @xmath164 became overwhelmingly more probable than any other unfolding order .
in none of the simulations at @xmath36 or @xmath35 , did the rupture of module @xmath160 occur as the final event .
_ unfolding within a module depended on proximity to the point of force application : _ at @xmath166 , @xmath37 is fairly brittle , and the rupture of contacts within a module occured nearly simultaneously . nevertheless , by carefully examining time - dependent contact maps such as those illustrated in fig .
[ fig:20 ] , we were able to determine ( 1 ) at both @xmath36 and @xmath35 , @xmath167 contacts were the first to rupture , and ( 2 ) only for module @xmath160 was this rupture event _ invariably _ followed by the loss of @xmath168 contacts .
when other modules ruptured , loss of @xmath167 contacts was occasionally followed by loss of the @xmath169 strand - pair contacts .
@xmath37 _ must achieve a larger _
@xmath7 _ to rupture at _ @xmath35 : figure [ fig:21]*a * illustrates the time dependence of @xmath170 at @xmath36 and at @xmath35 .
interestingly , the plot reveals that after the second rupture event , the @xmath171 increased more rapidly in the presence of crowders than in their absence .
this is likely a reflection of the fact that at @xmath3 modules @xmath159 and @xmath160 were the first two modules to unfold in 44% of the trajectories , while at @xmath35 these two modules were the first to unfold in only 19% of trajectories .
thus , @xmath171 increased more rapidly in the presence of crowders , because the @xmath7 of @xmath37 with two adjacent modules unfolded is larger than that with two unfolded but non - adjacent modules . interestingly , these differences are masked in the time - dependent increase of the end - to - end distance ( fig .
[ fig:21]*b * ) .
figure [ fig:21]*a * also reveals that the horizontal inflection points marking the third unfolding event occur at different times at @xmath172 and in the absence of crowders .
the difference between these two times was responsible for the difference in average unfolding forces of @xmath173 pn illustrated in fig .
[ fig:13 ] .
thus , it is clear that despite the highly non - equilibrium nature of this pulling experiment , depletion effects were substantial and it required a much greater force to reach @xmath111 in the presence of these crowders than in their absence .
we predict that systematic experiments will reveal that polyub molecules composed of larger numbers of modules will show greater increases in @xmath24 relative to @xmath174 than polyub molecules composed of fewer repeats .
the size of these differences should increase with decreasing loading rate .
finally , it may even be possible to observe differences in the @xmath133 as a function of unfolding event number ( as in fig .
[ fig:13 ] ) .
the increase in @xmath24 , at a fixed @xmath13 , for poly ub is likely to be even more significant for small crowding agents as shown in eq .
( [ eqn:6 ] ) ( see below for further discussion ) .
general theory , based on the concept of depletion effects ( see eqs .  ( [ eqn:8 ] ) and ( [ eqn:6 ] ) ) , shows that crowding should enhance the stability of proteins , and hence should result in higher forces to unfold proteins .
however , predicting the precise values of @xmath24 is difficult because of the interplay of a number of factors such as the size of the crowding agents and the number of amino acid residues in the protein . despite the complexity
a few qualitative conclusions can be obtained based on the observation that , when only excluded volume interactions are relevant , then the protein or polyprotein would prefer to be localized in a region devoid of crowding particles @xcite .
the size of such a region @xmath175 .
if @xmath176 then the crowding would have negligible effect on the unfolding forces .
the condition @xmath176 can be realized by using large crowding particles at a fixed @xmath1 . in the unfolded state ,
@xmath177 nm @xcite which for ub leads to @xmath178 nm .
thus , @xmath179 . these considerations suggest that the crowder with @xmath158 nm would have negligible effect on the unfolding force , which is in accord with the simulations .
on the other hand , @xmath180 when @xmath181 nm , and hence we expect that the smaller crowders would have measurable effect on the unfolding forces
. our simulations are in harmony with this prediction .
we expect that for the smaller crowding agent @xmath2 would scale with @xmath1 in a manner given by eq .
( [ eqn:6 ] ) . in general ,
appreciable effect of crowding on the unfolding forces can be observed only for large proteins or for polyproteins using relatively small crowding agents .
although we have only carried out simulations for ub and ( ub)@xmath0 at one non - zero @xmath13 , theoretical arguments can be used to predict the changes in @xmath24 as @xmath13 increases .
the expected changes in the force required to unfold a protein can be obtained using a generalization of the arguments of cheung et al.@xcite in the presence of crowding agents the protein is localized in a region that is largely devoid of the crowding particles @xcite .
the most probable size of the region is @xmath182 where @xmath4 is the size of the crowding agent .
if the structures in the dse are treated as a polymer with no residual structure then the increase in entropy of the dse upon confinement is @xmath183 where @xmath7 is the dimension of the unfolded state of the protein .
if native state stabilization is solely due to the entropic stabilization mechanism we expect : @xmath184 where @xmath22 is the critical force for unfolding the protein , @xmath185 is the gain in contour length at the unfolding transition , and @xmath186 relates @xmath7 to the number of amino acids through the relation @xmath187 ( @xmath188 varies between 2 - 4  ) .
a few comments regarding eq .
( [ eqn:6 ] ) are in order .
( 1 ) the @xmath13 dependence in eq . ( [ eqn:6 ] )
does not depend on the nature of the most probable region that is free of crowding particles .
as long as the confining region , which approximately mimics the excluded volume effects of the macromolecule , is characterized by a single length d , we expect eq .
( [ eqn:6 ] ) to be valid .
( 2 ) the additional assumption used in eq . ( [ eqn:6 ] )
is that @xmath189 , and hence there may be deviations due to finite size effects .
( 3 ) the equivalence between crowding and confinement breaks down at large @xmath13 values .
consequently , we do not expect eq .
( [ eqn:6 ] ) to fit the experimental data at all values of @xmath13 .
( 4 ) it follows from eq .
( [ eqn:6 ] ) that , for a given @xmath7 , small crowding agents are more effective in stabilizing proteins than large ones .
thus , the prediction based on eq .  [ eqn:6 ] is supported by our simulations . ( 5 ) from the variation of @xmath24 with @xmath13 ping et al .
@xcite suggest that @xmath190 .
however , the large errors in the measurements can not rule out the theoretical prediction in eq .
( [ eqn:6 ] ) .
we have successfully fit their experimental results using eq .
( [ eqn:6 ] ) ( fig .
[ fig:21]*c * ) .
additional quantitative experiments are required to validate the theoretical prediction .
it is difficult to map the concentrations in g / l used in the study of ping et al .
@xcite to an effective volume fraction because of uncertainties in the molecular weight of dextran used in the study .
hence , a quantitative comparison between theory and experiments is challenging .
a naive estimate may be obtained using the values reported by weiss et al .
given that the dextran used in the study is thought to have an average molecular weight of 40 kda and an estimated average hydrodynamic radius of 3.5 nm @xcite , we find that @xmath191
g / l corresponds to a volume fraction @xmath192 , which is very large .
nevertheless , @xmath13 must be large when @xmath191
alternatively , we estimated @xmath193 for dextran using @xmath194 where @xmath152 is the number of monomers in a 40 kda dextran is @xmath195 .
if the monomer size @xmath196 nm , then we find @xmath197 .
if we assume that @xmath191
g / l corresponds to @xmath198 and that @xmath199 ( eq .  ( [ eqn:6 ] ) ) , then at @xmath35 we would expect a nearly 18@xmath25 increase in @xmath133 . similarly , if @xmath191 g
/ l corresponds to @xmath200 then we would expect an increase in @xmath133 of approximately 21@xmath25 . in any case
, we can say that at physiologically relevant volume fractions ( @xmath13 @xmath201 [ 0.1,0.3 ] ) , the percent increase in @xmath133 is likely to be @xmath202 .
our simulations for @xmath203 nm predict an increase of 3 - 4% , which shows that a more detailed analysis is required to obtain an accurate value of @xmath4 for dextran before a quantitative comparison with experiments can be made .
the larger increase seen in experiments may also be a reflection of the use of ( ub)@xmath157 rather than a monomer .
regardless of the crowder size we find that the unfolding pathways are altered in the presence of crowding agents .
it is normally assumed that the rupture of secondary structure elements is irreversible if the applied force exceeds a threshold value . however , when unfolding experiments are carried out in the presence of crowding particles , that effectively localize the protein in a smaller region than when @xmath36 , reassociation between already ruptured secondary structures is facilitated as shown here .
thus , forced - unfolding can not be described using one dimensional free energy profiles with @xmath116 as the reaction coordinate @xcite .
we find that the average unfolding force for the final rupture event of the unfolding of @xmath37 occurred at much larger values in the presence of crowders than in their absence . with @xmath61 nm , which has practically no effect on the unfolding force of the monomer , and @xmath35
even with unfolding of two modules the interactions between the stretched modules and protein are small ( @xmath204 ) . only upon unfolding of the third ub
do crowding effects become relevant , which leads to an increase in @xmath24 .
our results suggest that @xmath205 should increase with the number of modules in the array and that it may be possible to detect differences in @xmath133 which are conditional on the unfolding event number .
we speculate that naturally occurring polyproteins that are subject to mechanical stress have evolved to take advantage of precisely such enhanced depletion effects .
this work was supported by a grant from the national science foundation ( che 05 - 14056 ) .
dlp is grateful for a ruth l. kirschstein postdoctoral fellowship from the national institute of general medical sciences ( f32gm077940 ) .
computational time and resources for this work were kindly provided by the national energy research scientific computing ( nersc ) center .
10    zhou ,  h.  x. ;  rivas ,  g. ;  minton ,  a.  p. _ annu .
biophys . _ * 2008 , * _ 37 , _ 375 - 397 .
minton ,  a.  p. _ biophys .
j. _ * 2005 , * _ 88 , _ 971 - 985 .
cheung ,  m.  s. ;  klimov ,  d. ;  thirumalai ,  d. _ proc .
* 2005 , * _ 102 , _ 4753 - 4758 .
zhou ,  h.  x. _ acc . chem .
* 2004 , * _ 37 , _ 123 - 130 .
asakura ,  s. ;  oosawa ,  f. _ j. chem . phys .
_ * 1954 , * _ 22 , _ 1255 - 1256 .
asakura ,  s. ;  oosawa ,  f. _ j. polym .
sci . _ * 1958 , * _ 33 , _ 183 - 192 .
vrij ,  a. _ pure appl .
chem . _ * 1976 , * _ 48 , _ 471 - 483 .
sasahara ,  k. ;  mcphie ,  p. ;  minton ,  a.  p. _ j. mol .
biol . _ * 2003 , * _ 326 , _ 1227 - 1237
.    stagg ,  l. ;  zhang ,  s.  q. ;  cheung ,  m.  s. ;  wittung - stafshede ,  p. _ proc .
* 2007 , * _ 104 , _ 18976 - 18981 .
homouz ,  d. ;  perham ,  m. ;  samiotakis ,  a. ;  cheung ,  m.  s. ;   wittung - stafshede ,  p. _ proc .
* 2008 , * _ 105 , _ 11754 - 11759 .
cheung ,  m.  s. ;  thirumalai ,  d. _ j. phys .
b _ * 2007 , * _ 111 , _ 8250 - 8257 .    ping ,  g. ;  yang ,  g. ;  yuan ,  j .- m . _
polymer _ * 2006 , * _ 47 , _ 2564 - 2570 .
yuan ,  j.  m. ;  chyan ,  c.  l. ;  zhou ,  h.  x. ;  chung ,  t.  y. ;  haibo ,  j. ;   ping ,  g. ;  yang ,  g. _ prot .
_ published online before print * sep 9 , 2008 * doi:10.1110 .
fisher ,  t. ;  marszalek ,  p. ;  fernandez ,  j. _ nature structural biology _ * 2000 , * _ 7 , _ 719 - 724 .
evans ,  e. ;  ritchie ,  k. _ biophys .
j. _ * 1997 , * _ 72 , _ 1541 - 1555 .
hyeon ,  c. ;  thirumalai ,  d. _ biophys . j. _ * 2006 , * _ 90 , _ 3410 - 3427 .
hyeon ,  c. ;  dima ,  r. ;  thirumalai ,  d. _ structure _ * 2006 , * _ 14 , _ 1633 - 1645 .
hyeon ,  c. ;  lorimer ,  g. ;  thirumalai ,  d. _ proc . natl .
* 2006 , * _ 103 , _ 18939 - 18944 .
chen ,  j. ;  dima ,  r. ;  thirumalai ,  d. _ j. mol .
biol . _ * 2007 , * _ 374 , _ 250 - 266 .
hyeon ,  c. ;  onuchic ,  j. _ proc .
natl . acad .
* 2007 , * _ 104 , _ 17382 - 17387 .
best ,  r.  b. ;  paci ,  e. ;  hummer ,  g. ;  dudko ,  o.  k. _ j. phys . chem .
b _ * 2008 , * _ 112 , _ 5968 - 5976 .
karanicolas ,  j. ;  brooks iii ,  c.  l. _ prot .
sci . _ * 2002 , * _ 11 , _ 2351 - 2361 .
li ,  m. ;  kouza ,  m. ;  hu ,  c. _ biophys . j. _ * 2007 , * _ 92 , _ 547 - 561 .
vijay - kumar ,  s. ;  bugg ,  c. ;  cook ,  w. _ j. mol .
_ * 1987 , * _ 194 , _ 531 - 544
.    kremer ,  k. ;  grest ,  g. _ j. chem .
* 1990 , * _ 92 , _ 5057 - 5086 .
pincus ,  d.  l. ;  cho ,  s.  s. ;  hyeon ,  c. ;  thirumalai ,  d. minimal models for proteins and rna : from folding to function . in _
progress in nucleic acid research and molecular biology _ ; conn ,  p.  m. ,  ed . ;
( in press ) .
allen ,  m.  p. ;  tildesley ,  d.  j. _ computer simulation of liquids ; _ oxford university press : 1987 .
ermak ,  d. ;  mccammon ,  j. _ j. chem . phys . _
* 1978 , * _ 69 , _ 1352 - 1360 .
veitshans ,  t. ;  klimov ,  d. ;  thirumalai ,  d. _ folding and design _ * 1996 , * _ 2 , _ 1 - 22 .
carrion - vazquez ,  m. ;  li ,  h. ;  lu ,  h. ;  marszalek ,  p. ;   oberhauser ,  a. ;  fernandez ,  j. _ nature structural biology _ * 2003 , * _ 10 , _ 738 - 743 .
bockelmann ,  u. ;  essevaz - roulet ,  b. ;  heslot ,  f. _ phys .
e _ * 1998 , *
_ 58 , _ 2386 - 2394 .
carrion - vazquez ,  m. ;  oberhauser ,  a. ;  fowler ,  s. ;  marszalek ,  p. ;   broedel ,  s. ;  clarke ,  j. ;  fernandez ,  j. _ proc .
usa _ * 1999 , * _ 96 , _ 3694 - 3699 .
hyeon ,  c. ;  thirumalai ,  d. _ j. phys .
matt . _ * 2007 , * _ 19 , _ 113101 .
guo ,  z. ;  thirumalai ,  d. _ biopolymers _ * 1995 , * _ 36 , _ 83 - 102 .
mickler ,  m. ;  dima ,  r. ;  dietz ,  h. ;  hyeon ,  c. ;  thirumalai ,  d. ;   rief ,  m. _ proc .
* 2007 , * _ 104 , _ 20268 - 20273 .
peng ,  q. ;  li ,  h. _ proc .
* 2008 , * _ 105 , _ 1885 - 1890 .
shaw ,  m.  r. ;  thirumalai ,  d. _ phys .
a _ * 1991 , * _ 44 , _ r4797-r4800 .
verma ,  r. ;  crocker ,  j. ;  lubensky ,  t. ;  yodh ,  a. _ phys .
lett . _ * 1998 , * _ 81 , _ 4004 - 4007 .    toan ,  n. ;  marenduzzo ,  d. ;  cook ,  p. ;  micheletti ,  c. _ phys .
_ * 2006 , * _ 97 , _ 178301(4 ) .
thirumalai ,  d. _ phys . rev .
a _ * 1988 , * _ 37 , _ 269 - 276 .
kohn ,  j.  e. ;  millet ,  i.  s. ;  jacob ,  j. ;  zagrovic ,  b. ;   dillon ,  t.  m. ;  cingel ,  n. ;  dothager ,  r.  s. ;  seifert ,  s. ;   thiyagarajan ,  p. ;  sosnick ,  t.  r. ;  hasan ,  m.  z. ;  pande ,  v.  s. ;   ruczinski ,  i. ;  doniach ,  s. ;  plaxco ,  k.  w. _ proc .
* 2005 , * _ 101 , _ 14475 - 14475 .
weiss ,  m. ;  elsner ,  m. ;  kartberg ,  f. ;  nilsson ,  t. _ biophys . j. _ * 2004 , * _ 87 , _ 3518 - 3524 .
humphrey ,  w. ;  dalke ,  a. ;  schulten ,  k. _ j. mol .
graph . _ * 1996 , * _ 14 , _ 33 - 38 .
.module unfolding order frequencies at @xmath36 and @xmath35 . [ cols="^,^,^",options="header " , ]     [ tbl : moduleorderfrequencies ]
fig .  [ fig:1 ] : ( * a * ) cartoon representation of the native structure of ubiquitin ( pdb accession i d 1ubq ) in the presence of spherical crowding agents .
the five beta - strands , labeled @xmath1291 through @xmath1295 , are colored in yellow .
the two alpha - helices ( @xmath301 and @xmath302 ) are shown in purple .
the n- and c - terminal beads are represented as spheres . in our simulations
the n - terminal bead was held fixed while the c - terminal bead was pulled via a tethered spring .
( * b * ) snapshots from an unfolding trajectory illustrating the main ubiquitin ( ub ) unfolding pathway ( brown - dashed arrows ) and an alternate unfolding pathway ( green - dotted arrow ) . in both pathways
the initial unfolding event corresponds to separation of the c - terminal strand @xmath1295 from the n - terminal strand @xmath1291 . along the main pathway , this is quickly followed by separation of @xmath1295 from @xmath1293 . the penultimate rupture event along the main pathway corresponds to disruption of the @xmath1293/@xmath1294 strand - pair , while the n - terminal @xmath1291/@xmath1292 strand pair is the last to break .
the trajectory illustrated here was generated at @xmath108 and @xmath113 pn / s .
an alternate pathway was observed at @xmath108 and @xmath126 pn / s . along this pathway separation of @xmath1295 from @xmath1291 is followed by separation of @xmath1291/@xmath1292 .
the final two rupture events correspond to those of the @xmath1295/@xmath1293 contacts and @xmath1293/@xmath1294 strand - pair respectively .
( figures generated with vmd @xcite )    fig .
[ fig:2 ] : force - extension curves ( fec s ) at two different loading rates , ( * a * ) @xmath206 pn / s , and ( * b * ) @xmath126 pn / s .
data from the simulation is presented as a red trace . for each trajectory a black arrow points to the unfolding force , @xmath111 .
@xmath116 corresponds to the extension of the molecule along the pulling ( i.e. , z- ) direction evaluated at @xmath111 .
@xmath112 ( dotted blue line ) is the contour length increment , and is a measure of the amount of chain released in an unfolding event .
we measured @xmath112 as @xmath207 , where @xmath208 nm is a nominal contour length of 28.5 nm .
stars in each subfigure mark the minimum force observed after an unfolding event , and chain conformations corresponding to the starred points in figures * a * and * b * are illustrated at the center of the figure . unfolding events at smaller @xmath34
resulted in larger molecular extensions before significant resistance was encountered .
( yellow arrows correspond to beta - strands and purple cylinders correspond to @xmath30-helices ) .
( figures generated with vmd @xcite ) .
( * c * ) @xmath133 vs.  @xmath34 evaluated at @xmath36 .
the red curve corresponds to a linear - least squares fit to the set of basis functions @xmath209 and demonstrates that @xmath210 .
( note that the abscissa is a log - scale ) .
each point is labeled mean @xmath115 standard error .
statistics at @xmath113 pn / s , @xmath85 pn / s , @xmath86 pn / s , and @xmath87 pn / s were calculated from 50 , 49 , 50 , and 16 trajectories respectively .
[ fig:12 ] : rupture events at @xmath108 ( * a * ) and at @xmath35 and @xmath61 nm ( * b * ) .
the figure illustrates that after the initial rupture event at step 0 , subsequent unfolding was affected by the crowding agent .
@xmath1291/@xmath1292 contacts and @xmath1293/@xmath1294 persisted longer and there were many more local refolding events in the presence of the crowders ( * b * ) than in their absence ( * a * ) .
( * c * ) snapshots from an unfolding trajectory at @xmath35 and @xmath126 pn / s illustrate the primary difference between unfolding at @xmath35 and at @xmath108 .
the brown - dotted arrow shows unfolding without any local - refolding events , while the sequence of black arrows illustrates local - rupture and -refolding events .
the hallmark of forced - unfolding in a crowded environment is repeated breaking and reforming of contacts after an initial rupture event .
( figures generated with vmd @xcite )    fig .
[ fig:9 ] : ( * a * ) examples of force - extension traces resulting from simulation in the presence of spherical crowding agents of diameter @xmath61 nm and obtained at @xmath136 pn / s .
( * b * ) examples of force - extension traces resulting from simulation in the presence of spherical crowding agents of diameter @xmath62 nm and obtained at @xmath136 pn / s .
both subfigures are labeled as in fig .
[ fig:2](*a * and * b * ) .
( * c * ) @xmath133 vs.  @xmath34 evaluated at @xmath35 .
black triangles and red circles correspond to spherical crowders of diameter @xmath62 nm and @xmath61 nm respectively .
each point is labeled mean @xmath115 standard error .
statistics for each point were calculated from 16 independent trajectories . only the @xmath62 nm had an appreciable effect on @xmath133 when compared to those obtained at identical @xmath34 and at @xmath108 ( see fig .  [ fig:2]*c * )    fig .
[ fig:14 ] : fec s for @xmath37 forced unfolding at @xmath108 ( * a * ) and at @xmath35 ( * b * ) .
trajectories were generated at @xmath88 pn / s and with crowders of diameter @xmath61 nm .
black arrows mark each trajectory s three unfolding events .
[ fig:13 ] : @xmath133 vs.  unfolding event number for the unfolding of @xmath37 in the presence of spherical crowders of diameter @xmath61 nm ( @xmath35 , black triangles ) and in their absence ( @xmath108 , red circles ) .
the individual modules of the polyub tandem were n - c - linked and the loading rate was @xmath211 pn / s . each point is labeled mean @xmath212 standard error . statistics for each point were calculated from 16 independent unfolding trajectories .
an unfolding event corresponded to the unfolding of an individual module .
note that although the crowders had little effect on @xmath133 for the first and second unfolding events , they had a substantial effect on the last unfolding event .
@xmath133 increased by @xmath173 pn for the last unfolding event .
[ fig:20 ] : illustration of the stochastic nature of module unfolding . in (
* a * ) module @xmath160 ( the most proximal to the applied force ) unfolds first , followed by module @xmath161 , and finally by rupture of module @xmath159 . in ( * b
* ) the order of module unfolding events is @xmath162 .
the two trajectories illustrated here were both collected at @xmath36 and @xmath88 pn / s .
table [ tbl : moduleorderfrequencies ] provides the frequencies at which the different possible orders were observed at @xmath36 and at @xmath35 .
the diameter of the crowders is @xmath158 nm .    fig .
[ fig:21 ] : ( * a * ) @xmath213 versus time at @xmath36 ( red ) and at @xmath35 ( black ) .
although the @xmath170 increases more rapidly with time after the second rupture event at @xmath35 than at @xmath36 , the inset reveals that a larger @xmath7 must be achieved to initiate the final rupture event in the presence of crowding particles with @xmath158 nm .
( see text for additional discussion ) .
( * b * ) @xmath214 versus time at @xmath36 ( red ) and at @xmath35 ( black ) .
@xmath215 can not discriminate between unfolding at the different volume fractions and hence is less suitable ( than @xmath7 ) as a potential reaction coordinate .
( * c * ) the experimental results of ping et al .
@xcite ( red circles ) for unfolding - force , @xmath133 as a function of dextran concentration , @xmath216 . the black line is a fit assuming @xmath217 and the blue assuming @xmath218 .
although both fits are consistent with the data , based on theoretical considerations ( see eq .
( [ eqn:6 ] ) ) we prefer the blue fit ( see text ) .
( standard deviations taken from ping et al .