Project description:Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. In this study, we refined the previously published mathematical model for bloodstream form parasites and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. Through this study, we conclude that lLevels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting acquisition of new control mechanisms during adaptation to mammalian parasitism. Ribosome profiling and mRNA libraries were constructed in triplicate from in vitro PCF and in duplicate from in vitro T. brucei Lister427, to understand global differntial gene transcription.
Project description:Tuberculosis (TB), the disease caused by Mycobacterium tuberculosis (Mtb), remains a major health problem with 10.6 million cases of the disease and 1.6 million deaths in 2021. It is well understood that pulmonary TB is due to Mtb growth in the lung but quantitative estimates of rates of Mtb replication and death in lungs of patients or animals such as monkeys or rabbits remain largely unknown. We performed experiments with rabbits infected with a novel, virulent clinical Mtb isolate of the Beijing lineage, HN878, carrying an unstable plasmid pBP10. In our in vitro experiments we found that pBP10 is more stable in HN878 strain than in a more commonly used laboratory-adapted Mtb strain H37Rv (the segregation coefficient being s = 0.10 in HN878 vs. s = 0.18 in H37Rv). Interestingly, the kinetics of plasmid-bearing bacteria in lungs of Mtb-infected rabbits did not follow an expected monotonic decline; the percent of plasmid-bearing cells increased between 28 and 56 days post-infection and remained stable between 84 and 112 days post-infection despite a large increase in bacterial numbers in the lung at late time points. Mathematical modeling suggested that such a non-monotonic change in the percent of plasmid-bearing cells can be explained if the lung Mtb population consists of several (at least 2) sub-populations with different replication/death kinetics: one major population expanding early and being controlled/eliminated, while another, a smaller population expanding at later times causing a counterintuitive increase in the percent of plasmid-bearing cells. Importantly, a model with one kinetically homogeneous Mtb population could not explain the data including when the model was run stochastically. Given that in rabbits HN878 strain forms well circumscribed granulomas, our results suggest independent bacterial dynamics in subsets of such granulomas. Our model predictions can be tested in future experiments in which HN878-pBP10 dynamics in individual granulomas is followed over time. Taken together, our new data and mathematical modeling-based analyses illustrate differences in Mtb dynamics in mice and rabbits confirming a perhaps somewhat obvious observation that "rabbits are not mice".
Project description:The distribution of nanomedicines inside solid tumors is often restricted to perivascular areas, leaving most distal tumor cells out of reach. This partly explains modest patient benefit of many nanomedicines compared to their free-form counterparts. The objective for this study is to develop a mathematical model to quantitatively analyze this phenomenon and the influencing factors to such perivascular distribution and seek for effective strategies to alleviate this. A spatial tumor distribution model was firstly constructed to mimic the geometrical structure of tumor vessels and the surrounding tumor cells. This tumor model was further integrated with a systemic pharmacokinetics model for nanoparticles. A variety of factors on the tumor spatial distributions of nanomedicines were considered in the model. With the model, we quantified the effect of these influencing factors on tumor delivery efficacy (ID %), the magnitude of heterogeneous distribution (H index), and the effect of enhanced permeability and retention (EPR). In particularly, we compared the spatial distributions of the nanoparticles and the free payloads insides tumors. The model predicted high degrees of distributional heterogeneity for both nanoparticles and free payloads. The degree of heterogeneity and the influencing factors for free payloads were markedly different from those for nanoparticles. We found that nanoparticle diffusion coefficient was the most effective factor in reducing the nanoparticle H index but exerted moderate influence on the free payloads H index. The most effective factor in reducing the H index of free payload was payload diffusion coefficient. The factors that improved free payload distribution were closely associated with higher drug efficacy. In contrast, the factors that improved nanoparticle spatial distributions did not always confer improved anti-tumor efficacy of the delivered drug. These findings highlight the importance of assessing the heterogeneous free payload distribution in tumors for the development of effective nanomedicines.
Project description:Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. In this study, we refined the previously published mathematical model for bloodstream form parasites and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. Through this study, we conclude that lLevels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting acquisition of new control mechanisms during adaptation to mammalian parasitism.
Project description:Human mesenchymal stem cell (hMSC) delivery has demonstrated promise in preclinical and clinical trials for myocardial infarction therapy; however, broad acceptance is hindered by limited understanding of hMSC-human cardiomyocyte (hCM) interactions. To better understand the electrophysiological consequences of direct heterocellular connections between hMSCs and hCMs, three original mathematical models were developed, representing an experimentally verified triad of hMSC families with distinct functional ion channel currents. The arrhythmogenic risk of such direct electrical interactions in the setting of healthy adult myocardium was predicted by coupling and fusing these hMSC models to the published ten Tusscher midcardial hCM model. Substantial variations in action potential waveform-such as decreased action potential duration (APD) and plateau height-were found when hCMs were coupled to the two hMSC models expressing functional delayed rectifier-like human ether à-go-go K+ channel 1 (hEAG1); the effects were exacerbated for fused hMSC-hCM hybrid cells. The third family of hMSCs (Type C), absent of hEAG1 activity, led to smaller single-cell action potential alterations during coupling and fusion, translating to longer tissue-level mean action potential wavelength. In a simulated 2-D monolayer of cardiac tissue, re-entry vulnerability with low (5%) hMSC insertion was approximately eight-fold lower with Type C hMSCs compared to hEAG1-functional hMSCs. A 20% decrease in APD dispersion by Type C hMSCs compared to hEAG1-active hMSCs supports the claim of reduced arrhythmogenic potential of this cell type with low hMSC insertion. However, at moderate (15%) and high (25%) hMSC insertion, the vulnerable window increased independent of hMSC type. In summary, this study provides novel electrophysiological models of hMSCs, predicts possible arrhythmogenic effects of hMSCs when directly coupled to healthy hCMs, and proposes that isolating a subset of hMSCs absent of hEAG1 activity may offer increased safety as a cell delivery cardiotherapy at low levels of hMSC-hCM coupling.
Project description:Limited access to human islets has prompted the development of human beta cell models. The human beta cell lines EndoC-βH1 and EndoC-βH2 are increasingly used by the research community. However, little is known of their electrophysiological and secretory properties. Here, we monitored parameters that constitute the glucose-triggering pathway of insulin release. Both cell lines respond to glucose (6 and 20 mM) with 2- to 3-fold stimulation of insulin secretion which correlated with an elevation of [Ca2+]i, membrane depolarisation and increased action potential firing. Similar to human primary beta cells, KATP channel activity is low at 1 mM glucose and is further reduced upon increasing glucose concentration; an effect that was mimicked by the KATP channel blocker tolbutamide. The upstroke of the action potentials reflects the activation of Ca2+ channels with some small contribution of TTX-sensitive Na+ channels. The repolarisation involves activation of voltage-gated Kv2.2 channels and large-conductance Ca2+-activated K+ channels. Exocytosis presented a similar kinetics to human primary beta cells. The ultrastructure of these cells shows insulin vesicles composed of an electron-dense core surrounded by a thin clear halo. We conclude that the EndoC-βH1 and -βH2 cells share many features of primary human β-cells and thus represent a useful experimental model.
Project description:How viral infections develop can change based on the number of viruses initially entering the body. The understanding of the impacts of infection doses remains incomplete, in part due to challenging constraints, and a lack of research. Gaining more insights is crucial regarding the measles virus (MV). The higher the MV infection dose, the earlier the peak of acute viremia, but the magnitude of the peak viremia remains almost constant. Measles is highly contagious, causes immunosuppression such as lymphopenia, and contributes substantially to childhood morbidity and mortality. This work investigated mechanisms underlying the observed wild-type measles infection dose responses in cynomolgus monkeys. We fitted longitudinal data on viremia using maximum likelihood estimation, and used the Akaike Information Criterion (AIC) to evaluate relevant biological hypotheses and their respective model parameterizations. The lowest AIC indicates a linear relationship between the infection dose, the initial viral load, and the initial number of activated MV-specific T cells. Early peak viremia is associated with high initial number of activated MV-specific T cells. Thus, when MV infection dose increases, the initial viremia and associated immune cell stimulation increase, and reduce the time it takes for T cell killing to be sufficient, thereby allowing dose-independent peaks for viremia, MV-specific T cells, and lymphocyte depletion. Together, these results suggest that the development of measles depends on virus-host interactions at the start and the efficiency of viral control by cellular immunity. These relationships are additional motivations for prevention, vaccination, and early treatment for measles.
Project description:Tuberculosis ( TB ), the disease caused by Mycobacterium tuberculosis ( Mtb ), remains a major health problem with 10.6 million cases of the disease and 1.6 million deaths in 2021. It is well understood that pulmonary TB is due to replication of Mtb in the lung but quantitative details of Mtb replication and death in lungs of patients and how these rates are related to the degree of lung pathology are unknown. We performed experiments with rabbits infected with a novel, virulent clinical Mtb isolate of the Beijing lineage, HN878, carrying an unstable plasmid pBP10. In our in vitro experiments we found that pBP10 is more stable in HN878 strain than in a more commonly used laboratory-adapted Mtb strain H37Rv (the segregation coefficient being s = 0.10 in HN878 vs. s = 0.18 in H37Rv). Interestingly, the kinetics of plasmid-bearing bacteria in lungs of Mtb-infected rabbits did not follow an expected monotonic decline; the percent of plasmid-bearing cells increased between 28 and 56 days post-infection and remained stable between 84 and 112 days post-infection despite a large increase in bacterial numbers in the lung at late time points. Mathematical modeling suggested that such a non-monotonic change in the percent of plasmid-bearing cells can be explained if the lung Mtb population consists of several (at least 2) sub-populations with different replication/death kinetics: one major population expanding early and being controlled/eliminated, while another, a smaller population expanding at later times causing a counterintuitive increase in the percent of plasmid-bearing cells. Given that HN878 forms well circumscribed granulomas in rabbits, our results suggest independent bacterial dynamics in subsets of such granulomas. Our model predictions can be tested in future experiments in which HN878-pBP10 dynamics in individual granulomas is followed over time.
Project description:Aneuploidy is the leading contributor to pregnancy loss, congenital anomalies, and in vitro fertilization (IVF) failure in humans. Although most aneuploid conceptions are thought to originate from meiotic division errors in the female germline, quantitative studies that link the observed phenotypes to underlying error mechanisms are lacking. In this study, we developed a mathematical modeling framework to quantify the contribution of different mechanisms of erroneous chromosome segregation to the production of aneuploid eggs. Our model considers the probabilities of all possible chromosome gain/loss outcomes that arise from meiotic errors, such as nondisjunction (NDJ) in meiosis I and meiosis II, and premature separation of sister chromatids (PSSC) and reverse segregation (RS) in meiosis I. To understand the contributions of different meiotic errors, we fit our model to aneuploidy data from 11,157 blastocyst-stage embryos. Our best-fitting model captures several known features of female meiosis, for instance, the maternal age effect on PSSC. More importantly, our model reveals previously undescribed patterns, including an increased frequency of meiosis II errors among eggs affected by errors in meiosis I. This observation suggests that the occurrence of NDJ in meiosis II is associated with the ploidy status of an egg. We further demonstrate that the model can be used to identify IVF patients who produce an extreme number of aneuploid embryos. The dynamic nature of our mathematical model makes it a powerful tool both for understanding the relative contributions of mechanisms of chromosome missegregation in human female meiosis and for predicting the outcomes of assisted reproduction.
Project description:In contrast to heart failure (HF) with reduced ejection fraction (HFrEF), effective interventions for HF with preserved ejection fraction (HFpEF) have proven elusive, in part because it is a heterogeneous syndrome with incompletely understood pathophysiology. This study utilized mathematical modeling to evaluate mechanisms distinguishing HFpEF and HFrEF. HF was defined as a state of chronically elevated left ventricle end diastolic pressure (LVEDP > 20mmHg). First, using a previously developed cardiorenal model, sensitivities of LVEDP to potential contributing mechanisms of HFpEF, including increased myocardial, arterial, or venous stiffness, slowed ventricular relaxation, reduced LV contractility, hypertension, or reduced venous capacitance, were evaluated. Elevated LV stiffness was identified as the most sensitive factor. Large LV stiffness increases alone, or milder increases combined with either decreased LV contractility, increased arterial stiffness, or hypertension, could increase LVEDP into the HF range without reducing EF. We then evaluated effects of these mechanisms on mechanical signals of cardiac outward remodeling, and tested the ability to maintain stable EF (as opposed to progressive EF decline) under two remodeling assumptions: LV passive stress-driven vs. strain-driven remodeling. While elevated LV stiffness increased LVEDP and LV wall stress, it mitigated wall strain rise for a given LVEDP. This suggests that if LV strain drives outward remodeling, a stiffer myocardium will experience less strain and less outward dilatation when additional factors such as impaired contractility, hypertension, or arterial stiffening exacerbate LVEDP, allowing EF to remain normal even at high filling pressures. Thus, HFpEF heterogeneity may result from a range of different pathologic mechanisms occurring in an already stiffened myocardium. Together, these simulations further support LV stiffening as a critical mechanism contributing to elevated cardiac filling pressures; support LV passive strain as the outward dilatation signal; offer an explanation for HFpEF heterogeneity; and provide a mechanistic explanation distinguishing between HFpEF and HFrEF.