Scroll-Wave Dynamics in Human Cardiac Tissue: Lessons from a Mathematical Model with Inhomogeneities and Fiber Architecture

Cardiac arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), are among the leading causes of death in the industrialized world. These are associated with the formation of spiral and scroll waves of electrical activation in cardiac tissue; single spiral and scroll waves are believed to be associated with VT whereas their turbulent analogs are associated with VF. Thus, the study of these waves is an important biophysical problem. We present a systematic study of the combined effects of muscle-fiber rotation and inhomogeneities on scroll-wave dynamics in the TNNP (ten Tusscher Noble Noble Panfilov) model for human cardiac tissue. In particular, we use the three-dimensional TNNP model with fiber rotation and consider both conduction and ionic inhomogeneities. We find that, in addition to displaying a sensitive dependence on the positions, sizes, and types of inhomogeneities, scroll-wave dynamics also depends delicately upon the degree of fiber rotation. We find that the tendency of scroll waves to anchor to cylindrical conduction inhomogeneities increases with the radius of the inhomogeneity. Furthermore, the filament of the scroll wave can exhibit drift or meandering, transmural bending, twisting, and break-up. If the scroll-wave filament exhibits weak meandering, then there is a fine balance between the anchoring of this wave at the inhomogeneity and a disruption of wave-pinning by fiber rotation. If this filament displays strong meandering, then again the anchoring is suppressed by fiber rotation; also, the scroll wave can be eliminated from most of the layers only to be regenerated by a seed wave. Ionic inhomogeneities can also lead to an anchoring of the scroll wave; scroll waves can now enter the region inside an ionic inhomogeneity and can display a coexistence of spatiotemporal chaos and quasi-periodic behavior in different parts of the simulation domain. We discuss the experimental implications of our study.

A Ten-microRNA Expression Signature Predicts Survival in Glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n = 222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio HR] = 2.4; 95% CI = 1.4-3.8; p < 0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR = 1.7; 95% CI = 1.1-2.8; p = 0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR = 2.0; 95% CI = 1.4-2.8; p < 0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR = 1.120; 95% CI = 1.04-1.20; p = 0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.

Angelman Syndrome Protein UBE3A Interacts with Primary Microcephaly Protein ASPM, Localizes to Centrosomes and Regulates Chromosome Segregation

Many proteins associated with the phenotype microcephaly have been localized to the centrosome or linked to it functionally. All the seven autosomal recessive primary microcephaly (MCPH) proteins localize at the centrosome. Microcephalic osteodysplastic primordial dwarfism type II protein PCNT and Seckel syndrome (also characterized by severe microcephaly) protein ATR are also centrosomal proteins. All of the above findings show the importance of centrosomal proteins as the key players in neurogenesis and brain development. However, the exact mechanism as to how the loss-of-function of these proteins leads to microcephaly remains to be elucidated. To gain insight into the function of the most commonly mutated MCPH gene ASPM, we used the yeast two-hybrid technique to screen a human fetal brain cDNA library with an ASPM bait. The analysis identified Angelman syndrome gene product UBE3A as an ASPM interactor. Like ASPM, UBE3A also localizes to the centrosome. The identification of UBE3A as an ASPM interactor is not surprising as more than 80% of Angelman syndrome patients have microcephaly. However, unlike in MCPH, microcephaly is postnatal in Angelman syndrome patients. Our results show that UBE3A is a cell cycle regulated protein and its level peaks in mitosis. The shRNA knockdown of UBE3A in HEK293 cells led to many mitotic abnormalities including chromosome missegregation, abnormal cytokinesis and apoptosis. Thus our study links Angelman syndrome protein UBE3A to ASPM, centrosome and mitosis for the first time. We suggest that a defective chromosome segregation mechanism is responsible for the development of microcephaly in Angelman syndrome.

Estimating the Threshold Surface Density of Gp120-CCR5 Complexes Necessary for HIV-1 Envelope-Mediated Cell-Cell Fusion

Reduced expression of CCR5 on target CD4(+) cells lowers their susceptibility to infection by R5-tropic HIV-1, potentially preventing transmission of infection and delaying disease progression. Binding of the HIV-1 envelope (Env) protein gp120 with CCR5 is essential for the entry of R5 viruses into target cells. The threshold surface density of gp120-CCR5 complexes that enables HIV-1 entry remains poorly estimated. We constructed a mathematical model that mimics Env-mediated cell-cell fusion assays, where target CD4(+)CCR5(+) cells are exposed to effector cells expressing Env in the presence of a coreceptor antagonist and the fraction of target cells fused with effector cells is measured. Our model employs a reaction network-based approach to describe protein interactions that precede viral entry coupled with the ternary complex model to quantify the allosteric interactions of the coreceptor antagonist and predicts the fraction of target cells fused. By fitting model predictions to published data of cell-cell fusion in the presence of the CCR5 antagonist vicriviroc, we estimated the threshold surface density of gp120-CCR5 complexes for cell-cell fusion as similar to 20 mu m(-2). Model predictions with this threshold captured data from independent cell-cell fusion assays in the presence of vicriviroc and rapamycin, a drug that modulates CCR5 expression, as well as assays in the presence of maraviroc, another CCR5 antagonist, using sixteen different Env clones derived from transmitted or early founder viruses. Our estimate of the threshold surface density of gp120-CCR5 complexes necessary for HIV-1 entry thus appears robust and may have implications for optimizing treatment with coreceptor antagonists, understanding the non-pathogenic infection of non-human primates, and designing vaccines that suppress the availability of target CD4(+)CCR5(+) cells.

Gene Expression Profiling of Preovulatory Follicle in the Buffalo Cow: Effects of Increased IGF-I Concentration on Periovulatory Events

The preovulatory follicle in response to gonadotropin surge undergoes dramatic biochemical, and morphological changes orchestrated by expression changes in hundreds of genes. Employing well characterized bovine preovulatory follicle model, granulosa cells (GCs) and follicle wall were collected from the preovulatory follicle before, 1, 10 and 22 h post peak LH surge. Microarray analysis performed on GCs revealed that 450 and 111 genes were differentially expressed at 1 and 22 h post peak LH surge, respectively. For validation, qPCR and immunocytochemistry analyses were carried out for some of the differentially expressed genes. Expression analysis of many of these genes showed distinct expression patterns in GCs and the follicle wall. To study molecular functions and genetic networks, microarray data was analyzed using Ingenuity Pathway Analysis which revealed majority of the differentially expressed genes to cluster within processes like steroidogenesis, cell survival and cell differentiation. In the ovarian follicle, IGF-I is established to be an important regulator of the above mentioned molecular functions. Thus, further experiments were conducted to verify the effects of increased intrafollicular IGF-I levels on the expression of genes associated with the above mentioned processes. For this purpose, buffalo cows were administered with exogenous bGH to transiently increase circulating and intrafollicular concentrations of IGF-I. The results indicated that increased intrafollicular concentrations of IGF-I caused changes in expression of genes associated with steroidogenesis (StAR, SRF) and apoptosis (BCL-2, FKHR, PAWR). These results taken together suggest that onset of gonadotropin surge triggers activation of various biological pathways and that the effects of growth factors and peptides on gonadotropin actions could be examined during preovulatory follicle development.

Insights into the Regulatory Characteristics of the Mycobacterial Dephosphocoenzyme A Kinase: Implications for the Universal CoA Biosynthesis Pathway

Being vastly different from the human counterpart, we suggest that the last enzyme of the Mycobacterium tuberculosis Coenzyme A biosynthetic pathway, dephosphocoenzyme A kinase (CoaE) could be a good anti-tubercular target. Here we describe detailed investigations into the regulatory features of the enzyme, affected via two mechanisms. Enzymatic activity is regulated by CTP which strongly binds the enzyme at a site overlapping that of the leading substrate, dephosphocoenzyme A (DCoA), thereby obscuring the binding site and limiting catalysis. The organism has evolved a second layer of regulation by employing a dynamic equilibrium between the trimeric and monomeric forms of CoaE as a means of regulating the effective concentration of active enzyme. We show that the monomer is the active form of the enzyme and the interplay between the regulator, CTP and the substrate, DCoA, affects enzymatic activity. Detailed kinetic data have been corroborated by size exclusion chromatography, dynamic light scattering, glutaraldehyde crosslinking, limited proteolysis and fluorescence investigations on the enzyme all of which corroborate the effects of the ligands on the enzyme oligomeric status and activity. Cysteine mutagenesis and the effects of reducing agents on mycobacterial CoaE oligomerization further validate that the latter is not cysteine-mediated or reduction-sensitive. These studies thus shed light on the novel regulatory features employed to regulate metabolite flow through the last step of a critical biosynthetic pathway by keeping the latter catalytically dormant till the need arises, the transition to the active form affected by a delicate crosstalk between an essential cellular metabolite (CTP) and the precursor to the pathway end-product (DCoA).

A Transcript Cleavage Factor of Mycobacterium tuberculosis Important for Its Survival

After initiation of transcription, a number of proteins participate during elongation and termination modifying the properties of the RNA polymerase (RNAP). Gre factors are one such group conserved across bacteria. They regulate transcription by projecting their N-terminal coiled-coil domain into the active center of RNAP through the secondary channel and stimulating hydrolysis of the newly synthesized RNA in backtracked elongation complexes. Rv1080c is a putative gre factor (MtbGre) in the genome of Mycobacterium tuberculosis. The protein enhanced the efficiency of promoter clearance by lowering abortive transcription and also rescued arrested and paused elongation complexes on the GC rich mycobacterial template. Although MtbGre is similar in domain organization and shares key residues for catalysis and RNAP interaction with the Gre factors of Escherichia coli, it could not complement an E. coli gre deficient strain. Moreover, MtbGre failed to rescue E. coli RNAP stalled elongation complexes, indicating the importance of specific protein-protein interactions for transcript cleavage. Decrease in the level of MtbGre reduced the bacterial survival by several fold indicating its essential role in mycobacteria. Another Gre homolog, Rv3788 was not functional in transcript cleavage activity indicating that a single Gre is sufficient for efficient transcription of the M. tuberculosis genome.

DNA analysis indicates that Asian elephants are native to Borneo and are therefore a high priority for conservation

The origin of Borneo's elephants is controversial. Two competing hypotheses argue that they are either indigenous, tracing back to the Pleistocene, or were introduced, descending from elephants imported in the 16th-18th centuries. Taxonomically, they have either been classified as a unique subspecies or placed under the Indian or Sumatran subspecies. If shown to be a unique indigenous population, this would extend the natural species range of the Asian elephant by 1300 km, and therefore Borneo elephants would have much greater conservation importance than if they were a feral population. We compared DNA of Borneo elephants to that of elephants from across the range of the Asian elephant, using a fragment of mitochondrial DNA, including part of the hypervariable d-loop, and five autosomal microsatellite loci. We find that Borneo's elephants are genetically distinct, with molecular divergence indicative of a Pleistocene colonisation of Borneo and subsequent isolation. We reject the hypothesis that Borneo's elephants were introduced. The genetic divergence of Borneo elephants warrants their recognition as a separate evolutionary significant unit. Thus, interbreeding Borneo elephants with those from other populations would be contraindicated in ex situ conservation, and their genetic distinctiveness makes them one of the highest priority populations for Asian elephant conservation.

Nature's Swiss Army Knives: Ovipositor Structure Mirrors Ecology in a Multitrophic Fig Wasp Community

Background: Resource partitioning is facilitated by adaptations along niche dimensions that range from morphology to behaviour. The exploitation of hidden resources may require specially adapted morphological or sensory tools for resource location and utilisation. Differences in tool diversity and complexity can determine not only how many species can utilize these hidden resources but also how they do so. Methodology and Principal Findings: The sclerotisation, gross morphology and ultrastructure of the ovipositors of a seven-member community of parasitic wasps comprising of gallers and parasitoids developing within the globular syconia (closed inflorescences) of Ficus racemosa (Moraceae) was investigated. These wasps also differ in their parasitism mode (external versus internal oviposition) and their timing of oviposition into the expanding syconium during its development. The number and diversity of sensilla, as well as ovipositor teeth, increased from internally ovipositing to externally ovipositing species and from gallers to parasitoids. The extent of sclerotisation of the ovipositor tip matched the force required to penetrate the syconium at the time of oviposition of each species. The internally ovipositing pollinator had only one type of sensillum and a single notch on the ovipositor tip. Externally ovipositing species had multiple sensilla types and teeth on their ovipositors. Chemosensilla were most concentrated at ovipositor tips while mechanoreceptors were more widely distributed, facilitating the precise location of hidden hosts in these wasps which lack larval host-seeking behaviour. Ovipositor traits of one parasitoid differed from those of its syntopic galler congeners and clustered with those of parasitoids within a different wasp subfamily. Thus ovipositor tools can show lability based on adaptive necessity, and are not constrained by phylogeny. Conclusions/Significance: Ovipositor structure mirrored the increasingly complex trophic ecology and requirements for host accessibility in this parasite community. Ovipositor structure could be a useful surrogate for predicting the biology of parasites in other communities.

Modulation of Catalytic Activity in Multi-Domain Protein Tyrosine Phosphatases

Signaling mechanisms involving protein tyrosine phosphatases govern several cellular and developmental processes. These enzymes are regulated by several mechanisms which include variation in the catalytic turnover rate based on redox stimuli, subcellular localization or protein-protein interactions. In the case of Receptor Protein Tyrosine Phosphatases (RPTPs) containing two PTP domains, phosphatase activity is localized in their membrane-proximal (D1) domains, while the membrane-distal (D2) domain is believed to play a modulatory role. Here we report our analysis of the influence of the D2 domain on the catalytic activity and substrate specificity of the D1 domain using two Drosophila melanogaster RPTPs as a model system. Biochemical studies reveal contrasting roles for the D2 domain of Drosophila Leukocyte antigen Related (DLAR) and Protein Tyrosine Phosphatase on Drosophila chromosome band 99A (PTP99A). While D2 lowers the catalytic activity of the D1 domain in DLAR, the D2 domain of PTP99A leads to an increase in the catalytic activity of its D1 domain. Substrate specificity, on the other hand, is cumulative, whereby the individual specificities of the D1 and D2 domains contribute to the substrate specificity of these two-domain enzymes. Molecular dynamics simulations on structural models of DLAR and PTP99A reveal a conformational rationale for the experimental observations. These studies reveal that concerted structural changes mediate inter-domain communication resulting in either inhibitory or activating effects of the membrane distal PTP domain on the catalytic activity of the membrane proximal PTP domain.

Clinical Proteomics of the Neglected Human Malarial Parasite Plasmodium vivax

Recent reports highlight the severity and the morbidity of disease caused by the long neglected malaria parasite Plasmodium vivax. Due to inherent difficulties in the laboratory-propagation of P. vivax, the biology of this parasite has not been adequately explored. While the proteome of P. falciparum, the causative agent of cerebral malaria, has been extensively explored from several sources, there is limited information on the proteome of P. vivax. We have, for the first time, examined the proteome of P. vivax isolated directly from patients without adaptation to laboratory conditions. We have identified 153 proteins from clinical P. vivax, majority of which do not show homology to any previously known gene products. We also report 29 new proteins that were found to be expressed in P. vivax for the first time. In addition, several proteins previously implicated as anti-malarial targets, were also found in our analysis. Most importantly, we found several unique proteins expressed by P. vivax. This study is an important step in providing insight into physiology of the parasite under clinical settings.

Structural Annotation of Mycobacterium tuberculosis Proteome

Of the similar to 4000 ORFs identified through the genome sequence of Mycobacterium tuberculosis (TB) H37Rv, experimentally determined structures are available for 312. Since knowledge of protein structures is essential to obtain a high-resolution understanding of the underlying biology, we seek to obtain a structural annotation for the genome, using computational methods. Structural models were obtained and validated for similar to 2877 ORFs, covering similar to 70% of the genome. Functional annotation of each protein was based on fold-based functional assignments and a novel binding site based ligand association. New algorithms for binding site detection and genome scale binding site comparison at the structural level, recently reported from the laboratory, were utilized. Besides these, the annotation covers detection of various sequence and sub-structural motifs and quaternary structure predictions based on the corresponding templates. The study provides an opportunity to obtain a global perspective of the fold distribution in the genome. The annotation indicates that cellular metabolism can be achieved with only 219 folds. New insights about the folds that predominate in the genome, as well as the fold-combinations that make up multi-domain proteins are also obtained. 1728 binding pockets have been associated with ligands through binding site identification and sub-structure similarity analyses. The resource (http://proline.physics.iisc.ernet.in/Tbstructuralannotation), being one of the first to be based on structure-derived functional annotations at a genome scale, is expected to be useful for better understanding of TB and for application in drug discovery. The reported annotation pipeline is fairly generic and can be applied to other genomes as well.

UDP-glucose 4, 6-dehydratase Activity Plays an Important Role in Maintaining Cell Wall Integrity and Virulence of Candida albicans

Candida albicans, a human fungal pathogen, undergoes morphogenetic changes that are associated with virulence. We report here that GAL102 in C. albicans encodes a homolog of dTDP-glucose 4,6-dehydratase, an enzyme that affects cell wall properties as well as virulence of many pathogenic bacteria. We found that GAL102 deletion leads to greater sensitivity to antifungal drugs and cell wall destabilizing agents like Calcofluor white and Congo red. The mutant also formed biofilms consisting mainly of hyphal cells that show less turgor. The NMR analysis of cell wall mannans of gal102 deletion strain revealed that a major constituent of mannan is missing and the phosphomannan component known to affect virulence is greatly reduced. We also observed that there was a substantial reduction in the expression of genes involved in biofilm formation but increase in the expression of genes encoding glycosylphosphatidylinositol-anchored proteins in the mutant. These, along with altered mannosylation of cell wall proteins together might be responsible for multiple phenotypes displayed by the mutant. Finally, the mutant was unable to grow in the presence of resident peritoneal macrophages and elicited a weak pro-inflammatory cytokine response in vitro. Similarly, this mutant elicited a poor serum pro-inflammatory cytokine response as judged by IFN gamma and TNF alpha levels and showed reduced virulence in a mouse model of systemic candidiasis. Importantly, an Ala substitution for a conserved Lys residue in the active site motif YXXXK, that abrogates the enzyme activity also showed reduced virulence and increased filamentation similar to the gal102 deletion strain. Since inactivating the enzyme encoded by GAL102 makes the cells sensitive to antifungal drugs and reduces its virulence, it can serve as a potential drug target in combination therapies for C. albicans and related pathogens.

Chimeras of Escherichia coli and Mycobacterium tuberculosis Single-Stranded DNA Binding Proteins: Characterization and Function in Escherichia coli

Single stranded DNA binding proteins (SSBs) are vital for the survival of organisms. Studies on SSBs from the prototype, Escherichia coli (EcoSSB) and, an important human pathogen, Mycobacterium tuberculosis (MtuSSB) had shown that despite significant variations in their quaternary structures, the DNA binding and oligomerization properties of the two are similar. Here, we used the X-ray crystal structure data of the two SSBs to design a series of chimeric proteins (m beta 1, m beta 1'beta 2, m beta 1-beta 5, m beta 1-beta 6 and m beta 4-beta 5) by transplanting beta 1, beta 1'beta 2, beta 1-beta 5, beta 1-beta 6 and beta 4-beta 5 regions, respectively of the N-terminal (DNA binding) domain of MtuSSB for the corresponding sequences in EcoSSB. In addition, m beta 1'beta 2(ESWR) SSB was generated by mutating the MtuSSB specific `PRIY' sequence in the beta 2 strand of m beta 1'beta 2 SSB to EcoSSB specific `ESWR' sequence. Biochemical characterization revealed that except for m beta 1 SSB, all chimeras and a control construct lacking the C-terminal domain (Delta C SSB) bound DNA in modes corresponding to limited and unlimited modes of binding. However, the DNA on MtuSSB may follow a different path than the EcoSSB. Structural probing by protease digestion revealed that unlike other SSBs used, m beta 1 SSB was also hypersensitive to chymotrypsin treatment. Further, to check for their biological activities, we developed a sensitive assay, and observed that m beta 1-beta 6, MtuSSB, m beta 1'beta 2 and m beta 1-beta 5 SSBs complemented E. coli Delta ssb in a dose dependent manner. Complementation by the m beta 1-beta 5 SSB was poor. In contrast, m beta 1'beta 2(ESWR) SSB complemented E. coli as well as EcoSSB. The inefficiently functioning SSBs resulted in an elongated cell/filamentation phenotype of E. coli. Taken together, our observations suggest that specific interactions within the DNA binding domain of the homotetrameric SSBs are crucial for their biological function.

Retention Time Variability as a Mechanism for Animal Mediated Long-Distance Dispersal

Long-distance dispersal (LDD) events, although rare for most plant species, can strongly influence population and community dynamics. Animals function as a key biotic vector of seeds and thus, a mechanistic and quantitative understanding of how individual animal behaviors scale to dispersal patterns at different spatial scales is a question of critical importance from both basic and applied perspectives. Using a diffusion-theory based analytical approach for a wide range of animal movement and seed transportation patterns, we show that the scale (a measure of local dispersal) of the seed dispersal kernel increases with the organisms' rate of movement and mean seed retention time. We reveal that variations in seed retention time is a key determinant of various measures of LDD such as kurtosis (or shape) of the kernel, thinkness of tails and the absolute number of seeds falling beyond a threshold distance. Using empirical data sets of frugivores, we illustrate the importance of variability in retention times for predicting the key disperser species that influence LDD. Our study makes testable predictions linking animal movement behaviors and gut retention times to dispersal patterns and, more generally, highlights the potential importance of animal behavioral variability for the LDD of seeds.