# Abstract

Experimental studies on enzyme evolution show that only a small fraction of all possible mutation trajectories are accessible to evolution. However, these experiments deal with individual enzymes and explore a tiny part of the fitness landscape. We report an exhaustive analysis of fitness landscapes constructed with an off-lattice model of protein folding where fitness is equated with robustness to misfolding. This model mimics the essential features of the interactions between amino acids, is consistent with the key paradigms of protein folding and reproduces the universal distribution of evolutionary rates among orthologous proteins. We introduce mean path divergence as a quantitative measure of the degree to which the starting and ending points determine the path of evolution in fitness landscapes. Global measures of landscape roughness are good predictors of path divergence in all studied landscapes: the mean path divergence is greater in smooth landscapes than in rough ones. The model-derived and experimental landscapes are significantly smoother than random landscapes and resemble additive landscapes perturbed with moderate amounts of noise; thus, these landscapes are substantially robust to mutation. The model landscapes show a deficit of suboptimal peaks even compared with noisy additive landscapes with similar overall roughness. We suggest that smoothness and the substantial deficit of peaks in the fitness landscapes of protein evolution are fundamental consequences of the physics of protein folding.

# Author Summary

Is evolution deterministic, hence predictable, or stochastic, that is unpredictable? What would happen if one could "replay the tape of evolution:" will the outcomes of evolution be completely different or evolution is so constrained that history will be repeated? Arguably, these questions are among the most intriguing and most difficult in evolutionary biology. In other words, the predictability of evolution depends on the fraction of the trajectories on fitness landscapes that are accessible for evolutionary exploration. Because direct experimental investigation of fitness landscapes is technically challenging, the available studies only explore a minuscule portion of the landscape for individual enzymes. We therefore sought to investigate the topography of fitness landscapes within the framework of a previously developed model of protein folding and evolution where fitness is equated with robustness to misfolding. We show that model-derived and experimental landscapes are significantly smoother than random landscapes and resemble moderately perturbed additive landscapes; thus, these landscapes are substantially robust to mutation. The model landscapes show a deficit of suboptimal peaks even compared with noisy additive landscapes with similar overall roughness. Thus, the smoothness and substantial deficit of peaks in fitness landscapes of protein evolution could be fundamental consequences of the physics of protein folding.

# Introduction

One of the most intriguing questions in evolutionary biology is: to what extent evolution is deterministic and to what extent it is stochastic and hence unpredictable? In other words, what happens if "the tape of evolution is replayed:" are we going to see completely different outcomes or the constraints are so strong that history will be repeated ? If evolution is envisaged as movement of a population across a fitness landscape, the question can be reworded more specifically: among the numerous trajectories connecting any two points on the landscape, what fraction is accessible to evolution? Until recently, these remained purely theoretical questions as experimental study of fitness landscapes in the actual sequence space was impractical, due both to the technical difficulty of producing and assaying numerous expressed sequence variants and to the more fundamental problem of defining an adequate quantitative measure of fitness. However, recent experimental studies of fitness landscapes could potentially shed light on the problem of evolutionary path predictability.

The most thoroughly characterized feature of empirical fitness landscapes is the structure near a peak. In experiments that examine the peak structure, a high fitness sequence is typically subjected to either random mutations or an exhaustive set of mutations at a small number of important sites. The resulting library of mutants is then assayed to measure a proxy of fitness . Significant sign epistasis (a situation in which the fitness effect of a particular mutation can be either positive or negative depending on the genetic context) near the peak has been typically observed whereas local deviations from the additive model have been found to be uncorrelated with the genetic context and derived from a nearly normal distribution . Because these studies characterize only a small region of the landscape, they cannot be used to address the question of path predictability.

Another broad class of experiments probes the evolutionary trajectories from low to high fitness. Usually, in such experiments, a random peptide is subjected to repeated rounds of random mutagenesis and purifying selection . During this process fitness grows with each generation and eventually stagnates at a suboptimal plateau. The characteristics of the fitness growth as well as the dependence of the plateau height on the library size can be used to classify landscapes . A quantitative comparison to the $NK$ model of random epistatic landscapes ($N$ is the number of sites in an evolving sequence and $K$ is the number of sites that affect the fitness contribution of a particular site through epistatic interactions) can even yield quantitative estimates of $N$ and $K$ . The directed evolution studies explore the evolutionarily accessible portion of the landscape and could in principle be used to shed light on the question of path predictability. However, the inaccessible regions of the landscape remain unexplored and the volume of data at this point is insufficient to obtain quantitative conclusions regarding path predictability.

A different type of landscapes has been explored in various microarray experiments where protein-DNA(RNA) binding affinity serves as the proxy for fitness . These experiments produce vast, densely sampled landscapes. A comparison with a sophisticated Landscape State Machine model of a correlated fitness landscapes yields estimates of the model parameters . The DNA binding landscapes, in principle, contain the information required for the analysis of path statistics, and could be a valuable resource for advancing the understanding of evolutionary path predictability.

Empirical studies that exhaustively sample a region of the fitness landscape allow one to actually assess the accessibility of the entire set of theoretically possible evolutionary trajectories in a particular (small) area of the fitness landscape. For example, all mutational paths between two states of an enzyme, e.g., the transition from an antibiotic-sensitive to an antibiotic resistant form of $\beta$-lactamase or the transition between different specificities of sesquiterpene synthase have been explored. The results of these experiments, which out of necessity explore only short mutational paths of $<10$ amino acid replacements, suggest that there is a substantial deterministic component to protein evolution: only a small fraction of the possible paths are accessible for evolution .

Theoretical analyses of available empirical fitness data reveal a tight link between genetic and molecular interactions which are responsible for the landscape ruggedness and ubiquitous sign epistasis . The emerging quantitative analysis of fitness landscapes can shed light on some of the most fundamental aspects of evolution but the interpretation of the currently available experimental results requires utmost caution as only a minuscule part of the sequence space can be explored, and that only for a few more or less arbitrarily selected experimental systems.

Here we focus on the question of the predictability of mutational paths which is intimately tied to the ruggedness/smoothness of the fitness landscapes. The study of random landscapes of low dimensionality revealed an intuitively plausible negative correlation between the roughness of a landscape and the availability of pathways of monotonic fitness . In the same study, Carneiro and Hartl showed that experimentally characterized landscapes are significantly smoother than their permuted counterparts and exhibit greater peak accessibility .

To gain insights into the structure of the fitness landscapes of protein evolution and in particular the accessibility of mutational paths we used a previously developed simple model of protein folding and evolution . The key assumption of this model, which is based on the concept of misfolding-driven evolution of proteins , is that the fitness of model proteins is determined solely by the number of misfolded copies that are produced before the required abundance of the correctly folded protein is reached. We have previously shown that this model accurately reproduces the shape of the universal distribution of the evolutionary rates among orthologous protein-coding genes along with the dependencies of the evolutionary rate on protein abundance and effective population size . These results appear to suggest that our folding model (described in detail the Methods section) is sufficiently rich to reproduce some of the salient aspects of evolution. The model is also simple enough to allow exhaustive exploration of the fitness landscapes, which prompted us to directly address the problem of evolutionary path predictability.

We build on the efforts of Carneiro and Hartl who examined the statistics of evolutionary trajectories. Although counting monotonic fitness paths reveals important features of the landscapes, we argue that reliable retrodiction of the evolutionary past is possible (i.e., evolution is quasi-deterministic) only when the available monotonic paths are similar to each other in a quantifiable way. We therefore propose a measure of path divergence to quantify the difference between the available monotonic paths. Our aims are to investigate the structure of the fitness landscapes of protein evolution and to elucidate the connection between the roughness of landscapes and the predictability of mutational trajectories. We analyze three classes of fitness landscapes: landscapes in which fitness is derived from the folding robustness of model polymers; additive random landscapes perturbed by noise; and experimental landscapes derived from the combinatorial mutation analysis of drug resistance and enzymatic activity. We show that all three classes of landscapes are markedly smoother than their randomly permuted counterparts and all exhibit a similar qualitative connection between roughness and path predictability. However, at the same level of path predictability, the folding landscapes have substantially fewer fitness peaks. Given that the statistical properties of the model landscapes can be directly traced to the constraints imposed by the energetics and kinetics of a folding heterpolymer, we hypothesize that the relative smoothness and the suppression of suboptimal peaks in fitness landscapes of protein evolution are fundamental consequences of protein folding physics.

# Results

## Quantitative characterization of fitness landscapes

Carneiro and Hartl compared small random landscapes to several empirical fitness landscapes using deviation from additivity as a measure of roughness . They found that empirical landscapes were significantly smoother than their random counterparts and that the degree of smoothness was correlated with the number of monotonic paths to the main summit. Deviation from additivity of a landscape is computed by fitting an additive model in which the fitness of each sequence is different from the peak fitness by the sum of contributions of the substitutions that differentiate it from the peak sequence. The negative fitness contributions of the substitutions to the peak fitness are adjusted to minimize the sum $S$ of squares of the differences between the actual fitnesses in the landscape and the fitnesses predicted by the additive model. Deviation from additivity is defined as $\sqrt{S/L}$, where $L$ is the number of points in the landscape.

Because roughness of a multidimensional landscape with variable degree connectivity is not an intuitive concept, we introduce three additional quantitative measures to probe alternative facets of the concept of roughness. First, local roughness is the root mean squared difference between the fitness of a point and its neighbors, averaged over the entire landscape. As defined, local roughness conflates the measures of roughness and "steepness." For example, a globally smooth landscape, in which fitness depends only on the distance from the peak, will have a non-zero local roughness. However, because there is a large number of directions that change the distance from the peak by one, the local roughness of a globally smooth landscape will be vanishingly small. In addition, our landscapes tend to be globally flat–so that the average decrease in fitness due to a single mutation step away from the main peak is much smaller than the local fitness variability–everywhere except a small region around the main peak (see Fig. ). Therefore, the landscape-average local roughness in our case is a true measure of the local fitness variability.

Second, the fraction of peaks is the number of points with no fitter neighbors divided by the total number of points in the landscape. A strictly additive landscape has a single peak whereas the peak fraction in landscapes derived from the folding model as well as the corresponding randomized landscapes depends on the method of landscape construction, alphabet size and sequence length.

Third, the roughness of a landscape can be assessed by identifying its tree component. The tree component is the set of all nodes with no more than one neighbor of higher fitness. Thus, the tree component includes peaks and plateaus. Monotonic fitness paths along the tree component form a single or several disjoint tree structures without loops. In the limit of high selection pressure, a mutational trajectory that finds itself on the tree component has a single path to the nearest peak or plateau, i.e. evolution on the tree component is completely deterministic. We use the mean distance to the tree component, i.e. the distance to the tree component averaged over the landscape, as a measure of roughness. In a fully additive landscape, only the peak sequence and its immediate neighbors belong to the tree component and therefore the mean distance to the tree component is a measure of the diameter of an additive landscape (which, for example, could be defined as the maximum pairwise distance between points on the landscape). Kauffman and Levin have shown that in a large class of correlated random landscapes, the mean distance to the tree component grows only logarithmically with the number of points in the landscape .

We utilize two quantitative measures of the predictability of evolutionary trajectories. First is fraction of monotonic paths to the main peak $F_m$ which is computed by counting the number $n_i$ of simple (without reverse substitutions or multiple substitutions at the same site) monotonic paths to the main peak from each point $i$ on the landscape, dividing it by the total number of simple paths $h_i!$ (where $h_i$ is the Hamming distance from point $i$ to the peak), and averaging over the landscape via $$\label{eq:7}
  F_m = \frac{1}{L}\sum_i \frac{n_i}{h_i!},$$ where $L$ is the number of points in the landscape and the sum excludes the main peak. The monotonic path fraction $F_m$ measures the scarcity of accessible evolutionary paths when selection is strong. When the monotonic path fraction is small, evolution is more constrained.

Second, the mean path divergence, is a fine-grained measure of evolutionary (un)predictability. We first define the divergence $d(p_1,p_2)$ of a pair of paths $p_1$ and $p_2$, as the average of the shortest Hamming distances from each point on one path to the other path. Suppose that we have a way of generating stochastic evolutionary paths. The outcome of a large number of evolutionary dynamics simulations is a collection of paths with their associated probabilities of occurrence. In general, the probability of occurrence of an evolutionary path is proportional to the product of fixation probabilities of its constituent mutation steps. Given a bundle of paths with the same starting and ending points, we define its mean path divergence to be $$\label{eq:1}
  \bar d = \sum_{p_1 \neq p_2}  d(p_1, p_2) \, O(p_1) \,
  O(p_2),$$ where $O(p)$ is the probability of occurrence of path $p$ in the ensemble. In other words, if two paths were drawn from the bundle at random with probabilities proportional to $O(p)$, their expected divergence would be $\bar d$. Alternatively, if we were to fix one path to be the most likely path in the bundle and to select the second path at random with probability proportional to $O(p)$, the divergence would be proportional to $\bar d$ as well.

In an additive landscape, the mutational trajectory is maximally ambiguous. As every substitution that brings the sequence closer to the peak increases fitness, substitutions can occur in any order and all shortest mutational trajectories to the peak–without reverse substitutions or multiple substitutions at the same site–are monotonic in fitness. In the strong selection limit of our model defined below, all monotonic trajectories have roughly the same probability of occurrence, so the mutational path cannot be predicted.

The mean path divergence is a better measure of the predictability of evolutionary trajectories than the number or fraction of accessible paths. Even when only a small fraction of paths are monotonic in fitness, these paths could potentially be quite different, perhaps randomly scattered over the landscape. In such a case, prediction of the evolutionary trajectory would be inaccurate despite the scarcity of accessible paths which will be reflected in a high value of path divergence.

Equation () introduces the mean path divergence of a bundle of paths with the same starting and ending points. The landscape-wide mean path divergence is measured by constructing representative path bundles with all possible \[start, peak\] pairs including suboptimal peaks as trajectory termination points. Path divergence is averaged over all bundles with the starting and ending points separated by the same Hamming distance. To construct the path bundles, we employed a low mutation rate model in which the attempted substitutions are either eliminated or fixed in the population before the next mutation attempt occurs.

We invoke the misfolding-cost hypothesis to assign a fitness to a sequence that folds with probability $P$ to a particular structure. To produce an abundance $A$ of correctly folded copies, an average of $A(1 - P)/P$ of misfolded copies are produced. The "fitness" of a sequence should be a monotonically decreasing function of the cost incurred by the misfolded proteins. Previously we showed that qualitative conclusions drawn from the average population dynamics on the fitness landscape did not depend on the precise functional relationship between the number of misfolded copies and fitness . We use simply the negative of the number of misfolded copies and assign a fitness $w = -A/P$, to a sequence whose probability of folding to the reference structure is $P$. Because the exact population dynamics model is not important, we use diploid population dynamics in the low mutation rate limit. Therefore, the probability of fixation of a mutant $j$ in the background of $i$ is given by $$\label{eq:2}
  \pi(i\rightarrow j) = \frac{1 - e^{-2(w_j - w_i)}}{1 - e^{-4
      N_e (w_j - w_i)}},$$ where $N_e$ is the effective population size which in all simulations was fixed at $N_e = 10,000$. The required abundance $A$ is a measure of the strength of selection. In the limit of large $A$, the probability of fixation of a beneficial mutation is unity whereas neutral and deleterious mutations are never fixed. In this limit, all uphill steps in the fitness landscape are equally likely and all monotonic uphill trajectories have equal evolutionary significance.

In the analysis that follows, we study the association between landscape roughness and path predictability for the folding landscapes and their randomized (also referred to as permuted or scrambled) versions. In the scrambled landscapes, the topology (i.e. connectivity) of the landscape is preserved but the fitness values are randomly shuffled. We also compare the roughness and path predictability characteristics of the model and the experimental landscapes for $\beta$-lactamase and sesquiterpene synthase to those for noisy additive landscapes with a continuously tunable amount of roughness.

## Evolutionary path predictability in fitness landscapes

### Deviation from additivity, local roughness, peak fraction, and monotonic paths

We first establish that the folding and the experimental landscapes are significantly different from their randomly permuted counterparts. The deviation from additivity of the folding landscapes is typically several standard deviations below the mean of their scrambled counterparts. Although the additivity hypothesis accounts for less than 40% of the fitness variability (computed by comparing the sum of the squares of the fitnesses in the landscape to the sum of the squares of the residuals of the additive fitness model fit) in all but one of the folding landscapes, the deviation from additivity of the permuted landscapes is substantially greater (Fig. A). The experimental landscapes follow the same pattern, in agreement with the earlier findings of Carneiro and Hartl . Furthermore, both in the folding and in the experimental landscapes, the fraction of monotonic paths to the main peak is several standard deviations greater than in the respective scrambled landscapes (Fig. B). An even more striking disparity exists between the fraction of peaks in the folding landscapes and their permuted versions: the folding landscapes contain at least an order of magnitude fewer peaks than their scrambled counterparts; the experimental landscapes resemble the folding landscapes more closely than their own randomized versions (Fig. C).

To further characterize the deviation of the folding and experimental landscapes from their permuted counterparts, each landscape metric was measured and the mean and standard deviation were computed among 100 randomly permuted landscapes. We then compute the Z-score (deviation from the mean measured in the units of the standard deviation) of the original non-permuted landscape compared to the ensemble of the permuted landscapes. This Z-score shows how much more correlated the original landscape is, as measured by the chosen characteristic, compared to its scrambled counterparts (Figure ). Notably, despite the considerable scatter of the Z-score values for the folding landscapes, they all showed extremely large difference (mean Z-score greater than 20 standard deviations) from the scrambled landscapes for all measures, with the sole exception of the monotonic path fraction (Figure ). The two experimental landscapes also significantly differed from the scrambled landscapes albeit less so than the folding landscapes, again with the exception of the monotonic path fraction in which case the two classes of landscapes had similar Z-scores (Figure ).

Aside from the significant correlation (Pearson $\rho = -0.68$) between peak fraction and mean distance to the tree component, there was little or no correlation between the four measures of landscape roughness (Fig. ). Roughness of landscapes of high and variable dimensionality is impossible to capture by a single quantity. Therefore, the different measures seam to reveal distinct aspects of landscape architecture. The strong negative correlation between the peak fraction and mean distance to the tree component is due to the fact that each peak spawns a distinct subset of the tree component. The higher the density of peaks on the landscape, the larger fraction of the landscape that is covered by the tree component. Therefore the average distance to the tree component declines with the increasing density of peaks.

### Path divergence

Starting from a random non-peak sequence in the landscape, we introduced random mutations and accepted or rejected them according to equation () until the trajectory arrived at a fitness peak. This procedure was repeated a large number of times, and path bundles were constructed for all pairs of starting and ending sequences. Then the mean path divergence was computed for each path bundle using equation () and averaged over all bundles for which starting and ending points were separated by the same Hamming distance. When selection is weak, all mutations which do not result in a sequence with zero folding probability are accepted. Thus, evolution is a random walk on the landscape and the statistical properties of evolutionary trajectories are fully determined by the topology of the landscape (i.e. the connectivity of each node). Conversely, in the strong selection limit, only mutations that increase fitness are fixed. The mean path divergence varies smoothly between the two limits (Fig. ) and saturates at high selection pressure. In our analysis, we focus on the strong selection limit plateau. In the weak selection limit, the diversity of trajectories stems solely from the number of neighbors of each point; by contrast, in the strong selection limit, the statistics of the monotonic trajectories depend on the roughness of the landscape. Thus, the weak selection limit probes only the topology of the landscape whereas the strong selection limit also exposes its topography which appears to be critical for assessing predictability of evolution.

### Predictors and correlates of path divergence and monotonic path fraction

All four measures of landscape roughness can serve as predictors of path divergence and monotonic path fraction to some degree (Fig. ), in agreement with the notion that each of these measures reflects salient properties of fitness landscapes. The properties of the folding and empirical landscapes are consistent with those of additive landscapes that were perturbed by a moderate amount of noise (see Methods for details). A striking exception is the dearth of peaks and monotonic paths in folding landscapes all other characteristics being similar. Deviation from additivity and fraction of peaks are negatively correlated with path divergence. This relationship captures the intuitive notion that in rough landscapes there are fewer accessible evolutionary paths than in smooth landscapes, and furthermore, in rough landscapes, even those paths that are accessible show the tendency to aggregate within small areas on the landscape. Indeed, in both the folding model-derived landscapes and the experimental landscapes, the mean path divergence for all Hamming distances between the starting and ending points was dramatically greater than in scrambled landscapes (Fig. ). Interpreting these findings in terms closer to biology, the fitness landscapes derived from the model as well as experimental landscapes show greater robustness to mutations than random landscapes: a random mutation in a model-derived or experimental fitness landscape is more likely than expected for random landscapes to have no deleterious effect, leading to another monotonic path to the peak. Consequently, evolution on the model-derived and experimental landscapes is less predictable (deterministic) than it would be on uncorrelated random landscapes.

In contrast to deviation from additivity, the mean distance to the tree component is positively correlated with path divergence. When the tree component comprises a large fraction of the landscape, the mean distance to the nearest tree branch is small. Consequently, the path divergence is reduced as the paths that reach the tree component do not deviate from each other from that point onward. By the same token, when the tree component is large, there are fewer monotonic paths.

The origin of the positive correlation between the local roughness and path divergence (Fig. ) is less obvious. Paradoxically, greater noise results in lower mean local roughness of noisy additive landscapes. The lowering of the overall mean fitness with noise and, more importantly, the flattening of the mean fitness dependence on the distance from the peak (Fig. ) appear to provide an explanation for this counter-intuitive result. Indeed we found that in noisy additive landscapes there is a characteristic fitness value of approximately 0.2 above which roughness increases with increasing noise and below which roughness declines with increasing noise. Given that roughly 75% of the points on the landscape have fitnesses below 0.2, the landscape-averaged local roughness declines with increasing noise amplitude.

# Discussion

Here we examined the fraction of monotonic paths and introduced mean path divergence as quantitative measures of the degree to which the starting and ending points determine the path of evolution on fitness landscapes. The lower the mean path divergence value, the more deterministic (and predictable) evolution is. Global measures of landscape roughness correlate with path divergence in the three analyzed classes of fitness landscapes: additive landscapes perturbed by noise, landscapes derived from our protein folding model and two small empirical landscapes. The folding landscapes are substantially smoother than their permuted counterparts. As a result, although in all analyzed landscapes only a small fraction of the theoretically possible evolutionary trajectories is accessible, this fraction is much greater in the folding and experimental landscapes than it is in randomized landscapes. In addition, the mean path divergence in the randomized landscapes is significantly smaller than in the original landscapes. Thus, the model and empirical landscapes possess similar global architectures with many more diverged monotonic paths to the high peaks than uncorrelated landscapes with the same distribution of fitness values. Consequently, evolution in fitness landscapes is substantially more robust to random mutations and less deterministic (less predictable) than expected by chance. These findings are compatible with the concept that might appear counter-intuitive but is buttressed by results of population genetic modeling, namely, that robustness of evolving biological systems promotes their evolvability . Additionally, the folding landscapes exhibit a substantial deficit of peaks compared to perturbed additive landscapes and experimental landscapes, a property that translates into a substantially greater fraction of paths leading to the main peak.

When it comes to the interpretation of the properties of fitness landscapes described here, an inevitable and important question is whether the folding model employed here is sufficiently complex and realistic to yield biologically relevant information. In selecting the complexity of our folding model, we attempted to construct the simplest model which exhibits 1) a rich spectrum of low energy conformations across the sequence space, and 2) a non-trivial distribution of substitutions effects on the low energy conformations. An important choice is whether the location of monomers is confined to a lattice or can be varied continuously. When the configuration space is continuous, the distribution of energy barriers between energetically optimal conformations can extend to zero. Therefore, the subtlety of distinctions between conformations can lead to a richer structure of the fitness landscape. We chose not increase the complexity of the model further and treated monomers as point-like particles in a chain where the distance between nearest neighbors is fixed but the angle between successive links in the chain in unrestricted. Our level of abstraction is therefore somewhere between lattice models and all-atom descriptions of proteins .

Another important choice is the number of the model monomer types. Again, we opted for an intermediate level of abstraction and chose four types of monomers: hydrophobic, hydrophilic, and charged. This choice drastically reduces the size of the sequence space while retaining some of the substitution complexity whereby hydrophilic and charged monomers can be swapped under some conditions without radically altering the native state. The intermediate level of abstraction in our approach has its pros and cons. Although the model reproduces key features of protein folding such as the existence of the hydrophobic folding nucleus and two-stage folding kinetics , compact conformations certainly do not represent proteins. Rather, we might think of our monomers as representing structurally grouped regions several (perhaps up to a dozen) amino-acids in length. Compact conformations in the model might therefore be analogous to tertiary structures of proteins. Representing sequence space with only four monomer types and treating mutations without reference to the underlying DNA or genetic code does not accurately reflect the natural mutation process. However, our goal was to isolate the features of fitness landscapes which could be traced directly to the constraints imposed by the heteropolymer folding kinetics and energetics. We therefore used a simple sequence space and a homogeneous mutation model to avoid compounding the fitness landscape structure by the complexity derived from the mutation process.

Most importantly, our folding model has been shown to reproduce the observed universal distribution of the evolutionary rates of protein-coding genes as well as the dependencies of the evolutionary rate on protein abundance and effective population sizes . Therefore, despite its simplicity, the behavior of this model might reflect important aspects of protein evolution. In particular, the conclusions drawn from the analysis of the model landscapes exhaustively explored here could also apply to the fitness landscapes of protein evolution. In the previous work, we concluded that the universal distribution of evolutionary rates and other features of protein evolution follow from the fundamental physics of protein folding . The results presented here suggest that the (relative) smoothness and a substantial deficit of peaks in the fitness landscapes of protein evolution that lead to mutational robustness and the ensuing evolvability could similarly follow from the fact that proteins are heteropolymers that have to fold in three dimensions to perform their functions.

The experimental landscapes considered here are decidedly incomplete. Due to experimental limitations, only the analysis of binary substitutions at a handful of sites is feasible at this time. The incompleteness of the empirical landscapes analyzed in this work could be the cause of the observed lack of peak suppression. This proposition will be put to test by the study of larger parts of experimental landscapes that are becoming increasingly available.

# Materials and Methods

## Folding model

The goal of this study is to explore the relationship between roughness and path divergence in realistic fitness landscapes. Our polymer folding model provides a simple way of constructing such landscapes. The model has been described in detail previously .

In brief, the model polymer is a flexible chain of monomers in which the the nearest neighbors interact via a stiff harmonic spring potential with rest length $a = 1$. The angles between the successive links in the chain are unrestricted. There are four types of monomers: hydrophobic <span class="sans-serif">H</span>, hydrophilic <span class="sans-serif">P</span>, and charged <span class="sans-serif">+</span> and <span class="sans-serif">–</span>. Next nearest neighbors $i$ and $j$ in the chain and beyond interact via a pairwise potential $$\label{eq:3}
  U_{ij}(r_{ij}) = \frac{A_{ij}}{r_{ij}^{12}} - \frac{C_{ij}}{r_{ij}^6} +
    \frac{q_i q_j e^{-D \, r_{ij}}}{r_{ij}},$$ where $r_{ij}$ is the distance between monomers $i$ and $j$, $q_i$ is the monomer's charge, $D$ is the Debye-Hückel screening length, and $A_{ij}$ and $C_{ij}$ depend on the pair in question. The interaction parameters are chosen to mimic the essential features of the amino-acid interactions. To emulate the effects of solvent, we assign a stronger attraction to the <span class="sans-serif">HH</span> pair than to the <span class="sans-serif">PP</span>, <span class="sans-serif">++</span>, and <span class="sans-serif">– –</span> pairs. There is also a long range repulsion between <span class="sans-serif">H</span> and <span class="sans-serif">P</span> and even stronger repulsion between <span class="sans-serif">H</span> and the charged monomers. The values of the parameters are $q_{\pm} = \pm 2$, Debye-Hückel screening length $D = 3$. The Lennard-Jones coefficients $A_{ij}$ and $C_{ij}$ are $$\begin{aligned}
  \label{eq:6}
  && A_{HH} = 4, \ A_{HP} = A_{H+} = 2, \ A_{PP} = A_{P+} = A_{++} =
  1, \nonumber \\
  && C_{HH} = 8, \ C_{HP} = -1, \ C_{H+} = -3, \ C_{PP} = C_{P+} =
  C_{++} = 2.
\end{aligned}$$ Note that a $+$ can be substituted by a $-$ in the subscripts and the coefficients are symmetric with respect to the interchange of the indices.

The energy of the chain is $$\label{eq:4}
  E = \sum_{|i - j| > 1} U_{ij} + \frac{b \, T}{2} \sum_{i = 1}^{N-1}
  (r_{i,i+1} - a)^2,$$ where the first term is the sum of the pairwise energies given by Eq. () over non-nearest neighbor pairs, and the second term reflects the springs connecting nearest neighbors. The spring constant is proportional to temperature $T$. The parameters are fixed for all simulation runs at $b = 300$, and the quench temperature $T =
1$. To mimic the observed tendency of the $N$ and $C$ termini to be in close proximity, we fixed the endpoint monomers of the model sequences to be of $+$ and $-$ types.

Dynamics of folding are simulated via overdamped Brownian kinetics which are appropriate when inertial and hydrodynamic effects are not important. Units are chosen so that each component $\alpha$ of the $i$'th monomer's coordinates $x_{\alpha i}$ is updated according to $$\label{eq:5}
  x_{\alpha i}(t + \Delta t) = x_{\alpha i}(t) - \frac{\Delta
    t}{T} \, \frac{\partial E}{\partial x_{\alpha i}}(t) +
  W_{\alpha i}(t),$$ where $\Delta t$ is the time step and $W_{\alpha i}(t)$ is a random variable with zero mean, variance $2 \Delta t$, uncorrelated with $W$ for other times, monomers and spatial directions.

## Native structure ensemble and correct folding probability

The "native structure" of a particular sequence is represented by an equilibrium ensemble of conformations. The ensemble is constructed by identifying the typical folded conformation and measuring the characteristic RMSD $D$ due to thermal fluctuations in the folded state. Three thousand quenches are then performed and the resulting folded conformations are accumulated. The equilibrium ensemble that represents the native structure is defined as the largest cluster of quenched conformations within RMSD distance $D$ from each other. Thus, each conformation in the ensemble differs from any other by an amount comparable to the differences introduced by thermal fluctuations alone.

The concept of the native structure ensemble allows us to compute the probability that a sequence folds to a particular structure in a natural, physically plausible fashion. Given a native structure ensemble we assess its conformation space density by computing the distance $d_i$ between each member $i$ of the ensemble and its closest neighbor. Given the set $\{d_i\}$ of these shortest distances we compute the median $Q$ and the median absolute deviation (MAD) $V$. A new conformation is deemed to belong to the ensemble if the shortest distance from this conformation to the members of the ensemble is smaller than $R = Q + 3V$.

Given a native structure ensemble of some sequence $s_1$ we compute the probability $P$ that sequence $s_2$ (which could be $s_1$ itself) folds to the this structure by accumulating $M = 100$ equilibrated quenched conformations of $s_2$ and using the above criterion to determine the fraction $P$ that belong to the native structure ensemble of $s_1$. Because $M = 100$ sample conformations are computed, the smallest measurable $P$ is $1/M = 0.01$. The sample size used to measure $P,$ dictated by the computational demands of the model, introduces a random component to the model fitness landscapes. As we report below, model landscapes turn out to be substantially smoother than random. Therefore the underlying global structure of the model landscapes appears to survive the modest amount of randomness introduced by the relatively small sample size used for measuring $P$.

## Search for compact robust folders

Robust folders (sequences with a high probability of correct folding) tend to have large linear regions stretched by repulsive Coulomb interactions. Because the linear regions have no contacts with other monomers, we focused our attention on compact conformations with a high monomer contact density. Substitutions in these higher complexity conformations were more likely to exhibit non-trivial effects. To find compact robust folders in the vast available sequence space of $23$-mers (the sequences are of length $N=25$ but the endpoint monomer types are fixed) with $4$ monomer types, we implemented a simulated annealing search which optimized the correct folding probability $P$ divided by the cube of the native conformation's radius of gyration. The search produced over 800 sequences with $P > 0.5$ and at least two distinct regions of the polymer in mutual contact.

## Assembly of the folding fitness landscapes

We examined each single substitution mutant of a robustly folding sequence and computed the folding probability $P$ to the structure of the original sequence. All mutants with $P > 0$ were added to the landscape and if $P\ge 0.1$ their mutants were also examined. This process is repeated until all mutants of the last sequence under consideration have $P < 0.1$.

From our study of complete landscapes we estimate that on average for each sequence with $P > 0$ which is included into the landscape, roughly 6 others with $P = 0$ need to be examined. Since each quench and equilibration takes about 2–4 seconds, landscape construction takes roughly 30 minutes to an hour per included sequence. Thus landscapes larger than 10,000 sequences take months to compile.

At the time of submission, 39 complete landscapes have been constructed, the largest comprising 12969 sequences.

## Additive landscapes perturbed by noise

The organization of the folding fitness landscapes and experimental landscapes were compared with perfectly additive landscapes perturbed by noise constructed as follows. Each substitution to the peak fitness sequence was assigned a negative fitness differential drawn at random from an exponential distribution with parameter $\lambda = 3$. The sum over the fitness differentials of a particular set of substitution was modified by either additive of multiplicative noise . Additive noise is drawn from a Gaussian distribution with zero mean and standard deviation $\nu$ which was varied between $0.05$ and $0.5$. The multiplicative perturbation is achieved by multiplying the fitness by a number drawn from a uniform distribution $[0,1)$ raised to a positive power $\mu$ varied between $0.1$ and $10.$ When $\mu$ is small, multiplicative factors are close to unity and the perturbation is small as well. If the perturbed fitness was positive, the mutant was included into the landscape. The noise amplitude was varied to obtain a family of landscapes of continuously varying roughness. Only the data for the additive landscapes with multiplicative noise were included in this manuscript. Landscapes perturbed by other types of noise exhibited essentially the same qualitative behavior.

## Experimental landscapes

The studies on experimental fitness landscapes typically involve constructing a library of all possible combinations of binary mutations at a small number of sites. The first study included in the present analysis measured the minimum inhibitory concentrations (MIC) of an antibiotic for a complete spectrum of mutants with modified TEM $\beta$-lactamases; the transition from the antibiotic-sensitive to the antibiotic-resistant form requires five mutation, so the landscape encompassed 120 mutational trajectories between the most distant points on the landscape (or 32 sequences) . The logarithm of MIC was used as the proxy for fitness. In the second study, catalytic activity of 419 sesquiterpene synthase mutants that differed by at most 9 substitutions was measured . We used the catalytic specificity (propensity for producing a particular reaction product rather than a broad spectrum of products) of the mutant enzymes as the proxy for fitness. Before performing the analysis, the fitnesses in the experimental landscapes are mapped onto the $[0.01,1)$ interval to enable meaningful quantitative comparisons of the roughness measures.

# Figure Legends


abstract: This paper presents a survey of energy efficiency of Medium Access Control (MAC) protocols for Wireless Body Area Sensor Networks (WBASNs). We highlight the features of MAC protocols along with their advantages and limitations in context of WBASNs. Comparison of Low Power Listening (LPL), Scheduled Contention and Time Division Multiple Access (TDMA) is also elaborated. MAC protocols with respect to different approaches and techniques which are used for energy minimization, traffic control mechanisms for collision avoidance are discussed.We also present a survey of path loss models for In-body, On-body and Off-body communications in WBASNs and analytically discuss that path loss is maximum in In-body communication because of low energy levels to take care of tissues and organs located inside the body. Survey of Power model for WBANs of CSMA/CA and beacon mode is also presented.
author: S. Hayat, N. Javaid, Z. A. Khan$^{\S}$, A. Shareef, A. Mahmood, S. H. Bouk
 Department of Electrical Engineering, COMSATS
Institute of Information Technology, Islamabad, Pakistan
$^{\S}$Faculty of Engineering, Dalhousie University, Halifax, Canada
title: Energy Efficient MAC Protocols

# Introduction

of wireless, medical and computer networking technology has merged into an emerging horizon of science and technology called Wireless Body Area Networks (WBANs). However, applications of WBANs are not limited to medical field only. Miniaturization and connectivity are notable parameters of this field. WBANs consist of three levels; first level is low power sensors or nodes which are battery powered and need to be operated for a long time without repairing and maintenance. These nodes may be placed on the body, around the body or implanted in the body. Second level is called master node, gateway or coordinator which controls its child nodes; its power requirements may be less strengthened than nodes due to its applications and flexibility. Third level is the local or metropolitan or internet network that serves for monitoring purposes.

Energy efficiency or effective power consumption of a system is one of the basic requirements for WBANs because of limited power of batteries. The most suitable layer for discussing energy and power issues is MAC Layer. The basic way of saving power or enhancing energy efficiency is to minimize the energy wastage. There are several sources of energy wastage including packet collisions, over hearing, idle listening, control packet overhead, etc. Major source of energy inefficiency among the above listed sources is packet collision for WBANs. Fig. 1. best explains that how a node's battery is consumed, in the process of communication.

Collision avoidance for energy efficiency, minimum latency, high throughput, and communication reliability, are basic requirements in the design of MAC protocol. The fundamental way of saving power or enhancing energy efficiency is to minimize the energy wastage. Simulations are performed in MATLAB for different scenarios to compute path loss. Results show that path loss is maximum in In-body communication, as compared to On-body and Off-body communication because human body is composed of tissues and organs in which communication is difficult and thus results in high path loss. On-body and Off-body also results some variations when the source and destination sensors or nodes are placed Line of Sight (LoS) and Non Line of Sight (NLoS).

In this paper, we therefore, provide a survey of energy efficient MAC protocols for WBANs. First, we elaborate the protocol features and then their advantages and limitations are discussed. Sources that contribute to the energy inefficiency in a particular protocol is also identified. Moreover, comparisons of MAC protocols in the context of WBANs are tabulated in detail.

# Related Work

Gopalan *et al.* \[1\] survey MAC protocols for WBANs along with the comparison of four protocols i.e., Energy Efficient MAC, MedMac, Low Duty Cycle MAC, and Body MAC. Some key requirements and sources of energy wastage are also discussed. They also discussed some open research issues in this survey. Still a lot of work has to be done in data link layer, network layer and cross layer design.

In \[2\], Shahjahan kutty *et al.* discuss the design challenges for MAC protocols for WBANs. They classify data traffic for WBANs into three categories: energy minimization techniques, frame structures and network architecture. However, the comparison of protocols is not provided by them.

Sana Ullah *et al.* in \[3\] provide relatively a comprehensive study of MAC protocols for WBANs. Comparison of the low power listening, scheduled contention and Time Division Multiple Access (TDMA) is provided. MAC requirements, frame structures and comparison of different protocols and their trade-offs are discussed in detail.

# Energy Minimization Techniques in MAC Protocols for WBANs

Low power mechanisms play an important role in performance enhancement of MAC protocol for WBANs. In this section, different approaches and techniques that provide energy efficiency in MAC protocols for WBANs are discussed and compared.

Energy efficiency is an important issue because the power of nodes in WBANs is limited and long duration of operation is expected. The key concept for low power consumption is to minimize the energy consumption in the following sources: sensing, data processing and communication.

Most of the energy wastage is caused during communication process because of the collision of packets, idle listening, over hearing, over-emitting, control packet overhead and traffic fluctuations. Idle listening can be reduced through duty cycling. To reduce energy waste in order to increase network's life time and to enhance the performance of MAC protocol, different wake-up mechanisms are used.

There are three main approaches adopted for the energy saving mechanisms in MAC protocols for WBANs, which are: Low Power Listening (LPL), Scheduled Contention, and TDMA.

## Low Power Listening

LPL procedure is that "node awakes for a very short period to check activity of channel". If the channel is not idle then the node remains in active state to receive data and other nodes go back to sleeping mode. This is also termed as channel polling \[3\]. This procedure is performed regularly without any synchronization among the nodes. A long preamble is used by the sender to check polling of the receiver. LPL is sensitive to traffic rates which results in degradation of performance in the scenario of highly varying traffic rates. However, it can be optimized effectively for already known periodic traffic rates. Wise-MAC \[3\] is one of the MAC protocols which is based on LPL. This protocol reduce Idle listening using non-persistent CSMA and preamble sampling technique.

## Scheduled Contention

Scheduled Contention is the combination of the scheduling and contention based mechanisms to effectively cope with the scalability and collision problems. In contention based protocols, contending nodes try to access the channel for data transmission therefore, ability of collision of packet is greatly increased. Example of contention based MAC protocol is Carrier Sense Multiple Access/Collision Avoidance (CSMA/CA) in which Clear Channel Assessment (CCA) is performed by the nodes before transmitting data.

Scheduling or contention free means that each node has the schedule of transmission in the form of bandwidth or time slot assignment. TDMA, CDMA and FDMA schemes are some examples of scheduling mechanisms. However, CDMA and FDMA are not suitable for WBANs because of high computational overhead and frequency limitations, respectively.

TDMA is the most suitable scheduling scheme, even though it requires extra power consumption due to its sensitivity for synchronization. The scheduled contention is the combination of scheduling and contention based mechanisms. In scheduled contention, a common schedule is adopted by all the nodes to transmit data. This schedule is exchanged periodically among the nodes to make communication adaptive, flexible and scalable.

Sensor MAC (S-MAC) is one of a MAC protocol based on the scheduled contention. In this protocol, low duty mode is set as default mode for all the nodes which assures the coordinated sleeping among neighboring nodes. The energy wastage due to collision, overhearing, idle listening etc. is minimized because the node is turned on only for transmission of data and remains in sleep mode, otherwise.

## Time Division Multiple Access

In TDMA mechanism, a super frame consists of a fixed number of time slots is used. Time slots are allocated to the sensor nodes by a central node and is known as Master Node (MN), Cluster Head (CH), coordinator or Base Station Transceiver (BST). Traffic rate is one of the key parameter used by the coordinator to allocate time for each contending node. The scheme is power efficient because a node gets time slot for transmission of data and remains in sleep mode for rest of the time. However, the synchronization requirements may degrade performance in terms of power consumption. Therefore, it is highly sensitive to clock drift, which may result in limited throughput. Preamble-Based TDMA (PB-TDMA) protocol is one of the TDMA based protocol. Other examples include Body-MAC (B-MAC) \[5\], MedMAC \[3\] etc.

These techniques are briefly compared in Table.I.

```latex
\begin{table*}[t]\begin{tabular} { | m {2cm}| m {5cm}| m {4cm}| m{5cm}| }
       \multicolumn{4}{c}{Table.1. Comparison between LPL, Schedule Contention, and TDMA} \\
    \hline
    \centering
   \textbf{Energy Saving Mechanisms} & \textbf{ LPL} & \textbf{Scheduled Contention} & \textbf{TDMA} \\ \hline
     Adaptability to traffic and delay & Scalable and adaptive to traffic load and low delay & Better delay performance due to sleep schedules & Better end-to-end reliability, smaller delays, high reliability \\ \hline
    Transmission latency and throughput & Flexible, high throughput, tolerable latency, and low power consumption  & High transmission latency, loosely synchronized, low throughput & Good for energy efficiency, prolonged network's lifetime, load balancing \\ \hline
    Synchronous/ Asynchronous & Asynchronous &  Synchronous  & Synchronous-Fine grained time synchronization  \\ \hline
   Traffic heterogeneity requirements &  Low duty cycle nodes do not accommodate aperiodic traffic. Very hard to satisfy the WBANs traffic heterogeneity requirements & Low duty cycle nodes do not require frequent synchronization of schedules. Hard to satisfy the WBANs traffic heterogeneity requirements & Low duty cycle nodes do not require frequent synchronization at the beginning of each superframe. Easy to satisfy the WBANs traffic heterogeneity requirements
    \\ \hline
   Sensitivity &  Sensitive to tuning for neighborhood size and traffic rate & Sensitive to clock drift & Very sensitive to clock drift \\ \hline
    Performance with respect to traffic rates&  Poor performance when traffic rates changes & With the increase in traffic, performance is improved & Throughput and number of active nodes are limited  \\ \hline
    Cost incurred by sender and receiver &  Receiver and polling efficiency is gained at much greater cost of senders & Similar cost incurred by sender and receiver & Require clustering \\ \hline
     Extravagant & It does not listen for full contention period as a result it is less expensive & Listening for full contention period  & Low duty cycle  \\ \hline
     Scalability and adaptability & Challenging to adapt LPL directly to new radios like IEEE 802.15.4 & Scalable, adaptive, and flexible & Limited scalability and adaptability to changes on number of nodes   \\ \hline

     \end{tabular}
     %\caption{COMPARISON BETWEEN LPL, SCHEDULE CONTENTION, AND TDMA}
   \end{table*}
```

# Energy Efficient MAC Protocols

In this section, we briefly discuss the energy efficient MAC protocols for WBAN.

## Okundu MAC Protocol

An energy efficient MAC protocol for single hop WBANs is proposed by Okundu *et al.* in \[4\]. This protocol consists of three main processes: link establishment, wakeup service, and alarm process. Basic energy saving mechanism of this protocol consists of central control of wakeup/sleep time and Wakeup Fall-back Time (WFT) processes. WFT mechanism is used to avoid collision due to continuous time slot. This mechanism states that, if a slave node wants to communicate with a MN and it fails in its task due to MN's other activities, then it goes back to sleep mode for a specific time computed by WFT. However, data is continuously being buffered during the sleep time.

To minimize time slot collision, the concept of WFT has been introduced. This concept helps every slave node to maintain a guaranteed time slot even if it fails to communicate with the MN. In this protocol, problems like idle listening and over-hearing can be reduced because of central management of traffic.

In one cluster, only $8$ slave nodes can be connected to MN which restricts inclusion of other slave nodes. In link establishment, wakeup service, and alarm processes, communication is initiated by the MN. Another main problem is that, only one slave node can join network at a time.

## MedMac Protocol

N. F. Timmons *et al.* in \[5\] propose a TDMA-Based MAC protocol for WBANs called MedMAC. This protocol consists of two schemes for the power saving: Adaptive Guard Band Algorithm (AGBA) and Drift Adjustment Factor (DAF). AGBA along with time stamp is used for synchronization among coordinator and other nodes. This synchronization is introduced using Guard Band (GB) between time slots to allow the node to sleep for many beacon periods. DAF is used to minimize bandwidth. GB is calculated by AGBA and shows the worst cases. However, practically gaps may be different between time slots depending upon application scenarios. DAF adjusts GB according to practical situation and avoids overlapping between consecutive slots.

MedMac outperforms IEEE 802.15.4 for Class 0 (lower data rate applications such as health monitoring and fitness) and Class 1 (medium data rate medical applications such as EEG). Energy waste due to collision is reduced by introducing Guaranteed Time Slot (GTS). Each device has exclusive use of a channel for a fixed time slot, therefore, synchronization overhead is also reduced.

This protocol works efficiently for low data rate applications, and work inefficiently for high data rate applications. However, In-body and On-body applications of WBAN are usually of higher data rate.

## Low Duty Cycle MAC Protocol

Low Duty Cycle MAC protocol for WBANs is designed in \[6\]. In this protocol, analog to digital conversion is performed by slave nodes while the other complex tasks such as digital signal processing is carried out at MN. MNs are supposed to be less power than slave nodes.

This protocol introduces the concept of Guard Time ($Tg$) to avoid overlapping between consecutive time slots. After T frames a Network Control (NC) packet is used for general network information. Power saving is achieved by using effective TDMA strategy.

This protocol is energy efficient because it sends data in short bursts. By using TDMA strategy, this protocol effectively overcomes the collision problem. It allows monitoring patient's condition and can reduce the work load on medical staff, while keeping minimum power usage.

TDMA strategy is used in WBANs, and it is found that TDMA is more suitable for static type of networks with a limited number of sensors generating data at a fixed rate therefore, this protocol may not respond well in a dynamic topology.

## B-MAC Protocol

B-MAC protocol achieves energy efficiency by using three bandwidth management schemes: Burst, Periodic and Adjust Bandwidth.

Burst bandwidth consists of temporary period of the bandwidth, which includes several MAC frames and recycled by the gateway (coordinator). Bandwidth is reduced to half if it does not fully utilized by the nodes, which is also informed about reduction of bandwidth. Periodic bandwidth is a provision for a node to have access to the channel exclusively within a portion of each MAC frame or few MAC frames. It is also allocated by the gateway based on node's QoS requirements and current availability of the bandwidth \[7\]. Adjust bandwidth defines the amount of bandwidth to be added to or reduced from previous Periodic Bandwidth \[7\].

Nodes can enter into sleep mode and wake up only when they have to receive and transmit any data to the gateway, because the nodes and the gateway are synchronized in time. The time slot allocation in Contention Free Period (CFP) is collision free, which improves packet transmission and thus, saves energy.

The protocol uses CSMA/CA in the uplink frame of Contention Access Period (CAP) period, which is not reliable scheme due to its unreliable CCA and collision issues.

## Ta-MAC Protocol

Traffic aware MAC (Ta-MAC) protocol utilizes traffic information to enable low-power communication. It introduces two wakeup mechanisms: a traffic-based wakeup mechanism, and a wakeup radio mechanism. Former mechanism accommodates normal traffic by exploiting traffic patterns of nodes whereas, later mechanism accommodates emergency and on-demand traffic by using a wakeup radio signal.

In the traffic-based wakeup mechanism, the operation of each node is based on traffic patterns. The initial traffic pattern is defined by the coordinator and can be changed later. The traffic patterns of all nodes are organized into a table called traffic-based wakeup table. In wakeup radio mechanism, a separate control channel is used to send a wakeup radio signal. The coordinator and the member node send wakeup radio signal in on-demand and emergency case.

In Ta-MAC, a node wakes up, whenever it has a packet to send/receive. Since the traffic patterns are pre-defined and known to the coordinator, it does not have to wait for resource allocation information/beacon. As a result, delay is minimized comparitive to other MAC protocols. This protocol accommodates normal, emergency and on-demand traffic in a reliable manner. To achieve energy efficiency in MAC protocol, the central coordination and resource allocation is based upon the traffic patterns of the nodes.

As, in this protocol, the traffic pattern are defined by the coordinator, in a static topology. Therefore, it does not work efficient in dynamic topology (in dynamic topology, traffic patterns are changed frequently).

```latex
\begin{tiny}
\begin{table}
\begin{tabular} { | m{2.5cm}| m{2.5cm}| m{2.5cm}| }
    \multicolumn{3}{c}{Table. 2. Qualitative Comparison of MAC Protocols} \\
    \hline
    \centering
   \textbf{Protocols} & \textbf{Advantages} & \textbf{Disadvantages} \\ \hline
    Okundu MAC  & Minimize time slot collision, reduce idle listening and overhearing   & Only 8 slave nodes can be communicated to MN   \\ \hline
    MedMAC &  Energy waste due to collision is reduced  & Do not support high data rate applications  \\ \hline
    Low Duty Cycle & Collision problem is reduced, allows  patients' monitoring & Not suitable for dynamic type of networks \\ \hline
    B-MAC & Improves packet transmission hence saves energy & Uses CSMA/CA in the uplink frame of CAP period, which is not a reliable scheme  \\ \hline
    Ta-MAC & Accommodates normal, emergency and on-demand traffic, energy efficient, reasonable delay & Not suitable for dynamic topologies  \\ \hline
     S-MAC & High latency and time synchronization overhead may be prevented due to sleep schedules & Low throughput, overhearing and collision may cause if packet is not destined to listening node  \\ \hline
      T-MAC & Packets are sent in burst and with low latency which collectively gives better result under variable load & Suffers from sleeping problems  \\ \hline
    H-MAC & Improves BSN's energy efficiency and reduces extra energy cost & Does not support sporadic events and posseseslow spectral/bandwidth efficiency   \\ \hline
     DTDMA & Reduce packet dropping rate, less energy consumption & Does not support emergency and on-demand traffic
    \\ \hline
     \end{tabular}
    % \caption{QUALITATIVE COMPARISON OF MAC PROTOCOLS}
   \end{table}
   \end{tiny}
```

## S-MAC Protocol

S-MAC \[8\] protocol is proposed for WBASNs. This protocol uses fixed duty cycles to solve idle listening problem. Nodes wakeup after a specific time, as assigned by coordinator, sends data and goes back to sleep mode again. As, all the nodes are synchronized, therefore, collision can also be easily avoided. S-MAC gives considerably low latency. In this protocol, time synchronization overhead may be prevented due to sleep schedules.

Fluctuating traffics are not supported and no priority is given to the emergency traffic scenarios by S-MAC. Therefore, it is not a reliable for WBANs. Overhearing and collision may occur if the packet is not destined to the listening node.

## T-Mac Protocol

Mihai *et al.* \[9\] suggested Time-out MAC (T-MAC) for WBASNs. It uses flexible duty cycles for increasing energy efficiency. In T-MAC, the node wakes up after time slot assignment, sends pending messages. If there is no activation event for Time Interval (TA), the node goes back to sleep mode again. If a node sends Route To Send (RTS) and does not receive Clear To Send (CTS), then sends RTS two more times before going to sleep. To solve early sleep problem, it uses future RTS for taking priority on full buffer.

In T-MAC, packets are sent in burst, as a result delay is minimized. It also outperforms other MAC protocols under variable load. The main disadvantage in this protocol is that it suffers from sleeping problems.

## H-MAC Protocol

Heartbeat Driven MAC (H-MAC) uses heart beat rhythm information for synchronization of nodes. This avoids the use of external clock and thus reducing the power consumption. Also guaranteed time slot (GTS) provision to each node helps to avoid collision.

H-MAC aims to improve BSNs energy efficiency by exploiting heartbeat rhythm information, instead of using periodic synchronization beacons to perform time synchronization \[3\].

Although, H-MAC protocol reduces extra energy cost of synchronization, however, it does not support sporadic events. Since TDMA slots are dedicated and are not traffic adaptive, H-MAC protocol encounters low spectral/bandwidth efficiency in case of low traffic. The heartbeat rhythm information varies depending on patient's condition. It may not reveal valid information for synchronization all the time \[3\].

## DTDMA Protocol

Reservation based dynamic TDMA (DTDMA) protocol uses slotted ALOHA in CAP field of super frame to reduce collisions and to enhance power efficiency.

Through the adaptive allocation of the slots in a DTDMA frame, WBAN's coordinator adjusts the duty cycle adaptively with traffic load. Comparing with IEEE 802.15.4 MAC protocol, DTDMA provides more dependability in terms of lower packet dropping rate and low energy consumption especially for an end device of WBAN \[3\]. It does not support emergency and on-demand traffic. Furthermore, DTDMA protocol has several limitations when considered for the Medical Implant Communication Service (MICS) band. The MICS band has ten sub-channels and each sub-channel has $300$ Kbps bandwidth. DTDMA protocol can operate on one sub-channel, however, cannot operate on ten sub-channels simultaneously \[3\].

The main purposes of a MAC protocol are to provide energy efficiency, network stability, bandwidth utilization and reduce packet collision.

The energy minimization techniques and mechanism in MAC protocols are summarized in Table. 3.

```latex
\begin{tabular}{| m {3cm}| m{4cm}| }
    \multicolumn{2}{c}{Table. 3. Energy Minimization Techniques and Mechanisms} \\
    \hline
    \textbf{Protocol} & \textbf{Energy Efficiency Mechanism}  \\ \hline

    Okundu MAC & Wake up Fall back Time (WFT)  \\ \hline
    MedMAC  & TDMA, Adaptive Guard Band Algorithm (AGBA) and Drift Adjustment Factor (DAF)\\ \hline
    Low Duty Cycle & TDMA, concept of Guard Time ($Tg$)\\ \hline
    B-MAC  & TDMA, Bandwidth mechanism \\ \hline
    Ta-MAC  & Central coordination according  to traffic patterns of the nodes \\ \hline
    S-MAC  & Scheduled based, organized in slots and operation based on schedules \\ \hline
    T-MAC  & Have slots and operation is based on schedules \\ \hline
    H-MAC  & Heartbeat Rhythm information is used for synchronization \\ \hline
    DTDMA & TDMA based, use of slotted aloha in CAP field\\ \hline

    \end{tabular}
```

# Performance Trade-offs made by MAC Protocols

In this section, we discuss the performance of the MAC protocols they achieve and price they pay. In other words, trade-offs, the MAC protocols have to make.

## Okundu MAC Protocol

Network's scalability is mainly application dependent, e.g., ECG can support upto maximum of 8 slave nodes because of 8 percent duty cycle. However, in practice this is 6 to allow for possible retransmissions. Therefore, we have a trade-off, for retransmission, slave nodes attached to the MN are reduced to attain scalability of network.

## MedMac Protocol

The low data rate applications of Class $0$ medical devices include monitoring of respiration system, temperature of human body, pulse monitoring etc. Power consumed by respiration transceiver is slightly high in MedMAC protocol with respect to other protocols, while temperature and pulse node show much low power consumption, as compared to other protocols. MedMAC trade-offs power consumption of respiration for less power of other two applications.

## Low Duty Cycle MAC Protocol

The number of extra slots needed for protocol robustness is dependent on Packet Error Rate (PER) and Packet Loss Ratio (PLR). When PER is high, it will increase PLR. However, PLR, may be reduced by using extra slots in the time frame. Therefore, this protocol can trade-offs extra slots for less PLR.

## B-MAC Protocol

B-MAC trade-offs idle listening for a reduced time to transmit and reception of data. As, we know that reducing duty cycle increases sleep time which in turn reduces idle listening. Another trade-off is between idle listening and packet length, because this overhead dominates the energy consumption.

## Ta-MAC Protocol

Ta-MAC uses two wakeup mechanisms one for handling data traffic and other for emergency traffic. By using these two mechanisms this protocol outperforms all other protocols in terms of power consumption because problems like idle listening, collision and overhearing are reduced. However, by sending frequent control messages to the nodes increases node's overhead, which is a trade-off. The initial traffic patterns of all the nodes are defined by the coordinator, as a result delay is also slightly increased.

## S-MAC Protocol

For transmission and reception of data in S-MAC, an extremely low duty cycle is used. When throughput increases SAC's duty cycle also increases , which further increases the overhead of SYNChronization (SYNC) period, as a result, power consumption is increase linearly. S-MAC can trade-offs throughput for energy, also it can trade-offs energy for latency.

## T-Mac Protocol

T-MAC uses adaptive duty cycle, implemented as a time out after the last event. At lower transmission rates, throughput increases because probability of packet loss is much less than received packet, however, the latency is increased between source and destination node.

## H-MAC Protocol

In H-MAC a Guard Band is introduced in time slots to avoid collision by overlapping of data, however, when time slots are completely aligned then there will be no data transmission in Guard Band, therefore, it reduces bandwidth utilization. The coordinator of BSN then uses this GB for synchronization, by sending re-synchronization control packets, hence achieving energy efficiency. Thus making a trade-offs between energy efficiency and bandwidth utilization efficiency.

## DTDMA Protocol

DTDMA is a TDMA based protocol which uses time slots for data transmission and as a result low power is consumed. However, TDMA requires synchronization between nodes and the coordinator, as a result, overhead is increased. This overhead is a trade-off for energy.

The trade-offs of each selected protocol are summarized and given in Table. IV.

```latex
\begin{tabular}{| m {3cm}| m{4cm}| }
    \multicolumn{2}{c}{Table. 4. Trade-offs Made by MAC Protocols} \\
    \hline
    \textbf{Protocol} & \textbf{Trade-Offs}  \\ \hline

    Okundu MAC &Trade-offs number of slave nodes attached to the MN are reduced for scalability of network\\ \hline
    MedMAC  &  Trade-offs idle listening for a reduced time to transmit and reception of data\\ \hline
    Low Duty Cycle & Can trade-off extra slots for less PLR\\  \hline
    B-MAC  &   Trade-off is between idle listening and packet length\\\hline
    Ta-MAC  &  Trade-off delay for low power consumption\\\hline
    S-MAC  &  Can trade-off energy for latency.\\\hline
    T-MAC  &  Can trade-off latency for high throughput\\\hline
    H-MAC  &  Trade-offs between energy efficiency and  bandwidth utilization efficiency\\\hline
    DTDMA &  Trade-offs overhead for a low power consumption\\\hline

    \end{tabular}
```

# MAC Frame Structure

MAC frame structure consists of control portion or control packet and data portion. Control portion is responsible for the management and control messages (beacon period, request period, topology management period) to control and manage dynamic topology and varying data rate traffic. Data portion consist of two sub parts: CAP and Contention Free Period (CFP). CAP consists of CSMA/CA while the nodes contend in CAP transmit MAC control packets. Similarly, small size data packets can also be transmitted in CAP.

In \[7\], the allocation of time slots is controlled by the coordinator. The coordinator arrange the duration of control and data packet on the basis of current traffic of topology that is why the slots are allocated to CFP and are collision free. In each frame, bandwidth allocation in CFP can be changed.

In \[5\], the GB is used to maintain synchronization among devices even if a node is sleeping for many beacon periods.

In \[4\]\[9\]\[10\]\[11\], MAC protocols use slotted ALOHA in its frame structure to divide a slot into 4 equal mini slots. In \[6\], $Tg$ is introduced in its frame structure to reduce overlapping between the two following nodes.

# Technique for Collision Avoidance for Traffic Control

The main schemes of MAC protocol for WBANs are divided into two groups: contention based i.e., CSMA and contention free i.e., TDMA. Most of the traffic is interrelated in WBANs, therefore, contention based solutions are not suitable for them. For example, if a patient is suffering from fever, the body temperature increases which increases blood pressure, hence, the sensor sensing temperature variation and the sensor that senses blood pressure variation, both become active. Along with them other respiration sensors also become active at the same time and try to access the channel/coordinator. However, in this situation, collision occurs in CSMA. In TDMA, each node communicate to MN according to the assigned pattern by the coordinator. As a result, collision in data traffic is low, as compared to CSMA.

## Okundu MAC

This protocol controls traffic using centrally controlled wakeup/sleep time. Slots are assigned to sensors change every time when coordinator detects any change in traffic pattern. Assignment of different time slots, decreases collision between the nodes . It makes the system to handle fluctuating traffic. The sensor nodes establishe link with the coordinator after listening to the Radio Frequency (RF)-channel for a fixed time period. MN sends request to the sensor node for information by setting and communicating the next wakeup time after establishing the link.

| Table. 5. Comparison between IEEE 802.15.4 MAC and Original IEEE 802.15.4 |  |
|:--:|:--:|
| **IEEE 802.15.4 MAC** | **Original IEEE 802.15.4** |
| Low power consumption | High power consumption |
| Higher data rate | Low data rate |
| Higher flexibility | Low flexibility |
| TDMA based | Contention based |
| Collision Free | Greater collisions |
| Sleep mode | Idle listening |

## MedMAC

MedMac reduces the collision by using AGBA. AGBA allows the sensor nodes to sleep for a GB time period between each time slot. Each node has specific time slot to communicate with master node/coordinator, which means there is no collision. Thus minimizes the synchronization overhead.

## Ta-MAC

Ta-MAC protocol uses two channel access mechanisms for traffic control i.e., traffic based wakeup mechanism for normal traffic, and wakeup radio mechanism for on-demand and emergency traffic. In traffic-based wakeup mechanism, all nodes have traffic patterns that are assigned by the coordinator. The initial patterns are defined and updated by the coordinator. The traffic patterns of all nodes are synchronized and arranged in a specific table, known as $Traffic$ $Based$ $Wakeup$ $Table$. Node's ID and its respective traffic patterns are stored in this table.

Normally, all the nodes become active/wakeup according to their traffic patterns. If two or more nodes have same wakeup pattern then the node with high priority is treated first by the coordinator, as shown in Fig. 3. By assigning these patterns, load at the coordinator is minimized, and chances of collision is also reduced.

## B-MAC

B-MAC uses downlink and uplink schemes along with sleeping mode for data traffic control. Downlink is only used by MN, therefore, traffic and data load on downlink are reduced. Uplink is divided into CAP and CFP. MN allocates time slots to CFP according to data traffic which makes CFP collision free. In case, when nodes have no data to transmit or receive then they go to sleeping mode.

## Low Duty Cycle

Low duty cycle MAC protocol is based on TDMA. In TDMA, time slots are assigned to the sensor nodes by the coordinator. To avoid collision between the data traffic, the concept of $Tg$ is introduced. Use of $Tg$ between every consecutive slots prevents the transmission overlaps and controls data traffic.

## IEEE 802.15.4 MAC

The basic requirement of QoS is to minimize delay and maximize the probability of successful transmission. CFP scheme is used to control data traffic to guarantee the QoS. If a node wants to send data, first it listens for the network beacon. After node finds the beacon that is sent by the coordinator, the node synchronizes to the super frame structure.

IEEE 802.15.4 supports up to $250$ Kbps data rate with possible coverage of $10$ meters. This data rate is not enough to support the required rates of WBANs that is up to $10$ Mbps. According to IEEE 802.15.4, packets are transmitted in the contention period, which may result longer delays in real time critical applications. When traffic is increased, the nodes compete for the contention based slots, resulting in long delays and the actual size of the network is almost doubled \[7\].

In \[10\], to satisfy the requirements of WBANs including QoS, network scalability, support for multiple PHY's and multiple application traffics, IEEE 802.15.4 MAC is proposed. It is the modified version of original IEEE 802.15.4. QoS means to decrease the packet latency and increase the probability of successful transmission of data packets without collision and loss of data. In original IEEE 802.15.4, GTS mechanism is provided to support the emergency data. GTS is very effective for data transfer, however, inherently the limit of GTS in a super frames is seven. As a result, it cannot support more than seven devices simultaneously in CFP. Whereas, in IEEE 802.15.4 MAC, the coordinator may allocate more than seven GTS simultaneously to the sensor devices.

IEEE 802.15.4 MAC and 802.15.4 original is compared briefly in Table. 5.

# Conclusion

We present a survey of different MAC protocols with respect to energy efficiency and their advantages and disadvantages in WBANs. Low power listening, scheduled contention and TDMA are also compared. It is observed that TDMA is more power efficient, however, suffers with synchronization sensitivity. Techniques for collision avoidance of different MAC protocols are also comparatively analyzed. Path loss model for In-body, On-body and Off-body communication in WBANs is also described. Because human body is composed of tissues and organs in which communication is difficult and thus results in high path loss. On-body and Off-body also show some variations in results when the source and destination sensors or nodes are LoS and NLoS.