Assessing the degree of hepatic ischemia-reperfusion injury using PBPK modeling of sodium fluorescein disposition in ex vivo machine-perfused livers.

Ischemia-reperfusion injury (IRI) is an intrinsic risk associated with liver transplantation. Ex vivo hepatic machine perfusion (MP) is an emerging organ preservation technique that can mitigate IRI, especially in livers subjected to prolonged warm ischemia time (WIT). However, a method to quantify the biological response to WIT during MP has not been established. Previous studies used physiologically-based pharmacokinetic (PBPK) modeling to demonstrate that a decrease in hepatic transport and biliary excretion of the tracer molecule sodium fluorescein (SF) could correlate with increasing WIT in situ. Furthermore, these studies proposed intracellular sequestration of the hepatocyte canalicular membrane transporter multi-drug resistance-associated protein 2 (MRP2) leading to decreased MRP2 activity (maximal transport velocity; Vmax) as the potential mechanism for decreased biliary SF excretion. We adapted an extant PBPK model to account for ex vivo hepatic MP and fit a 6-parameter version of this model to control time course measurements of SF in MP perfusate and bile. We then identified parameters whose values were likely insensitive to changes in WIT and fixed them to generate a reduced model with only 3 unknown parameters. Finally, we fit the reduced model to each individual biological replicate SF time course with differing WIT and found the mean estimated value for each parameter and compared them using a one-way ANOVA. We demonstrated that there was a significant decrease in the estimated value of Vmax for MRP2 at 30 min WIT. These studies provide the foundation for future studies investigating real-time assessment of liver viability during ex vivo MP.

Replication efficiencies of human cytomegalovirus-infected epithelial cells are dependent on source of virus production.

Human cytomegalovirus (HCMV) is a prevalent betaherpesvirus, and infection can lead to a range of symptomatology from mononucleosis to sepsis in immunocompromised individuals. HCMV is also the leading viral cause of congenital birth defects. Lytic replication is supported by many cell types with different kinetics and efficiencies leading to a plethora of pathologies. The goal of these studies was to elucidate HCMV replication efficiencies for viruses produced on different cell types upon infection of epithelial cells by combining experimental approaches with data-driven computational modeling. HCMV was generated from a common genetic background of TB40-BAC4, propagated on fibroblasts (TB40Fb) or epithelial cells (TB40Epi), and used to infect epithelial cells. We quantified cell-associated viral genomes (vDNA), protein levels (pUL44, pp28), and cell-free titers over time for each virus at different multiplicities of infection. We combined experimental quantification with data-driven simulations and determined that parameters describing vDNA synthesis were similar between sources. We found that pUL44 accumulation was higher in TB40Fb than TB40Epi. In contrast, pp28 accumulation was higher in TB40Epi which coincided with a significant increase in titer for TB40Epi over TB40Fb. These differences were most evident during live-cell imaging, which revealed syncytia-like formation during infection by TB40Epi. Simulations of the late lytic replication cycle yielded a larger synthesis constant for pp28 in TB40Epi along with increase in virus output despite similar rates of genome synthesis. By combining experimental and computational modeling approaches, our studies demonstrate that the cellular source of propagated virus impacts viral replication efficiency in target cell types.

Computational modeling of protracted HCMV replication using genome substrates and protein temporal profiles.

Human cytomegalovirus (HCMV) is a major cause of illness in immunocompromised individuals. The HCMV lytic cycle contributes to the clinical manifestations of infection. The lytic cycle occurs over ∼96 h in diverse cell types and consists of viral DNA (vDNA) genome replication and temporally distinct expression of hundreds of viral proteins. Given its complexity, understanding this elaborate system can be facilitated by the introduction of mechanistic computational modeling of temporal relationships. Therefore, we developed a multiplicity of infection (MOI)-dependent mechanistic computational model that simulates vDNA kinetics and late lytic replication based on in-house experimental data. The predictive capabilities were established by comparison to post hoc experimental data. Computational analysis of combinatorial regulatory mechanisms suggests increasing rates of protein degradation in association with increasing vDNA levels. The model framework also allows expansion to account for additional mechanisms regulating the processes. Simulating vDNA kinetics and the late lytic cycle for a wide range of MOIs yielded several unique observations. These include the presence of saturation behavior at high MOIs, inefficient replication at low MOIs, and a precise range of MOIs in which virus is maximized within a cell type, being 0.382 IU to 0.688 IU per fibroblast. The predicted saturation kinetics at high MOIs are likely related to the physical limitations of cellular machinery, while inefficient replication at low MOIs may indicate a minimum input material required to facilitate infection. In summary, we have developed and demonstrated the utility of a data-driven and expandable computational model simulating lytic HCMV infection.

Molecular dissection of cone photoreceptor-enriched genes encoding transmembrane and secretory proteins.

Cone photoreceptors mediate color perception and daylight vision through intricate synaptic circuitry. In most mammalian retina, cones are greatly outnumbered by rods and exhibit inter-dependence for functional maintenance and survival. Currently, we have limited understanding of cone-specific molecular components that mediate response to extrinsic signaling factors or are involved in communication with rods and other retinal cells. To fulfill this gap, we compared the recently-published transcriptomes of developing S-cone-like photoreceptors from the Nrl-/- mouse retina with those of rods and identified candidate genes responsible for cone cell functions and communication. We generated an in silico expression profile of 823 genes that encode candidate transmembrane and secretory proteins and are up-regulated in Nrl-/- cone photoreceptors compared to wild type cones. In situ hybridization analysis validated high expression of seven of the selected candidate genes in the Nrl-/- retina. To examine their relevance to cone function, we performed in vivo knockdown of Epha10 in the Nrl-/- retina and demonstrated aberrant morphology and mislocalization of the photoreceptor cell bodies. Thus, the receptor tyrosine kinase Ephrin type-A receptor 10 appears to influence cone morphogenesis. Our studies reveal novel cone-enriched genes involved in interaction of cones with other retinal cell types and provide a framework for examining molecular pathways associated with intercellular communication.