Secondary Logo

Journal Logo

Prevention Research

The Impact of Preexposure Prophylaxis Among Men Who Have Sex With Men

An Individual-Based Model

Kasaie, Parastu PhD*; Pennington, Jeff BS; Shah, Maunank S. MD, PhD; Berry, Stephen A. MD; German, Danielle PhD*; Flynn, Colin P. ScM§; Beyrer, Chris MD, PhD; Dowdy, David W. MD, PhD

Author Information
JAIDS Journal of Acquired Immune Deficiency Syndromes: June 1, 2017 - Volume 75 - Issue 2 - p 175-183
doi: 10.1097/QAI.0000000000001354

Abstract

BACKGROUND

Preexposure prophylaxis (PrEP) is increasingly recognized as an important component of any comprehensive strategy to end local HIV epidemics. Several clinical trials have established the effectiveness and safety of PrEP,1–3 and the United States Food and Drug Administration (FDA) has approved the provision of tenofovir/emtricitabine for prevention of HIV among high-risk populations.4 Since then, many cities have begun implementation of PrEP, using a variety of venues and strategies.5 Despite these promising developments, questions remain regarding the likely population-level impact of PrEP, as implemented in various settings and with varying degrees of targeting to populations at highest risk of HIV infection.2,6,7

In the absence of cluster-randomized trials (which would encounter ethical problems given the known efficacy of PrEP), mathematical models can be a useful way to project the potential impact of PrEP on HIV incidence and prevalence at the community level. In particular, agent-based simulations have the ability to incorporate several sources of individual-level heterogeneity, such as sexual mixing patterns, behavioral rules, and treatment adherence, as relevant to transmission and control of the HIV epidemic.8–12 In this work, we present an agent-based simulation model of PrEP delivery among men who have sex with men (MSM) to assess the population-level impact on HIV incidence and prevalence. Using Baltimore City as an illustrative case, we study various PrEP strategies targeting subpopulations at high-risk of HIV infection and compare the epidemiological impact of these strategies under various assumptions regarding the population-level coverage of PrEP and the degree of protection afforded by PrEP at the individual level (for example, reflecting different levels of adherence).

METHODS

Our agent-based simulation model of the HIV epidemic among MSM in Baltimore City is structured as a collection of different modules that govern various aspects of population demographics, partnerships, HIV natural history, and the cascade of care (Fig. 1). Each “agent” represents a single MSM, characterized by age, race, and place of residence. The simulation modules are programmed at the individual agent level but calibrated against available population-level epidemiological data. A natural history module characterizes the progression of HIV among infected individuals according to disease stage (acute, early, and late). Each stage is associated with a different per-act risk of HIV transmission, and progression from early to late disease can be prevented (and/or reversed) by provision of antiretroviral therapy (ART). A continuum of care module estimates individuals' probabilities of HIV testing, linkage to care, disengagement/reengagement, and ART provision/viral suppression. A sexual network and transmission module creates and modifies each individual's sexual network as a series of stable (long-term) and casual (short-term) partnerships, modeling HIV transmission as a per-act probability among serodiscordant partners according to frequency of unprotected sex acts, HIV stage of the infected partner, and ART/PrEP use. The probability of engaging in new partnerships reflects an individual's age, level of sexual activity, and current partnership status. New partnerships are modeled as assortative according to age, race, and location of residence [within Community Statistical Areas (CSAs), and within contiguous “CSA groups” that are similar according to race and socioeconomic status]. Finally, the demographic module accounts for aging, death, and birth. The model is updated using a time step of 1 week and followed for 5 years. Further details regarding the structure of simulation modules are provided in section 1 of the Appendices http://links.lww.com/QAI/A989.

FIGURE 1
FIGURE 1:
Overview of simulation modules. This figure illustrates the schematic simulation logic for modeling the demographic and partnership networks among MSM (panel A), as well as HIV natural history and the cascade of care (panel B). Panel A: The population is structured into a collection of Community Statistical Areas (CSAs) and CSA neighborhood groups, which are in turn based on geographical proximity, and the level of similarity in income and racial structure (represented schematically by same colors). Partnership domains are determined via discrete choice within an individual's own CSA of residence, a random neighboring CSA, or a non-neighbor CSA. Once the partnership domain is established, individuals follow a search mechanism based on a combination of race- and age-dependent mixing patterns to select their future partner from the pool of eligible people in that domain. Panel B: HIV infection is modeled as a gradual decline in CD4 count and via 3 main states in the absence of HIV treatment (acute/chronic/late infection). The cascade of care models the processes for screening infected individuals, linking to care, retaining in care and starting ART.

Parameterization/Calibration

Simulation parameters are categorized as either fixed (known) or variable (unknown). Fixed parameters are those for which strong individual-level evidence from the literature exists (Table 1); these parameters are set to specific values throughout the modeling process (though varied in sensitivity analysis). Variable parameters are those for which an individual-level estimate is unknown or unavailable. These parameters are calibrated in an iterative stepwise procedure designed to create a simulated population that matches selected population-level targets. Using Baltimore City as an illustrative case, demographic parameters (such as the birth rates within CSAs) are calibrated to replicate the average MSM population size and demographic structure available from census data13; partnership parameters (such as the individual likelihood of forming new partnerships) are tuned to replicate the reported frequency and type of sexual partnerships available from the National HIV Behavioral Surveillance study in Baltimore (BESURE)14; the parameters used for modeling HIV care (such as probability of HIV testing) are tuned to replicate the proportion of the infected population estimated to be aware of their HIV status, linked and virally suppressed in this population15; and the per-week probability of HIV transmission over each active contact was tuned to replicate the reported level of HIV prevalence among Baltimore City's MSM in year 2012.15 For more details regarding the calibration procedure and a list of model parameters and values see section 2 of the Appendices http://links.lww.com/QAI/A989.

TABLE 1
TABLE 1:
List of Fixed Model Parameters*

Simulated Population and Experiments

Baltimore City is a moderate-sized city (2015 population 620,96113) whose HIV epidemic encapsulates some of the dynamics of HIV transmission seen in many US cities, including heterogeneity in transmission, disparity in prevalence by age and race, and a relatively weak continuum of care. In 2012, Baltimore City had 469 reported new HIV diagnoses (89.9 infections per 100,000), 85.7% of which occurred among non-Hispanic blacks (referred to as Black) and 55.6% in MSM.15 Among young (<30) MSM in Baltimore in 1996–2000, estimated annual HIV incidence was 11.0% (95% CI: 5.5% to 19.7%) in Blacks, versus 0.6% in non-Hispanic whites.16 Rates by ZIP code likewise vary by a factor of 10.15 Only 31% of all PLHIV in Baltimore are estimated to have a suppressed viral load, with prominent disparities by age, race, and HIV risk factors.15 Using the current estimate of Baltimore City's male population who are 15 years or older in age and the estimated percentage of adult MSM in each racial group (7.5% of white males and 5.8% of black males17), we estimate the size of Baltimore's MSM population at approximately ∼15,000 men in Baltimore City (See section 1 of the Appendices http://links.lww.com/QAI/A989).

Corresponding to known risk factors for increased risk of HIV infection in this population, we further consider a set of hypothetical scenarios for annual provision of PrEP to eligible MSM, where eligibility criteria are based on the frequency of reported partnerships, age, and race as follows:

We define emblematic scenarios of PrEP administration based on eligibility criteria as follows:

  • S1: Universal administration of PrEP: All HIV-negative MSM are eligible for PrEP.
  • S2: Focused administration for MSM with >1 partner: All HIV-negative MSM with more than 1 partner in the last year are eligible.
  • S3: Focused administration for MSM with >5 partners: All HIV-negative MSM with more than 5 partners in the last year are eligible.
  • S4: Focused administration for young MSM: All HIV-negative MSM younger than 30 years old are eligible.
  • S5: Focused administration for young Black MSM: All HIV-negative Black MSM younger than 30 years old are eligible.

Although these scenarios do not precisely reflect the CDC guidelines for implementation of PrEP, they may more closely describe the likely implementation of these guidelines in actual clinical or public health practice (S2). Scenarios S3–S5 were designed using categories in available data to target sufficiently large subpopulations at increased risk of HIV infection (eg, those with a large number of partnerships18 or young age19).

In each scenario, we assume that a proportion of eligible individuals (population coverage) receive PrEP at the beginning of each year and continue to take PrEP at a level of adherence capable of providing protection against HIV infection on a specified percentage of days (protection), for at least 1 year. At the end of each year, those currently on PrEP who are no longer eligible to receive PrEP (eg, because of new HIV infection or reduction in condomless sex) are discontinued. Those discontinuing PrEP can restart PrEP in subsequent years, if still eligible. We model PrEP delivery for 5 consecutive years (year 1–5), and map the percent reduction in HIV incidence and prevalence compared with a baseline scenario without PrEP, at various levels of individual-level HIV protection and population coverage. Given the variation in the size of the PrEP-eligible population under each scenario (for example, about 11,000 MSM are eligible for PrEP in scenario 1, versus 3300 in scenario 5), we control the total number of people on treatment during each year and use this number as an alternative measure of PrEP coverage across all scenarios.

Simulation Procedure

Each simulation run was first carried out over a long “burn-in” period of 200 years until equilibrium was reached (baseline), targeting Baltimore's estimate of 3300 MSM living with HIV in 2014.15 At baseline, we then carried each simulation forward for an additional 5 years under each PrEP scenario described above, repeating the analysis for different levels of individual-level HIV protection provided by PrEP [range (0%–100%)] and population coverage [in terms of the number of people receiving PrEP, as well as the proportion of the eligible population receiving PrEP, range (0%–100%)] in each PrEP scenario. We estimated the required number of independent simulations to provide precision of at most 1% around the main output (HIV prevalence) and repeated the analysis across all (∼100) independent models. For more information, see section 2.1 of the Appendices http://links.lww.com/QAI/A989.

Sensitivity Analysis

One-way sensitivity analysis was performed by varying all simulation parameters (fixed and variable) to ±25% of their original values. Since small variations in parameters governing transmission had a disproportionate impact on HIV incidence, we evaluated those parameters in a way that changed baseline HIV incidence by <±25%. In each sensitivity analysis, we compared the annual impact of PrEP on HIV incidence after parameter variation relative to its impact at baseline (original parameter setting), using scenario 1 and maintaining 2000 randomly–selected-MSM with more than 1 partner on PrEP each year, assuming adherence sufficient to provide protection on 60% of days.

RESULTS

At baseline, the median age in the simulated population was 39.04 (95% uncertainty range 39–39.08) years, and 54.4% (54.3%–54.5%) of MSM were Black. This is the 95% uncertainty range around the mean of 100 simulations (see section 2.2.1 of the appendices http://links.lww.com/QAI/A989 for a description of methods). Figure 2 provides a summary of the model's primary calibration targets, illustrating that the model population is similar to that of available data. In any given year, 22.5% (22.4%–22.6%) of MSM in our model are only engaged in stable partnerships, whereas 40.9% (40.8%–41%) have both stable and casual partnerships. The simulated median number of annual partners, by age group, was 5.2 (ages 15–24), 4.2 (ages 25–34), 3.6 (ages 35–44) and 3.0 (ages 45 and above).

FIGURE 2
FIGURE 2:
Model calibration to epidemiological data. Shown are the mean values of 200 simulations (in green) compared against empirical data (in orange). The black bars represent the range of observations around each simulated measure. The data bars represent available point estimates on each measure including measures of partnership frequency (panels A and B) reported through unpublished data from the latest survey of Baltimore's MSM in 2014,14 as well as available estimates of HIV prevalence, racial disparities, and the HIV continuum of care among Baltimore's MSM in year 2012 (Fig. 2, panels C through E).15

After calibrating to an HIV prevalence of 22% among MSM15 (Fig. 2, panel C), our model projects an HIV incidence of 1390 cases per 100,000 person years (1374–1406) in the absence of PrEP. Young Black MSM between ages of 15–29 account for 18.5% (17.9%–19.1%) of the population but 36% (35.5%–36.6%) of incident HIV infections. Black MSM account for 80% (75%–85%) of HIV incidence in this model, and they have 30%–50% lower levels of being diagnosed, in care on ART, and virally suppressed, relative to white MSM. For more information, see section 3.1 in the Appendices http://links.lww.com/QAI/A989.

Impact of PrEP as a Function of Coverage and Individual-Level Protection

The impact of PrEP largely reflected the number of individuals in the population who were protected at any given time and was thus highly sensitive to both population-level coverage and the proportion of days protected at the individual level (a proxy for adherence). At 40% coverage and 40% individual-level protection, community-wide provision of PrEP to individuals with more than 1 partner in the last 12 months (scenario 2) is estimated to reduce HIV incidence by 9.5% (8%–11%) after 1 year (Fig. 3, panel A, bottom center), and by 10.6% (9%–12%) after 5 years (Fig. 3 panel B, bottom center). If, however, coverage and individual-level protection could both be increased to 80% (for example, through a mass community campaign using an injectable formulation), 5-year impact on HIV transmission could be increased to 43% (42%–44%). The impact of PrEP on HIV incidence is achieved rapidly on PrEP scale-up (Fig. 3, left panel) and lost with similar speed when PrEP protection ends, as PrEP has little effect on HIV prevalence over 5 years (Fig. 3, right panel). The impact of improving coverage (as a proportion of the eligible population covered) is similar to that of improving individual-level HIV protection. For example, starting from a baseline of 40% coverage and 40% protection, achieving a 2-fold increase in coverage would reduce incidence by 19.4% (18%–21%), whereas a 2-fold increase in individual-level protection would reduce incidence by 24% (23%–25%). Simultaneous reductions in both coverage and individual-level HIV protection generate approximately additive reductions in population-level impact (Fig. 3, left and center panels). For more information, see section 3.2 in the Appendices http://links.lww.com/QAI/A989.

FIGURE 3
FIGURE 3:
Percent reduction in HIV incidence after delivering PrEP to MSM reporting more than one partner in a year (scenario 2) at different individual-level HIV protection and population-level coverage. Contour lines show the percent reduction in HIV incidence at various levels of population coverage (on the x-axis, with the percent of the eligible population covered on top and the numbers of MSM on PrEP in parentheses) and individual-level protection (on the y-axis). Panel A shows impact after 1 year of PrEP delivery, Panel B shows the same impact after sustaining the same level of PrEP for 5 years, and Panel C shows the impact 1 year after stopping PrEP (at the end of year 5). All impact estimates are shown relative to a baseline of no PrEP delivery.

Comparing the Impact of PrEP Targeted at Various Risk Groups

Given that some scenarios have large eligible pools with relatively less targeting (eg, scenario 1), whereas others have smaller but more targeted eligible pools (eg, scenario 5), we compared the scenarios as a function of the number of people on PrEP in each year and individual-level protection (Fig. 4). Per person prescribed PrEP, targeting PrEP to a demographic group at high risk for HIV infection (young Black MSM, scenario 5) results in the greatest reduction in HIV incidence. For example, under the assumption of full individual-level protection, providing PrEP to 1500 young Black MSM has equivalent projected impact [27% (25.5%–28%) reduction in incidence] to covering 3000 randomly selected eligible MSM in the general population [30% (29%–31%) reduction in incidence]. Although targeting PrEP to those at higher sexual risk (>1 partner in the previous 12 months: scenario 2, green line in Fig. 4) modestly improved impact over untargeted delivery (scenario 1, black line), additional targeting to those at very high sexual risk (>5 partners per year: scenario 3, blue line) or young age alone (scenario 4, yellow line) did not further increase PrEP impact. For example, maintaining 2000 randomly selected uninfected MSM on PrEP was projected to reduce HIV incidence at the end of 5 years by 20% (18.3%–21.4%) in scenario 1 (middle of x-axis, Fig. 4). A similar campaign targeting MSM with >1 annual partner augmented PrEP impact 1.2-fold, versus 1.3-fold if targeted to those with >5 partners or to those below age 30. By contrast, a campaign delivering PrEP to 2000 young Black MSM (scenario 5, red line) achieved 1.8 times greater reduction in community-wide HIV incidence compared with the baseline scenario. Although the absolute impact of targeted delivery declined proportionally with decreasing individual-level protection, the relative impact of targeting was preserved, regardless of individual-level protection levels.

FIGURE 4
FIGURE 4:
Comparative impact of targeted PrEP scenarios on HIV incidence at various level of HIV protection. Each panel maps the comparative impact of PrEP scenarios on HIV incidence at a different level of achievable protection under PrEP, ranging from 20% (top left) to 100% (bottom left). The y-axis of each panel shows the percent reduction in incidence under each PrEP scenario at the end of the fifth year of implementation, relative to a baseline with no PrEP. The x-axis represents the number of MSM receiving PrEP in each of 5 scenarios (described in the inset and the manuscript text). Shaded areas represent the 95% uncertainty range around the simulated means. The end of each line represents the total eligible population in each scenario except for scenario 1 (eg, Scenario 3 is continued until the estimated 2000 MSM in Baltimore with more than 5 partners are all covered). For a given number of MSM on PrEP, all targeted scenarios (scenarios 2–4) show modest improvements over random selection (scenario 1), and scenario 5 shows the greatest impact of all strategies evaluated.

In contrast to the rapid impact of PrEP on HIV incidence that quickly plateaus (as shown in Fig. 3), impact on HIV prevalence is slower but builds over time. For example, delivering PrEP to MSM with more than 1 partner in a year (scenario 2) was projected to reduce HIV prevalence by only 15% (12%–17%) after 5 years, but if this program could be sustained over 20 years, the corresponding reduction in HIV prevalence was projected to reach 54% (51%–56%). For more information, see section 3.3 in the Appendices http://links.lww.com/QAI/A989.

Sensitivity Analysis

The relative impact of PrEP was most sensitive to increased frequency of stable and casual partnerships, HIV viral load during the acute and chronic states, and overall probability of HIV transmission on contact (Appendix Figure 16 http://links.lww.com/QAI/A989), but there was no parameter for which a 25% variation (or variation sufficient to change HIV incidence by 25%) resulted in a sustainable change (>25%) in the projected impact of PrEP after 10 years. The results suggest that the efficacy of PrEP for prevention of HIV is likely to be similar in a variety of other epidemiological settings. A complete description of sensitivity analysis and the results are provided in section 4 of the Appendices http://links.lww.com/QAI/A989.

DISCUSSION

This agent-based simulation of HIV transmission among MSM illustrates that the impact of PrEP depends both on population-level coverage and individual-level protection (as a proxy for adherence). If delivered with exceptionally high population-level coverage and individual-level protection, PrEP can have a major and immediate impact on the rate of HIV incidence in this population. For example, we estimated that a delivery strategy that covers 80% of the HIV-negative MSM population with >1 sex partner in the past year, taken with sufficient adherence to provide protection on 80% of days, can reduce city-wide HIV incidence among MSM by almost half within 1 year. However, under more realistic assumptions regarding coverage and adherence, the impact of PrEP alone is much more modest—emphasizing the importance of including PrEP with other interventions if local HIV epidemics are to be brought under control. The impact of PrEP can be augmented (on a per-person basis) by targeting demographic groups at high risk of HIV infection (for example, young Black MSM in Baltimore). The scope of this impact, however, is limited to the size of the high-risk population in question and the ability to provide PrEP to that population with high coverage and adherence using available resources. Providing PrEP to populations with lower risk (but still eligible) does further enhance impact on incidence—and may be logistically more feasible (for example, providing PrEP to eligible individuals already connected to healthcare). Decisions about targeting PrEP delivery must, therefore, balance expected increases in efficiency from targeting demographic groups at high transmission risk, achievable levels of uptake and adherence in these groups, and logistical/resource considerations that may favor other strategies.

Our results are consistent with those of other studies of the population-level impact of PrEP. For example, these results agree with the argument of van der Straten et al that low population uptake and suboptimal adherence can undermine the impact of PrEP at the population level.20 Other modeling studies have forecasted the possibility of substantial reductions in HIV incidence in the United States, South Africa, and India21–24 and have suggested higher cost-effectiveness ratios for programs targeting younger and higher-risk populations.21,22 Our results extend these earlier findings by incorporating heterogeneities in the transmission and control of HIV that can serve as the basis for targeted PrEP delivery. As a result, we are able to quantify, for example, the additional impact that might be achieved by focusing PrEP interventions among young Black MSM in an urban US city characterized by racial disparities in HIV. Another modeling study incorporating CDC guidelines for implementation of PrEP among MSM in the United States estimated that a PrEP campaign with coverage of 40% and adherence of 62% could reduce HIV incidence by 33% over a decade,25 almost identical to our corresponding scenario 1 estimate of a 32% (Uncertainty Range (UR): 12.2%–44.8%) reduction under similar coverage and adherence.

Our results have a number of practical implications for the implementation of PrEP programs at the local level. First, if PrEP is taken with relatively low adherence (for example, sufficient to provide protection against HIV on only 20%–40% of days), very high levels of population coverage must be achieved to attain substantive reductions in HIV incidence. Similarly, if population coverage is capped at a given threshold (because of drug supply or inability to reach certain subpopulations), then adherence must be exceptionally high. Our results suggest that at low level of coverage and adherence, a small number of settings may still experience an increase in HIV incidence after implementation of PrEP because of random variation and the low number of HIV infections that occur on a population level in any given year (ie, 95% uncertainty range being less than zero). This model also suggests that targeting of race/age-based demographic groups at high risk for HIV infection may substantially increase the efficiency of PrEP delivery and may—at least in this simulated population—have greater impact than targeting groups with high sexual frequency or low age alone. These findings highlight that HIV transmission largely occurs in specific assortative sexual networks (which may correlate with certain demographic characteristics such as race), rather than simply occurring among those who have the largest number of sex partners.

As with any modeling study, our results are limited by certain simplifying assumptions used in design and calibration of the simulation model. In the absence of spatial sexual-network data, we adopted a simplified approach in modeling geographical mixing based on location of residence. The simplified definition of stable/casual contacts in this model further limits our ability to represent MSM with very high rates of partner turnover such as in sex work or group sex settings. Additional simplifying assumptions that may also impact our ability to precisely project HIV transmission dynamics include homogeneous sexual position preference, frequency of sexual acts in each type of partnership, and the likelihood of condom use (conditional on partnership type). For purposes of maintaining an appropriately simple model structure, we limited our analysis to MSM; future studies could extend this work to other population groups including high-risk heterosexuals, injection drug users and female sex workers. Additional efforts could also seek to characterize migration patterns across cities, behavioral changes while on PrEP, differences in PrEP uptake by race and age, interplay between PrEP and ART, and the role of other preventive interventions (eg, condoms and male circumcision) in modulating the impact of PrEP.

In summary, the advent of PrEP presents a major opportunity to make substantial progress toward ending HIV epidemics at the local level. This simulation shows that, although the projected epidemic impact of PrEP at currently feasible levels of implementation is modest, PrEP can rapidly and substantially reduce HIV transmission if ambitious targets for coverage and individual-level protection can be achieved through additional innovation in delivery and/or formulation. This impact is largely proportional to the number of individuals at high risk for transmission who successfully take PrEP, and it can be enhanced by targeting demographic groups at high risk for HIV infection. Although our results are based on a single setting (Baltimore), we feel they may be generalizable to epidemics among MSM in many settings. To maximize the impact of PrEP as an individual- and population-level HIV prevention strategy, it is ultimately necessary to consider PrEP as part of a comprehensive package that also includes ART for HIV-positive individuals and other behavioral intervention options. Nonetheless, our model illustrates that this novel biomedical intervention can play an important role in altering the face of the HIV epidemic in high-risk populations at the city level.

REFERENCES

1. Grant RM, Lama JR, Anderson PL, et al. Preexposure chemoprophylaxis for HIV prevention in men who have sex with men. N Engl J Med. 2010;363:2587–2599.
2. Baeten JM, Donnell D, Ndase P, et al. Antiretroviral prophylaxis for HIV prevention in heterosexual men and women. N Engl J Med. 2012;367:399–410.
3. Thigpen MC, Kebaabetswe PM, Paxton LA, et al. Antiretroviral preexposure prophylaxis for heterosexual HIV transmission in Botswana. N Engl J Med. 2012;367:423–434.
4. U.S. Food and Drug Administration. FDA Approves First Medication to Reduce HIV Risk. Silver Spring, MD: U.S. Food and Drug Administration; 2016.
5. Mayer KH, Hosek S, Cohen S, et al. Antiretroviral pre-exposure prophylaxis implementation in the United States: a work in progress. J Int AIDS Soc. 2015;18:54. doi:.
6. Underhill K, Morrow KM, Colleran CM, et al. Access to healthcare, HIV/STI testing, and preferred pre-exposure prophylaxis providers among men who have sex with men and men who engage in street-based sex work in the US. PLoS One. 2014;9:e112425.
7. Abdool Karim SS, Baxter C. Antiretroviral prophylaxis for the prevention of HIV infection: future implementation challenges. HIV Ther. 2009;3:3–6.
8. Marshall BD, Paczkowski MM, Seemann L, et al. A complex systems approach to evaluate HIV prevention in metropolitan areas: preliminary implications for combination intervention strategies. PLoS One. 2012;7:e44833.
9. Beyrer C, Baral SD, van Griensven F, et al. Global epidemiology of HIV infection in men who have sex with men. Lancet. 2012;380:367–377.
10. Alam SJ, Meyer R, Norling E. A model for HIV Spread in a South African Village. Berlin, Germany: Springer Berlin Heidelberg; 2009:33–45.
11. Goodreau SM, Carnegie NB, Vittinghoff E, et al. Can male circumcision have an impact on the HIV epidemic in men who have sex with men? PLoS One. 2014;9:e102960.
12. Hallett TB, Gregson S, Dube S, et al. The impact of monitoring HIV patients prior to treatment in resource-poor settings: insights from mathematical modelling. PLoS Med. 2008;5:e53.
13. Alliance BNI. Vital Signs 12. 2012. Available at: http://bniajfi.org/vital_signs/. Accessed January 14, 2016.
14. Maryland Department of Health and Mental Hygiene. BESURE Study 2004–2014: Baltimore Site of National HIV Behavioral Surveillance (NHBS). Baltimore, MD: Maryland Department of Health and Mental Hygiene; 2015.
15. Maryland Department of Health and Mental Hygiene, MDHMH. 2012 Baltimore City Annual HIV Epidemiological Profile. Baltimore, MD; 2012.
16. Sifakis F, Hylton JB, Flynn C, et al. Racial disparities in HIV incidence among young men who have sex with men: the Baltimore young men's survey. J Acquir Immune Defic Syndr. 2007;46:343–348.
17. Lieb S, Fallon SJ, Friedman SR, et al. Statewide estimation of racial/ethnic populations of men who have sex with men in the U.S. Public Health Rep. 2011;126:60–72.
18. Maulsby C, Sifakis F, German D, et al. Partner characteristics and undiagnosed HIV seropositivity among men who have sex with men only (MSMO) and men who have sex with men and women (MSMW) in Baltimore. AIDS Behav. 2012;16:543–553.
19. Maulsby C, Jain K, Sifakis F, et al. Individual-level and partner-level predictors of newly diagnosed HIV infection among black and white men who have sex with men in Baltimore, MD. AIDS Behav. 2015;19:907–917.
20. van der Straten A, Van Damme L, Haberer JE, et al. Unraveling the divergent results of pre-exposure prophylaxis trials for HIV prevention. AIDS. 2012;26:F13–F19.
21. Paltiel AD, Freedberg KA, Scott CA, et al. HIV preexposure prophylaxis in the United States: impact on lifetime infection risk, clinical outcomes, and cost-effectiveness. Clin Infect Dis. 2009;48:806–815.
22. Abbas UL, Anderson RM, Mellors JW. Potential impact of antiretroviral chemoprophylaxis on HIV-1 transmission in resource-limited settings. PLoS One. 2007;2:e875.
23. Vissers DCJ, Voeten HACM, Nagelkerke NJD, et al. The impact of pre-exposure prophylaxis (PrEP) on HIV epidemics in Africa and India: a simulation study. PLoS One. 2008;3:e2077.
24. Desai K, Sansom SL, Ackers ML, et al. Modeling the impact of HIV chemoprophylaxis strategies among men who have sex with men in the United States: HIV infections prevented and cost-effectiveness. AIDS. 2008;22:1829–1839.
25. Jenness S, Goodreau S. Impact of the centers for disease control's HIV preexposure prophylaxis guidelines for men who have sex with men in the United States. J Infect Dis. 2016;214:1800–1807.
26. Baggaley RF, Ferguson NM, Garnett GP. The epidemiological impact of antiretroviral use predicted by mathematical models: a review. Emerg Themes Epidemiol. 2005;2:9.
    27. Salomon JA, Hogan DR, Stover J, et al. Integrating HIV prevention and treatment: from slogans to impact. PLoS Med. 2005;2:e16.
      28. Currie S, Rogstad KE, Piyadigamage A, et al. Time taken to undetectable viral load, following the initiation of HAART. Int J STD AIDS. 2009;20:265–266.
      29. Wit FW, Blanckenberg DH, Brinkman K, et al. Safety of long-term interruption of successful antiretroviral therapy: the ATHENA cohort study. AIDS. 2005;19:345–348.
      30. Maggiolo F, Ripamonti D, Gregis G, et al. Effect of prolonged discontinuation of successful antiretroviral therapy on CD4 T cells: a controlled, prospective trial. AIDS. 2004;18:439–446.
      31. Ortiz GM, Wellons M, Brancato J, et al. Structured antiretroviral treatment interruptions in chronically HIV-1-infected subjects. Proc Natl Acad Sci. 2001;98:3288–13293.
        32. El-Sadr WM, Lundgren JD, Neaton JD, et al. CD4+ count-guided interruption of antiretroviral treatment. New Engl J Med. 2006;355:2283–2296.
        33. The INSIGHT. START Study Group. Initiation of antiretroviral therapy in early asymptomatic HIV infection. N Engl J Med. 2015;373:795–807.
        34. Kitahata MM, Gange SJ, Abraham AG, et al. Effect of early versus deferred antiretroviral therapy for HIV on survival. N Engl J Med. 2009;360:1815–1826.
        35. Palella FJ Jr, Delaney KM, Moorman AC, et al. Declining morbidity and mortality among patients with advanced human immunodeficiency virus infection. N Engl J Med. 1998;338:853–860.
        36. Rana AI, Liu T, Gillani FS, et al. Multiple gaps in care common among newly diagnosed HIV patients. AIDS Care. 2015;27:687–697.
          Keywords:

          HIV; pre-exposure prophylaxis; prevention; computer simulation; homosexuality; male

          Supplemental Digital Content

          Copyright © 2017 Wolters Kluwer Health, Inc. All rights reserved.