Drug tolerance and persistence in bacteria, fungi and cancer cells: Role of non-genetic heterogeneity.

A common feature of bacterial, fungal and cancer cell populations upon treatment is the presence of tolerant and persistent cells able to survive, and sometimes grow, even in the presence of usually inhibitory or lethal drug concentrations, driven by non-genetic differences among individual cells in a population. Here we review and compare data obtained on drug survival in bacteria, fungi and cancer cells to unravel common characteristics and cellular pathways, and to point their singularities. This comparative work also allows to cross-fertilize ideas across fields. We particularly focus on the role of gene expression variability in the emergence of cell-cell non-genetic heterogeneity because it represents a possible common basic molecular process at the origin of most persistence phenomena and could be monitored and tuned to help improve therapeutic interventions.

Cell-state transitions and density-dependent interactions together explain the dynamics of spontaneous epithelial-mesenchymal heterogeneity.

Cancer cell populations comprise phenotypes distributed among the epithelial-mesenchymal (E-M) spectrum. However, it remains unclear which population-level processes give rise to the observed experimental distribution and dynamical changes in E-M heterogeneity, including (1) differential growth, (2) cell-state switching, and (3) population density-dependent growth or state-transition rates. Here, we analyze the necessity of these three processes in explaining the dynamics of E-M population distributions as observed in PMC42-LA and HCC38 breast cancer cells. We find that, while cell-state transition is necessary to reproduce experimental observations of dynamical changes in E-M fractions, including density-dependent growth interactions (cooperation or suppression) better explains the data. Further, our models predict that treatment of HCC38 cells with transforming growth factor ß (TGF-ß) signaling and Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/3) inhibitors enhances the rate of mesenchymal-epithelial transition (MET) instead of lowering that of E-M transition (EMT). Overall, our study identifies the population-level processes shaping the dynamics of spontaneous E-M heterogeneity in breast cancer cells.

Prevalence of Epstein Barr Virus and Herpes Simplex Virus Among Human Papillomavirus Negative Oral Cancer Patients: A Cross-Sectional Study from South India.

The high incidence of oral carcinomas is due to its multifactorial etiology and the presence of various risk factors. Human Papillomavirus (HPV) has a proven role in the pathogenesis of oral carcinomas, but in the recent times there has been an increasing incidence of oral cancers who are negative for HPV infection. Also, these patients are non-smokers and non-drinkers so it could be speculated that these oral cancers are due to some other etiological factor probably of other viral infections. Therefore, this study examined the prevalence of Epstein Barr Virus (EBV) and Herpes Simplex Virus (HSV) among oral cancer patients. This cross-sectional study was conducted from January 2019 to June 2020. Biopsy samples from 47 newly diagnosed untreated patients with oral malignancies were collected along with their demographic and clinicopathological information. DNA extracted from the biopsies was processed for nested PCR for the detection of EBV and HSV. All the samples tested negative for HPV and HSV infection. Nested PCR detected 29 cases (70.7%) to be positive for EBV. The non-cancerous adjacent tissues also were negative for HPV, EBV and HSV. The prevalence of EBV was found to be more in males (62.1%) and the highest number of cases was of the left buccal mucosa compromising 34% of the total cases. From the present study it can be concluded that EBV but not HSV infection is associated with an increased risk of developing oral cancers. Although, 70.7% of the patients were found to be positive for EBV whether the viral infection played any role in the driving the malignancy needs to be further elucidated.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity.

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Epigenetic memory acquired during long-term EMT induction governs the recovery to the epithelial state.

Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFß need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFß withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.

Stochastic population dynamics of cancer stemness and adaptive response to therapies.

Intratumoral heterogeneity can exist along multiple axes: Cancer stem cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states, and a spectrum of epithelial-hybrid-mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple stable gene expression patterns ('attractors'). These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

Rapid, sensitive, and low-cost detection of Escherichia coli bacteria in contaminated water samples using a phage-based assay.

Inadequate drinking water quality is among the major causes of preventable mortality, predominantly in young children. Identifying contaminated water sources remains a significant challenge, especially where resources are limited. The current methods for measuring Escherichia coli (E. coli), the WHO preferred indicator for measuring fecal contamination of water, involve overnight incubation and require specialized training. In 2016, UNICEF released a Target Product Profile (TPP) to incentivize product innovations to detect low levels of viable E. coli in water samples in the field in less than 6 h. Driven by this challenge, we developed a phage-based assay to detect and semi-quantify E. coli. We formulated a phage cocktail containing a total of 8 phages selected against an extensive bacterial strain library and recombined with the sensitive NanoLuc luciferase reporter. The assay was optimized to be processed in a microfluidic chip designed in-house and was tested against locally sourced sewage samples and on drinking water sources in Nairobi, Kenya. With this assay, combined with the microfluidic chip platform, we propose a complete automated solution to detect and semi-quantify E. coli at less than 10 MPN/100 mL in 5.5 h by minimally trained personnel.

Population Dynamics of Epithelial-Mesenchymal Heterogeneity in Cancer Cells.

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAMhigh subpopulation, (b) re-establishment of parental distributions from the EpCAMhigh and EpCAMlow subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL-an EMT-inducing transcription factor-during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.

A student centric method for calculation of fatty acid energetics: Integrated formula and web tool.

Mitochondrial beta-oxidation is one of the most common modes of fatty acids' oxidation in most organisms, particularly mammals. Biochemistry undergraduate curriculum often contains the description of the process, with emphasis on ATP calculations for various types of fatty acids. During our decade long teaching experience in biochemistry, we observed the difficulty faced by students in calculating energetics of several fatty acids beyond palmitic acid. We developed a canonical formula by mathematical transformations and logical derivation to aid the calculation in a much simpler manner to ease both teaching and learning experience. ATP yield of even-numbered fatty acids may be calculated using a formula [(7C - 6 - 1.5 D) - 2(D-2)], andadenosine triphosphate (ATP) yield for odd-numbered fatty acids can be calculated using [(7C - 19 - 1.5 D) - 2(D-2)], where C is the number of carbon atoms in fatty acids, D is the number of double bonds. The unbold part of the formulae is limited to polyunsaturated fatty acids. Moreover, we integrated these formulae into an HTML based web-interface for handily calculations, which is likely to augment fatty acids oxidation learning-teaching processes easier. This tool has been recently tested in our classroom programs on biochemistry and received an excellent feedback from the learners. Also, the mathematical formula is ready for being incorporated into standard biochemistry textbooks. The webtool as an opensource biochemical calculator can be effectively used in classrooms by both instructors and students while solving comprehension based questions on lipid metabolism.

Differential seminal plasma proteome signatures of acute lymphoblastic leukemia survivors.

With advances in therapeutic methods, there is a high survival rate among leukemia patients, of an extent more than 80%. However, chemotherapeutic drugs used to treat these patients have adverse effects on their overall health profile including fertility. The primary aim of this study was to identify differentially expressed proteins in seminal plasma of acute lymphoblastic leukemia (ALL) survivors compared to age-matched healthy controls, which can provide molecular basis of idiopathic infertility in such survivors. Differential proteome profiling was performed by 2D-differential in-gel electrophoresis, protein spots were identified by mass spectrometry and selective differentially expressed proteins (DEPs) were validated by western blotting and ELISA method. Out of eight DEPs identified, five proteins (isocitrate dehydrogenase 1, semenogelin 1, lactoferrin, prolactin-inducible protein, and human serum albumin) were upregulated and three (pepsinogen, prostate specific antigen and prostatic acid phosphatase) were downregulated. Expression profiles of these proteins are suggestive of reduction in semen quality in ALL survivors and can further be explored to determine their fertility status.

Enhanced respiration prevents drug tolerance and drug resistance in Mycobacterium tuberculosis.

Persistence, manifested as drug tolerance, represents a significant obstacle to global tuberculosis control. The bactericidal drugs isoniazid and rifampicin kill greater than 99% of exponentially growing Mycobacterium tuberculosis (Mtb) cells, but the remaining cells are persisters, cells with decreased metabolic rate, refractory to killing by these drugs, and able to generate drug-resistant mutants. We discovered that the combination of cysteine or other small thiols with either isoniazid or rifampicin prevents the formation of drug-tolerant and drug-resistant cells in Mtb cultures. This effect was concentration- and time-dependent, relying on increased oxygen consumption that triggered enhanced production of reactive oxygen species. In infected murine macrophages, the addition of N-acetylcysteine to isoniazid treatment potentiated the killing of Mtb Furthermore, we demonstrate that the addition of small thiols to Mtb drug treatment shifted the menaquinol/menaquinone balance toward a reduced state that stimulates Mtb respiration and converts persister cells to metabolically active cells. This prevention of both persister cell formation and drug resistance leads ultimately to mycobacterial cell death. Strategies to enhance respiration and initiate oxidative damage should improve tuberculosis chemotherapies.

Characterization of DNA topoisomerase I from Mycobacterium tuberculosis: DNA cleavage and religation properties and inhibition of its activity.

Type I DNA topoisomerases from bacteria catalyse relaxation of negatively supercoiled DNA in a Mg(2+) dependent manner. Although topoisomerases of distinct classes have been subjected for anti-cancer and anti-infective drug development, bacterial type I enzymes are way behind in this regard. Our studies with Mycobacterium smegmatis topoisomerase I (MstopoI) revealed several of its distinct properties compared to the well studied Escherichia coli topoisomerase I (EctopoI) suggesting the possibility of targeting the mycobacterial enzyme for inhibitor development. Here, we describe Mycobacterium tuberculosis topoisomerase I (MttopoI) and compare its properties with MstopoI and EctopoI. The enzyme cleaves DNA at preferred sites in a pattern similar to its ortholog from M. smegmatis. Oligonucleotides containing the specific recognition sequence inhibited the activity of the enzyme in a manner similar to that of MstopoI. Substitution of the acidic residues, D111 and E115 which are involved in Mg(2+) co-ordination, to alanines affected the DNA relaxation activity. Unlike the wild type enzyme, D111A was dependent on Mg(2+) for DNA cleavage and both the mutants were compromised in religation. The monoclonal antibody (mAb), 2F3G4, developed against MstopoI inhibited the relaxation activity of MttopoI. These studies affirm the characteristics of MttopoI to be similar to MstopoI and set a stage to target it for the development of specific small molecule inhibitors.

Indispensable, functionally complementing N and C-terminal domains constitute site-specific topoisomerase I.

Mycobacterium smegmatis topoisomerase I differs from the typical type IA topoisomerase in many properties. The enzyme recognizes both single and double-stranded DNA with high affinity and makes sequence-specific contacts during DNA relaxation reaction. The enzyme has a conserved N-terminal domain and a highly varied C-terminal domain, which lacks the characteristic zinc binding motifs found in most of the type I eubacterial enzymes. The roles of the individual domains of the enzyme in the topoisomerase I catalyzed reactions were examined by comparing the properties of full-length topoisomerase I with those of truncated polypeptides lacking the conserved N-terminal or the divergent C-terminal region. The N-terminal larger fragment retained the site-specific binding, DNA cleavage and religation properties, hallmark characteristics of the full-length M.smegmatis topoisomerase I. In contrast, the non-conserved C-terminal fragment lacking the typical DNA binding motif, exhibited non-specific DNA binding behaviour. The two polypeptide fragments, on their own do not catalyze DNA relaxation reaction. The relaxation activity is restored when both the fragments are mixed in vitro reconstituting the enzyme function. These results along with the DNA interaction pattern of the proteins implicate an essential role for the C-terminal region in single-strand DNA passage between the two transesterification reactions catalyzed by the N-terminal domain.