Tarunendu Mapder

Pharmacokinetics and pharmacodynamics (PKPD) are key considerations in any study of molecular therapies. It is thus imperative to factor their effects into any in silico model of biological tissue involving such therapies. Furthermore, creating a standardized and flexible framework will benefit the community by increasing access to such modules and enhancing their communicability. PhysiCell is an open-source physics-based cell simulator, i.e., a platform for modeling biological tissue, that is quickly being adopted and utilized by the mathematical biology community. We present here PhysiPKPD, an open-source PhysiCell-based package that allows users to include PKPD in PhysiCell models. Availability & Implementation The source code for PhysiPKPD is located here: https://github.com/drbergman/PhysiPKPD.

We extend our established agent-based multiscale computational model of infection of lung tissue by SARS-CoV-2 to include pharmacokinetic and pharmacodynamic models of remdesivir. We model remdesivir treatment for COVID-19; however, our methods are general to other viral infections and antiviral therapies. We investigate the effects of drug potency, drug dosing frequency, treatment initiation delay, antiviral half-life, and variability in cellular uptake and metabolism of remdesivir and its active metabolite on treatment outcomes in a simulated patch of infected epithelial tissue. Non-spatial deterministic population models which treat all cells of a given class as identical can clarify how treatment dosage and timing influence treatment efficacy. However, they do not reveal how cell-to-cell variability affects treatment outcomes. Our simulations suggest that for a given treatment regime, including cell-to-cell variation in drug uptake, permeability and metabolism increase the likelihood of uncontrolled infection as the cells with the lowest internal levels of antiviral act as super-spreaders within the tissue. The model predicts substantial variability in infection outcomes between similar tissue patches for different treatment options. In models with cellular metabolic variability, antiviral doses have to be increased significantly (>50% depending on simulation parameters) to achieve the same treatment results as with the homogeneous cellular metabolism.

Marine sedimentation rate and bottom-water O ₂ concentration control the remineralization/sequestration of organic carbon across continental margins; but whether/how they shape microbiome architecture (the ultimate effector of all biogeochemical phenomena), across shelf/slope sediments, is unknown. Here we reveal distinct microbiome structures and functions, amidst comparable pore-fluid chemistries, along ~3 m sediment-horizons underlying the seasonal (shallow coastal) and perennial (deep sea) oxygen minimum zones (OMZs) of the Arabian Sea, situated across the western-Indian margin (water-depths: 31 m and, 530 and 580 m, respectively). Along the perennial- and seasonal-OMZ sediment-cores microbial communities were predominated by Gammaproteobacteria/Alphaproteobacteria and Euryarchaeota/Firmicutes respectively. As a perennial-OMZ signature, a cryptic methane production-consumption cycle was found to operate near the sediment-surface; overall diversity, as well as the relative abundances of simple-fatty-acids-requiring anaerobes (methanogens, anaerobic methane-oxidizers, sulfate-reducers and acetogens), peaked in the topmost sediment-layer and then declined via synchronized fluctuations until the sulfate-methane transition zone was reached. The entire microbiome profile was reverse in the seasonal-OMZ sediment-horizon. We discerned that in the perennial-OMZ sediments organic carbon deposited was higher in concentration, and marine components-rich, so it potentially degraded readily to simple fatty acids; lower sedimentation rate afforded higher O ₂ exposure time for organic matter degradation despite perennial hypoxia in the bottom-water; thus, the resultant abundance of reduced metabolites sustained multiple inter-competing microbial processes in the upper sediment-layers. Remarkably, the whole geomicrobial scenario was opposite in the sediments of the seasonal/shallow-water OMZ. Our findings create a microbiological baseline for understanding carbon-sulfur cycling across distinct marine depositional settings and water-colum n oxygenation regimes.

The ecology of aerobic microorganisms is never explored in marine oxygen minimum zone (OMZ) sediments. Here we reveal aerobic bacterial communities along ∼3 m sediment-horizons of the eastern Arabian Sea OMZ. Sulfide-containing sediment-cores retrieved from 530 mbsl (meters beneath the sea-level) and 580 mbsl were explored at 15–30 cm intervals, using metagenomics, pure-culture-isolation, genomics and metatranscriptomics. Genes for aerobic respiration, and oxidation of methane/ammonia/alcohols/thiosulfate/sulfite/organosulfur-compounds, were detected in the metagenomes from all 25 sediment-samples explored. Most probable numbers for aerobic chemolithoautotrophs and chemoorganoheterotrophs at individual sample-sites were up to 1.1 × 107 (g sediment)-1. The sediment-sample collected from 275 cmbsf (centimeters beneath the seafloor) of the 530-mbsl-core yielded many such obligately aerobic isolates belonging to Cereibacter, Guyparkeria, Halomonas, Methylophaga, Pseudomonas and Sulfitobacter which died upon anaerobic incubation, despite being provided with all possible electron acceptors and fermentative substrates. High percentages of metatranscriptomic reads from the 275 cmbsf sediment-sample, and metagenomic reads from all 25 sediment-samples, matched the isolates’ genomic sequences including those for aerobic metabolisms, genetic/environmental information processing and cell division, thereby illustrating the bacteria's in-situ activity, and ubiquity across the sediment-horizons, respectively. The findings hold critical implications for organic carbon sequestration/remineralization, and inorganic compounds oxidation, within the sediment realm of global marine OMZs.

To explore the potential role of tetrathionate in the sedimentary sulfur cycle, population ecology of microorganisms capable of metabolizing this polythionate was revealed at 15–30 cm resolution along two, ∼3 m long, cores collected from 530 and 580 m below the sea level, off India's west coast, within the oxygen minimum zone (OMZ) of the Arabian Sea. Metagenome analysis along the cores revealed widespread occurrence of genes involved in the formation, oxidation, and reduction of tetrathionate; high diversity and relative abundance were also detected for bacteria that are known to render these metabolisms in vitro. Results of slurry culture of the sediment samples in thiosulfate- or tetrathionate-containing microbial growth media, data obtained via pure-culture isolation, and finally metatranscriptome analyses corroborated the in situ functionality of the tetrathionate-forming, tetrathionate-oxidizing, and tetrathionate-reducing microorganisms. Ion chromatography of pore waters revealed the presence of up to 11.1 µM thiosulfate in the two cores, whereas tetrathionate remained undetected in spectroscopic assay based on its reaction with cyanide. While thiosulfate oxidation by chemolithotrophic bacteria prevalent in situ is the apparent source of tetrathionate in this ecosystem, high biochemical and geochemical reactivity of this polythionate could be instrumental in its cryptic status in the sulfur cycle. Potential abiotic origin of tetrathionate in the sediment horizon explored could neither be ruled out nor confirmed from the geochemical information available. On the other hand, tetrathionate potentially present in the system can be either oxidized to sulfate or reduced back to thiosulfate/sulfide via chemolithotrophic oxidation and respiration by native bacterial populations, respectively. Up to 2.01 mM sulfide present in the sediment cores may also reduce tetrathionate abiotically to thiosulfate and elemental sulfur. However, in the absence of measured data for O₂ or other oxyanions having possibilities of serving as electron acceptors, the biogeochemical modalities of the oxidative half of the tetrathionate cycle remained unresolved.

The host-vector shuttle and the bottleneck in dengue transmission is a significant aspect with regard to the study of dengue outbreaks. As mosquitoes require 100–1000 times more virus to become infected than human, the transmission of dengue virus from human to mosquito is a vulnerability that can be targeted to improve disease control. In order to capture the heterogeneity in the infectiousness of an infected patient population towards the mosquito population, we calibrate a population of host-to-vector virus transmission models based on an experimentally quantified infected fraction of a mosquito population. Once the population of models is well-calibrated, we deploy a population of controls that helps to inhibit the human-to-mosquito transmission of the dengue virus indirectly by reducing the viral load in the patient body fluid. We use an optimal bang-bang control on the administration of the defective virus (transmissible interfering particles (TIPs)) to symptomatic patients in the course of their febrile period and observe the dynamics in successful reduction of dengue spread into mosquitoes.