knotted proteins provide an ideal ground for examining how amino acid interactions ( which are local ) can favor their folding into a native state of non - trivial topology ( which is a global property ) .
some of the mechanisms that can aid knot formation are investigated here by comparing coarse - grained folding simulations of two enzymes that are structurally similar , and yet have natively knotted and unknotted states , respectively . in folding simulations that exclusively promote the formation of native contacts , neither protein forms knots .
strikingly , when sequence - dependent non - native interactions between amino acids are introduced , one observes knotting events but only for the natively - knotted protein .
the results support the importance of non - native interactions in favoring or disfavoring knotting events in the early stages of folding .
a much debated problem in the thermodynamics and kinetics of protein folding @xcite is the process of knot formation in proteins  @xcite .
the topological entanglement found in naturally - occurring proteins presents several differences from that of compact or collapsed flexible polymers  @xcite .
firstly , the overall percentage of knotted native states in the protein data bank ( pdb ) is lower than for globular flexible chains with the same degree of polymerization  @xcite .
in addition , the free energy landscape of proteins , unlike that of homopolymers or random heteropolymers , is sufficiently smooth to ensure that the same knot in the same protein location is observed in folded structures  @xcite .
these distinctive features are arguably the consequence of concomitant effects , including the propensity of polypeptides to form secondary structure elements , which enhance their local order compared to a typical collapsed polymer chain  @xcite , as well as functionally - oriented evolutionary mechanisms @xcite .    to gain insight into such mechanisms as well as to highlight general physico - chemical mechanisms favoring knot formation , an increasing number of experimental and numerical studies of knotted proteins have been carried out @xcite .
notably , the latest experimental results have clarified that native knots can form spontaneously and efficiently from an unknotted initial state . in particular , yeates and coworkers
@xcite successfully designed a protein which can refold into the natively - knotted structure , albeit much more slowly than an unknotted counterpart .
finally , the recent study of mallam and jackson showed that newly translated , knot free , yibk molecules fold spontaneously to the native trefoil - knotted state , and that chaperones can significantly speed up the folding process  @xcite .
numerical simulations have been a valuable complement to experiments , particularly regarding the characterization of the pathways leading to the self - tying of knotted proteins . in particular , folding simulations based on simplified protein representations and/or force fields , have indicated two main mechanisms leading to knot formation : the threading of one of the termini through a loop  @xcite , or by slipping a pseudoknot through a loop  @xcite .
the two mechanisms are not mutually exclusive , as reported by sulkowska _
@xcite for protein mj0366 from _
methanocaldococcus jannaschii_. importantly , in coarse - grained ( c@xmath0-trace ) simulations where folding was promoted by exclusively favoring native contacts it was seen that the yield of folding trajectories was low and knots would form at late folding stages . specifically , for the deeply - knotted protein yibk , only 1 - 2% of the trajectories reached the native state and the native knot was formed when about @xmath1 of the native contacts had already been established .
arguably , the coarse - grained nature of the models contributes to the observed low folding yield .
however , even when employing atomistic representations , purely native - centric models tend to favor a rather late onset of knotting  @xcite . for example , a relatively late stage of knot formation @xmath2 was also estimated from the analysis of the free - energy landscape of protein mj0366 , despite being half as long as yibk and with a shallower knot  @xcite .
these findings are aptly complemented by the early observation of wallin _ et al .
_ @xcite that the low yield of purely native - centric models can be dramatically enhanced by promoting the attraction of specific protein regions that are not in contact in the native state . in particular , it was seen that the formation of the native knot in yibk was particularly facilitated by introducing an _ ad hoc _ non - native attraction between a loop and the c - terminus that , by threading the former , produced a sizeable population of knots even at early folding stages , @xmath3% @xcite .    these results , together with recent experimental ones @xcite ,
pose the question of whether natively - knotted proteins can form non - trivial entanglements during early folding stages and whether non - native interactions can provide a general mechanisms for such self - tying events .
clarifying such aspects is important to advance the understanding of some of the general mechanisms aiding knot formation .
it must however , be borne in mind that there could co - exist independent pathways where knots are established at different stages of the folding process and that active mechanisms , such as interactions with chaperones , could be involved _ in vivo _
@xcite .
in this study we address these questions by simulating the early folding process of two transcarbamylase proteins that are structurally very similar and yet their native states are differently knotted  @xcite . specifically , one of them is a trefoil - knotted n - acetylornithine carbamoyltransferase ( aotcase ) , while the other is an unknotted ornithine carbamoyltransferase ( otcase ) , see fig .
[ figure_1 ] .
their evolutionary relatedness has posed the interesting question of understanding the source of their different native topology , with particular regards to the role of specific loop regions whose `` virtual '' excision or addition alters the native topology @xcite .    here , for both the otcase and aotcase ,
more than a hundred folding simulations are carried out using a coarse - grained model and two different energy functions .
we first evolve a fully extended configuration and steer it towards the native state by exclusively promoting native contact interactions . in a second set of simulations we enrich the energy function by adding quasi - chemical and electrostatic non - native interactions . in both cases ,
the stochastic evolution is followed up to the formation of about 35% of the native contacts and the knotted topology of the partially - folded structures is monitored .
we find that knot occurrence is negligible for both the aotcase and otcase when the energy function favoring only native contacts is used . a dramatic difference in knotting propensity
is instead seen when the quasi - chemical non - native interactions are added . in this case , the level of self - entanglement remains negligible for the otcase but is greatly enhanced for the aotcase , which is natively knotted .
notably , all proper and improper knots formed in these early folding stages have the correct ( native ) trefoil topology and chirality .
the analysis of the repeated knotting / unknotting events observed in the simulated trajectories indicates that knotting usually results from the threading of the c - terminal through loops present in the loose protein globule .
the threading events are favored by the effective ( non - native ) attraction of the c - terminus to other protein regions .
additional folding simulations carried out for `` in silico mutants '' of the aotcases , provide further support for the effect of the c - terminus chemical composition on the knotting propensities .
the present investigation of knot formation in proteins that are only partially folded is focused on two transcarbamylases : aotcase , pdb entry 2g68 , and otcase , pdb entry 1pvv .
the monomeric units of these two proteins have nearly the same length ( 332 and 313 amino acids , respectively ) and are evolutionarily related .
this is established from their significant sequence homology ( same cath code , 3.40.50.1370 @xcite and sequence identity equal to @xmath440% ) . despite their statistically - significant structural alignability ( mistral p - value @xmath5
@xcite ) they have a different knotted topology @xcite , see fig .  [ figure_1 ] .
specifically , the otcase is unknotted , while the aotcase contains a right - handed trefoil knot .
the knotted region , established using the algorithm of ref .
@xcite , corresponds to the k172-g255 segment and is at a distance of about 80 amino acids from the c terminus , which is the nearest one in sequence . because the biological units of the two proteins have a different oligomeric state it is assumed that , consistently with other multimeric knotted proteins , the folding of the monomeric units precedes their assembly .
we shall therefore limit considerations to the monomeric units .
for the purpose of our study , both proteins are described with a simplified structural model , where each amino acid is represented by one interaction center , coinciding with the c@xmath0 atom .
the folding dynamics is simulated using the stochastic monte carlo ( mc ) approach introduced and validated in ref .
the mc evolution entails a brief relaxation from the initial fully - extended state using both local and non - local moves . after this stage , indicated with a shaded region in fig.[figure_2 ] and
related ones , the mc dynamics proceeds exclusively through local moves . the moves amplitude is sufficiently small that no chain crossings occurs , as verified _ a posteriori .
the stochastic scheme is analogous to the kink - jump dynamics  @xcite which , by virtue of the local character of the moves , is known to provide a physically - viable description of biopolymers kinetics in thermal equilibrium  @xcite . _
the folding simulations are carried out using two alternative energy functions .
the first one promotes exclusively the formation of native contacts between pairs of amino acids . following refs .
@xcite , the strength of the attractive interaction between any given native amino acid pair is derived from the strength of their hydrogen bonding in the native state and is expressed in thermal units , @xmath6 , at the nominal monte carlo temperature of @xmath7k .
the second energy function includes electrostatic and non - native pairwise interactions in addition to the native - centric potential .
the relative strength of the non - native interactions is set according to the quasi - chemical potentials of miyazawa and jernigan ( mj ) @xcite that reflect the statistical propensity of amino acid pairs to be in contact in proteins native states .
following , again , ref.s
@xcite , the average strength of the added non - native interactions is set to one tenth of the native one .
no change to the mc temperature was made using this second potential because previous studies have shown that the addition of non - native interactions modifies the effective temperature of the system by less than @xmath8  @xcite .
consistently with this fact we have verified that the native states of both aotcase and otcase remain stable when evolved with the mc scheme using the second type of potential .
in fact , the asymptotic value for average fraction of native contacts was @xmath9% .
the fact that this value is smaller than the asymptotic one , @xmath10% , of the first type of potential indicates that the strength non - native interactions while weak enough to maintain stable the native state , can nevertheless compete with the native attractive interactions .
the addition of the quasi - chemical interactions are therefore expected to be capable of accounting for sequence - specific non - native contact propensities in the partly - folded state .
the sequence - dependent character of the non - native interaction differentiates the present approach from the early one of wallin _ et al . _
@xcite where the knot - promoting effect of non - native interactions was probed by systematically introducing attractive interactions between various pairs of segments of yibk .    for each of the two proteins and for each of the two energy functions
, we generated @xmath4150 mc trajectories each consisting of @xmath11 mc moves per amino acid ( corresponding to 3000 units of mc time reported in fig .
[ figure_2 ] and fig .
[ figure_3 ] ) , for a total of 75,000 cpu hours .
the acceptance rate of the local mc moves for both proteins was nearly constant throughout the simulation and approximately equal to 50% .    the entanglement of the structures sampled during the mc evolution was established using a combination of two unrelated knot detection schemes , which are described in the materials and methods .
this choice was made to maximize the robustness of the criterion used to establish the knotted state of an open chain - which is mathematically properly defined only after the chain termini are joined by a segment or arc thus giving a circular chain ( chain closure operation ) .
as described in the materials and methods section we used many alternative closures for each chain and adopted a majority rule to single out chains that have non trivial self - entanglements .
these chains can correspond to proper , fully developed knots ( fully accommodated in the knotted structures ) as well as improper ones ( partly accommodated in the closing arcs ) .      in order to illustrate how the folding process differs for the two potential energy functions , we show in fig .
[ figure_2 ] the evolution of the gyration radius ( upper panel ) and the fraction of native contacts ( lower panel ) of the knotted aotcase , obtained by employing these two potentials .
it is seen that the effect of the non - native and electrostatic interactions results in the higher compactness of the protein globule , since the gyration radius in this case lies systematically below the one calculated with the purely native - centric model . at any given stage of the mc evolution ,
the fraction of established native contacts is systematically lower for non - native interactions thus indicating that the latter introduce frustration that slows down the folding dynamics .
analogous results hold for the unknotted otcase , see top panel of fig .
[ figure_3 ] .    the data in fig .
[ figure_3 ] represent how the average radius of gyration and fraction of formed native contacts evolve in the course of the folding simulations of the two transcarbamylases .
it is seen that , over the duration of the simulation , the average fraction of formed native contacts grows to @xmath12% on average ( the maximum value in all trajectories was @xmath13% .
overall , the probability of formation of native contacts decreases with the sequence separation of the amino acid pair , so that native @xmath14-helical contacts are substantially more probable than inter - strand ones @xcite .    from the steady increase of the fraction of native contacts , @xmath15 ,
it is extrapolated that folding completion would occur on time - scales at least ten times larger than considered here . from the same figure
it emerges that the radius of gyration gradually decreases to about 25 @xmath16 which is 30% larger than the native one , consistently with the less compact character of the partially folded states .    to characterize the overall propensity of the two proteins to form knots while still largely unfolded ,
we extracted configurations at regular intervals of the simulations and analyzed their topological state .
the knotting propensity of aotcase , which is natively knotted , is illustrated in the upper panel of fig .
[ figure_4 ] which portrays the fraction of configurations that are properly or improperly knotted during a given time - window of the mc dynamics .
it is seen that when the purely native - centric potential is used , the knotting propensity is always negligible .
however , when non - native interactions are added , the number of knotted configurations increases to a definite fraction of the total : at the end of the simulations , when only 35% of the native contacts are formed , the fraction of conformations that are properly or improperly knotted is about 1% .
notice that this value is comparable with the yield of purely native - centric coarse - grained models where pathways with a late - stage formation of knots are typically observed @xcite . within the limitations of coarse - grained approaches , the results indicate that sequence - dependent non - native interactions can produce a detectable fraction of configurations with the correct native entanglement already at early folding stages .
the analysis was repeated for the natively - unknotted otcase and the results are shown in the lower panel of fig .
[ figure_4 ] . the contrast with the case of the atocase is striking .
in fact , the knotting probability is negligible not only for the purely - native case , but remains so even when non - native interactions are introduced .
the results indicate that the sequence - dependent non - native interactions  promoted by the quasi - chemical potential
increase dramatically the incidence of knots in the partially folded state of the aotcase compared to the otcase .
equivalently , the model calculations indicate that the two related natively - knotted and unknotted proteins have different knotting propensities already in the partially folded state , and that this difference can be ascribed to sequence - specific non - native interactions .      to further elucidate the mechanisms responsible for these differences we monitored several parameters in the course of the simulations of the two transcarbamylases .
more specifically , we searched for systematic differences in the interactions that the termini of these proteins establish with other parts of the peptide chain and for preferential locations of fully - developed knots .    in this respect ,
it is important to point out that , in the course of the simulations for the aotcase with non - native interactions , the process of knot formation is not irreversible . as illustrated in fig .
[ figure_5 ] , trefoil knots can be formed and untied during each simulation . in particular
, fully - developed knots can persists for up to one tenth of the duration of our simulations .
notably , only knots with the correct ( i.e. native ) chirality are observed .    by analysing the configurations preceding and following the knotting / unknotting events it is found that the change in topological state is typically caused by the threading of the helical c - terminus through loops formed by various protein regions .    among the independent knotting events we identified six for which the stochastic closure returned a non - trivial topology in at least 80% of the cases or the portion accommodating the knot had a depth of at least 20 amino acids .
the analysis of the loop threading events leading to such persistent fully - developed knots showed that they involved the c terminus interacting with amino acids 55 - 65 , 90 - 110 , or 125 - 155 . only in one case
it was observed the threading of the n - terminal through a loop involving segments 252 - 277 .
representative configurations are shown in the left panel of fig .
[ figure_6 ] where it is clear that the helical c - terminus of the aotcase typically points towards the rest of the protein , and the stiffness of the helix facilitates the threading of various loop regions .    at variance with the above situation , for both potential energy functions the otcase c
- terminus is typically exposed to the solvent and pointing away from the rest of the protein , as in the example shown in the right panel of fig .
[ figure_6 ] .
the terminus also shows a lower propensity to form @xmath14-helices .
an analogous situation is found for the aotcase in the presence of only native interactions .    the effect is quantitatively illustrated in fig .
[ figure_7 ] , which portrays the average strength of the non - native interaction potential energy between the amino acids in the c - terminal @xmath14-helix of the aotcase or the otcase and all the other residues in the chain .
it is seen that the average non - native attraction is consistently stronger in the case of the knotted protein by about 15% , which is a sizable amount considering that the average is taken over all sampled structures , irrespective of their compactness and knotted topology .
the same figure also illustrates that , if the amino acids in the c - terminal @xmath14-helix of protein aotcase are replaced by neutral hydrophilic residues ( glu , gln , asn ) , then the average non - native interaction energy becomes slightly repulsive , explaining the tendency of the @xmath14-helix to point outwards .
this provides an additional argument in support of the picture where the hydrophobicity character of the residues plays a role in the folding of knotted proteins  @xcite .    finally , in fig .
[ figure_8 ] we plot the average electrostatic potential energy between the @xmath14-helix and the rest of the chain and show that , overall , the electrostatic interaction is about two orders of magnitude smaller than that of the quasi - chemical potential .    to further clarify the role of the c - terminal @xmath14-helix in the knotting of the partially - folded states of the aotcase , we have carried out another set of simulations on a mutant protein .
the last 25 residues of the otcase ( which are involved in the c - terminal @xmath14-helix ) were substituted with the same number of residues that form the c - terminal @xmath14-helix in the knotted aotcase .
apart from such a substitution , all other attributes , namely the length of the chain and the native contact map , were identical to the unknotted protein .
the outcome of the simulation is summarized in fig .
[ figure_9 ] , from which we can conclude that the mutant protein forms a substantial fraction of knotted configurations in the presence of non - native interactions in its early stage of the folding .      a simplified coarse - grained protein model was used for a comparative study of early folding stages of two evolutionarily related transcarbamylase proteins : an aotcase ( pdb i d 2g68 ) and an otcase ( pdb i d 1pvv ) .
the two proteins are well - alignable in sequence and structure and yet possess differently knotted native states : aotcase is trefoil - knotted while otcase is unknotted .
the role of sequence - dependent non - native interactions in promoting the correct native topology in the early folding stages was investigated by using two different energy functions in the simulations : one with a purely native - centric potential energy ( favoring only native contact interaction ) and one with the additional contribution of quasi - chemical and electrostatic non - native interactions .
we found that in the absence of quasi - chemical interactions , neither protein shows an appreciable propensity to self - entangle in the early folding stages .
however , once non - native interactions are introduced , the natively - knotted aotcase does show a strongly enhanced propensity to form proper and improper knots with the correct native topology and chirality .
by contrast , the knot enhancement effect is completely absent in the natively - unknotted otcase .
inspection of the ensemble of folding trajectories of the two proteins suggest that knotting in the partly unfolded aotcase mostly results from the approach of the hydrophobic c - terminal @xmath14-helix to other regions of the protein which are eventually threaded through . because of the non - compact character of the partly - unfolded structures , the @xmath14-helix in the c - terminal can also retract from the threaded regions , so that various events of formation / disruption of proper and improper knots can be observed in a given trajectory .
the influence of the c - terminal region on the different knotting propensity of the two carbamylases is further supported by the folding simulations for an `` in silico mutant '' of the otcase obtained by replacing the c - terminal sequence with the one of the knotted aotcase . in fact , the folding simulations based on the quasi - chemical non - native interactions yielded a portion of knotted structures similar to the natively - knotted aotcase , from which the c terminus was taken .
we recall that , the seminal study of wallin _ et al . _
@xcite , which was based on a coarse - grained model of yibk with _ ad hoc _ non - native interactions , had suggested the relevance of non - native interactions for steering the early formation of the native knotted topology ( by the threading of a loop by the c - terminal ) .
our results , which are consistent with these conclusions , provide additional elements to the picture by taking advantage of a coarse - grained framework where the folding of two carbamylases with similar structures , and yet different knotted topology are compared on equal footing . in particular our results indicate that non - native interaction propensities that are encoded in the primary sequence , can favor or disfavor the formation of knots already at early folding stages .
the necessarily limited scope of the coarse - grained approach used here does not clarify whether additional pathways leading to a late - stage formation of knots can be present in transcarbamylases .
it would be most interesting to address this standing issue in future studies within the present , or alternative comparative schemes .
we have adopted the coarse - grained model developed in ref.s  @xcite , in which the effective degrees of freedom are the amino acid residues , represented by a spherical bead located at the position of the corresponding c@xmath17 atom .
the potential energy of the model consists of bonded and non - bonded terms .
the bonded part of the of the potential consists of stretching potentials of the pseudo - bond and pseudo - angle , as well as potential for pseudo - torsions : @xmath18    @xmath19 the bond - stretching potential has the form @xmath20 where @xmath21 is the distance between the residues @xmath22 and @xmath23 , while @xmath24 is the equilibrium length of the c@xmath17-c@xmath17 pseudo - bond .
@xmath19 the double - well pseudo - angle potential is given by @xmath25 ,     \label{equation_angle}\ ] ] where @xmath26 is the pseudo - angle formed by the residues @xmath22 , @xmath23 and @xmath27 , while @xmath28 and @xmath29 are the equilibrium values of the helical and the extended pseudo - angles , respectively .
@xmath19 the torsion - angle potential for pseudo - torsions is of the form @xmath30
v_\text{n } , \label{equation_torsion}\ ] ] with @xmath31 being dihedral angle between the planes ( @xmath22,@xmath23,@xmath27 ) and ( @xmath23,@xmath27,@xmath32 ) .
the constants @xmath33 and @xmath34 are depending only on the type of the middle two residues @xmath23 and @xmath27 and are adopted from karanicolas and brooks  @xcite .
all the values of the spring and energy constants appearing in the expressions of the bonded energy terms in eqs .
( [ equation_bond ] ) , ( [ equation_angle ] ) and ( [ equation_torsion ] ) can be found in ref .
@xcite .    the non - bonded part of the potential consists of three terms , comprising native , non - native and electrostatic interactions , respectively : @xmath35    @xmath19 @xmath36 denotes the g@xmath37-type potential developed by karanicolas and brooks  @xcite .
within such an approach , the native contact map is defined on the basis of the network of the hydrogen bonds in the native state , as well as on the degree of the proximity of the backbone atoms side - chains .
namely , two residues are defined to be in the native contact if the hydrogen bond between them is stronger than -0.5 kcal / mol , or if any of their non - hydrogen side - chain atoms are within the distance of 4.5 @xmath16 , in the native state .
we do not consider contacts between residues with a distance in sequence smaller than 3 amino - acids .
these interactions are described by the following functional form : @xmath38,\ ] ] where @xmath39 is the native - state separation of residues @xmath22 and @xmath23 .
the strength of the interaction @xmath40 is chosen to be that of their hydrogen bond , for the residues that are hydrogen - bonded in the native state , while for the side - chain interacting residues it is a value proportional to the corresponding mj contact potential  @xcite and it was suitably renormalized in order to match the hydrogen bond native contact energy scale .
the residues in @xmath41-sheets or hairpins are often in contact via multiple hydrogen bonds , so in order to stabilize these structures within the g@xmath37-model , the additional network of weaker hydrogen bonds ( having the strength @xmath42 ) involving the neighboring pairs @xmath43 , @xmath44 , @xmath45 and @xmath46 was introduced around each residue pair @xmath47 that was found to interact via either two hydrogen bonds or a hydrogen bond and a side - chain contact .
@xmath19 depending on the amino acid type of the residues @xmath22 and @xmath23 , the non - native interactions in the model developed by kim and hummer in ref .
@xcite can be both attractive and repulsive .
repulsive interactions are applied between amino acid pairs that interact less favorably with each other than with the solvent and _ vice versa_. for a pair of residues that experiences an effective attractive interaction with a strength @xmath48 , the non - native interaction potential is given by @xmath49\,,\ ] ] with @xmath50 , where @xmath51 and @xmath52 are van der waals radii of the residues @xmath22 and @xmath23 .    for pairs of residues that effectively repel each other , so that @xmath53 , the non - native potential energy function takes the following form : @xmath54 + 2 \varepsilon_{ij } , & r_{ij } <
\sqrt[6]2\,\sigma_{ij}\\                                     -4 \varepsilon_{ij } \bigg [                                        \big(\frac{\sigma_{ij}}{r_{ij}}\big)^{12}-                                        \big(\frac{\sigma_{ij}}{r_{ij}}\big)^{6 }                                      \bigg ] , & r_{ij } \geq \sqrt[6]2\,\sigma_{ij }                                   \end{array }                                   \label{veff }                                   \right .\ ] ]    the effective lennard - jones interaction strength @xmath40 between residues @xmath22 and @xmath23 in the coarse - grained model of kim and hummer  @xcite is defined as @xmath55 the coefficients @xmath56 are negative and coincide with the entries of the mj matrix @xcite , while @xmath57 is an offset parameter .
hence , eq .  ( [ veff ] ) defines a statistical knowledge - based potential which measures the preference of residue - residue interactions relative to residue - solvent interactions .
the parameter @xmath58 scales the strength of the lennard - jones interaction compared to the physical electrostatic interactions .
the two free parameters @xmath58 and @xmath57 are fitted in order to correctly reproduce the binding affinity of the broad set of experimentally well - characterized protein complexes  @xcite . as an illustration ,
the effective interaction strength between the neutral hydrophilic residues glu and gln is @xmath59 kcal / mol , so that effective interaction is repulsive , while for the hydrophobic residues ile and leu is @xmath60 kcal / mol , so that the effective interaction is attractive .
@xmath19 @xmath61 is the long - range electrostatic interaction between residues @xmath22 and @xmath23 and it is modeled by debye - h@xmath62ckel - type of the potential @xmath63 } { r_{ij}}\,,\ ] ] where @xmath64 and @xmath65 are the electrostatic charges of residues @xmath22 and @xmath23 , @xmath66 is the debye screening length , @xmath67 is dielectric vacuum constant and @xmath68 is the relative dielectric constant of water in near - ambient conditions .
the non - native interactions in the force field of our coarse - grained model introduce frustration and make the potential energy surface quite rugged . as a consequence , even within such a simple model , molecular dynamics ( md ) simulations of the folding for chains of several hundreds amino acids are very expensive . within the simulation time intervals which were accessible to our computer resources
, such md trajectories did not allow to monitor significant changes in the chain conformations .    to cope with this problem , we have simulated mc dynamics , using an algorithm which combines different types of moves , namely :    @xmath19 local crankshaft moves @xcite , that consist of the rotation of a randomly selected single bead around the axis defined by its nearest neighbors .
the angle of the rotation was randomly selected in the interval @xmath69 , +    @xmath19 local end - point moves , in which the last 10 residues on both terminals are rotated rigidly with respect to the rest of the chain by up to @xmath70 around a random axis passing through the most interior bead of the end segment .
@xmath19 local cartesian moves , that involve the displacement of the coordinates of a single randomly selected bead in the chain , within a sphere of radius @xmath71 , +    @xmath19 global pivot moves @xcite , where one amino acid is picked at random and the chain portion involving all amino acids with smaller ( or alternatively larger ) sequence index are rotated by up to @xmath70 around a random axis passing through the picked amino acid .
the moves were accepted or rejected according to the standard metropolis criterion .
mc algorithms based on local crankshaft and end moves are commonly employed in the polymer physics  @xcite to study _ dynamic properties _ , since it is was shown that they can mimic the intrinsic dynamics of a polymer in solution  @xcite at a much lower computational cost of md simulations  @xcite .
however , we emphasize that since our computational scheme is based on the mc evolution , we can not establish a quantitative mapping between the simulation time and physical time .
we have generated 150 independent mc trajectories for each system , starting from the same stretched coil configuration .
in the early stage of the folding , these types of moves were attempted with the constant probabilities of 0.8 , 0.1 and 0.1 , respectively .
such an algorithm generates a rapid collapse of the chain from a fully stretched configuration to one in which the gyration radius is reduced from about 70 @xmath16 to about 30 @xmath16 in about @xmath72 mc steps per particle , that corresponds to 100 units of mc time ( see fig .  [ figure_3 ] ) . during such collapse ,
no knot was observed . at the end of this stage , the acceptance rate of
the global pivot moves typically dropped below 1% .
after this happened , we switched off the pivot moves , except from those involving few residues near the the two terminal of the chain . from this point on ,
the conformational changes of the chain are driven by the local crank - shaft and cartesian moves with an overall acceptance of about 50% .
when computing the fraction native contacts @xmath15 along the calculated mc trajectories , we adopt a criterion according to which two residues with index difference @xmath73 are said to be in contact if the distance of their @xmath74 is less than @xmath75 .
the conformations visited during the monte carlo dynamics were topologically classified by computing the alexander determinants after suitable closure into a ring @xcite .    for robustness
, two alternative closure schemes are used : the minimally - interfering closure  @xcite and a modified version of the stochastic one  @xcite .
the minimally - interfering closure is used first because it is very computationally effective , in that it entails a single , optimally chosen closure . in case of positive knot detection
we further validate the non - trivial entanglement by performing 100 closures where each terminus is prolonged far out of the protein along a stochastically chosen direction , and the end of the prolonged segments are closed by an arc ( that does not intersect the protein ) .
the stochastic exit directions are picked uniformly among those that are not back - turning .
specifically , they must form an angle of more than 90@xmath76 with the oriented segment going from each terminus to the @xmath74 at a sequence distance of 10 .
if the majority of the stochastic closures return non - trivial alexander determinants , than the conformation is non - trivially entangled .
such conformations can correspond to both proper , fully developed knots , and improper ones .
the two can be distinguished using knot localization criteria  @xcite : proper knots are entirely accommodated within the original protein chain , while improper ones span the exit segments .
the knot type was determined using the scheme of ref.@xcite , which is based on the knotfind algorithm .
we thank p. carloni , a.v . finkelstein and l. tubiana for valuable discussions and suggestions .
pf is also a member of the interdisciplinary laboratory for computational science ( lisc ) , a joint venture between university of trento and fondazione bruno kessler .
t is supported by the _ provincia autonoma di trento _ , through the aurorascience project .
simulations were performed on the aurora supercomputer at the lisc .
99 finkelstein av , ptitsyn ob ( 2002 ) , protein physics : a course of lectures , academic press .
creighton te ( 1992 ) , proteins : structures and molecular properties , w. h. freeman .
mansfield ml ( 1994 ) are there knots in proteins ?
nature struct biol 1 : 213 - 214 . mansfield ml ( 1997 )
fit to be tied .
nature struct biol 4 : 166 - 167 .
taylor wr ( 2000 ) a deeply knotted protein structure and how it might fold .
nature 406 : 916 - 919 .
nureki o , shirouzu m , hashimoto k , ishitani r , terada t _
( 2002 ) an enzyme with a deep trefoil knot for the active - site architecture .
acta crystallog d 58 : 1129 - 1137 .
michel g , sauve v , larocque r , li y , matte a , cygler m ( 2002 ) the structure of the rlmb 23s rrna methyltransferase reveals a new methyltransferase fold with a unique knot .
structure 10 : 1303 - 1315 .
taylor wr , lin k ( 2003 ) a tangled problem .
nature 421 : 25 .
lua rc , grosberg ay ( 2006 ) statistics of knots , geometry of conformations , and evolution of proteins .
plos comput biol 2 : doi : 10.1371/journal.pcbi.0020045 .
delbruck m ( 1962 ) knotting problems in biology
. mathematical problems in the biological sciences , edited by r e bellman , proceedings of symposia in applied mathematics , american mathematical society , providence , rhode island 14 : 55 - 63 .
frish hl , wasserman e ( 1961 ) chemical topology .
journal of american physical society 83 : 3789 - 3795 .
vologodskii av , lukashin av , frank - kamenetskii md , anshelevich vv ( 1974 ) the knot problem in statistical mechanics of polymer chains .
sov phys - jetp 39 : 1059 - 1063 .
sumners dw , whittington sg ( 1988 ) knots in self - avoiding walks .
journal of physics a - mathematical and general 21 : 1689 - 1694 .
janse van rensburg ej , sumners dw , wasserman e , whittington sg ( 1992 ) entanglement complexity of self - avoiding walks . journal of physics a - mathematical and general 25 : 6557 - 6566 .
micheletti c , marenduzzo d , orlandini e , sumners dw ( 2006 ) knotting of random ring polymers in confined spaces .
j chem phys 124 : 064903 .
michels jpg , wiegel fw ( 1986 ) on the topology of a polymer ring .
proc r soc lond a403 : 269 - 284
. micheletti c , marenduzzo d , orlandini e ( 2011 ) polymers with spatial or topological constraints : theoretical and computational results .
phys rep 504 : 1 .
marenduzzo d , orlandini e , stasiak a , sumners dw , tubiana l and micheletti c ( 2009 ) dna  dna interactions in bacteriophage capsids are responsible for the observed dna knotting proc natl acad sci usa 106 : 22269 - 22274 .
taylor wr ( 2007 ) protein knots and fold complexity : some new twists .
comput biol chem 31 : 151 - 162 .
mallam al , jackson se ( 2012 ) knot formation in newly translated proteins is spontaneous and accelerated by chaperonins .
nat chem biol 8 : 147 - 53 .
king np , jacobitz aw , sawaya mr , goldschmidt l , yeates to ( 2008 ) structure and folding of a designed knotted protein .
proc natl acad sci u s a 107 : 20732 - 7 .
yeates to , norcross ts , king np ( 2007 ) knotted and topologically complex proteins as models for studying folding and stability .
current opinion in chemical biology 11 : 595 - 603 .
mallam al ( 2009 )
how does a knotted protein fold ?
febs j 276 : 365 - 375 .
virnau p , mallam a , jackson s ( 2011 ) structures and folding pathways of topologically knotted proteins . j phys : condens matter 23 : 33101 - 33118 .
noel jk , sulkowska ji , onuchic j ( 2010 ) slipknotting upon native - like loop formation in a trefoil knot protein .
proc natl acad
sci usa 107 : 15403 - 15408 .
virnau p , mirny la , kardar m ( 2006 ) intricate knots in proteins : function and evolution .
plos comput biol 2 : doi : 10.1371/journal.pcbi.0020122 .
potestio r , micheletti c , orland h ( 2010 ) knotted vs. unknotted proteins : evidence of knot - promoting loops .
plos comput biol 6 : doi : 10.1371/journal.pcbi.1000864 .
shakhnovich e ( 2011 ) protein folding : to knot or not to knot ?
nat mater 10 : 84 - 6 .
mallam al , jackson se ( 2005 ) folding studies on a knotted protein .
j mol biol 346 : 1409 - 1421 .
mallam al , jackson se ( 2006 ) probing nature s knots : the folding pathway of a knotted homodimeric protein .
j mol biol 359 : 1420 - 1436 .
mallam al , jackson se ( 2006 ) a comparison of the folding of two knotted proteins : ybea and yibk .
j mol biol 366 : 650 - 665 .
mallam al , jackson se ( 2008 ) knotted fusion proteins reveal unexpected possibilities in protein folding .
molecular cell 30 : 642 - 648 .
mallam al , morris er , jackson se ( 2009 ) exploring knotting mechanisms in protein folding .
proc natl acad sci usa 105 : 18740 - 18745 .
mallam al , rogers jm , jackson se ( 2010 ) experimental detection of knotted conformations in denaturated proteins .
proc natl acad sci usa 107 : 8189 - 8194 .
wallin s , zeldovich kb , shakhnovich ei ( 2007 ) the folding mechanics of a knotted protein .
j mol biol 368 : 884 - 893 .
sulkowska ji , sulkowski p , onuchic j ( 2008 ) dodging the crisis of folding proteins with knots .
proc natl acad
sci usa 106 : 3119 - 3124 .
faisca pfn , travasso rdm , charters t , nunes a , cieplak m ( 2010 ) the folding of knotted proteins : insights from lattice simulations .
phys biol 7 : 16009 - 16020 .
orengo ca , michie ad , jones s , jones dt , swindells mb , et al .
( 1997 ) cath - a hierarchical classification of protein domain structures .
structure 5 : 1093 - 1109 .
micheletti c , orland h ( 2009 ) mistral : a tool for energy - based multiple structural alignment of proteins .
bioinformatics 25 : 2663 - 2669 .
kolesov g , virnau p , kardar m , mirny la ( 2007 ) protein knot server : detection of knots in their structures .
nucleic acids res 35 : w425-w428 ( doi : 10.1093/nar / gkm312 ) .
quake sr ( 1995 ) fast monte carlo algorithms for knotted polymers .
phys rev e 52 : 1176 - 1180 .
shuler ke , orwoll ra , stockmayer wh .
stochastic models for chain dynamics , in advances in chemical physics : stochastic processes in chemical physics , volume 15 ( ed k. e. shuler ) , john wiley & sons , inc . ,
hoboken , nj , usa . de pablo jj , escobedo fa ( 2007 ) monte carlo methods for polymeric systems , in advances in chemical physics : monte carlo methods in chemical physics , volume 105 ( eds i. prigogine and s. a. rice ) , john wiley & sons , inc . , hoboken , nj , usa . best rb , hummer g ( 2010 ) coordinate - dependent diffusion in protein folding .
proc natl acad sci usa 107:1088 - 1093 .
kim yc , hummer g ( 2008 ) coarse - grained models for simulations of multi - protein complexes : application to ubiquitin binding .
j mol biol 375 : 1416 - 1433 .
miyazawa s , jernigan rl ( 1996 ) residue - residue potentials with a favorable contact pair term and an unfavorable high packing density term for simulation and threading .
j mol biol 256 : 623 - 644 .
a beccara s , krbi t , covino r , and faccioli p ( 2012 ) dominant folding pathways of a ww domain .
proc natl acad sci usa 109 : 2330 - 2335 .
ivankov dn , garbuzynskiy so , alm e , plaxco kw , baker d , finkelstein av ( 2003 ) contact order revisited : influence of protein size on the folding rate .
protein sci 12 : 2057 - 2062 .
karanicolas j , brooks cl ( 2002 ) the origins of the asymmetry in the folding transition states of protein l and protein g. protein sci 11 : 2351 - 2361 .
kremer k , binder k ( 1988 ) monte carlo simulation of lattice models for macromolecules .
comput phys rep 7 : 259 - 310 .
madras n , sokal ae ( 1988 ) the pivot algorithm : a highly efficient monte carlo method for the self - avoiding walk .
j stat phys 50 : 109 - 186 .
jorgensen wl , tirado - rives j ( 1996 ) monte carlo vs molecular dynamics for conformational sampling .
j phys chem 100 : 14508 - 14513 .
tubiana l , orlandini e , micheletti c ( 2011 ) profiling the arc entanglement of compact ring polymers : a comparison of different arc - closure schemes with applications to knot localization .
prog theor phys 191 : 192 - 204 .
k. millett , a. dobay and a. stasiak ( 2005 ) linear random knots and their scaling behaviour .
macromolecules 38 : 601 - 606 .
tubiana l , orlandini e , micheletti c ( 2011 ) multiscale entanglement in ring polymers under spherical confinement .
phys rev lett 107 : 188302 .
traces of aotcase ; the n and c termini are colored in red and blue , respectively .
bottom panel : monte carlo evolution of the fraction of native contacts of aotcase .
the overall fraction is shown in red , while the the fraction of formed native contacts involved in @xmath14-helices and @xmath41-sheets are shown in green and blue , respectively.,width=453 ]            -helix and the rest of the chain , in the knotted aotcase and unknotted otcase . * the upper curve shows the same quantity for a mutant of the natively knotted protein in which all residues in the c - terminal @xmath14-helix are replaced by hydrophilic residues.,width=453 ]


cognitive radio ( cr ) is a promising approach to efficiently utilize the scarce rf spectrum resource . in this paradigm ,
knowledge about spectrum occupancy in time , frequency , and space that is both accurate and timely is crucial in allowing cr networks to opportunistically use the spectrum , and avoid interference to a primary user ( pu ) .
in particular , information about pu location will enable several key capabilities in cr networks including improved spatio - temporal sensing , intelligent location - aware routing , as well as aiding spectrum policy enforcement .
the localization problem in cr networks is in general different from localization in other applications such as wireless sensor networks ( wsn ) and global positioning system ( gps ) , in which the target to be localized cooperates with the localization devices .
in contrast , a pu does not communicate directly with the crs during the localization process .
this scenario is referred to as non - cooperative localization and does not allow the use of conventional techniques such as time - of - arrival ( toa ) and time - delay - of - arrival ( tdoa ) ranging .
several techniques have been proposed in the literature for estimating the position of active pus in a cr network .
these schemes include both range - based @xcite and range - free @xcite algorithms .
range - free positioning schemes , such as centroid localization @xcite , have attracted a lot of interest because of their simplicity and robustness to changes in wireless propagation properties such as path loss .
this characteristic makes them suitable candidates for systems requiring coarse - grained , but reliable and cost - effective techniques . on the other hand , range - based techniques such as multi - lateration @xcite ,
offer better estimates , but require accurate information about the target s distance that can only be achieved either through cooperation with the target or a very precise knowledge of the path loss model .    in this paper
, we focus on analyzing the performance of a particular , low complexity , coarse - grained localization algorithm , referred to as weighted centroid localization ( wcl ) @xcite . in this technique , pu location is approximated as the weighted average of all secondary user positions within its transmission range . strictly speaking ,
wcl is not an entirely range - free technique because it requires additional information aside from simple connectivity , namely received signal strength ( rss ) measurements .
wcl and wcl - like methods are suitable for pu localization because of their low computational complexity .
they also do not require cooperation from the pu and only use the readily available rss information .
in addition , wcl always converges to a solution , as opposed to some range - based techniques , and do not require estimation of the path loss exponent beforehand .
recent papers on wcl have proposed improvements on the basic algorithm to enhance its performance for specific scenarios @xcite .
all previous wcl studies reported in the literature have been analyzed either via simulation or field experiments @xcite . to the best of our knowledge , there has been no analysis of the probability distribution of the localization error for even the simplest version of wcl in the presence of shadowing in both wsn and cr literature .
we present the first theoretical analysis of the error distribution of wcl in this paper .
the presented analytical framework models varying levels of shadowing , including both independent and correlated shadowing environments , and also takes inaccuracy of sensor node positioning into consideration . using this analysis in conjunction with numerical simulations we observe and quantify the robustness of wcl to variable channel conditions and node placements .
we also investigate variable node participation as a simple modification that can improve wcl performance .    in the latter part of this paper we propose and investigate a distributed cluster - based implementation of wcl in the context of cr networks .
previous work primarily assumes a centralized fusion node aggregating data from other nodes . in the context of wsn where wcl was first utilized @xcite ,
the primary motivation was node self - localization where a node applies wcl using measurements from beacon nodes in close proximity to it .
this approach naturally dictates a centralized system where the node to be localized acts as the fusion center and only requires minimal interaction with the beacons . in the case of pu localization in cr networks
, however , where the primary motivation is estimating the position of a non - cooperative node , there no longer exists any particular node that has all the requisite information to perform wcl . in this case a fusion node is arbitrarily chosen .
this approach does not scale very well in terms of communication overhead and transmit power , which limits its application to low node densities or small network size .
there is therefore a need to create distributed implementations of wcl for this technique to be practically applied in cr localization .
a standard approach for distributing control in sensor networks is clustering .
several algorithms exist in the literature for creating either disjoint @xcite or overlapping @xcite clusters . in @xcite , a distributed technique based on the concept of mean - shift @xcite
. however , the technique requires repeatedly forming a virtual cluster until it converges to the target location .
this introduces additional communication overhead and also requires several iterations to converge .
we present an implementation which saves total transmit power needed for localization process communication , and investigate its accuracy as well as its effect on per node computational complexity and total transmit power consumption .
the rest of the paper is organized as follows .
the system model is introduced in section [ sec : systemmodel ] .
the theoretical analysis of centralized wcl accuracy is presented in section [ sec : analysis ] , while a practical and distributed implementation of wcl is proposed and analyzed in section [ sec : distributed ] .
numerical results for various scenarios used to evaluate performance of wcl algorithms are presented in section [ sec : evaluation ] .
finally , the paper is concluded in section [ sec : conclusion ] .
this section presents the general framework for our analysis .
we first review the assumptions about the system setting and propagation model , and then present the basic wcl algorithm .
we adopt the pu - centric model @xcite in which the pu transmitter is assumed to be at the center area of a circle .
the rationale behind this choice is that in a practical scenario , sensor nodes that obtain enough received power ( i.e. above sensitivity threshold ) approximately form a circle around the primary transmitter .
thus , sensors within this circle are used to jointly estimate the position of the pu . an implicit assumption in this model
is that the pu is surrounded by enough nodes so that edge effects can be disregarded .
assume @xmath0 sensor nodes in a circle of fixed radius @xmath1 .
the number of sensors within the circle is jointly affected by several factors , such as transmit power of the pu ( which effectively determines @xmath1 ) , node density , node selection methodology , etc .
@xmath0 is a design parameter for wcl deployment , the choice of which will have significant influence on the localization accuracy , as will be shown in our analysis .
the pu located at the center has a coordinate defined as @xmath2^t=[0 , 0]^t$ ] . for the @xmath3th sensor node , we define its 2-dimensional location as @xmath4^t$ ] where @xmath5 .
sensors are placed independently in both dimensions .
we assume sensors obtain information about their own locations , which is modeled as imperfect .
the measured location of the @xmath3th sensor is distributed as @xmath6 , where @xmath7 defines standard deviation of the error in sensor node self - localization .
we further assume the self - localization errors are independent among sensors and are independent from received powers .
we adopt the following channel model for received signal strength used in wcl .
the received power of the @xmath3th node from the primary transmitter , @xmath8 , is given by @xmath9 in which the first two terms characterize path loss and the last term describes shadowing effect .
the path loss is characterized by three constants : reference power @xmath10 , reference distance @xmath11 and path loss exponent @xmath12 .
denote collection of shadowing variables as @xmath13 $ ] , then its distribution is characterized by @xmath14 .
we consider the following shadowing cases    _ = _ s^2 * i_n , & i.i.d case , + \{_}_ij = _
i-_j /x_c , & correlated case , [ shadowing ] *    where @xmath15 is the correlation distance @xcite within which the shadowing effects of nodes are correlated , and @xmath16 is the identity matrix with dimension @xmath0 .    the theoretical analysis presented in section [ sec : analysis ] is based on the centralized architecture , which will also give us an upper bound on the performance . in each localization period , all sensors involved in the wcl will send their rss measurements to the fusion center , where an estimate will be made about the coordinate of the primary transmitter . the considered algorithm is a modification of relative - span wcl proposed in @xcite . in this scheme ,
the weights are guaranteed to be non - negative in db scale .
the 2-dimensional estimated location of pu is formed as @xmath17 the weight factor of node @xmath3 is defined as @xmath18 , @xmath19 where @xmath20 is a constant calculated from average received power of the node on the border with some margin , to guarantee the probability that receive power of all nodes in the area below @xmath20 is sufficiently small , say 1% .
@xmath21 is the maximum received power among nodes , and @xmath22 is the span of received power .
the localization error is then given by @xmath23^t$ ] , where @xmath24 and @xmath25 are the 1-dimensional errors along the x and y - axis respectively .
finally , the performance of the wcl scheme is evaluated by distance error , which is given by @xmath26
this section presents the analysis of wcl error using probability theory , ultimately leading to analytical expressions for the error distribution , first of wcl applied to 1-dimensional coordinates , and finally the 2-dimensional wcl error given by ( [ error - dist ] ) .        first we derive the 1-dimensional localization error for the i.i.d . shadowing case .
we define a constant @xmath27 for each node @xmath3 as @xmath28 one can verify that @xmath29 .
the location estimation in one dimension , say , x - axis , is naturally given by @xmath30 based on our i.i.d assumption among all @xmath8 s , one can verify that @xmath31 for the statistical property of @xmath32 , since @xmath33 s are also i.i.d and are independent of @xmath8 s , the mean of @xmath32 is given by @xmath34 .
the variance of the @xmath3th term of @xmath32 is given by    rcl [ variance a ] [ ( p_i - p _ ) x_i ] & = & [ ( p_i - p_)^2 ] [ x_i^2 ] - [ p_i - p_]^2 [ x_i]^2 + & = & _ l^2 _ i^2 + _ s^2 _ i^2 + _ l^2 _ s^2 .    thus the variance of @xmath32 is derived as @xmath35 . from central limit theorem ,
we approximate @xmath32 as a single gaussian variable , i.e. @xmath36 , which is verified by numerical simulations that are not presented here due to space limit .
since @xmath32 and @xmath37 are correlated due to summation over the same set of variables , i.e. @xmath8 s , @xmath24 is a ratio of two correlated gaussian variables . in order to obtain statistical properties of @xmath24
, we first calculate correlation coefficient @xmath38 of @xmath32 and @xmath37 , which is given by @xmath39/(\sigma_a \sigma_b )        = \operatorname{\mathbb{e}}\left [ \sum_{i=1}^{n } \sum_{j=1}^{n } x_i s_i s_j \right]/(\sigma_a \sigma_b ) .
\label{correlation - coeff1}\ ] ] note that    = _ i _ s^2 , & @xmath40 , + 0 , & @xmath41 , [ correlation ]    thus the correlation coefficient is obtained as @xmath42 where @xmath43^t$ ] , and @xmath44 is the @xmath45th norm of the vector .    from the previous analysis we know that @xmath24 is the ratio of two dependent gaussian random variables with non - zero means for which the closed - form pdf is available in ( * ? ? ?
* eqn 7.14 ) . based on the pdf of @xmath24
, we calculate its mean @xmath46 and variance @xmath47 by numerical integration .
alternatively , we also calculate @xmath46 and @xmath47 through approximation results of @xcite , which assume @xmath24 follows normal distribution .
the approximation results are given by    rcl m__p & & ( m_a / m_b ) + _
b^2 m_a / m_b^3 - _ ab _ a _ b /m_b^2 + _ _ p^2 & & _ b^2 m_a^2/m_b^4 + _ a^2/m_b^2 - 2 _ ab _ a
_ b m_a /m_b^3 .
[ approx ]    the gaussian approximation ( [ approx ] ) is verified via numerical simulations .
specifically , the mean and variance calculated from integration of the closed - form pdf are compared with the approximated mean and variance , for different node numbers .
our results show that for @xmath48 , the approximation error for the means is essentially zero and less than 3% for the variances .    up to this stage , we are able to calculate the mean and variance of the estimation error in the x - axis as @xmath49 , @xmath50 .
the same analysis applies for calculating the mean and variance of the estimation error along the y - axis , namely , @xmath51 and @xmath52 .
note that the estimation errors in x - axis and y - axis are correlated , therefore the statistical properties of the 2-dimensional error , @xmath53 , can not be obtained directly .      in this case , since @xmath8 s in ( [ est x ] )
are correlated , thus the distribution of @xmath24 needs to be evaluated again . since @xmath37 in ( [ est x ] ) is the linear combination of correlated gaussian variables , it is still gaussian distributed @xcite , i.e. @xmath54 .
thus , @xmath55 , and @xmath56 is given by @xmath57 , where @xmath58 is the @xmath59th element of the covariance matrix of @xmath60^{t}$ ] , defined as @xmath61                   = \operatorname{\mathbb{e}}[s_i s_j ] = \sigma_s^2 \lambda_{ij } , \label{covariance_b}\ ] ] where @xmath62 .
for the distribution of @xmath32 , it follows directly that @xmath63 .
its variance is given by @xmath64 , where @xmath65 is the @xmath59th element of the covariance matrix of @xmath66^{t}$ ] , which can be represented by @xmath67[(p_j - p_{min})x_j- x_j \mu_j]\ } \ieeenonumber \\
= \operatorname{\mathbb{e}}[x_i s_i x_j s_j ] , \label{covariance_a}\ ] ] where @xmath68 $ ] is given by     = _ s^2 ( _ i^2 + _
l^2 ) , & @xmath40 , + _ i _ j _ s^2 _ ij , & @xmath41 .
[ correlation ]    based on central limit theorem , we approximate @xmath32 as a gaussian random variable , i.e. @xmath69 . the approximation is validated via simulations , and it applies for @xmath70 .    the correlation coefficient @xmath38 is given by @xmath71/(\sigma_a \sigma_b )                                 = ( \sigma_s^2 \sum_{i=1}^{n }                                     \sum_{j=1}^{n } \overline{x}_i \lambda_{ij } ) / ( \sigma_a \sigma_b ) .
\label{correlation coeff}\ ] ] substituting the expressions for @xmath72 and @xmath73 into ( [ correlation coeff ] ) , after some simple manipulation , @xmath74 where @xmath75 and @xmath76 is a @xmath77 vector of all @xmath78 s .
since the location estimate @xmath24 is formed by @xmath79 , we use the closed - form pdf given by ( * ? ? ?
* eqn 7.14 ) to calculate @xmath46 and @xmath47 by numerical integration .
gaussian approximation can also be used similarly to i.i.d shadowing case in sec .
[ subsubsec : iidshadowing ] .      in our previous analysis
, we obtained the statistical properties of 1-dimensional location estimation error for both i.i.d . and correlated shadowing case . from gaussian approximation of 1-dimensional localization error ( [ approx ] ) , we know that @xmath80 , where @xmath81^t$ ] , @xmath82,\ ] ] and @xmath83 is the correlation coefficient between 1-dimensional errors , the calculation of which is presented in the appendix .
the calculation of localization error @xmath53 involves computing 2-norm of @xmath84 . since the two variables in @xmath84 are dependent gaussian r.v.s with different variances , we perform the following transformations .
first we form two independent gaussian variables from @xmath84 , by the standard technique of decorrelation @xcite , @xmath85^t$ ] , where @xmath86 is the orthogonal eigenvalue decomposition matrix of @xmath87 , which satisfies @xmath88 , where @xmath89 is a diagonal matrix .
therefore , one can verify that ( * ? ? ?
* eqn 2.1 - 159 ) @xmath90 . then by the orthogonal - invariant property of 2-norm @xcite , we know @xmath91 .    in order to obtain the pdf of @xmath92 ,
let us first consider pdf of @xmath93 .
@xmath94 and @xmath95 are two chi - square variables with one degree of freedom .
then the pdf of @xmath96 is the convolution of these two chi - square pdfs , denoted as @xmath97 .
the pdf of the localization error @xmath98 is given by @xmath99 .
as a result , the closed - form pdf of @xmath53 , for given sensor node positions , is given by    rcl p_e_l(x ) & = & ( - ) + & & _ i=0^ _ l=0^ , x 0 [ pdf_2d ]    where @xmath100 , @xmath101 and @xmath102 , @xmath103 are means and variances of @xmath104 and @xmath105 , respectively , and @xmath106 and @xmath107 are gamma function and @xmath108th - order modified bessel function of the first kind , respectively . the mean and variance of @xmath53 , denoted as @xmath109 and @xmath110 , are calculated using their definitions , by numerical integration over the above pdf .
note that the above calculations are based on a particular fixed node placement .
if the nodes are randomly placed within the area according to some distribution , say uniform distribution , integrating @xmath111 times over all node coordinates gives the theoretical average performance .
based on the theoretical analysis , wcl requires large number of nodes in order to achieve a high level of accuracy .
this causes large communication overhead in terms of number of messages and transmit power consumption when implemented in a centralized manner .
we now present a practical distributed implementation of wcl to address this problem in the context of cr networks .
clustering is commonly used in implementation of distributed algorithms . the challenge in designing these algorithms is to reduce the amount of information exchanged among clusters .
the proposed algorithm accomplishes this by splitting the task into two phases .
first , the cluster with the highest average rss is selected .
this cluster then acts as the coordinator for the second phase which performs the actual wcl algorithm .
nodes , @xmath112 , in the cr network are assumed to be organized into clusters , @xmath113 , using the clustering algorithm presented in @xcite .
this particular algorithm is based on geographic boundaries and generates clusters which have a predefined shape ( hexagonal in this case ) .
the regularity in the shape of the cluster is beneficial for the first step of the distributed algorithm since wcl works best for a radially symmetric distribution of nodes .
however , the proposed algorithm would work for any arbitrary clustering algorithm as long as nodes are clustered based on proximity .
the pu is assumed to be randomly placed in the entire area and can thus be in any of these clusters .
each cluster selects one node , @xmath114 , to be the cluster head .
the cluster head maintains a table of all adjacent clusters , @xmath115 , and the location of its member nodes , @xmath116 .
once a pu is detected by a node , using an existing pu detection method such as energy detection @xcite , it transmits its rss measurement to the corresponding cluster head which then participates in the distributed wcl algorithm outlined in the next section . using a predetermined fusion rule , such as a majority vote ,
the cluster head makes a decision on the presence of the pu .
all clusters that decide that the pu is present , participate in the distributed wcl algorithm detailed in the next subsection .
these clusters form a set defined as @xmath117 . if in addition
, some prior knowledge of the pu transmit power is available , a situation typical of a cr network coexisting with a licensed incumbent network , then a predetermined threshold can be used to further reduce the number of clusters in @xmath117 .
the threshold can be chosen based on the average node spacing to get a certain desired average number of nodes .
we investigate the effect of this number on wcl accuracy in sec .
[ subsec : variableparticipation ] .      the first phase , which selects the head cluster , is shown in algorithm [ alg : selection ] .
the objective is to select the cluster with the highest average rss , @xmath118 .
this is done in a distributed manner by comparing the average rss between neighboring clusters . however , directly implementing this procedure would require a large amount of inter - cluster communication .    to solve this issue , a cluster only communicates its average rss to the neighboring cluster , @xmath119 , in the direction of the pu relative to the cluster s centroid .    although techniques for estimating the angle of arrival of a transmitter could be beneficial , these methods require cluster heads to have multiple antennas which increases complexity of sensors . in order to get an approximation of the pu direction
, we propose a metric , @xmath120 , shown in line 5 of algorithm [ alg : selection ] .
we get @xmath120 by subtracting the wcl result of all nodes within the same cluster from the geometric centroid of those nodes .
this vector approximates the rss gradient within the cluster .
it corresponds to an approximate direction of the pu since on average the rss in ( [ received power ] ) decreases as distance to pu increases .
therefore , the negative of the gradient is directed towards the point with maximum rss .
once the cluster with the highest average rss is selected , algorithm [ alg : mwcl ] is performed . in this step
, the final wcl calculation is done using nodes within a given radius , @xmath121 , of the selected cluster s strongest node ( sn ) , @xmath122 .
@xmath121 is adaptively chosen to be the minimum of the distance to the closest map edge , @xmath123 , and the cluster radius @xmath124 .
this reduces the border effect , which is the tendency of the location estimate to be biased towards the center when the pu is close to the map border , by ensuring a radially symmetric node distribution even when the pu is close to the network edge .
in order to improve the localization accuracy , additional nodes from adjacent clusters are also included in the wcl calculation . using the location of @xmath122 and @xmath121
, these adjacent clusters can reduce the amount of information exchanged by only transmitting the location and rss measurements of nodes that fall within this circle . in the case of large clusters , each cluster could also choose to calculate a weighted centroid using its own nodes and transmit the result of this calculation only . since wcl is essentially an averaging algorithm , this technique allows even algorithm [ alg : mwcl ] to be distributed .
@xmath125 @xmath126 @xmath127 \mathbf{l}(n_j ) } { \sum_{n_j\in c_i}(\operatorname{p}(n_j ) - \min(\operatorname{p}(n_j))}$ ] @xmath128 @xmath129 select @xmath130 proceed to algorithm [ alg : mwcl ]    @xmath131 @xmath132 @xmath133 @xmath134 \mathbf{l}(n_j ) } { \sum_{n_j \in c^*_{wcl}}\left[\operatorname{p}(n_j ) - \min(\operatorname{p}(n_j)\right]}$ ]    * function *  @xmath135 * return * @xmath136 ,   @xmath137 * end function *     the power efficiency and computational complexity of this distributed approach could further be improved by taking into account a practical mobility model of the pu .
there has been a lot of work such as @xcite , which discussed how to incorporate past observations of the target position and a mobility model to improve localization .
however , applying these methods introduces complexity to wcl and defeats the purpose of such a low complexity scheme .    intuitively , assuming the pu continuously transmits during successive wcl estimations
, a pu would most likely be in the same cluster as it was in the previous wcl calculation , or it could move to an adjacent cluster . by considering this
fact , algorithm [ alg : selection ] can further limit the clusters involved in choosing the head cluster , which is the most power consuming phase , by only taking the previous cluster head and its adjacent clusters as @xmath117 .
we analyze the total transmit power and computational complexity of cwcl and the proposed dwcl algorithm . in order to do so
, we first need to investigate the transmit power on a single link of length @xmath138
. the minimum transmit power is obtained from ( [ received power ] ) as @xmath139 , where @xmath140 is the minimum received power to correctly receive the message . for the computation complexity estimates , we assume the number of operations ( ops ) for addition , substraction and comparison is 1 ; for multiplication and division is 10 .      in cwcl ,
the total number of messages is the sum of reporting messages from nodes to the fusion center , thus , @xmath0 . for the transmit power consumption , we assume the fusion center is roughly at the center of a circle area of radius @xmath1 , and nodes are placed according to a uniform distribution . in this case
, one can verify that , @xmath141 , the distance from the @xmath3th node to the fusion center satisfies @xmath142 , where @xmath143 is the uniform distribution from @xmath32 to @xmath37 .
thus the average total power consumption is expressed as @xmath144 = \operatorname{\mathbb{e } } [ p_{r,\min } \sum_{i=1}^{n } ( d_i^{\gamma}10^{-s_i/10 } ) ] \approx n p_{r,\min } \operatorname{\mathbb{e}}[d_i^{\gamma } ] \operatorname{\mathbb{e}}[10^{-s_i/10}],\ ] ] where we use statistical average to replace sum of realizations .    note that @xmath145 = \operatorname{\mathbb{e}}[{(d_i^2)}^{\gamma/2}]$ ] , which is the @xmath146th moment of an uniform variable . in order to calculate @xmath147 $ ] , we need the distribution of @xmath148 , which is a function of central gaussian variable @xmath149 . thus we can derive the desired pdf as @xcite @xmath150    finally , the average total power consumption is given by @xmath151 = n p_{r,\min } \operatorname{\mathbb{e}}[d_i^{\gamma } ] \int_0^{\infty } x p(x ) d x.\ ] ]    from the definition of wcl estimate , number of ops of cwcl is @xmath152 .
assume that : @xmath153 is the average number of nodes in one cluster ; @xmath154 is the average number of clusters ; @xmath155 is the number of neighboring clusters for each cluster ; @xmath156 is the average percentage of clusters who decide they are the head cluster in step 7 of algorithm 1 ( clusters that have to perform step 8) .
the total number of messages for algorithm 1 is @xmath157 . to calculate transmit power , in step 1 )
we use the same procedure as cwcl ( with a smaller radius ) ; step 2 ) and 3 ) are communication among cluster heads , we use the average distance between centers of two clusters as the estimate of the average distance among cluster heads . in algorithm 2
, the total number of messages is @xmath158 and the calculation of power consumption for this process is the same as step 2 ) of algorithm 1 .    in terms of computational complexity , the number of ops for algorithm 1 is @xmath159 while the number of ops for algorithm 2 is @xmath160 .
in this section , we evaluate the performance of cwcl algorithm and its proposed distributed implementation . for cwcl algorithm , using the theoretical framework presented in section [ sec : analysis ] , we consider impact of sensor node positioning error , different levels of shadowing @xmath161 and different correlation distances @xmath162 .
the effect of node placement is investigated under random grid and uniform distribution .
prior work on wcl uses mean localization error as the performance metric of interest . in @xcite , a slightly modified performance metric where the mean error is normalized by the transmission range @xmath1 of pu is adopted .
simulations in these prior work have shown that the accuracy of wcl improves as the density of nodes is increased .
this result is intuitive since the node spacing decreases as the density of nodes is increased . in this work ,
we introduce a novel performance metric where the mean localization error @xmath163 is normalized with respect to the average node spacing @xmath164 . the motivation for choosing this metric is to eliminate the effect of geometry scaling which intuitively improves accuracy .        in order to isolate the effect of shadowing on wcl performance we consider a fixed grid placement of nodes with perfect positioning ( i.e. @xmath165 ) and fixed pu position . the normalized mean localization error for this scenario is shown in fig .  [ fig : fixgridsquare_shadow_mean ] . for given shadowing variance , as the number of nodes increases , there is a slight improvement in the accuracy , which is a consequence of averaged shadowing over many nodes .
note that this gain saturates at approximately 200 nodes or more .
we also observe that wcl is quite robust to shadowing , since for an additional @xmath166 of shadowing from @xmath167 to @xmath168 , normalized error increases only 5% .
therefore wcl is well suited for pu localization in severely shadowed environments .
the results for mean localization error under correlated shadowing are shown in fig .
[ fig : fixgridsquare_corr_mean ] . assuming fixed value of @xmath169 , and different levels of @xmath162 normalized to @xmath164
are presented . in
highly correlated shadowing such as @xmath170 , increasing the number of nodes always worsens the error . for the high node density case , e.g. 400 nodes , the average error is as high as 47% of @xmath164 .
we conclude that when all nodes are included in the wcl calculation , correlation has a degrading effect to accuracy .
further investigation shows that the correlation distance resulting in the worst performance is proportional to the number of nodes involved .      the impact of positioning accuracy of sensing nodes is presented in fig .
[ fig : imperfect location knowledge ] with a fixed grid placement of nodes in an i.i.d . shadowing environment with @xmath171 .
analysis shows that errors in sensor node positioning decrease the wcl accuracy but the impact is marginal .
for example , @xmath7 of @xmath172 m increases the error for about @xmath78% .    ) as a function of variance in the participating nodes positioning error , @xmath7 . for this scenario ,
nodes are placed on a grid with equal spacing ( fixed grid placement ) in an area with @xmath173 m .
pu is at the center of the area . ]
the previous subsections considered fixed node placement and pu location . realistically
, node placement can be random and the pu also appears at random locations relative to the sensor nodes .
next we consider two cases of random node placement : random grid and uniformly distributed placement .    in the random grid scenario , the sensor nodes are still on a grid , but the pu can appear uniformly within the center box .
a grid distribution of nodes can be practically implemented if nodes are explicitly deployed for the purpose of acting as a dedicated localization sensor network .
evenly placed beacons have been shown to provide better accuracy for wcl in wsns @xcite .
results for random grid placement in fig .
[ fig : randomgridanduniformrandomplacement ] show approximately 10% performance loss compared to fixed pu location scenarios .
the uniformly distributed placement has been used in @xcite to more accurately model scenarios where the sensor node positions are not predefined , such as in mobile ad - hoc networks .
the results for uniformly distributed placement are also shown in fig .
[ fig : randomgridanduniformrandomplacement ] .
it is observed that randomness in the node placement can increase the error as much as three times compared to the fixed node placement .    ) with varying correlation distance normalized to the node spacing .
results of two placements , random grid ( rg ) and uniform random ( ur ) , are shown in an area with @xmath173 m .
for both placements , @xmath174 .
]      the results given so far assumed an ideal circular coverage area for the pu .
works such as @xcite indicated that in practical scenarios , the circular model can not be guarranteed due to several non - idealities such as fading , shadowing , and interference in the system .
intuitively , the degree - of - irregularity ( _ doi _ ) parameter determines how closely the coverage area can be approximated by a circle . _
= 0 means the transmission range is perfectly circular .
the model used in our simulation is the _
model briefly described in @xcite where the transmission range is a correlated gaussian r.v .
centered at @xmath175 with variance @xmath176 .    in fig .
[ fig : varydoi ] the effect of _ doi _ is investigated in the presence of uncorrelated medium level shadowing ( @xmath177 ) .
the result shows that _ doi _ between 0% and 40% increases the error by 4% .
_ doi _ is not a significant impairment for wcl because nodes closest to the pu , where the effect of irregularity is negligible , are weighted more than nodes closer to the edge of the coverage area . however , a slight increase in the normalized mean error can be observed as more nodes are used at high doi ( 30% or 40% ) .
the denser the node distribution , the more nodes are included in wcl , but the more nodes there are , the more deviation from a perfect circular coverage coverage area , leading to lower accuracy .     with varying degree - of - irregularity ( doi ) . for this scenario ,
nodes are uniformly and independently distributed ( uniform random placement ) in an area with average @xmath175 . ]      from fig .
[ fig : randomgridanduniformrandomplacement ] we observe that in correlated shadowing increasing the number of nodes also increases the error .
one solution to this problem is to reduce the number of participating nodes used in wcl for this scenario .
this could be achieved by ignoring nodes farther away from the pu to reduce the number of nodes resulting in less degradation .    in fig .
[ fig : varyingparticipation ] we investigate varying levels of node participation for a fixed average node spacing , @xmath164 . in this case
the total number of nodes is 100 , and the density remains the same while the number of involving nodes is varied .
to select which nodes are included in the wcl calculation , nodes are sorted according to their rss values .
the nodes with high received power are included in wcl . as a result , increasing the participation for the uncorrelated case always improves the accuracy .
thus , 100% participation is optimal .
however , even with a little correlation , this number is greatly reduced and about 10% of the total number of nodes is optimal .
another selection scheme was previously proposed in @xcite where only nodes whose rss is within 15% of the stronges node s rss are used .
the localization accuracy of wcl is also compared to other techniques suitable for cr scenarios .
this includes lateration , centroid and strongest node ( sn ) .
lateration is a range - based technique that uses least - squares estimation using distance estimates to the target . the centroid algorithm , on the other hand
, is a range - free technique which calculates the average of all node positions in range of the pu .
finally , the sn technique simply selects the node with the highest rss and uses its location as an estimate of the pu position .
results in fig .
[ fig : comaprison with other techiques ] show that wcl performs better under higher @xmath161 compared to lateration .
although this seems counter - intuitive at first , the result can be explained by taking note that range - based techniques require an accurate estimate of the distance to the target to provide accurate results . in the case of non - interactive localization ,
where the only available information is the rss from the pu , shadowing induces a very large error in the range estimates . under low shadowing ( @xmath178 db ) , the lateration method is more accurate than wcl . at 0 db
shadowing ( not shown in the figure ) , the range - based technique will achieve perfect localization .
this result supports the claim that wcl is a viable method in the context of non - interactive localization because of its robustness .
the centroid algorithm performs slightly worse than wcl because it does not utilize rss information .
however , the performance difference becomes smaller at higher level of shadowing .    .
for this scenario , nodes are uniformly and independently distributed ( uniform random placement ) in an area with @xmath175 ( uncorrelated shadowing ) .
note that the centroid technique is unaffected by @xmath179 . ]
[ fig : dwcl compare ] presents simulation results that compare normalized mean error of cwcl , sn and dwcl under different shadowing environments .
sn is included in the comparison because it forms the basis of the head cluster selection for dwcl so it acts as an upperbound to the error .
we observe that the proposed dwcl performance is between cwcl and sn algorithms .
it is also interesting to note that the accuracy of dwcl degrades relative to cwcl and approaches the performance of sn when shadowing is high .
the reason behind this is that for cluster head selection , dwcl is dependent on the sn algorithm .
the accuracy of the sn degrades at large shadowing shown in fig .
[ fig : comaprison with other techiques ] .
[ fig : fixgrid_nodepower_theo ] compares the transmit power consumption for cwcl and dwcl with different number of clusters .
in the simulation we set @xmath180 which is the mean received power at the border of the area of @xmath175 , assuming @xmath181 dbm , @xmath182 m and @xmath183 in ( [ received power ] ) . we observe that dwcl with 16 clusters consumes about @xmath184dbm less power for each node compared with cwcl .
the analysis also shows that transmit power of each node is reduced as the number of node increases , since the average communication distance is getting smaller in dwcl .
[ fig : fixgrid_cpu_theo ] compares computational complexity estimated in terms of number of ops .
it shows that distributed algorithms require about three times more operations than cwcl , however total transmit power is reduced .
in this paper we present the first analytical framework for the performance of wcl for localization of pu transmitters in cr networks .
we derived the closed - form expression for the error distribution of wcl parameterized by node density , node placement , shadowing variance , correlation distance and errors in sensor nodes positions . using the theoretical results in conjunction with numerical simulations ,
the robustness of wcl against various physical conditions is investigated and quantified .
based on the above analysis , we provide guidelines in terms of number of nodes and node placement required to achieve certain localization accuracy using wcl and propose variable node participation for improved wcl performance . to facilitate the deployment of wcl in an energy efficient context
, a distributed cluster - based implementation of wcl is also proposed , which does not assume a fusion center that collects all rss information from nodes .
the distributed method achieves comparable accuracy with its centralized counterpart , and greatly reduces total power consumption .        for notational convenience
we define @xmath186 and form the vector @xmath187^{t}$ ] , so that @xmath188 with @xmath189 .
note that this definition includes both i.i.d and correlated shadowing case .
the first term in ( [ dimensional correlation ] ) then becomes @xmath190 = \operatorname{\mathbb{e}}\left[\frac{\mathbf{x}^{t } \mathbf{q } \mathbf{q}^{t } \mathbf{y}}{\mathbf{1}^{t } \mathbf{q } \mathbf{q}^{t }     \mathbf{1}}\right].\ ] ]    define three new variables @xmath191 , @xmath192 and @xmath193 , and their vector form @xmath194^{t}$ ] .
now the first term in ( [ dimensional correlation ] ) becomes @xmath195 $ ] .
as shown in sec .
ii , @xmath196 , and @xmath197 .
@xmath198 can be approximated as a single gaussian variable as shown in sec .
[ subsubsec : iidshadowing ] since it is the same as @xmath32 in ( [ est x ] ) .
we denote its distribution as @xmath199 , similarly @xmath200 .
the distribution of @xmath201 is given as @xmath202 , where @xmath203 and @xmath204 .
as a result @xmath205 , where @xmath206^t$ ] is the mean vector and @xmath207 is the covariance matrix . to derive @xmath207 we find @xmath208 as @xmath209 - m_1 m_2                           = \operatorname{\mathbb{e}}[(\sum_{i=1}^{n } x_i q_i)(\sum_{j=1}^{n } y_j q_j ) ] - m_1 m_2 = \overline{\textbf{x}}^{t } \mathbf{\omega } \overline{\textbf{y } } - m_1 m_2\ ] ] where in the last the equation we used the fact @xmath210 = \overline{x}_i \overline{y}_j \mathbf{\omega}_{ij}$ ] .
similarly we calculate all entries of @xmath207 .    in order to simplify the calculation of ( [ reform ] ) , we perform the following transformation to @xmath211 , @xmath212 , thus @xmath213
. clearly @xmath214 , where @xmath215 is the @xmath3th column of @xmath216
. then we have @xmath217 = \operatorname{\mathbb{e}}[\frac{(m_1 + \boldsymbol{\omega}_{1}^{t } \widetilde{\textbf{t}})(m_2 + \boldsymbol{\omega}_{2}^{t } \widetilde{\textbf{t}})}{(m_3 + \boldsymbol{\omega}_{3}^{t } \widetilde{\textbf{t}})^2}].\ ] ]    since the expectation in ( [ results ] ) involves multiplication and division of three correlated random variables , the calculation can not be done directly .
in order to reduce the number of variables involved , first we form @xmath218 $ ] , then perform qr factorization to reduce dimensionality , obtaining @xmath219 , in which @xmath86 is orthogonal and @xmath220 is upper triangular @xmath221.\ ] ]      note that since @xmath225 , denominator of ( [ results ] ) can further be expressed as @xmath226 \widetilde{\textbf{t}}\right)^{2 }      \doteq \left(m_3 + r_{11 } \widetilde{t}_{1}\right)^{2 } , \label{den2_new}\ ] ] where we use the fact that @xmath227 is the first column of @xmath220 ; and @xmath228^{t}$ ] . similarly ,
numerator of ( [ results ] ) can also be further expressed as      therefore , the expectation in ( [ reform ] ) becomes @xmath229     = \frac { r_{12 } r_{13 } } { r_{11}^2 } \operatorname{\mathbb{e}}\left[\frac{\widetilde{t}_{1}^2 + \left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right ) \widetilde{t}_{1 } + \left(\frac{m_1 m_2 + r_{22 } r_{23}}{r_{12 } r_{13}}\right)}{\left(\widetilde{t}_{1 } + \frac{m_3}{r_{11}}\right)^{2}}\right ] .
\label{exp1_new}\ ] ] after factorization ( [ exp1_new ] ) becomes @xmath230 = \frac { r_{12 } r_{13 } } { r_{11}^2 } \!\!\left(\!\!1 + \frac{\left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right)-2 \frac{m_3}{r_{11}}}{\frac{m_3}{r_{11 } } } + \frac{\left(\frac{m_1 m_2 + r_{22 } r_{23}}{r_{12 } r_{13}}\right)+\left({\frac{m_3}{r_{11}}}\right)^2-\left(\frac{m_1}{r_{12 } } + \frac{m_2}{r_{13}}\right ) \frac{m_3}{r_{11}}}{1+\left({\frac{m_3}{r_{11}}}\right)^2 } \right).\ ] ] substituting ( [ corr - coeff_new ] ) into ( [ dimensional correlation ] ) gives the result .
j.  wang , p.  urriza , y.  han , and d.  abri , `` performance analysis of weighted centroid algorithm for primary user localization in cognitive radio networks , '' in _ proc .
asilomar conference on signals , systems , and computers _ , pacific grove , ca , usa , nov . 710 , 2010 .
h.  chen , q.  shi , r.  tan , h.  poor , and k.  sezaki , `` mobile element assisted cooperative localization for wireless sensor networks with obstacles , '' _ ieee trans .
wireless commun .
_ , vol .  9 , no .  3 , pp . 956 963 , mar . 2010
s.  liu , y.  chen , w.  trappe , and l.  j. greenstein , `` non - interactive localization of cognitive radios based on dynamic signal strength mapping , '' in _ proc .
conf . on wireless on - demand network systems and services
_ , snowbird , ut , usa , feb . 24 , 2009 .
r.  behnke , j.  salzmann , r.  grossmann , d.  lieckfeldt , d.  timmermann , and k.  thurow , `` strategies to overcome border area effects of coarse grained localization , '' in _ proc .
6th workshop on positioning , navigation and communication _ ,
hannover , germany , eu , mar .
19 , 2009 .
q.  zhou , x.  li , and y.  xu , `` mean shift based collaborative localization with dynamically clustering for wireless sensor networks , '' in _ proc .
conf . on communications and mobile computing _ ,
vol .  2 ,
yunnan , china , jan . 68 , 2009 .
m.  k. simon , _ probability distributions involving gaussian random variables : a handbook for engineers , scientists and mathematicians_. 1em plus 0.5em minus 0.4emnew york , ny , usa : springer - verlag , 2006 .
r.  chandra , j.  padhye , a.  wolman , and b.  zill , `` a location - based management system for enterprise wireless lans , '' in _ proc .
symposium on networked systems design and implementation _ , cambridge , ma , usa , apr . 1113 , 2007 .