An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Translocation and encapsulation of siRNA inside carbon nanotubes

We report spontaneous translocation of small interfering RNA (siRNA) inside carbon nanotubes (CNTs) of various diameters and chirality using all atom molecular dynamics simulations with explicit solvent. We use umbrella sampling method to calculate the free energy landscape of the siRNA entry and translocation event. Free energy profiles show that siRNA gains free energy while translocating inside CNT, and barrier for siRNA exit from CNT ranges from 40 to 110 kcal/mol depending on CNT chirality and salt concentration. The translocation time tau decreases with the increase of CNT diameter with a critical diameter of 24 angstrom for the translocation. In contrast, double strand DNA of the same sequence does not translocate inside CNT due to large free energy barrier for the translocation. This study helps in understanding the nucleic acid transport through nanopores at microscopic level and may help designing carbon nanotube based sensor for siRNA. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4773302]

Chromosomal translocations among the healthy human population: implications in oncogenesis

Chromosomal translocations are characteristic features of many cancers, especially lymphoma and leukemia. However, recent reports suggest that many chromosomal translocations can be found in healthy individuals, although the significance of this observation is still not clear. In this review, we summarize recent studies on chromosomal translocations in healthy individuals carried out in different geographical areas of the world and discuss the relevance of the observation with respect to oncogenesis.

On finding a set of healthy individuals from a large population

In this paper, we explore fundamental limits on the number of tests required to identify a given number of ``healthy'' items from a large population containing a small number of ``defective'' items, in a nonadaptive group testing framework. Specifically, we derive mutual information-based upper bounds on the number of tests required to identify the required number of healthy items. Our results show that an impressive reduction in the number of tests is achievable compared to the conventional approach of using classical group testing to first identify the defective items and then pick the required number of healthy items from the complement set. For example, to identify L healthy items out of a population of N items containing K defective items, when the tests are reliable, our results show that O(K(L - 1)/(N - K)) measurements are sufficient. In contrast, the conventional approach requires O(K log(N/K)) measurements. We derive our results in a general sparse signal setup, and hence, they are applicable to other sparse signal-based applications such as compressive sensing also.

Flying between Sky Islands: The Effect of Naturally Fragmented Habitat on Butterfly Population Structure

High elevation montane areas are called ``sky islands'' when they occur as a series of high mountains separated by lowland valleys. Different climatic conditions at high elevations makes sky islands a specialized type of habitat, rendering them naturally fragmented compared to more continuous habitat at lower elevations. Species in sky islands face unsuitable climate in the intervening valleys when moving from one montane area to another. The high elevation shola-grassland mosaic in the Western Ghats of southern India form one such sky island complex. The fragmented patches make this area ideal to study the effect of the spatial orientation of suitable habitat patches on population genetic structure of species found in these areas. Past studies have suggested that sky islands tend to have genetically structured populations, possibly due to reduced gene flow between montane areas. To test this hypothesis, we adopted the comparative approach. Using Amplified Fragment Length Polymorphisms, we compared population genetic structures of two closely related, similar sized butterfly species: Heteropsis oculus, a high elevation shola-grassland specialist restricted to the southern Western Ghats, and Mycalesis patnia, found more continuously distributed in lower elevations. In all analyses, as per expectation the sky island specialist H. oculus exhibited a greater degree of population genetic structure than M. patnia, implying a difference in geneflow. This difference in geneflow in turn appears to be due to the natural fragmentation of the sky island complexes. Detailed analysis of a subset of H. oculus samples from one sky island complex (the Anamalais) showed a surprising genetic break. A possible reason for this break could be unsuitable conditions of higher temperature and lower rainfall in the intervening valley region. Thus, sky island species are not only restricted by lack of habitat continuity between montane areas, but also by the nature of the intervening habitat.

G-Quadruplex Structures Formed at the HOX11 Breakpoint Region Contribute to Its Fragility during t(10;14) Translocation in T-Cell Leukemia

The t(10;14) translocation involving the HOX11 gene is found in several T-cell leukemia patients. Previous efforts to determine the causes of HOX11 fragility were not successful. The role of non-B DNA structures is increasingly becoming an important cause of genomic instability. In the present study, bioinformatics analysis revealed two G-quadruplex-forming motifs at the HOX11 breakpoint cluster. Gel shift assays showed formation of both intra- and intermolecular G-quadruplexes, the latter being more predominant. The structure formation was dependent on four stretches of guanines, as revealed by mutagenesis. Circular dichroism analysis identified parallel conformations for both quadruplexes. The non-B DNA structure could block polymerization during replication on a plasmid, resulting in consistent K K+-dependent pause sites, which were abolished upon mutation of G-motifs, thereby demonstrating the role of the stretches of guanines even on double-stranded DNA. Extrachromosomal assays showed that the G-quadruplex motifs could block transcription, leading to reduced expression of green fluorescent protein (GFP) within cells. More importantly, sodium bisulfite modification assay showed the single-stranded character at regions I and II of HOX11 in the genome. Thus, our findings suggest the occurrence of G-quadruplex structures at the HOX11 breakpoint region, which could explain its fragility during the t(10;14) translocation.

A landscape-level assessment of Asian elephant habitat, its population and elephant-human conflict in the Anamalai hill ranges of southern Western Ghats, India

Spatial information at the landscape scale is extremely important for conservation planning, especially in the case of long-ranging vertebrates. The biodiversity-rich Anamalai hill ranges in the Western Ghats of southern India hold a viable population for the long-term conservation of the Asian elephant. Through rapid but extensive field surveys we mapped elephant habitat, corridors, vegetation and land-use patterns, estimated the elephant population density and structure, and assessed elephant-human conflict across this landscape. GIS and remote sensing analyses indicate that elephants are distributed among three blocks over a total area of about 4600 km(2). Approximately 92% remains contiguous because of four corridors; however, under 4000 km2 of this area may be effectively used by elephants. Nine landscape elements were identified, including five natural vegetation types, of which tropical moist deciduous forest is dominant. Population density assessed through the dung count method using line transects covering 275 km of walk across the effective elephant habitat of the landscape yielded a mean density of 1.1 (95% Cl = 0.99-1.2) elephant/km(2). Population structure from direct sighting of elephants showed that adult male elephants constitute just 2.9% and adult females 42.3% of the population with the rest being subadults (27.4%), juveniles (16%) and calves (11.4%). Sex ratios show an increasing skew toward females from juvenile (1:1.8) to sub-adult (1:2.4) and adult (1:14.7) indicating higher mortality of sub-adult and adult males that is most likely due to historical poaching for ivory. A rapid questionnaire survey and secondary data on elephant-human conflict from forest department records reveals that villages in and around the forest divisions on the eastern side of landscape experience higher levels of elephant-human conflict than those on the western side; this seems to relate to a greater degree of habitat fragmentation and percentage farmers cultivating annual crops in the east. We provide several recommendations that could help maintain population viability and reduce elephant-human conflict of the Anamalai elephant landscape. (C) 2013 Deutsche Gesellschaft far Saugetierkunde. Published by Elsevier GmbH. All rights reserved.

Variation at Diabetes- and Obesity-Associated Loci May Mirror Neutral Patterns of Human Population Diversity and Diabetes Prevalence in India

South Asian populations harbor a high degree of genetic diversity, due in part to demographic history. Two studies on genome-wide variation in Indian populations have shown that most Indian populations show varying degrees of admixture between ancestral north Indian and ancestral south Indian components. As a result of this structure, genetic variation in India appears to follow a geographic cline. Similarly, Indian populations seem to show detectable differences in diabetes and obesity prevalence between different geographic regions of the country. We tested the hypothesis that genetic variation at diabetes-and obesity-associated loci may be potentially related to different genetic ancestries. We genotyped 2977 individuals from 61 populations across India for 18 SNPs in genes implicated in T2D and obesity. We examined patterns of variation in allele frequency across different geographical gradients and considered state of origin and language affiliation. Our results show that most of the 18 SNPs show no significant correlation with latitude, the geographic cline reported in previous studies, or by language family. Exceptions include KCNQ1 with latitude and THADA and JAK1 with language, which suggests that genetic variation at previously ascertained diabetes-associated loci may only partly mirror geographic patterns of genome-wide diversity in Indian populations.

Temporal variability of forest communities: empirical estimates of population change in 4000 tree species

Long-term surveys of entire communities of species are needed to measure fluctuations in natural populations and elucidate the mechanisms driving population dynamics and community assembly. We analysed changes in abundance of over 4000 tree species in 12 forests across the world over periods of 6-28years. Abundance fluctuations in all forests are large and consistent with population dynamics models in which temporal environmental variance plays a central role. At some sites we identify clear environmental drivers, such as fire and drought, that could underlie these patterns, but at other sites there is a need for further research to identify drivers. In addition, cross-site comparisons showed that abundance fluctuations were smaller at species-rich sites, consistent with the idea that stable environmental conditions promote higher diversity. Much community ecology theory emphasises demographic variance and niche stabilisation; we encourage the development of theory in which temporal environmental variance plays a central role.

Role of the loop L-4,L-5 in allosteric regulation in mtHsp70s: in vivo significance of domain communication and its implications in protein translocation

Mitochondrial Hsp70 (mtHsp70) is essential for a vast repertoire of functions, including protein import, and requires effective interdomain communication for efficient partner-protein interactions. However, the in vivo functional significance of allosteric regulation in eukaryotes is poorly defined. Using integrated biochemical and yeast genetic approaches, we provide compelling evidence that a conserved substrate-binding domain (SBD) loop, L-4,L-5, plays a critical role in allosteric communication governing mtHsp70 chaperone functions across species. In yeast, a temperature-sensitive L-4,L-5 mutation (E467A) disrupts bidirectional domain communication, leading to compromised protein import and mitochondrial function. Loop L-4,L-5 functions synergistically with the linker in modulating the allosteric interface and conformational transitions between SBD and the nucleotide-binding domain (NBD), thus regulating interdomain communication. Second-site intragenic suppressors of E467A isolated within the SBD suppress domain communication defects by conformationally altering the allosteric interface, thereby restoring import and growth phenotypes. Strikingly, the suppressor mutations highlight that restoration of communication from NBD to SBD alone is the minimum essential requirement for effective in vivo function when primed at higher basal ATPase activity, mimicking the J-protein-bound state. Together these findings provide the first mechanistic insights into critical regions within the SBD of mtHsp70s regulating interdomain communication, thus highlighting its importance in protein translocation and mitochondrial biogenesis.

Effects of Withania somnifera and Tinospora cordifolia Extracts on the Side Population Phenotype of Human Epithelial Cancer Cells: Toward Targeting Multidrug Resistance in Cancer

Recent reports suggest the existence of a subpopulation of stem-like cancer cells, termed as cancer stem cells (CSCs), which bear functional and phenotypic resemblance with the adult, tissue-resident stem cells. Side population (SP) assay based on differential efflux of Hoechst 33342 has been effectively used for the isolation of CSCs. The drug resistance properties of SP cells are typically due to the increased expression of ABC transporters leading to drug efflux. Conventionally used chemotherapeutic drugs may often leads to an enrichment of SP, revealing their inability to target the drug-resistant SP and CSCs. Thus, identification of agents that can reduce the SP phenotype is currently in vogue in cancer therapeutics. Withania somnifera (WS) and Tinospora cordifolia (TC) have been used in Ayurveda for treating various diseases, including cancer. In the current study, we have investigated the effects of ethanolic (ET) extracts of WS and TC on the cancer SP phenotype. Interestingly, we found significant decrease in SP on treatment with TC-ET, but not with WS-ET. The SP-inhibitory TC-ET was further fractionated into petroleum ether (TC-PET), dichloromethane (TC-DCM), and n-butyl alcohol (TC-nBT) fractions using bioactivity-guided fractionation. Our data revealed that TC-PET and TC-DCM, but not TC-nBT, significantly inhibited SP in a dose-dependent manner. Furthermore, flow cytometry-based functional assays revealed that TC-PET and TC-DCM significantly inhibited ABC-B1 and ABC-G2 transporters and sensitized cancer cells toward chemotherapeutic drug-mediated cytotoxicity. Thus, the TC-PET and TC-DCM may harbor phytochemicals with the potential to reverse the drug-resistant phenotype, thus improving the efficacy of cancer chemotherapy.

A comparative study of a direct discretization and an operator-splitting solver for population balance systems

A direct discretization approach and an operator-splitting scheme are applied for the numerical simulation of a population balance system which models the synthesis of urea with a uni-variate population. The problem is formulated in axisymmetric form and the setup is chosen such that a steady state is reached. Both solvers are assessed with respect to the accuracy of the results, where experimental data are used for comparison, and the efficiency of the simulations. Depending on the goal of simulations, to track the evolution of the process accurately or to reach the steady state fast, recommendations for the choice of the solver are given. (C) 2015 Elsevier Ltd. All rights reserved.

Observation of coherent population trapping in a V-type two-electron system

We observe coherent population trapping (CPT) in a two-electron atom-Yb-174-using the S-1(0), F= 0 -> P-3(1), F `= 1 transition. CPT is not possible for such a transition according to one-electron theory because the magnetic sublevels form a V-type system, but in a two-electron atom like Yb, the interaction of the electrons transforms the level structure into a V-type system, which allows the formation of a dark state and hence the observation of CPT. Since the two levels involved are degenerate, we use a magnetic field to lift the degeneracy. The single fluorescence dip then splits into five dips-the central unshifted one corresponds to coherent population oscillation, while the outer four are due to CPT. The linewidth of the CPT resonance is about 300 kHz and is limited by the natural linewidth of the excited state, which is to be expected because the excited state is involved in the formation of the dark state.

REPERTOIRE OF PROTEIN KINASES ENCODED IN THE GENOME OF ZEBRAFISH SHOWS REMARKABLY LARGE POPULATION OF PIM KINASES

In recent times, zebrafish has garnered lot of popularity as model organism to study human cancers. Despite high evolutionary divergence from humans, zebrafish develops almost all types of human tumors when induced. However, mechanistic details of tumor formation have remained largely unknown. Present study is aimed at analysis of repertoire of kinases in zebrafish proteome to provide insights into various cellular components. Annotation using highly sensitive remote homology detection methods revealed ``substantial expansion'' of Ser/Thr/Tyr kinase family in zebrafish compared to humans, constituting over 3% of proteome. Subsequent classification of kinases into subfamilies revealed presence of large number of CAMK group of kinases, with massive representation of PIM kinases, important for cell cycle regulation and growth. Extensive sequence comparison between human and zebrafish PIM kinases revealed high conservation of functionally important residues with a few organism specific variations. There are about 300 PIM kinases in zebrafish kinome, while human genome codes for only about 500 kinases altogether. PIM kinases have been implicated in various human cancers and are currently being targeted to explore their therapeutic potentials. Hence, in depth analysis of PIM kinases in zebrafish has opened up new avenues of research to verify the model organism status of zebrafish.