Search for: All records

Abstract Despite the revolutionary impact of immune checkpoint inhibition on cancer therapy, the lack of response in a subset of patients, as well as the emergence of resistance, remain significant challenges. Here we explore the theoretical consequences of the existence of multiple states of immune cell exhaustion on response to checkpoint inhibition therapy. In particular, we consider the emerging understanding that T cells can exist in various states: fully functioning cytotoxic cells, reversibly exhausted cells with minimal cytotoxicity, and terminally exhausted cells. We hypothesize that inflammation augmented by drug activity triggers transitions between these phenotypes, which can lead to non-genetic resistance to checkpoint inhibitors. We introduce a conceptual mathematical model, coupled with a standard 2-compartment pharmacometric (PK) model, that incorporates these mechanisms. Simulations of the model reveal that, within this framework, the emergence of resistance to checkpoint inhibitors can be mitigated through altering the dose and the frequency of administration. Our analysis also reveals that standard PK metrics do not correlate with treatment outcome. However, we do find that levels of inflammation that we assume trigger the transition from the reversibly to terminally exhausted states play a critical role in therapeutic outcome. A simulation of a population that has different values of this transition threshold reveals that while the standard high-dose, low-frequency dosing strategy can be an effective therapeutic design for some, it is likely to fail a significant fraction of the population. Conversely, a metronomic-like strategy that distributes a fixed amount of drug over many doses given close together is predicted to be effective across the entire simulated population, even at a relatively low cumulative drug dose. We also demonstrate that these predictions hold if the transitions between different states of immune cell exhaustion are triggered by prolonged antigen exposure, an alternative mechanism that has been implicated in this process. Our theoretical analyses demonstrate the potential of mitigating resistance to checkpoint inhibitors via dose modulation. 

Abstract Mathematical models are increasingly being developed and calibrated in tandem with data collection, empowering scientists to intervene in real time based on quantitative model predictions. Well-designed experiments can help augment the predictive power of a mathematical model but the question of when to collect data to maximize its utility for a model is non-trivial. Here we define data as model-informative if it results in a unique parametrization, assessed through the lens of practical identifiability. The framework we propose identifies an optimal experimental design (how much data to collect and when to collect it) that ensures parameter identifiability (permitting confidence in model predictions), while minimizing experimental time and costs. We demonstrate the power of the method by applying it to a modified version of a classic site-of-action pharmacokinetic/pharmacodynamic model that describes distribution of a drug into the tumor microenvironment (TME), where its efficacy is dependent on the level of target occupancy in the TME. In this context, we identify a minimal set of time points when data needs to be collected that robustly ensures practical identifiability of model parameters. The proposed methodology can be applied broadly to any mathematical model, allowing for the identification of a minimally sufficient experimental design that collects the most informative data. 

PremisePhylogenetic relationships within major angiosperm clades are increasingly well resolved, but largely informed by plastid data. Areas of poor resolution persist within the Dipsacales, including placement ofHeptacodiumandZabelia, and relationships within the Caprifolieae and Linnaeeae, hindering our interpretation of morphological evolution. Here, we sampled a significant number of nuclear loci using a Hyb‐Seq approach and used these data to infer the Dipsacales phylogeny and estimate divergence times. MethodsSampling all major clades within the Dipsacales, we applied the Angiosperms353 probe set to 96 species. Data were filtered based on locus completeness and taxon recovery per locus, and trees were inferred using RAxML and ASTRAL. Plastid loci were assembled from off‐target reads, and 10 fossils were used to calibrate dated trees. ResultsVarying numbers of targeted loci and off‐target plastomes were recovered from most taxa. Nuclear and plastid data confidently placeHeptacodiumwith Caprifolieae, implying homoplasy in calyx morphology, ovary development, and fruit type. Placement ofZabelia, and relationships within the Caprifolieae and Linnaeeae, remain uncertain. Dipsacales diversification began earlier than suggested by previous angiosperm‐wide dating analyses, but many major splitting events date to the Eocene. ConclusionsThe Angiosperms353 probe set facilitated the assembly of a large, single‐copy nuclear dataset for the Dipsacales. Nevertheless, many relationships remain unresolved, and resolution was poor for woody clades with low rates of molecular evolution. We favor expanding the Angiosperms353 probe set to include more variable loci and loci of special interest, such as developmental genes, within particular clades. 

Abstract Conjugation is used by bacteria to propagate antimicrobial resistance (AMR) in the environment. Central to this process are widespread conjugative F-pili that establish the connection between donor and recipient cells, thereby facilitating the spread of IncF plasmids among enteropathogenic bacteria. Here, we show that the F-pilus is highly flexible but robust at the same time, properties that increase its resistance to thermochemical and mechanical stresses. By a combination of biophysical and molecular dynamics methods, we establish that the presence of phosphatidylglycerol molecules in the F-pilus contributes to the structural stability of the polymer. Moreover, this structural stability is important for successful delivery of DNA during conjugation and facilitates rapid formation of biofilms in harsh environmental conditions. Thus, our work highlights the importance of F-pilus structural adaptations for the efficient spread of AMR genes in a bacterial population and for the formation of biofilms that protect against the action of antibiotics. 

The properties of alkanes are consequential for understanding many chemical processes in nature and industry. We use molecular dynamics simulations with the Amber force field GAFF2 to examine the structure of pure liquids at each respective normal boiling point, spanning the 15 n-alkanes from methane to pentadecane. The densities predicted from the simulations are found to agree well with reported experimental values, with an average deviation of 1.9%. The enthalpies of vaporization have an average absolute deviation from experiment of 10.4%. Radial distribution functions show that short alkanes have distinct local structures that are found to converge with each other with increasing chain length. This provides a unique perspective on trends in the n-alkane series and will be useful for interpreting similarities and differences in the n-alkane series as well as the breakdown of ideal solution behavior in mixtures of these molecules. 

Escherichia coliexpress adhesion pili that mediate attachment to host cell surfaces and are exposed to body fluids in the urinary and gastrointestinal tracts. Pilin subunits are organized into helical polymers, with a tip adhesin for specific host binding. Pili can elastically unwind when exposed to fluid flow forces, reducing the adhesin load, thereby facilitating sustained attachment. Here we investigate biophysical and structural differences of pili commonly expressed on bacteria that inhabit the urinary and intestinal tracts. Optical tweezers measurements reveal that class 1a pili of uropathogenicE. coli(UPEC), as well as class 1b of enterotoxigenicE. coli(ETEC), undergo an additional conformational change beyond pilus unwinding, providing significantly more elasticity to their structure than ETEC class 5 pili. Examining structural and steered molecular dynamics simulation data, we find that this difference in class 1 pili subunit behavior originates from an α-helical motif that can unfold when exposed to force. A disulfide bond cross-linking β-strands in class 1 pili stabilizes subunits, allowing them to tolerate higher forces than class 5 pili that lack this covalent bond. We suggest that these extra contributions to pilus resiliency are relevant for the UPEC niche, since resident bacteria are exposed to stronger, more transient drag forces compared to those experienced by ETEC bacteria in the mucosa of the intestinal tract. Interestingly, class 1b ETEC pili include the same structural features seen in UPEC pili, while requiring lower unwinding forces that are more similar to those of class 5 ETEC pili.