Virulence-linked pathways in opportunistic pathogens are putative therapeutic targets that may be associated with less potential for resistance than targets in growth-essential pathways. However, efficacy of virulence-linked targets may be affected by the contribution of virulence-related genes to metabolism. We evaluate the complex interrelationships between growth and virulence-linked pathways using a genome-scale metabolic network reconstruction of Pseudomonas aeruginosa strain PA14 and an updated, expanded reconstruction of P. aeruginosa strain PAO1. The PA14 reconstruction accounts for the activity of 112 virulence-linked genes and virulence factor synthesis pathways that produce 17 unique compounds. We integrate eight published genome-scale mutant screens to validate gene essentiality predictions in rich media, contextualize intra-screen discrepancies and evaluate virulence-linked gene distribution across essentiality datasets. Computational screening further elucidates interconnectivity between inhibition of virulence factor synthesis and growth. Successful validation of selected gene perturbations using PA14 transposon mutants demonstrates the utility of model-driven screening of therapeutic targets.

Antibiotic regimens often include the sequential changing of drugs to limit the development and evolution of resistance of bacterial pathogens. It remains unclear how history of adaptation to one antibiotic can influence the resistance profiles when bacteria subsequently adapt to a different antibiotic. Here, we experimentally evolved Pseudomonas aeruginosa to six 2-drug sequences. We observed drug order-specific effects, whereby adaptation to the first drug can limit the rate of subsequent adaptation to the second drug, adaptation to the second drug can restore susceptibility to the first drug, or final resistance levels depend on the order of the 2-drug sequence. These findings demonstrate how resistance not only depends on the current drug regimen but also the history of past regimens. These orderspecific effects may allow for rational forecasting of the evolutionary dynamics of bacteria given knowledge of past adaptations and provide support for the need to consider the history of past drug exposure when designing strategies to mitigate resistance and combat bacterial infections. Author summaryBacteria readily adapt to their environments and can develop ways to survive and grow in the presence of antibiotics. While many studies have investigated how bacteria evolve to become resistant to single drugs, it is unclear how adaptation to other drugs and environments in the past affect the way bacteria adapt to new drugs and environments. In this study, we allowed bacteria in a laboratory setting to adapt to three different antibiotics. We first exposed wild-type susceptible bacteria to high concentrations of the three antibiotics individually and then exposed these populations to each of the other drugs. By tracking the levels of resistance to all three drugs in all of the treatments, we identified cases in which past adaptation to one treatment influenced subsequent evolutionary dynamics with regard to both phenotypes (levels of resistance) and genotypes (genes that became mutated). Additionally, by allowing bacterial isolates originating from human patients to adapt to the three drugs, we recapitulated a subset of the adaptation history-dependent evolutionary dynamics. Overall, this study sheds light on how adaptation history in PLOS Biology | https://doi.org/10.1371/journal

Antibiotic regimens often include the sequential changing of drugs to limit the development and evolution of resistance of bacterial pathogens. It remains unclear how history of adaptation to one antibiotic can influence the resistance profiles when bacteria subsequently adapt to a different antibiotic. Here, we experimentally evolved Pseudomonas aeruginosa to six 2-drug sequences. We observed drug order–specific effects, whereby adaptation to the first drug can limit the rate of subsequent adaptation to the second drug, adaptation to the second drug can restore susceptibility to the first drug, or final resistance levels depend on the order of the 2-drug sequence. These findings demonstrate how resistance not only depends on the current drug regimen but also the history of past regimens. These order-specific effects may allow for rational forecasting of the evolutionary dynamics of bacteria given knowledge of past adaptations and provide support for the need to consider the history of past drug exposure when designing strategies to mitigate resistance and combat bacterial infections.

Burkholderia cenocepacia and Burkholderia multivorans are opportunistic drug-resistant pathogens that account for the majority of Burkholderia cepacia complex infections in cystic fibrosis patients and also infect other immunocompromised individuals. While they share similar genetic compositions, B. cenocepacia and B. multivorans exhibit important differences in pathogenesis. We have developed reconciled genome-scale metabolic network reconstructions of B. cenocepacia J2315 and B. multivorans ATCC 17616 in parallel (designated iPY1537 and iJB1411, respectively) to compare metabolic abilities and contextualize genetic differences between species. The reconstructions capture the metabolic functions of the two species and give insight into similarities and differences in their virulence and growth capabilities. The two reconstructions have 1,437 reactions in common, and iPY1537 and iJB1411 have 67 and 36 metabolic reactions unique to each, respectively. After curating the extensive reservoir of metabolic genes in Burkholderia, we identified 6 genes essential to growth that are unique to iPY1513 and 13 genes uniquely essential to iJB1411. The reconstructions were refined and validated by comparing in silico growth predictions to in vitro growth capabilities of B. cenocepacia J2315, B. cenocepacia K56-2, and B. multivorans ATCC 17616 on 104 carbon sources. Overall, we identified functional pathways that indicate B. cenocepacia can produce a wider array of virulence factors compared to B. multivorans, which supports the clinical observation that B. cenocepacia is more virulent than B. multivorans. The reconciled reconstructions provide a framework for generating and testing hypotheses on the metabolic and virulence capabilities of these two related emerging pathogens. Multidrug-resistant pathogens are a severe health concern and can cause chronic infections in a variety of patient populations with limited recourse for treatment. Here, we investigate two multidrug-resistant species, Burkholderia cenocepacia and Burkholderia multivorans, of the Burkholderia cepacia complex (BCC) which are considered dangerous and difficult to treat in patients with cystic fibrosis (CF), chronic granulomatous disease, or compromised immune systems (1). With larger genomes (8.06 Mbp and 7.01 Mbp, respectively) than many other multidrug-resistant pathogens, they also contain an expanded reservoir of genes that may assist their ability to avoid clinical eradication (2-4). Because they are nosocomial, transmissible between patients, and also routinely acquired (and reacquired) from the environment, B. cenocepacia and B. multivorans are the two BCC species most commonly isolated from the sputum of CF patients (5-8). In patients with CF, pulmonary infection with BCC can contribute to the rapid deterioration of lung function known as cepacia syndrome, a necrotizing pneumonia that can lead to bacteremia, septicemia, and increased mortality (9, 10). A combination of high antibiotic resistance and decreased immune function in patients makes B. cen...

Vascular endothelial growth factor (VEGF) is a key regulator of angiogenesis – the growth of new microvessels from existing microvasculature. Angiogenesis is a complex process involving numerous molecular species, and to better understand it, a systems biology approach is necessary. In vivo preclinical experiments in the area of angiogenesis are typically performed in mouse models; this includes drug development targeting VEGF. Thus, to quantitatively interpret such experimental results, a computational model of VEGF distribution in the mouse can be beneficial. In this paper, we present an in silico model of VEGF distribution in mice, determine model parameters from existing experimental data, conduct sensitivity analysis, and test the validity of the model.The multiscale model is comprised of two compartments: blood and tissue. The model accounts for interactions between two major VEGF isoforms (VEGF120 and VEGF164) and their endothelial cell receptors VEGFR-1, VEGFR-2, and co-receptor neuropilin-1. Neuropilin-1 is also expressed on the surface of parenchymal cells. The model includes transcapillary macromolecular permeability, lymphatic transport, and macromolecular plasma clearance. Simulations predict that the concentration of unbound VEGF in the tissue is approximately 50-fold greater than in the blood. These concentrations are highly dependent on the VEGF secretion rate. Parameter estimation was performed to fit the simulation results to available experimental data, and permitted the estimation of VEGF secretion rate in healthy tissue, which is difficult to measure experimentally. The model can provide quantitative interpretation of preclinical animal data and may be used in conjunction with experimental studies in the development of pro- and anti-angiogenic agents. The model approximates the normal tissue as skeletal muscle and includes endothelial cells to represent the vasculature. As the VEGF system becomes better characterized in other tissues and cell types, the model can be expanded to include additional compartments and vascular elements.