A Decision Analytics Model to Optimize Investment in Interventions Targeting the PrEP Cascade of Care ===================================================================================================== * Samuel M. Jenness * Gregory Knowlton * Dawn K. Smith * Julia L. Marcus * Emeli J. Anderson * Aaron J. Siegler * Jeb Jones * Patrick S. Sullivan * Eva Enns ## ABSTRACT **Background** Gaps between recommended and actual levels of HIV preexposure prophylaxis (PrEP) use remain among men who have sex with men (MSM) in the United States. Interventions can address these gaps within a PrEP cascade framework. However, it is unknown how public health initiatives like the *Ending the HIV Epidemic* (EHE) plan should invest prevention funds into these interventions to maximize their population impact. **Methods** We used a stochastic network-based HIV transmission model for MSM in the Atlanta area paired with an economic budget optimization model. The model simulated MSM particpating in up to three real-world PrEP cascade interventions designed to improve initiation, adherence, or persistence. The primary outcome was infections averted over 10 years. The budget optimization model identified the investment combination under different budgets that maximized this outcome given intervention costs from a payer perspective. **Results** From the base 15% PrEP coverage level, the three interventions could increase coverage to 27%, resulting in 12.3% of infections averted over 10 years. Uptake of each intervention was interdependent: maximal use of the adherence and persistence interventions depended on new PrEP users generated by the initiation intervention. At lower budgets, investment entirely in the persistence intervention was optimal. As the budget increased, optimal investment strategy involved a mixture of the initiation and persistence interventions, but not the adherence intervention. However, if adherence intervention costs were halved, the optimal investment distribution was roughly equal across interventions at the highest budget level. **Conclusions** Investments into the PrEP cascade through initiatives like the EHE plan should account for the interactions of the interventions as they are collectively deployed. Optimal funding depends on measuring service capacity against user demand. Given current intervention efficacy estimates, the total population impact of each intervention may be improved with greater total budgets or reduced intervention costs. KEYWORDS * Mathematical model * Sexual networks * HIV * Preexposure prophylaxis * Men who have sex with men * Ending the HIV Epidemic * Budget optimization * Economic model ## INTRODUCTION Preexposure prophylaxis (PrEP) is highly effective at HIV prevention, but large gaps exist between actual and recommended levels of use [1]. Men who have sex with men (MSM) are a priority population for PrEP in the United States [2], but only 5–20% of MSM with indications currently use it [3]. The new *Ending the Epidemic* (EHE) plan seeks to reduce HIV incidence by 90% by 2030 [4]. Key to the EHE plan is the “Prevent” pillar, which includes efforts to close the gaps between recommended and actual levels of PrEP use among MSM in high-burden areas like Atlanta. Opportunities to improve PrEP care have been characterized using a prevention cascade framework [5]. Steps in the PrEP cascade can be classified into three categories: initiation (including awareness, access, prescription), adherence, and persistence (i.e., retention in PrEP care). Several studies have been designed to evaluate behavioral and clinical interventions that target these steps in the PrEP cascade for MSM [6–8]. Achieving the EHE goals with PrEP will likely require a multi-pronged approach that improves all steps in the PrEP cascade. Understanding how combinations of initiation, adherence, and persistence interventions scale and interact is needed to direct EHE and local investment priorities. Transmission modeling of the PrEP cascade suggests that improving initiation and persistence are needed to further advance HIV prevention [9,10]. However, modeling studies to date have represented cascade improvements non-mechanistically: models project what happens when rates of PrEP initiation increase and rates of discontinuation decrease, but not the mechanisms by which those changes could actually occur. These modifications to the PrEP cascade may require investments in activities that address the challenges of attaining maximal engagement in PrEP care. These activities involve health system change (e.g., telemedicine), and targeted interventions using both traditional (e.g., counseling) and novel (e.g., mobile apps) approaches. Because these activities require financial resources, health authorities must consider how to invest EHE and other public health funding to balance the epidemiological benefits of these interventions against their costs. In this study, we used a decision analytic modeling approach that paired a stochastic network-based model of HIV transmission dynamics among MSM with a budget optimization model. These models simulated HIV transmission under varying levels of engagement in three categories of interventions (initiation, adherence, persistence) reflecting three real-world tools with unique impacts on the PrEP cascade. Our research objective was to determine the optimal investment of a fixed budget across these three categories of interventions to maximize the expected number of infections averted over 10 years. ## METHODS ### Study Design Our network-based model of HIV transmission dynamics was built with the EpiModel platform [11]. Building on our previous PrEP modeling applications [12], this study implemented the mechanisms for increasing PrEP coverage through engagement in specific interventions, and then integrated a budget optimization analysis. Full methodological details are provided in an Appendix [**LINK**]. The model represented main, casual, and one-time sexual partnerships for Black, Hispanic, and White/Other MSM, aged 15 to 65, in the Atlanta area. The starting network size in the model was 10000 MSM, which stochastically fluctuated over time with arrival (sexual debut) and departure (mortality or sexual cessation) (**Appendix S5**). ### HIV Transmission and Progression The epidemic model consisted of five main components: 1) statistical network models to generate dynamic sexual partnerships; 2) statistical models to predict behavior within partnerships; 3) simulation of HIV transmission in partnerships; 4) simulation of HIV progression and other natural history features; and 5) simulation of prevention and treatment service engagement. To fit the network models, we used data from ARTnet, an egocentric network study conducted in 2017–2019 of US MSM [13] (**Appendix S2**). Parameters were weighted by census-based race/ethnicity and age distributions to account for sampling biases. Predictors of partnership formation included partnership type, degree distributions by partnership types, heterogeneity in degree and assortative mixing by race/ethnicity, age, and mixing by sexual position (**Appendix S3**). Partnership durations were modeled with dissolution rates stratified by partnership type and age. Models were then fit to ARTnet to predict coital frequency and condom use probability as a function of race/ethnicity, age, diagnosed HIV status, and partnership type and duration (**Appendix S4**). In each weekly time step, MSM could be screened for HIV. MSM who screened positive entered the HIV care cascade (**Appendix S7**). MSM linked to care and initiating antiretroviral therapy (ART) had reduced vial load (VL) and increased longevity. MSM on ART could cycle off and back on ART based on rates calibrated to local surveillance data [14]. Lower VL with sustained ART use was associated with a reduced probability of HIV transmission per act. Other factors modifying the HIV transmission probability included PrEP use, condom use, sexual position, and circumcision (**Appendix S8**). ### Baseline PrEP Cascade HIV screening was also the entry point to PrEP [5]. MSM screening negative entered the PrEP cascade (initiation, adherence, and persistence) if they met indications based on CDC guidelines [2]. Baseline PrEP care was parameterized to reflect PrEP engagement observed in the Atlanta area (**Table 1**). Indicated MSM initiated PrEP at a probability generating an overall coverage level (percent of indicated MSM on PrEP) of 15% [15]. Heterogeneous PrEP adherence was represented as two levels (high versus medium/low), which each corresponded to different relative reductions in HIV acquisition risk [16,17]. Spontaneous PrEP discontinuation was based on data for the proportion of MSM who were retained in PrEP care at 6 months [18]. MSM also stopped PrEP if they no longer exhibited PrEP indications [2]. View this table: [Table 1.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/T1) Table 1. Epidemiological and Economic Model Parameters ### PrEP Cascade Interventions With the calibrated baseline model, we then simulated three categories of PrEP interventions at different levels of capacity. MSM could join one or more interventions, subject to capacity constraints. Engagement in each of the three interventions was which was associated with improvements to in PrEP initiation, adherence, or persistence, respectively. MSM not on an intervention continued status quo PrEP care in that domain. The initiation intervention was modeled after **Healthmindr**, currently in a clinical trial [7]. This intervention features a mobile phone app that delivers information and resources about PrEP, including individualized risk assessments and PrEP provider locations. The goal of the Healthmindr intervention is to lower PrEP barriers to increase the rate of initiation for PrEP-indicated MSM. In our model, MSM who were indicated for but not using PrEP were eligible to start Healthmindr, with a range of app initiation probabilities that translated to different levels of Healthmindr capacity (**Table 1**). MSM on Healthmindr then had an additional PrEP initiation probability (elevating the baseline probability) that corresponded to the intervention efficacy, which was based on power analyses for the study protocol [7]. MSM dropped off the intervention when indications lapsed, but they could reenter Healthmindr upon new indications. The adherence intervention was modeled after **Life-Steps for PrEP**, completed in 2017 [8]. Life-Steps used cognitive behavioral counseling delivered by nurses to improve medication adherence for MSM currently using PrEP. The trial found that PrEP drug levels were significantly higher in the intervention arm. In our model, MSM were eligible for Life-Steps at PrEP initiation. The intervention had a fixed number of new users each week, which we varied up to the maximum clinical service capacity. This intervention had the effect of shifting the probability of assignment into the high adherence category from 60% (baseline) to 100% for those enrolled. This was a one-time intervention with effects lasting the duration of a PrEP episode (intervention impact was not retained for future PrEP episodes). The persistence intervention was modeled after **ePrEP** and **PrEP@Home** interventions, currently in trials [6,19]. Both feature home-based PrEP care accomplished through a mobile phone app, specimen self-collection, and video-based clinical consultations (PrEP@Home adds annual in-person clinical assessment, which we did not represent). Primary intervention outcomes measure persistence on PrEP after initiating the intervention. In our model, MSM were eligible for the persistence intervention at the point of PrEP initiation. The intervention had a fixed capacity of prevalent (active) users, which we varied up to the maximum clinical service capacity expected if the intervention were fully scaled-up. Intervention participants had a lower rate of spontaneous PrEP discontinuation compared to non-intervention PrEP users. This lower rate was informed by power analyses in the study protocol [6]. PrEP users initiating ePrEP/PrEP@Home were assumed to stay on the intervention for the full duration of their PrEP episode. Thus, MSM stopped the intervention only if they discontinued PrEP, either spontaneously or through lapsed indications. ### Epidemic Model Calibration and Simulation We calibrated the epidemic model with a Bayesian approach that defined prior distributions for uncertain parameters and fit the model to empirical surveillance-based estimates for the target population (Appendix S9). After calibration, we randomly drew from a multivariate distribution of capacity parameters for the three interventions (**Table 1**) using Latin hypercube sampling. This generated 1000 capacity parameter sets. For each set, we then simulated the epidemic model 250 times for 10 years in weekly time steps. The primary outcome was the median cumulative 10-year incidence over the 250 simulations. For each of the 1000 capacity combinations, the number and percent of infections averted compared median cumulative incidence in each sampled parameter set against that in the base scenario (no PrEP cascade interventions). ### Budget Optimization Model A nonlinear programming (NLP) model was used to solve for the optimal allocation of the budget for these three interventions, maximizing infections averted. The NLP was defined by a nonlinear objective function, a linear constraint related to the total budget, and bounds for the minimum and maximum allowable intervention capacities. The optimization can be represented mathematically as: ![Formula][1] where *f*() is the objective function, *c**I,A,R* are the unit costs associated with the interventions for the entire 10-year time horizon, and *B* is the 10-year budget allocated to improving the PrEP cascade. The objective function is the expected infections averted over the 10-year time horizon as a function of the clinical capacities of the initiation (*I*), adherence (*A*), and persistence (*R*) interventions. In terms of capacity measures, *I* is the percentage of newly indicated MSM reached by the initiation intervention, *A* is the number of new PrEP users allowed to start the adherence intervention, and *R* is the maximum number of PrEP users engaged in the persistence intervention at any one time. We estimated the objective function through a generalized additive model (GAM) fit to the simulated infections averted over the considered ranges of intervention capacities (**Appendix S10**). For all three interventions, we assumed that the intervention costs increased linearly with capacity. The budget optimization was conducted from a payer-perspective, so only costs that would be paid by a hypothetical federal EHE or local health jurisdiction budget allocated to PrEP improvement were included. For each intervention, we estimated the resources needed to administer the intervention per person or person-time and then estimated unit costs of each resource type, using Atlanta-specific estimates when available (**Table 1**). For the initiation intervention, costs included the advertising needed to attract new PrEP users. These intervention costs also included 10% of the HIV/STI screening resources and laboratory costs incurred by users, which was the estimated proportion of these costs that would be covered directly by the budget based on the local uninsured rate [20]. Adherence intervention costs included personnel and physical space to administer the intervention. ## RESULTS **Figure 1** shows the impact of scaling up the capacity of the initiation and persistence interventions on PrEP coverage (**Panel A**) and HIV incidence (**Panel B**). In **Panel A**, scaling up the initiation and persistence interventions to their maximum considered capacities resulted in PrEP use among 27% of indicated MSM, an increase of 12 percentage points from the baseline PrEP coverage of 15%. Higher capacity on both the initiation and persistence interventions increased overall PrEP coverage, but PrEP coverage improvements from scaling up capacity in the persistence intervention were constrained by the capacity of the initiation intervention. PrEP coverage as a function of persistence intervention capacity plateaued when there was capacity to provide nearly all prevalent PrEP users with the persistence intervention. Additional persistence intervention capacity provided no benefit beyond this point. Before this plateau was reached, however, the marginal effect of scaling up the persistence intervention was similar across different levels of initiation intervention scale-up. **Panel B** shows how the percent infections averted changes as a function of these different combinations of initiation and persistence capacities. When both interventions were scaled-up to their maximum capacities, the total effect was a 9.3% reduction in cumulative incidence. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/11/2020.12.10.20247270/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/F1) Figure 1. Panel A (left) shows the relationship between capacity on initiation (probability of starting initiation intervention) and capacity on persistence interventions on PrEP coverage. Each dot represents one parameter scenario set (medians across 250 individual simulations). Panel B (right) shows the relationship of the initiation and persistence capacity on the percent of cumulative infections averted over 10 years, with lighter colors indicating more infections inverted due to the two interventions in tandem. **Figure 2** demonstrates the joint interaction of the three interventions. Maximizing capacity of all three interventions averted 12.3% of HIV infections compared to the base scenario. The marginal increase in infections averted with increasing persistence intervention capacity was independent of the adherence intervention capacity, except when PrEP initiation was low. When the initiation intervention capacity was low, a unit increase in persistence capacity provided roughly the same benefit for both medium and high adherence intervention capacities. However, for medium-to-high levels of initiation intervention capacity, increasing adherence intervention capacity had the impact of shifting the overall persistence intervention response curve higher. Overall, the more person-time on PrEP generated through initiation and persistence interventions, the more valuable it was to make that person-time highly adherent. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/11/2020.12.10.20247270/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/F2) Figure 2. Three panels show the impact of variations in three interventions together on 10-year cumulative infections averted. Subpanels show low (minimum capacity), medium (median capacity), and high (maximum capacity) for the initiation intervention, where the individual lines in each panel show low, medium, and high capacities of the adherence intervention (same definitions as intervention capacity). Thick lines represent point estimates and bands represent 95% confidence intervals around those estimates. **Figure 3** and **Table 2** present the results of the model of investment strategies across 10-year budget levels that maximize infections averted. The optimal strategy involved a mixture of initiation and persistence interventions investment across all budget levels. At no budget level was it considered optimal to invest into the adherence intervention. At lower budgets (less than $1.45 million), investment entirely in the persistence intervention was optimal. As the budget increased, a mixture of investments in the initiation and persistence interventions achieved the greatest benefits. As **Panel B** shows, higher levels of investment in the initiation intervention did not displace the budget allocated to the persistence intervention. Instead, the optimal budget allocation was to expand the persistence intervention capacity to fully meet the demand from MSM on PrEP who could utilize it; once that demand was satisfied, any remaining funds were allocated to increase initiation. At the highest 10-year budget level ($6 million), the optimal allocation was to invest 26% of the budget into the persistence intervention, while investing the remaining budget into initiation. At this budget level and investment strategy, PrEP coverage would increase to 24% (from 15% in the base model), averting 76.5 (6.4%) HIV infections over 10 years. The intervention cost per HIV infection averted increased with the total budget, ranging from $25,048 at the lowest budget to $69,281 for the highest total budget considered, indicating diminishing returns as the budget was increased. View this table: [Table 2.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/T2) Table 2. Allocation to Three PrEP Cascade Interventions (Initiation, Adherence, Persistence) across Different 10-Year Budgets to Optimize Cumulative HIV Incidence Reduction over 10 Years (Primary Budget Scenarios) ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/11/2020.12.10.20247270/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/F3) Figure 3. Results of the primary optimization analysis on budgetary allocation to each of the three interventions over different total 10-year budget sizes. Panel A (left) shows the fraction of the budget allocated to each intervention category, and Panel B (right) shows the overall funds allocated. Initiation and persistence interventions dominate the adherence intervention across all considered budget sizes. Despite having an epidemiological benefit on HIV incidence (**Figure 2**), the adherence intervention was not selected for investment in the primary budget optimization model. The reason was its higher relative costs per infection averted compared to the initiation and persistence interventions. This was driven by the higher unit costs associated with nurse personnel time for behavioral counseling in that intervention. In a sensitivity analysis (**Figure 4** and **Supplemental Table 14**), we explored how reductions in the costs for this intervention (e.g., by having the counseling conducted by non-nurse staff) would impact the budget optimization. When the adherence intervention unit costs were reduced by approximately half (55%), the adherence intervention was a favorable investment. The persistence intervention still absorbed all of the funding at the lowest budget levels. However, unlike with the primary optimization model, when the persistence capacity was met, the adherence intervention was first selected for investment over the initiation intervention. The adherence intervention then exhibited a diminishing return at middle budget sizes (approximately $3.9 million for 10 years), at which point some investment in the initiation intervention was optimal. At the highest budget levels, the distribution of spending for the initiation, adherence, and persistence interventions was roughly equal (35%, 40%, and 25%, respectively). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/11/2020.12.10.20247270/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/12/11/2020.12.10.20247270/F4) Figure 4. Results of the sensitivity optimization analysis on budgetary allocation to each of the three interventions over different total 10-year budget sizes, assuming that the unit costs of the adherence intervention are reduced by approximately half. Panel A (left) shows the fraction of the budget allocated to each intervention category, and Panel B (right) shows the overall funds allocated. ## DISCUSSION This study investigated how to optimally allocate funds targeting gaps in the PrEP cascade across three categories of interventions. This was the first decision analytics model using a budget optimization approach for PrEP investment. We found that with increasing EHE budgets, the persistence intervention should first be funded until its capacity fully satisfies user demand. At that point, the initiation intervention should be funded to increase overall PrEP coverage, and to further grow demand for the persistence intervention. Investments in the adherence intervention were not optimal based on current estimated unit costs; if those costs were reduced by half, t the optimal allocation would fund three interventions at roughly equal levels. This study highlights that EHE and related investments into the PrEP cascade should account for the interacting impacts of different interventions as they are deployed. Despite progress in delivery of PrEP to populations in need, persistent gaps in the PrEP cascade limit its optimal use [21]. Our model estimated the optimal distribution of investments in interventions to close these gaps, given the epidemiology of an EHE target jurisdiction [1]. EHE investments have the effect of changing the current baseline conditions in meaningful ways. Each level of increasing investment in these interventions effectively shifts the baseline to a better level. Further investment increases the overall infections averted but does so at a higher cost per infection averted. Our study demonstrates the importance of understanding the mechanisms for how to close these PrEP cascade gaps. An important implementation science question is why weaker than projected prevention benefits are seen as PrEP is scaled up in practice [22]. One analysis found only modest declines in HIV diagnoses in states with higher PrEP prescription rates [23]. While this may be partially due to diagnoses being an imperfect proxy for incidence, PrEP prescriptions do not necessarily translate into optimal levels of PrEP use [24]. Retaining MSM in PrEP care continues to be a major challenge [25], and MSM are at much higher HIV risk immediately after discontinuing PrEP than while on PrEP [26]. Our model suggests that spending on initiation is optimal only once PrEP persistence gaps are closed. The benefits of generating new PrEP prescriptions (e.g., with the initiation intervention) will not translate into maximum prevention benefits unless those new users are retained on PrEP during their sexual risk period. Although the adherence intervention also had an epidemiological impact, reducing its unit costs would be required for selection into an optimal investment mix. This might be possible by having peer educators rather than nurses deliver the counseling sessions, but the effectiveness of that approach has not been evaluated. Our representation of the PrEP cascade mechanisms drives the methodological innovation of this study. Many HIV economic models have used cost-effectiveness analysis (CEA) approaches to balance the financial costs against the epidemiological benefits of PrEP [27,28]. However, CEAs primarily focus on the efficiency of PrEP itself and assume it can be provided and scaled up to high levels of coverage among MSM; CEAs of PrEP have not, to date, accounted for the resources needed to achieve those levels of scale-up and to support high levels of adherence and persistence. This is an illustration of a common bias in CEAs to focus on the clinical and pharmaceutical aspects of HIV prevention and treatment interventions, while ignoring the cost of the supportive activities needed to ensure the success of these biomedical interventions. Furthermore, CEAs can inform whether an HIV prevention strategy, such as PrEP, achieve health benefits at a reasonable level of health sector efficiency. However, CEAs do not inform whether a given intervention should be adopted by a given payer, especially when that payer is subject to substantially smaller budgets than the entirety of US health care spending. Our budget optimization approach, in contrast, takes a direct payer perspective and can therefore specifically inform how the CDC or local health jurisdictions should spend EHE funds on services to amplify PrEP benefits in order to maximize HIV prevention benefits. As the price of PrEP medication decreases through generic TDF/FTC [29] or subsidies from the *Ready, Set, PrEP* program, the costs of ancillary clinical services (e.g., routine STI screening in PrEP care [30]) and interventions to close the gaps in the PrEP cascade will constitute a larger portion of PrEP-related spending. Thus, financial considerations for PrEP programs will shift from financing medication access to funding initiatives that support optimal PrEP use. ### Limitations The primary limitation of our study is the assumptions about the individual-level efficacy of the three modeled interventions. For Healthmindr and ePrEP/PrEP@Home, efficacy has not been precisely determined because their clinical trials are incomplete. However, our efficacy estimates were drawn from power analyses in the trial protocols, which themselves were based on pilot data on the expected efficacy. We expect our results would apply to other interventions with demonstrated efficacy similar our assumed values. Second, we made our best effort to quantify the unit costs for these interventions that would likely be paid by EHE or local funds if the interventions were deployed. This required some cost assumptions (e.g., online advertising costs needed to generate a new PrEP user). Formal costing studies are also underway for the ongoing trials that could change these selected costs. Finally, our results may not be transportable to settings with different baseline HIV incidence or PrEP coverage. Our target population of Atlanta MSM was selected based on its reflection of many EHE target jurisdictions precisely where these PrEP cascade interventions are needed. ### Conclusions Closing the gaps between actual and recommended PrEP use will depend on deploying interventions that address the steps of the PrEP cascade. Several such interventions have been developed, but models for how to collectively invest in them within HIV prevention initiatives like the EHE are currently required needed. Our study highlights the importance of understanding the mechanistic synergies of these interventions when optimizing these investments. ## Supporting information Supplemental Appendix [[supplements/247270_file02.pdf]](pending:yes) ## Data Availability The model code and parameters are all open and availability at the Github repository linked below. [https://github.com/EpiModel/PrEP-Optimize](https://github.com/EpiModel/PrEP-Optimize) ## Footnotes * **FUNDING** This work was supported by Centers for Disease Control and Prevention cooperative agreement number U38 PS004646 and National Institutes of Health grant R01 AI138783. Dr. Marcus is supported in part by National Institutes of Health grant K01 AI122853. * **DISCLAIMER** The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the official positions of the Centers for Disease Control and Prevention. * **DISCLOSURES** The authors declare no conflicts of interest. * Received December 10, 2020. * Revision received December 10, 2020. * Accepted December 11, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. Siegler AJ, Mouhanna F, Giler RM, Weiss K, Pembleton E, Guest J, et al. The prevalence of pre-exposure prophylaxis use and the pre-exposure prophylaxis-to-need ratio in the fourth quarter of 2017, United States. Ann Epidemiol 2018; 28:841–849. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F11%2F2020.12.10.20247270.atom) 2. Centers for Disease Control and Prevention. Preexposure prophylaxis for the prevention of HIV infection in the United States–2017 Update: A Clinical Practice Guideline.; 2017. [https://www.cdc.gov/hiv/pdf/risk/prep/cdc-hiv-prep-guidelines-2017.pdf](https://www.cdc.gov/hiv/pdf/risk/prep/cdc-hiv-prep-guidelines-2017.pdf) 3. Sullivan PS, Sanchez TH, Zlotorzynska M, Chandler CJ, Sineath RC, Kahle E, et al. National trends in HIV pre-exposure prophylaxis awareness, willingness and use among United States men who have sex with men recruited online, 2013 through 2017. J Int AIDS Soc 2020; 23:e25461. 4. Fauci AS, Redfield RR, Sigounas G, Weahkee MD, Giroir BP. Ending the HIV Epidemic. JAMA 2019; 321:844. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2019.1343&link_type=DOI) 5. Horn T, Sherwood J, Remien RH, Nash D, Auerbach JD. Towards an integrated primary and secondary HIV prevention continuum for the United States: A cyclical process model. J Int AIDS Soc 2016; 19:21263. 6. Siegler AJ, Brock JB, Hurt CB, Ahlschlager L, Dominguez K, Kelley CF, et al. An Electronic Pre-Exposure Prophylaxis Initiation and Maintenance Home Care System for Nonurban Young Men Who Have Sex With Men: Protocol for a Randomized Controlled Trial. JMIR Res Protoc 2019; 8:e13982. 7. Jones J, Dominguez K, Stephenson R, Stekler JD, Castel AD, Mena LA, et al. A Theoretically Based Mobile App to Increase Pre-Exposure Prophylaxis Uptake Among Men Who Have Sex With Men: Protocol for a Randomized Controlled Trial. JMIR Res Protoc 2020; 9:e16231. 8. Mayer KH, Safren SA, Elsesser SA, Psaros C, Tinsley JP, Marzinke M, et al. Optimizing Pre-Exposure Antiretroviral Prophylaxis Adherence in Men Who Have Sex with Men: Results of a Pilot Randomized Controlled Trial of “Life-Steps for PrEP.” AIDS Behav 2017; 21:1350–1360. 9. Jenness SM, Maloney KM, Smith DK, Hoover KW, Goodreau SM, Rosenberg ES, et al. Addressing Gaps in HIV Preexposure Prophylaxis Care to Reduce Racial Disparities in HIV Incidence in the United States. Am J Epidemiol 2019; 188:743–752. 10. Chan PA, Goedel WC, Nunn AS, Sowemimo-Coker G, Galárraga O, Prosperi M, et al. Potential Impact of Interventions to Enhance Retention in Care During Real-World HIV Pre-Exposure Prophylaxis Implementation. AIDS Patient Care STDS 2019; 33:434–439. 11. Jenness SM, Goodreau SM, Morris M. EpiModel: An R Package for Mathematical Modeling of Infectious Disease over Networks. J Stat Softw 2018; 84:1–47. 12. Jenness SM, Johnson JA, Hoover KW, Smith DK, Delaney KP. Modeling an integrated HIV prevention and care continuum to achieve the ending the HIV epidemic goals. AIDS Published Online First: 8 September 2020. doi:10.1097/QAD.0000000000002681 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAD.0000000000002681&link_type=DOI) 13. Weiss KM, Goodreau SM, Morris M, Prasad P, Ramaraju R, Sanchez T, et al. Egocentric sexual networks of men who have sex with men in the United States: Results from the ARTnet study. Epidemics 2020; 30:100386. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2020.100386&link_type=DOI) 14. Georgia Department of Public Health, HIV/AIDS Epidemiology Section. HIV Surveillance Summary, Georgia 2017.; 2019. [https://dph.georgia.gov/data-fact-sheet-summaries](https://dph.georgia.gov/data-fact-sheet-summaries) (accessed 20 Apr2020). 15. Weiss KM, Prasad P, Ramaraju R, Zlotorzynska M, Jenness SM. Estimated Number of Men who have Sex with Men with Indications for HIV Pre-Exposure Prophylaxis in a National Sexual Network Study. Journal of Acquired Immune Deficiency Syndromes doi:10.1097/QAI.0000000000002300 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAI.0000000000002300&link_type=DOI) 16. Liu AY, Cohen SE, Vittinghoff E, Anderson PL, Doblecki-Lewis S, Bacon O, et al. Preexposure Prophylaxis for HIV Infection Integrated With Municipal- and Community-Based Sexual Health Services. JAMA Intern Med 2016; 176:75–84. 17. Grant RM, Anderson PL, McMahan V, Liu A, Amico KR, Mehrotra M, et al. Uptake of pre-exposure prophylaxis, sexual practices, and HIV incidence in men and transgender women who have sex with men: a cohort study. Lancet Infect Dis 2014; 14:820–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(14)70847-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25065857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F11%2F2020.12.10.20247270.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000341342700025&link_type=ISI) 18. Chan PA, Mena L, Patel R, Oldenburg CE, Beauchamps L, Perez-Brumer AG, et al. Retention in care outcomes for HIV pre-exposure prophylaxis implementation programmes among men who have sex with men in three US cities. J Int AIDS Soc 2016; 19:20903. 19. Siegler AJ, Mayer KH, Liu AY, Patel RR, Ahlschlager LM, Kraft CS, et al. Developing and Assessing the Feasibility of a Home-based Preexposure Prophylaxis Monitoring and Support Program. Clin Infect Dis 2019; 68:501–504. 20. Smith DK, Van Handel M, Huggins R. Estimated Coverage to Address Financial Barriers to HIV Preexposure Prophylaxis Among Persons With Indications for Its Use, United States, 2015. J Acquir Immune Defic Syndr 2017; 76:465–472. 21. Kelley CF, Kahle E, Siegler A, Sanchez T, Del Rio C, Sullivan PS, et al. Applying a PrEP Continuum of Care for Men who Have Sex with Men in Atlanta, GA. Clin Infect Dis 2015; 61:1590–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ664&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26270691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F11%2F2020.12.10.20247270.atom) 22. Jenness SM, Goodreau SM, Rosenberg E, Beylerian EN, Hoover KW, Smith DK, et al. Impact of the Centers for Disease Control’s HIV preexposure prophylaxis guidelines for men who have sex with men in the United States. J Inf Dis 2016; 214:1800–1807. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiw223&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27418048&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F11%2F2020.12.10.20247270.atom) 23. Smith DK, Sullivan PS, Cadwell B, Waller LA, Siddiqi A, Mera-Giler R, et al. Evidence of an Association of Increases in Pre-exposure Prophylaxis Coverage With Decreases in Human Immunodeficiency Virus Diagnosis Rates in the United States, 2012-2016. Clin Infect Dis Published Online First: 25 February 2020. doi:10.1093/cid/ciz1229 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciz1229&link_type=DOI) 24. Marcus JL, Paltiel AD, Walensky RP. Has Pre-exposure Prophylaxis Made a Difference at a Population Level? Jury Is Still Out. Clin Infect Dis Published Online First: 25 February 2020. doi:10.1093/cid/ciz1236 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciz1236&link_type=DOI) 25. Coy KC, Hazen RJ, Kirkham HS, Delpino A, Siegler AJ. Persistence on HIV preexposure prophylaxis medication over a 2-year period among a national sample of 7148 PrEP users, United States, 2015 to 2017. J Int AIDS Soc 2019; 22:e25252. 26. Spinelli MA, Laborde N, Kinley P, Whitacre R, Scott HM, Walker N, et al. Missed opportunities to prevent HIV infections among pre-exposure prophylaxis users: a population-based mixed methods study, San Francisco, United States. J Int AIDS Soc 2020; 23:e25472. 27. Nosyk B, Zang X, Krebs E, Enns B, Min JE, Behrends CN, et al. Ending the HIV epidemic in the USA: an economic modelling study in six cities. The Lancet HIV Published Online First: 5 March 2020. doi:10.1016/S2352-3018(20)30033-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2352-3018(20)30033-3&link_type=DOI) 28. Juusola JL, Brandeau ML, Owens DK, Bendavid E. The cost-effectiveness of preexposure prophylaxis for HIV prevention in the United States in men who have sex with men. Ann Intern Med 2012; 156:541–550. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/0003-4819-156-8-201204170-00004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22508731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F11%2F2020.12.10.20247270.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303151800013&link_type=ISI) 29. Walensky RP, Horn T, McCann NC, Freedberg KA, Paltiel AD. Comparative Pricing of Branded Tenofovir Alafenamide-Emtricitabine Relative to Generic Tenofovir Disoproxil Fumarate-Emtricitabine for HIV Preexposure Prophylaxis: A Cost-Effectiveness Analysis. Ann Intern Med 2020; 172:583–590. 30. Chandra C, Weiss KM, Kelley CF, Marcus JL, Jenness SM. Gaps in Sexually Transmitted Infection Screening among Men who Have Sex with Men in PrEP Care in the United States. Clin Infect Dis Published Online First: 23 July 2020. doi:10.1093/cid/ciaa1033 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa1033&link_type=DOI) [1]: /embed/graphic-2.gif