Volume 9 Supplement 1

## The OptAIDS project: towards global halting of HIV/AIDS

- Research
- Open Access

# Can we spend our way out of the AIDS epidemic? A world halting AIDS model

- Robert J Smith
^{1}Email author, - Jing Li
^{2}, - Richard Gordon
^{3}and - Jane M Heffernan
^{4}

**9 (Suppl 1)**:S15

https://doi.org/10.1186/1471-2458-9-S1-S15

© Smith? et al; licensee BioMed Central Ltd. 2009

**Published:**18 November 2009

## Abstract

### Background

There has been a sudden increase in the amount of money donors are willing to spend on the worldwide HIV/AIDS epidemic. Present plans are to hold most of the money in reserve and spend it slowly. However, rapid spending may be the best strategy for halting this disease.

### Methods

We develop a mathematical model that predicts eradication or persistence of HIV/AIDS on a world scale. Dividing the world into regions (continents, countries etc), we develop a linear differential equation model of infectives which has the same eradication properties as more complex models.

### Results

We show that, even if HIV/AIDS can be eradicated in each region independently, travel/immigration of infectives could still sustain the epidemic. We use a continent-level example to demonstrate that eradication is possible if preventive intervention methods (such as condoms or education) reduced the infection rate to two fifths of what it is currently. We show that, for HIV/AIDS to be eradicated within five years, the total cost would be ≈ $63 billion, which is within the existing $60 billion (plus interest) amount raised by the donor community. However, if this action is spread over a twenty year period, as currently planned, then eradication is no longer possible, due to population growth, and the costs would exceed $90 billion.

### Conclusion

Eradication of AIDS is feasible, using the tools that we have currently to hand, but action needs to occur immediately. If not, then HIV/AIDS will race beyond our ability to afford it.

## Keywords

- Susceptible Individual
- Emigration Rate
- Cumulative Cost
- Trivial Equilibrium
- Metapopulation Model

## Introduction

There has been a sudden increase in the amount of money donors are willing to spend on the worldwide HIV/AIDS epidemic. Present plans are to hold most of the money in reserve and spend it slowly. However, rapid spending may be the best strategy for halting this disease. The purpose of this paper is to estimate how and when to best spend the available money to achieve halting HIV/AIDS as soon as possible.

In the year 2000, leaders from around the globe derived a set of goals aimed at making the world safer, healthier and more equitable. One of these goals was to halt and begin to reverse the global HIV/AIDS epidemic by 2015. To realise this goal, in 2006, world leaders pledged to move towards universal access to HIV prevention, treatment, care and support by 2010. In the past decade, there has been a marked decrease in the number of AIDS-related deaths worldwide and the percentage of individuals living with HIV has begun to plateau. Indeed, a number of highly infected countries have seen declines in the number of new HIV infections [1]. This success is owed largely to a six-fold increase in financing for HIV activities in low- and middle-income countries enabling better access to HIV treatments and care. However, the picture is not uniformly positive and some countries continue to see an increase in infections. Also, in spite of the improved access to HIV treatment and care, the majority of those who need antiretroviral therapies are not currently receiving them [1]. As a result, the epidemic is outpacing the rate at which these drugs are being delivered. In fact, the number of new infections is significantly greater than the increase in the number of people on antiretrovirals. Eradication of HIV/AIDS is therefore not achievable using the current funding model. In the absence of additional funding, we contend that a significant improvement in results can be achieved by optimising implementation, i.e. how donated funds are spent on drug delivery, educational programs and the distribution of HIV/AIDS prevention measures.

While estimates on spending have been made for many scenarios [2–4], these are usually projected in decades or generations. World events could overtake the HIV/AIDS epidemic on this time scale, or even be precipitated by it if it continues to grow [5–9]. Furthermore, donor fatigue [10] could return, or other crises (tsunamis, wars, etc) could divert presently committed funds. Is there a quick fix that could avoid passing the problem on to the next generation?

Consider, for example, Bill and Melinda Gates, who have made stopping AIDS the top priority of their foundation [11], with $60 billion available [12]. As an extreme case, for instance, ignoring the ability of economies to absorb so much cash [13], the Gates Foundation might spend all their money in two years, meeting the worldwide needs of $27 billion per year for all health needs of all developing countries [14] (including HIV/AIDS), with the hope that infrastructure would then be in place and the rest of the donors would pick up from there. The projected spending rate of the Gates Foundation is $3 billion per year [15], which is roughly comparable to the expected annual interest. It is being spread out over 20 years, more than doubling the duration of the epidemic, which has already lasted over 25 years. Something in between a jump start and a trickle would seem much more optimal than either of these extremes.

Furthermore, if we fail to spend at the optimal rate, the epidemic may blow up. Currently, the percentage of people living with HIV has stabilised; however, the overall number of people living with HIV has steadily increased [1]. The result is that the epidemic could become totally out of control or cost immensely more to bring under control. Many modeling studies have shown a bigger bang for the buck if the money is spent earlier in an epidemic [16, 17], but none have investigated the optimal course of spending to bring an epidemic under control as rapidly as possible and avoid a growing gap [18].

A number of views have been expressed on how to tackle HIV/AIDS. One approach involves getting back to the fundamentals of epidemiology [19] and avoiding AIDS exceptionalism [20], i.e. treating HIV/AIDS just like other epidemics. Most people advocate a mixture of prevention and treatment options [17, 21]. One worldwide single (but not magic) bullet approach that has been proposed would use the fact that highly active antiretroviral therapy (HAART) reduces infectivity:

"This optimistic population-based model shows that, in 45 years, HIV prevalence could be reduced by more than 70 times from more than 7 cases per 1000 people to less than 0.1 case per 1000. The number of HIV- infected people could be reduced from 38 million to less than 1 million. The cost of therapy would be about US $7 billion per year, with costs declining from $15 to $1 billion. Such a programme would be expected to cost $338 billion over 45 years" [22].

The average cost over the 45 years would be $7.5 billion per year.

Condoms offer another single-mode approach. Above an estimable and possibly attainable threshold, condoms could halt an HIV/AIDS epidemic [23]. If distributed to 1.5 billion men (half the male population of the world), at 100 condoms per person per year [24], and a cost of $0.02 each [25], total cost would be $3 billion per year. (Female condoms have their place, but are about ten times more expensive each.) For comparison, the USA spends $12 billion per year domestically on HIV/AIDS, plus $3 billion per year on HIV/AIDS research [26].

Condoms have the advantage of reducing other STIs (sexually transmitted infections), some of which enhance HIV infectivity. The threshold for decline of the epidemic is the product of compliance in use and condom effectiveness, and there is room for technical improvement of the latter [24]. The problems to be tackled are not simply behavioural, but also involve design, manufacturing, marketing and subsidisation [27]. Even in the behavioural realm, the emphasis on risk factors instead of collective behaviour may be counterproductive [28, 29].

The HIV/AIDS epidemic is often spoken of in terms of reducing the spread [16], or achieving sustainable financing [30]. Decades-long investments have been made in vaccines [31] and more recently in microbicides [32]. These approaches are slow in coming, and do not take into account the greatly increased amount of cash now available, at least for the short term. The donor community has created the potential of stepping ahead of the exploding HIV/AIDS epidemic by pooling money faster than the disease is spreading. Our aim is to use mathematical modelling to find ways to spend the available money as fast, efficiently and effectively as possible to bring this epidemic to a full halt, now, while this transient advantage is still available.

A large number of mathematical models have been used to describe the epidemiology of HIV [33–46]. These range from the simple to the complicated, and have used a variety of techniques to analyse and determine conclusions. One such technique is that of metapopulation models: dividing space into distinct regions, within which disease transmission occurs (often at different rates) and with "travel" between regions allowing different regions to be connected. Such models have been used to describe HIV in a spatial context without resorting to partial differential equations [47–53].

This paper presents a simple, linear, metapopulation model that describes the world regions at a number of different scales (continent, country, etc). We use this model to show that interventions can lead to eradication, but such interventions must be applied globally, not locally, and that time is of the essence if our interventions are to succeed.

## An eradication model

In its simplest form, HIV infection can be described by the classical SI (susceptible-infected) model. The dynamics are not linear, depending, at a minimum, on the behaviour of susceptible (uninfected) individuals, as well as their interactions with those who are infected. While the SI model does not capture all the dynamics of HIV, it has been used to model the disease throughout its entire history [54–60].

Here, Λ is the rate of appearance of new susceptibles, *β* is the infection rate, *μ* is the background death rate and *γ* is the death rate due to disease.

*λ*= -

*μ*, -

*μ*-

*γ*. It follows that the threshold for eradication [61] in the SI model is

The SI model, while it captures the dynamics of HIV spread in a population, is difficult to manipulate in a metapopulation scenario. We propose instead a linear approximation of this model. The approximation of the epidemic as linear in the number of infectives is clearly a coarse one. However, although this model structure does not capture the transient dynamics of infection and interaction, it nevertheless serves as a predictor for eradication.

This is because the major stability analysis carried out on such models - finding the basic reproductive ratio - consists of linearisation of the nonlinear models. Thus, a linear version of the model has the advantage of having the same invasion threshold, although it has the disadvantage of only describing local stability of the disease-free equilibrium.

will always overestimate the epidemic. We shall refer to this as the I-only model.

This model has only the trivial equilibrium *I* = 0. If the trivial equilibrium is stable, then all trajectories will approach it. If the trivial equilibrium is unstable, then solutions will increase without bound.

which is the same eradication threshold as the SI model (as expected). It follows that there will be eradication in the I-only model if and only if there is eradication in the SI model.

*R*

_{0,SI}< 1 (Figure 1) and

*R*

_{0,}

_{ SI }> 1 (Figure 2). Parameters used in the simulations were Λ = 20 people·years

^{-1},

*μ*= years

^{-1},

*γ*= years

^{-1}, with

*β*= 0.00007 people

^{-1}years

^{-1}(Figure 1) and

*β*= 0.0002 people

^{-1}years

^{-1}(Figure 2). Despite the fact that the transient dynamics are vastly different in the two models, the linear approximation has the same eradication threshold as the more accurate SI model and, furthermore, always overestimates the epidemic. It follows that, for eradication purposes, the simple, linear model should determine whether our control methods will be sufficient.

We now extend the *SI* model to a metapopulation model with *p* regions. For 1 ≤ *i* ≤ *p*, we can write

*S*_{
i
}= susceptible individuals in the *i*th region;

*I*_{
i
}= infected individuals in the *i*th region;

Λ_{
i
}= the rate of appearance of new susceptible individuals in the *i*th region;

*μ*_{
i
}= background death rate in the *i*th region;

*γ*_{
i
}= death rate due to disease in the *i*th region;

*n*_{
ij
}= migration rate of susceptible individuals from *j*th region to *i*th region;

*m*_{
ij
}= migration rate of infected individuals from *j*th region to *i*th region.

If the system (2) is well-posed, we can find a bound for the population of susceptibles in each patch. We can then extend the same idea to obtain the linear system of *I*-equations only. (See Appendix for details.)

Thus, the total population of susceptible individuals in the *i*th region is limited. So we can write

*π*_{
i
}= *β* × (total population in the *i*th region without infection)

*d*_{
i
}= *μ*_{
i
}+ *γ*_{
i
},

*d*

_{ i }represents the total death rate in the

*i*th region. Thus, for

*i*= 1,...,

*p*, the equations of the

*I*-only model are

## Analysis

### A two-region example

*p*= 2, which can be presented as follows:

This is the eradication threshold condition for this simple case *p* = 2. The case for *p* = 2 can be further simplified under different assumption of the travel (or migration) between the two regions. Here we list three possibilities.

#### Case 1

*i*= 1, 2,

#### Case 2

#### Case 3

which is also the eradication condition for the disease in this case.

From this, one can derive conditions for the travel rate of susceptible individuals which can either help a disease which is otherwise dying out to persist locally or cause an otherwise partially persistent disease to persist globally in the two regions (see Appendix for the details). It follows that the disease must be reduced within each region and the effects of travel between regions must be accounted for. This explains why HIV eradication must be considered as a world problem, not just a problem for individual countries, or continents, to tackle independently.

#### Finding the eradication threshold condition *T*_{0}

**K**. The only equilibrium for the

*I*-only model (3) is the trivial, disease-free equilibrium (0,...,0). The linear stability of the trivial equilibrium for (3) is determined by

*s*(

**K**), the stability modulus of

**K**, which is defined as the maximal real part of all eigenvalues of the matrix

**K**. Thus, we can define the overall eradication threshold as

The trivial equilibrium is globally asymptotically stable if *T*_{0} < 1 (for this linear system (3), the local stability is equivalent to global stability) and unstable if *T*_{0} > 1.

Noting that the right-hand side of *I*-equations in (2) is never bigger than the right-hand side of (3) and then by the standard comparison theorem (see [63, 64]), one can show that *T*_{0} is also the condition to guarantee the stability of the disease-free equilibrium of system (2). This implies that HIV can be eradicated under the condition that *T*_{0} < 1, but it can not be eradicated under the condition that *T*_{0} > 1.

## A continent-level example

Birth rates and death rates by continent [65].

Population | Births/Popn | Birth Rate | Deaths/Popn | Death Rate | |
---|---|---|---|---|---|

AF | 954879489 | 33203380 | 0.0348 | 12625314 | 0.0132 |

AS | 4043347897 | 76536061 | 0.0189 | 27947177 | 0.0069 |

EU | 729546003 | 7321919 | 0.0100 | 8437988 | 0.0116 |

NA | 527814776 | 8712554 | 0.0165 | 3798366 | 0.0072 |

OC | 33970173 | 545421 | 0.0161 | 244250 | 0.0072 |

SA | 383907961 | 6773241 | 0.0177 | 2354551 | 0.0061 |

Foreign born. Number of foreign born by continent of origin (vertical) and continent of destination (horizontal) [66].

Des/Ori | AF | AS | EU | NA | OC | SA |
---|---|---|---|---|---|---|

AF | 15973 | 16987 | 4263 | 0 | 0 | |

AS | 286806 | 298431 | 13062 | 1621 | 32500 | |

EU | 5312095 | 3566013 | 1167954 | 252860 | 1337972 | |

NA | 1176374 | 9690228 | 8573379 | 353095 | 4466748 | |

OC | 221003 | 1380652 | 2470078 | 136804 | 76873 | |

SA | 16595 | 218415 | 1087422 | 197314 | 3324 |

*T*

_{0}. The off-diagonal entries of our matrix are the off-diagonal entries from Table 2. The diagonal entries of our matrix are calculated according to the formula

where Λ_{
i
}is the influx of new susceptibles, *μ*_{
i
}is the death rate (Table 1), *β* is the infection rate, *γ* is the AIDS-specific death rate, *η*_{
ij
}is the number of individuals born in continent *j* residing on continent *i* per year (Table 2) and *P*_{
i
}is the population of continent *i* (Table 1).

*β*= 0.0002 people

^{-1}years

^{-1}and

*γ*= 0.1 years

^{-1}, reflecting the fact that the typical time course of the disease is ten years. We calculate the influx of new susceptibles according to the formula

where *p*_{
i
}is the continent's birth rate and *c* is an infected individual's annual number of (uninfected) sexual partners. We set *c* = 20; however, the factor
modifies an individual's annual number of sexual partners according to the growth (or decline) of the continent they are in. Thus, for Europe, the annual number of sexual partners is 17, whereas for Oceania the number is 44. These numbers are likely overestimates, but our intention is to be conservative in our approach.

If we now intervene and reduce *β* to 0.0002 x 2/5, then the overall threshold is
. It follows that reducing the infection rate to 0.00008 will result in eradication.

*β*by three fifths is, of course, a mammoth task. Furthermore, this model implicitly assumes that such a reduction is instantaneous. To compare this to slowly reducing

*β*over a period of decades, we assumed that

*β*would be reduced to two fifths of its current value in twenty years' time, but that reductions would occur evenly throughout that time. However, we also included a population increase of 3% per year. The results are shown in Figure 4. In this case, the population "overruns" our goal. The overall

*T*

_{0}is not reduced below one, because the population in twenty years' time will be too high for a 3/5 reduction in the infection rate to lead to eradication.

## Costs

Costs of HIV/AIDS intervention methods.

Intervention | Cost (USD) | Reference |
---|---|---|

HIV testing in low-income countries | 200-4000 per test | [74] |

Combination antiretroviral therapy in high-income countries | 10000-15000 per patient year | [75] |

Combination antiretroviral therapy in low-income countries | 350-4000 per patient year | |

Costs of monitoring viral load and cell count | 25-100 per test | [76] |

Health education | 500-3000 per patient | |

Patient out-of-pocket expenses to receive medical care in a low-income area | 60-250 | [77] |

Adult male cirumcision | 40-100 | [82] |

Male condom | 0.02 per condom | [25] |

Along with the cost of delivering drugs to a patient, it is also important that we determine what benchmark of viral load and cell count we will use for initiation of therapy. It is obvious that the costs of delivering therapy and monitoring will increase if these are started at an earlier time in HIV infection. For example, it is common practice to initiate therapy in the developed world when CD4 counts are less than 350 cells/mm^{3}. However, in the developing world, it is advised that treatment start when CD4 counts are less than 200 cell/mm^{3}. In a study by Vijayaraghavan *et al.* [67], it was found that treating patients with HIV according to developed- versus developing-world guidelines is highly cost-effective and may result in substantial long-term savings. We choose to initiate individuals earlier in life for our cost analysis. This means that we will include a cost of therapy in the upper range of that listed in Table 3. However, we will only choose this cost from those for developing countries since the donor funds are applied mainly to the epidemic in the developing world.

Note that we assume that the costs of delivering therapy and health care to HIV infected individuals will stay relatively constant over a short period of time. Since the donor money will be spent over a span of a few years rather than decades, this assumption can be expected to give realistic results. If we were to determine the cost of such interventons over a long period of time, we must also include variables describing a decrease in the costs of therapy and testing as drug manufacturing becomes more efficient, the discovery and development of new drugs which will affect the cost of therapy, the increase in education and other medical expenses (since administrative costs and the payment of healthcare workers will increase over time), etc.

To effectively determine the cost of treatment, we also must include the possibility of the evolution of drug resistance [37]. It is difficult to put a dollar figure on the evolution of drug resistance, since it may occur during different periods of infection in different patients. In our cost analysis, we choose costs of treatment and testing in the upper range of that reported (see Table 3), since it is assumed that the evolution of resistance will increase the cost of treatment and require more tests to monitor a patient's viral load and cell count.

Referring back to model (3), it is evident that condoms and education will reduce *π*_{
i
}, while treatment will reduce *π*_{
i
}, but will also decrease *d*_{
i
}. Education will include condom awareness, but may include other factors, such as abstinence and monogamy [68]. Travel/immigration restrictions will reduce *m*_{
ij
}and *m*_{
ji
}, but we do not expect that these will be realistic; past attempts to restrict movement based on travel or immigration have simply driven the epidemic underground [69].

Note that we assume a population growth of 3% per year. This corresponds to the maximum population growth rate reported by the CIA World Factbook. We use the maximum growth rate so that cost estimates give an upper bound of what may be realistic for donor spending. We also assume that HIV-positive cases will continue to increase at a rate of 3% per year, if no interventions are undertaken.

where *C* is the cumulative cost, *n* is the proportion of men who must ultimately receive condoms, *m* is the fraction of infected individuals who receive treatment, *r* is the timescale, ϵ is the (fixed) cost of distribution and education, and *θ* is the cost of treatment (including the cost of paid health workers, testing etc). Thus, in the first term of (17), there are 3 billion men, whose numbers increase by 3% every year, given 100 condoms each, at a cost of $0.02, with a fixed condom education campaign of $60 million. In the second term of (17), we assume that one tenth of infected individuals require treatment; thus, currently there are 3.3 million individuals requiring treatment, whose numbers increase by 3% every year, at an average cost of $2500 per patient per year.

assuming no fixed costs of condom education.

Thus, treatment is significantly more expensive than prevention (as expected).

Thus, we could reach our eradication threshold within five years, provide treatment to everyone who needs it and stay within our $60 billion budget (plus interest). See Figure 5B.

## Discussion

We have provided a method for predicting treatment and prevention levels necessary to eradicate HIV/AIDS, based on population and immigration data at the continental level. This method is easily applicable to a finer scale, such as country-level data. The mathematical model is linear, which allows it to be easily generalised. We have also developed a formula to estimate the cost of some of these interventions.

It should be stressed that our model does not attempt to quantify the prevalence of the disease or the time course of infection. The model is a predictor of eradication only and thus should not be used in other contexts. However, it provides us with insight into the degree of intervention measures that are needed to eradicate the disease: an across-the-board reduction in the infection rate by a factor of three fifths or more, and antiretroviral treatment for everyone who requires it.

These are steep requirements that will cost a lot of money. However, we have demonstrated that spending only some of the money now and saving the rest for the future will result in a loss of ability to eradicate the disease, using existing estimates. We have clear, specific goals and they need to be acted upon immediately.

Our model also demonstrates the folly of focusing locally. Even if the basic reproductive ratio in all continents were reduced below one, the disease could still be sustained by travel or immigration (see Appendix). We note that travellers/immigrants do not present a high risk factor at present, but they may be a significant obstacle when the disease is close to eradication. This is why we must consider HIV/AIDS from a global perspective.

It should be noted that, in addition to our model's obvious shortcomings, our parameter estimates also have limitations. In particular, the data we use currently assumes that the disease-specific death rate is negligible, compared to the background death rate. While this may be approximately true in some western countries, it is clearly not the case in Africa [70]. However, while the numbers used here are illustrative, our calculations of *T*_{0} are based on analytical thresholds. Thus, the sensitivity of the results to our parameter choices is easily measured: variations in the data in Tables 1 and 2, as well as formulas (15) and (16) allows similar calculations to be performed at other scales than continental (e.g. country-level).

Another assumption that should be noted is that we assume that travellers/immigrants receive no intervention help, such as education or treatment. This may be likely if foreign nationals receive no health care and have limited access to education, but will not be universally true. However, in each case we have overestimated the severity of the disease, in order to provide a strict criterion: if our overall *T*_{0} - calculated from the linear model with the overestimated parameters - is less than one, then we will have eradication. Our immigration/emigration data did not include tourism, especially sex tourism [71] which potentially adds a huge effect to the impact of travellers.

Although the cost of a condom is relatively cheap compared to other intervention measures, it is necessary to acknowledge the fact that a condom taken does not mean a condom used. Condom promotion requires a well- developed program infrastructure with excellent logistic competence of its staff which may be difficult to achieve. Many areas, especially rural and remote villages, are not accessible; for example, the healthcare system covers less than 50% of the population in many sub-Saharan African countries. In addition, the effectiveness of condoms in HIV prevention is much less than 100%, perhaps mostly due to condom failure [24, 72]. In addition, there are civil-war zones, such as eastern Congo, where services cannot be provided.

This further limits the lack of accessibility to reach out universally to all parts of the populations. However, model (3) has the property that decreasing *T*_{0} will still be beneficial. Thus, even if intervention methods are imperfect, a large-scale attempt at eradication will still be enormously effective.

Similarly, although we have estimated treatment costs, we have used aggregate numbers that do not explicitly take into account many of the specific costs associated with treatment, such as healthcare worker training. It has been shown that 100% HAART coverage in sub-Saharan Africa could be achieved in 10 years, if the number of healthcare workers were doubled every year [73], but this becomes expensive quite quickly.

There are, of course, many other intervention methods that could be examined, such as circumcision, STD control etc. These methods could be applied together with the proposed interventions; e.g. condom promotion could be combined with syphilis treatment without adding lots of additional costs. As we have learned from family planning programs, the more methods are offered to a client and the greater the choice is, the more likely the client is using one. This is likely to apply for HIV/AIDS control as well; more methods are advantageous and more likely to lead to positive results. Recently, UNAIDS has also made this proposition to promote as many methods as are available and proven effective.

It is currently advised that countries with different epidemic patterns will require different national strategies for implementing effective HIV prevention programs. Countries classified as low-level, concentrated, generalised or hyperendemic should concentrate their efforts on different national strategies.

In our model, we have applied the same strategy to all regions and each region includes countries with one or more of these classifications. We note, however, that the education component could be specific to each continent, or to each country.

Furthermore, although model (3) predicts eradication if *T*_{0} < 1, the threshold also has the property that lower values are beneficial. Thus, even if actual eradication cannot be achieved (as may be the case, in reality), intervention methods such as those mentioned here that reduce *T*_{0} will nevertheless be beneficial.

Thus, for example, if some countries make significant inroads against the disease, it will still have a global benefit even if eradication does not occur.

We view this model as a starting point. It provides a criterion for eradication, but we believe that a more accurate model, one that reflected the time course of the disease in all regions, is achievable. However, such a model requires accurate input from all countries; ideally, it could be updated as more data becomes available, or as new intervention methods are developed. Disease eradication has occurred on a global scale before, as happened with smallpox and is currently underway for polio. While the mechanics of organising the international community in this way are beyond the scope of this paper, given the scale of the problem and the immediacy of the required solution, we believe that this model can provide us with an easy-to-grasp way of understanding what needs to be done and when to do it. For the former, our strategy involves harnessing all existing intervention techniques in significant strength. For the latter, the time is now.

## Appendix

### Well-posedness of system (2)

The following theorem assures that the system is well-posed.

#### Theorem 1

*If* (*S*_{10},⋯, *S*_{p0}, *I*_{10},⋯, *I*_{p0}) ∈ *D*_{
L
}, *then*, *for any L* >*L**, *the set D*_{
L
}*is positively invariant for solutions of (2) with (18).*

#### Proof

**S**(

*t*) = (

*S*

_{1}(

*t*),⋯,

*S*

_{ p }(

*t*))

^{ T },

**I**(

*t*) = (

*I*

_{1}(

*t*),⋯,

*I*

_{ p }(

*t*))

^{ T },

**Λ**= (Λ

_{1},⋯,Λ

_{ p })

^{ T },

*i*= 1,...,

*p*. Noting that the off-diagonal matrix elements of

**A**(

*t*) are nonnegative, we conclude that the entries of the matrix are all nonnegative. Indeed, let and rewrite

**A**(

*t*) as

**E**denotes the

*p*×

*p*identity matrix. Then all entries of (

*t*) are nonnegative, and hence so are entries of . We also have

implying **I**(*t*) ≥ 0 for *t* ≥ 0.

Finally, we show that *S*_{
i
}(*t*) and *I*_{
i
}(*t*) are bounded for *t* ≥ 0 and *i* = 1,..., *p*. Let *N*(*t*) = *S*_{1}(*t*) + ⋯ + *S*_{
p
}(*t*) + *I*_{1}(*t*) + ⋯ + *I*_{
p
}(*t*). By the nonnegativity of *S*_{
i
}(*t*) and *I*_{
i
}(*t*), for *i* = 1,..., *p*, we

This implies that *N*(*t*) is bounded with an upper bound
; hence, so are *S*_{
i
}(*t*) and *I*_{
i
}(*t*)

for *t* ≥ 0 and *i* = 1,..., *p*. This completes the proof. □

### Comparison of Case 1 and 2 in Section 3.1

It is obvious from the formulas in (13) that
and
reflect the influence of travel of susceptible individuals between the two patches, and hence may be called the *travel-modified basic reproductive ratios* for Region 1 and Region 2, respectively. The following observations are direct consequences of (12)-(13).

Their verifications are straightforward and are thus omitted.

then > 1 and < 1. By symmetry, the conditions parallel to the above can lead to < 1 and > 1. Here, and in the following, we omit such parallel conditions.

then > 1 but > 1.

The biological meanings of **(A1)-(A2)** can be obtained in terms of the threshold condition meaning of
and
. For example, **(A1)** implies that travel of the susceptible individuals can help an otherwise dying out disease persist locally. Roughly speaking, a larger inflow of susceptible individuals favours the persistence of the disease in the region. **(A2)** shows that appropriate travel rates may cause an otherwise partially persistent disease to persist globally in both regions.

## Declarations

### Acknowledgements

We thank Shoshana Magnet, Andrea Brandon, Anthony Wilson, Avril Newman, Matt Davison, Frederick Chen and Beni Sahai for technical discussions. We are also grateful to two anonymous reviewers for comments that greatly improved the manuscript. RJS is supported by an NSERC Discovery Grant, an Early Researcher Award and funding from MITACS. JMH is supported by an NSERC Discovery Grant and funding from MITACS. Handling editor for this manuscript was Connell McCluskey.

This article has been published as part of *BMC Public Health* Volume 9 Supplement 1, 2009: The OptAIDS project: towards global halting of HIV/AIDS. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2458/9?issue=S1.

## Authors’ Affiliations

## References

- UNAIDS: Report on the global AIDS epidemic. 2008, [http://www.unaids.org/en/KnowledgeCentre/HIVData/GlobalReport/2008/2008_Global_report.asp]Google Scholar
- Piot P: AIDS: from crisis management to sustained strategic response. Lancet. 2006, 368 (9534): 526-530. 10.1016/S0140-6736(06)69161-7.View ArticlePubMedGoogle Scholar
- Stover J, Bertozzi S, Gutierrez JP, Walker N, Stanecki KA, Greener R, Gouws E, Hankins C, Garnett GP, Salomon JA, Boerma JT, De Lay P, Ghys PD: The global impact of scaling up HIV/AIDS prevention programs in low- and middle- income countries. Science. 2006, 311 (5766): 1474-1476. 10.1126/science.1121176.View ArticlePubMedGoogle Scholar
- UNAIDS: AIDS in Africa: Three Scenarios to 2025. 2005, Geneva: Joint United Nations Programme on HIV/AIDS, [http://www.unaids.org/unaids_resources/images/AIDSScenarios/AIDS-scenarios-2025_report_en.pdf]Google Scholar
- Barnett T, Whiteside A: AIDS in the Twenty-First Century: Disease and Globalization. 2006, Basingstoke, Hampshire: Palgrave Macmillan, [http://www.unaids.org/unaids_resources/images/AIDSScenarios/AIDS-scenarios-2025_report_en.pdf]2Google Scholar
- Bloom DE, Rosenfield A: A moment in time: AIDS and business. AIDS Pat Care and STDs. 2000, 14 (9): 509-517. 10.1089/108729100438908.View ArticleGoogle Scholar
- Bloom DE, Mahal A: Does the AIDS epidemic threaten economic growth?. J Econometrics. 1997, 1-296.Google Scholar
- Haaker M: The Macroeconomics of HIV/AIDS. 2004, Washington DC: International Monetary FundGoogle Scholar
- Tandon A: Macroeconomic Impact of HIV/AIDS in the Asian and Pacific Region. 2005, Manila: Asian Development BankGoogle Scholar
- Ewing T: AIDS programme faces donor fatigue. Nature. 1990, 346 (6285): 595-PubMedGoogle Scholar
- Gates W, Gates M: Keynote Address, Opening Session, XVI International AIDS Conference, August 13 2006, Toronto. 2006, [http://www.kaisernetwork.org/health_cast/hcast_index.cfm?display=detail&hc=1791]Google Scholar
- Okie S: Global health-the Gates-Buffett effect. N Engl J Med. 2006, 355 (11): 1084-1088. 10.1056/NEJMp068186.View ArticlePubMedGoogle Scholar
- Steinbrook R: Message from Toronto-deliver AIDS treatment and prevention. N Engl J Med. 2006, 355 (11): 1081-1084. 10.1056/NEJMp068167.View ArticlePubMedGoogle Scholar
- Shiffman J: Donor funding priorities for communicable disease control in the developing world. Health Policy Plan. 2006, 21 (6): 411-420. 10.1093/heapol/czl028.View ArticlePubMedGoogle Scholar
- McMahon R: Global Health and the Gates Foundation [Council on Foreign Relations.]. 2006, [http://www.cfr.org/publication/11004/global_health_and_the_gates_foundation.html]Google Scholar
- Rowley JT, Anderson RM, Ng TW: Reducing the spread of HIV infection in sub-Saharan Africa: some demographic and economic implications. AIDS. 1990, 4: 47-56. 10.1097/00002030-199001000-00007.View ArticlePubMedGoogle Scholar
- Stover J: Influence of mathematical modeling of HIV and AIDS on policies and programs in the developing world. Sex Transm Dis. 2000, 27 (10): 572-578. 10.1097/00007435-200011000-00005.View ArticlePubMedGoogle Scholar
- Kates J: Financing the response to HIV/AIDS in low and middle income countries: funding for HIV/AIDS from the G7 and the European Commission. 2005, [http://www.sarpn.org.za/documents/d0001507/index.php]Google Scholar
- Clark CF: AIDS and the Arrows of Pestilence. 1994, Golden, Colorado: Fulcrum PublishingGoogle Scholar
- Bayer R, Fairchild AL: Changing the paradigm for HIV testing-the end of exceptionalism. N Engl J Med. 2006, 355 (7): 647-649. 10.1056/NEJMp068153.View ArticlePubMedGoogle Scholar
- Schwartländer B, Stover J, Walker N, Bollinger L, Gutierrez JP, McGreevey W, Opuni M, Forsythe S, Kumaranayake L, Watts C, Bertozzi S: Resource needs for HIV/AIDS. Science. 2001, 292 (5526): 2434-2436. 10.1126/science.1062876.View ArticlePubMedGoogle Scholar
- Montaner JS, Hogg R, Wood E, Kerr T, Tyndall M, Levy AR, Harrigan PR: The case for expanding access to highly active antiretroviral therapy to curb the growth of the HIV epidemic. Lancet. 2006, 368 (9534): 531-536. 10.1016/S0140-6736(06)69162-9.View ArticlePubMedGoogle Scholar
- Moghadas SM, Gumel AB, McLeod RG, Gordon R: Could condoms stop the AIDS epidemic?. J Theor Med. 2004, 5 (3-4): 171-181. 10.1080/10273660412331315147.View ArticleGoogle Scholar
- Gordon R: A critical review of the physics and statistics of condoms and their role in individual versus societal survival of the AIDS epidemic. J Sex Marital Ther. 1989, 15: 5-30.View ArticlePubMedGoogle Scholar
- Dugger CW: U.S. jobs shape condoms' role in foreign aid, New York Times (October 29). 2006, [http://www.nytimes.com/2006/10/29/world/29condoms.html]Google Scholar
- UNAIDS: Joint United Nations Program on HIV/AIDS. 2004, [http://www.unaids.org/en]Google Scholar
- Cohen DA, Farley TA: Social marketing of condoms is great, but we need more free condoms. Lancet. 2004, 364 (9428): 13-14. 10.1016/S0140-6736(04)16611-7.View ArticlePubMedGoogle Scholar
- Cohen DA, Farley TA, Mason K, Ridgeway G: The collectivity of sexual behaviour. Int J STD AIDS. 2006, 17 (3): 151-156. 10.1258/095646206775809213.View ArticlePubMedGoogle Scholar
- Trucco E: Studies in imitative behavior: A generalization of the Rashevsky model; Its mathematical properties. Bull Math Biol. 1954, 16 (4): 279-316.Google Scholar
- United Nations General Assembly, Sixtieth session, Agenda item 45. Follow-up to the outcome of the twenty-sixth special session: implementation of the Declaration of Commitment on HIV/AIDS. Scaling up HIV prevention, treatment, care and support. 2006, [http://data.unaids.org/pub/InformationNote/2006/20060324_HLM_GA_A60737_en.pdf]
- Kremer M, Snyder CM: Why is There No AIDS Vaccine? [Working Paper]. 2004, Cambridge, MA: National Bureau of Economic ResearchGoogle Scholar
- Smith RJ, Bodine EN, Wilson DP, Blower SM: Evaluating the potential impact of vaginal microbicides to reduce the risk of acquiring HIV in female sex workers. AIDS. 2005, 19: 413-421. 10.1097/01.aids.0000161771.44276.92.View ArticlePubMedGoogle Scholar
- Hethcote HW, Van Ark JW: Modeling HIV Transmission and AIDS in the United States. Lect Notes in Biomath 95. 1992, Berlin-Heidelberg-New York: Springer-VerlagGoogle Scholar
- Huang W, Cooke KL, Castillo-Chavez C: Stability and bifurcation for a multiple group model for the dynamics of HIV transmission. SIAM J Appl Math. 1992, 52: 835-854. 10.1137/0152047.View ArticleGoogle Scholar
- Mode CJ, Sleeman CK: Stochastic Processes in Epidemiology. HIV/AIDS, Other Infectious Diseases and Computers. 2000, Singapore, New Jersey: World ScientificGoogle Scholar
- Korenromp EL, van Vliet C, Grosskurth H, Gavyole A, Ploeg van der CPB, Fransen L, Hayes RJ, Habbema JDF: Model-based evaluation of single-round mass treatment of sexually transmitted diseases for HIV control in a rural African population. AIDS. 2000, 14: 573-93. 10.1097/00002030-200003310-00013.View ArticlePubMedGoogle Scholar
- Blower SM, Aschenbach AN, Gershengorn HB, Kahn JO: Predicting the unpredictable: Transmission of drug-resistant HIV. Nat Med. 2001, 7: 1016-1020. 10.1038/nm0901-1016.View ArticlePubMedGoogle Scholar
- Boily MC, Bastos FI, Desai K, Mâsse B: Changes in the transmission dynamics of the HIV epidemic after the wide-scale use of antiretroviral therapy could explain increases in sexually transmitted infections-results from mathematical models. Sex Transm Dis. 2004, 31: 100-112. 10.1097/01.OLQ.0000112721.21285.A2.View ArticlePubMedGoogle Scholar
- Rapatski BL, Suppe F, Yorke JA: HIV epidemics driven by late disease stage transmission. J Acquir Immune Defic Syndr. 2005, 38: 241-253.PubMedGoogle Scholar
- Vickerman P, Terris-Prestholt F, Delany S, Kumaranayake L, Rees H, Watts C: Are targeted HIV prevention activities still cost-effective in high prevalence settings? Results from an STI treatment intervention for sex workers in Hillbrow, South Africa. Sex Transm Dis. 2006, 33 (10 Suppl): S122-32.View ArticlePubMedGoogle Scholar
- Coffee M, Lurie MN, Garnett GP: Modelling the impact of migration on the HIV epidemic in South Africa. AIDS. 2007, 21: 343-350. 10.1097/QAD.0b013e328011dac9.View ArticlePubMedGoogle Scholar
- Nagelkerke NJD, Moses S, de Vlas SJ, Bailey RC: Modelling the public health impact of male circumcision for HIV prevention in high prevalence areas in Africa. BMC Infect Dis. 2007, 7: 16-10.1186/1471-2334-7-16.PubMed CentralView ArticlePubMedGoogle Scholar
- Klein E, Laxminarayan R, Smith D, Gilligan C: Economic incentives and mathematical models of disease. Environment and Development Economics. 2007, 12: 707-732. 10.1017/S1355770X0700383X.View ArticleGoogle Scholar
- Garira Mukandavire: Effects of public health educational campaigns and the role of sex workers on the spread of HIV/AIDS among heterosexuals. Theor Pop Biol. 2007, 72 (3): 346-365. 10.1016/j.tpb.2007.07.002.View ArticleGoogle Scholar
- Xu Y, Xu GQ: Modeling of HIV/AIDS Transmission and Its Analysis. DCDIS Proceedings of the International Conference on Life System Modeling and Simulation. 2007, 488-492.Google Scholar
- Gran JM, Wasmuth L, Amundsen EJ, Lindqvist BH, Aalen OO: Growth rates in epidemic models: Application to a model for HIV/AIDS progression. Statistics in Medicine. 2008, 27 (23): 4817-4834. 10.1002/sim.3219.View ArticlePubMedGoogle Scholar
- Lloyd AL, May RM: Spatial Heterogeneity in Epidemic Models. J Theor Biol. 1996, 179: 1-11. 10.1006/jtbi.1996.0042.View ArticlePubMedGoogle Scholar
- Arino J, Davis JR, Hartley D, Jordan R, Miller JM, van den Driessche P: A multi-species epidemic model with spatial dynamics. Math Med Biol. 2005, 22 (2): 129-42. 10.1093/imammb/dqi003.View ArticlePubMedGoogle Scholar
- Lloyd AL, Jansen VAA: Spatiotemporal dynamics of epidemics: synchrony in metapopulation models. Math Biosci. 2004, 188 (1-2): 1-16. 10.1016/j.mbs.2003.09.003.View ArticlePubMedGoogle Scholar
- Liebovitch LS, Schwartz IB: Migration induced epidemics: dynamics of flux-based multipatch models. Physics Letters A. 2004, 332: 256-267. 10.1016/j.physleta.2004.09.071.View ArticleGoogle Scholar
- Arino J, van den Driessche P: A multi-city epidemic model. Math Pop Stud. 2003, 10 (3): 175-193. 10.1080/08898480306720.View ArticleGoogle Scholar
- Fitzgibbon WE, Langlais M, Morgan JJ: A mathematical model for indirectly transmitted diseases. Math Biosci. 2007, 206 (2): 233-248. 10.1016/j.mbs.2005.07.005.View ArticlePubMedGoogle Scholar
- Hsieh YH, van den Driessche P, Wang L: Impact of Travel Between Patches for Spatial Spread of Disease. Bull Math Biol. 2007, 69: 1355-1375. 10.1007/s11538-006-9169-6.View ArticlePubMedGoogle Scholar
- Jacquez JA, Simon CP, Koopman J, Sattenspiel L, Perry T: Modeling and analyzing HIV transmission: the effect of contact patterns. Math Biosci. 1988, 92: 119-199. 10.1016/0025-5564(88)90031-4.View ArticleGoogle Scholar
- Castillo-Chavez CC, Cooke K, Huang W, Levin SA: On the role of long incubation periods in the dynamics of AIDS I: single population models. J Math Biol. 1989, 27: 373-398. 10.1007/BF00290636.View ArticlePubMedGoogle Scholar
- Busenberg S, Cooke KL, Thieme HR: Demographic change and persistence of HIV/AIDS in a heterosexual population. SIAM J Appl Math. 1991, 51: 1030-1051. 10.1137/0151052.View ArticleGoogle Scholar
- Thieme HR, Castillo-Chavez CC: How May Infection-Age-Dependent Infectivity Affect the Dynamics of HIV/AIDS?. SIAM J App Math. 1993, 53 (5): 1447-1479. 10.1137/0153068.View ArticleGoogle Scholar
- Brauer F, van den Driessche P: Models for transmission of disease with immigration of infectives. Math Biosci. 2001, 171 (2): 143-154. 10.1016/S0025-5564(01)00057-8.View ArticlePubMedGoogle Scholar
- McCluskey CC: A model of HIV/AIDS with staged progression and amelioration. Math Biosci. 2003, 181: 1-16. 10.1016/S0025-5564(02)00149-9.View ArticlePubMedGoogle Scholar
- Gumel AB, McCluskey CC, van den Driessche P: Mathematical Study of a Staged-Progression HIV Model with Imperfect Vaccine. Bull Math Biol. 2006, 68: 2105-2128. 10.1007/s11538-006-9095-7.View ArticlePubMedGoogle Scholar
- Heffernan JM, Smith RJ, Wahl LM: Perspectives on the Basic Reproductive Ratio. J R Soc Interface. 2005, 2 (4): 281-293. 10.1098/rsif.2005.0042.PubMed CentralView ArticlePubMedGoogle Scholar
- van den Driessche P, Watmough J: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002, 180: 29-48. 10.1016/S0025-5564(02)00108-6.View ArticlePubMedGoogle Scholar
- Smith HL: Monotone Dynamical Systems. An Introduction to the Theory of Competitive and Cooperative Systems. Math Surveys Monogr. 1995, 41: 1-174. Amer. Math. Soc., Providence, R.I.View ArticleGoogle Scholar
- Smith HL, Waltman P: The Theory of the Chemostat. 1995, Cambrige, UK: Cambrige University PressView ArticleGoogle Scholar
- Central Intelligence Agency World Factbook. [https://www.cia.gov/library/publications/the-world-factbook]
- The Global Migrant Origin Database. [http://www.migrationdrc.org/research/typesofmigration/global_migrant_origin_database.html]
- Vijayaraghavan A, Efrusy MB, Mazonson PD, Ebrahim O, Sanne IM, Santas CC: Cost-effectiveness of alternative strategies for initiating and monitoring highly active antiretroviral therapy in the developing world. J Acquir Imme Defic Syndr. 2007, 46: 91-100.Google Scholar
- Gordon R, Tyson JE: Fatal sexually transmitted diseases (FSTDs), such as AIDS, select for the evolution of monogamy and provide a model for background extinction. J Biol Systems. 1993, 1: 1-26. 10.1142/S0218339093000021.View ArticleGoogle Scholar
- Burris S: Driving the epidemic underground? A new look at law and the social risk of HIV testing. AIDS Public Policy Journal. 1997, 12 (2): 66-78.PubMedGoogle Scholar
- Lopez AD, Mathers CD: Measuring the global burden of disease and epidemiological transitions: 2002-2030. Ann Trop Med Parasitol. 2006, 100 (5-6): 481-499. 10.1179/136485906X97417.View ArticlePubMedGoogle Scholar
- Abdullah AS, Ebrahim SH, Fielding R, Morisky DE: Sexually transmitted infections in travelers: implications for prevention and control. Clin Infect Dis. 2004, 39: 533-538. 10.1086/422721.View ArticlePubMedGoogle Scholar
- White ND, Hill DM, Bodemeier S: Male condoms that break in use do so mostly by a "blunt puncture" mechanism. Contraception. 2008, 77: 360-365. 10.1016/j.contraception.2008.01.014.View ArticlePubMedGoogle Scholar
- Barnighausen T, Bloom DE, Humair S: Human Resources for Treating HIV/AIDS: Needs, Capacities, and Gaps. AIDS Pat Care and STDs. 2007, 21 (11): 799-812. 10.1089/apc.2007.0193.View ArticleGoogle Scholar
- Balakrishnan P, Solomon S, Kumarasamy N, Mayer KH: Low-cost monitoring of HIV infected individuals on highly active antiretroviral therapy (HAART) in developing countries. Indian J Med Res. 2005, 121: 345-355.PubMedGoogle Scholar
- Yazdanpanah Y: Costs associated with combination therapy in HIV-infected patients. J Antimicrob Chemother. 2004, 53: 558-561. 10.1093/jac/dkh142.View ArticlePubMedGoogle Scholar
- Kumarasamy N, Solomon S, Chaguturu SK, Cecelia AJ, Vallabhaneni S, Flanigan TP, Mayer KH: The changing natural history of HIV disease: Before and after the introduction of generic antiretroviral therapy in southern India. Clin Infect Dis. 2005, 41: 1525-1528. 10.1086/497267.View ArticlePubMedGoogle Scholar
- Meng X, Anderson AF, Hou X, Wang Y, Sun L, Zhang X, Li Z, Qui B, Lang Y, Zhang L: A pilot project for the effective delivery of HAART in rural China. AIDS Pat Care and STDs. 2006, 20 (3): 213-219. 10.1089/apc.2006.20.213.View ArticleGoogle Scholar
- Dunne M: Antiretroviral drug development: the challenge of cost and access. AIDS. 2007, 21 (suppl 4): S73-S79. 10.1097/01.aids.0000279709.47298.8a.View ArticlePubMedGoogle Scholar
- Hounton SH, Akonde A, Zannou DM, Bashi J, Meda N, Newlands D: Costing universal access of highly active antiretroviral therapy in Benin. AIDS Care. 2008, 20 (5): 582-587. 10.1080/09540120701868303.View ArticlePubMedGoogle Scholar
- Beaston-Blaakman A, Shepard DS, Stone N, Shevitz H: Cost-effectiveness of clinical interventions for AIDS wasting. AIDS Care. 2007, 19 (8): 996-1001. 10.1080/09540120701335238.View ArticlePubMedGoogle Scholar
- Shevitz AH, Wilson IB, McDermott AY, Spiegelman D, Skinner SC, Antonsson K, Layne JE, Beaston-Blaakman A, Shepard DS, Gorbach SL: A comparison of the clinical and cost-effectiveness of 3 intervention strategies for AIDS wasting. J Acquir Immune Defic Syndr. 2005, 38 (4): 399-406. 10.1097/01.qai.0000152647.89008.2b.View ArticlePubMedGoogle Scholar
- Fieno JV: Costing adult male circumcision in high HIV prevalence, low circumcision rate countries. AIDS Care. 2008, 20 (5): 515-520. 10.1080/09540120701867149.View ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.