Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control chemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation.

Chlamydial infertility in women: Diagnosis, epidemiology and immune response

This project primarily focused on the development of a novel, non-invasive peptide-based ELISA to diagnose C. trachomatis-related infertility in women. In this study, an in silico approach was applied to identify and design novel peptide antigens that are specific to women with C. trachomatis-related infertility. The serological assay was developed such that it could potentially replace invasive techniques for early infertility investigation. The thesis further investigated the sero-prevalence of chlamydial infertility in women from a low-resource, high prevalence setting. In addition, the thesis also identified immune markers that were associated with the C. trachomatis-related infertility in women.

Examination of EXOSC4 as a new prognostic marker and a novel therapeutic avenue in lung adenocarcinoma

Although prevention and early detection of the disease greatly improved over the past few years, lung cancer remains the leading cause of cancer deaths. In order to be able to treat a larger population, we are in urgent need for novel treatments. While it is known that DNA repair genes play a major role in the oncogenic transformation, they also represent a weakness of cancers that constitute a therapeutic opportunity. To identify novel DNA repair genes implicated in Lung cancers, we conducted an in silico investigation to identify genes co-regulated with two DNA repair factors, BRCA2 and hSSB1. This approach allowed for the identification of EXOSC4, a component of the RNA Exosome machinery, as a potential factor involved in the maintenance of genome stability and that is deregulated in lung cancer.

Defective in Cuticular Ridges (DCR) of Arabidopsis thaliana, a Gene Associated with Surface Cutin Formation, Encodes a Soluble Diacylglycerol Acyltransferase

A key step in the triacylglycerol (TAG) biosynthetic pathway is the final acylation of diacylglycerol (DAG) by DAG acyltransferase. In silico analysis has revealed that the DCR (defective in cuticular ridges) (At5g23940) gene has a typical HX4D acyltransferase motif at the N-terminal end and a lipid binding motif VX(2)GF at the middle of the sequence. To understand the biochemical function, the gene was overexpressed in Escherichia coli, and the purified recombinant protein was found to acylate DAG specifically in an acyl-CoA-dependent manner. Overexpression of At5g23940 in a Saccharomyces cerevisiae quadruple mutant deficient in DAG acyltransferases resulted in TAG accumulation. At5g23940 rescued the growth of this quadruple mutant in the oleate-containing medium, whereas empty vector control did not. Lipid particles were localized in the cytosol of At5g23940-transformed quadruple mutant cells, as observed by oil red O staining. There was an incorporation of 16-hydroxyhexadecanoic acid into TAG in At5g23940-transformed cells of quadruple mutant. Here we report a soluble acyl-CoA-dependent DAG acyltransferase from Arabidopsis thaliana. Taken together, these data suggest that a broad specific DAG acyltransferase may be involved in the cutin as well as in the TAG biosynthesis by supplying hydroxy fatty acid.

Analysis of Conformational Variation in Macromolecular Structural Models

Experimental conditions or the presence of interacting components can lead to variations in the structural models of macromolecules. However, the role of these factors in conformational selection is often omitted by in silico methods to extract dynamic information from protein structural models. Structures of small peptides, considered building blocks for larger macromolecular structural models, can substantially differ in the context of a larger protein. This limitation is more evident in the case of modeling large multi-subunit macromolecular complexes using structures of the individual protein components. Here we report an analysis of variations in structural models of proteins with high sequence similarity. These models were analyzed for sequence features of the protein, the role of scaffolding segments including interacting proteins or affinity tags and the chemical components in the experimental conditions. Conformational features in these structural models could be rationalized by conformational selection events, perhaps induced by experimental conditions. This analysis was performed on a non-redundant dataset of protein structures from different SCOP classes. The sequence-conformation correlations that we note here suggest additional features that could be incorporated by in silico methods to extract dynamic information from protein structural models.

Computational interrogation of cis-regulatory elements of genes that are common targets of luteotropin and luteolysin in the primate corpus luteum

The rapid recent increase in microarray-based gene expression studies in the corpus luteum (CL) utilizing macaque models gathered increasing volume of data in publically accessible microarray expression databases. Examining gene pathways in different functional states of CL may help to understand the factors that control luteal function and hence human fertility. Co-regulation of genes in microarray experiments may imply common transcriptional regulation by sequence-specific DNA-binding transcriptional factors. We have computationally analyzed the transcription factor binding sites (TFBS) in a previously reported macaque luteal microarray gene set (n = 15) that are common targets of luteotropin (luteinizing hormone (LH) and human chorionic gonadotropin (hCG)) and luteolysin (prostaglandin (PG) F-2 alpha). This in silico approach can reveal transcriptional networks that control these important genes which are representative of the interplay between luteotropic and luteolytic factors in the control of luteal function. Our computational analyses revealed 6 matrix families whose binding sites are significantly over-represented in promoters of these genes. The roles of these factors are discussed, which might help to understand the transcriptional regulatory network in the control of luteal function. These factors might be promising experimental targets for investigation of human luteal insufficiency. (C) 2012 Elsevier B.V. All rights reserved.

Structural Integrity of the Greek Key Motif in beta gamma-Crystallins Is Vital for Central Eye Lens Transparency

Background: We highlight an unrecognized physiological role for the Greek key motif, an evolutionarily conserved super-secondary structural topology of the beta gamma-crystallins. These proteins constitute the bulk of the human eye lens, packed at very high concentrations in a compact, globular, short-range order, generating transparency. Congenital cataract (affecting 400,000 newborns yearly worldwide), associated with 54 mutations in beta gamma-crystallins, occurs in two major phenotypes nuclear cataract, which blocks the central visual axis, hampering the development of the growing eye and demanding earliest intervention, and the milder peripheral progressive cataract where surgery can wait. In order to understand this phenotypic dichotomy at the molecular level, we have studied the structural and aggregation features of representative mutations. Methods: Wild type and several representative mutant proteins were cloned, expressed and purified and their secondary and tertiary structural details, as well as structural stability, were compared in solution, using spectroscopy. Their tendencies to aggregate in vitro and in cellulo were also compared. In addition, we analyzed their structural differences by molecular modeling in silico. Results: Based on their properties, mutants are seen to fall into two classes. Mutants A36P, L45PL54P, R140X, and G165fs display lowered solubility and structural stability, expose several buried residues to the surface, aggregate in vitro and in cellulo, and disturb/distort the Greek key motif. And they are associated with nuclear cataract. In contrast, mutants P24T and R77S, associated with peripheral cataract, behave quite similar to the wild type molecule, and do not affect the Greek key topology. Conclusion: When a mutation distorts even one of the four Greek key motifs, the protein readily self-aggregates and precipitates, consistent with the phenotype of nuclear cataract, while mutations not affecting the motif display `native state aggregation', leading to peripheral cataract, thus offering a protein structural rationale for the cataract phenotypic dichotomy ``distort motif, lose central vision''.

Morphological Effects of G-Quadruplex Stabilization Using a Small Molecule in Zebrafish

Zebrafish (Danio rerio) embryos are transparent and advantageous for studying early developmental changes due to ex utero development, making them an appropriate model for studying gene expression changes as a result of molecular targeting. Zebrafish embryos were injected with a previously reported G-quadruplex selective ligand, and the phenotypic changes were recorded. We report marked discrepancies in the development of intersegmental vessels. In silico analysis determined that the putative G-quadruplex motif occur in the upstream promoter region of the Cdh5 (N-cadherin) gene. A real-time polymerase chain reaction-based investigation indicated that in zebrafish, CDH-2 (ZN-cad) was significantly downregulated in the ligand-treated embryos. Biophysical characterization of the interaction of the ligand with the G-quadruplex motif found in this promoter yielded strong binding and stabilization of the G-quadruplex with this ligand. Hence, we report for the first time the phenotypic impact of G-quadruplex targeting with a ligand in a vertebrate organism. This study has unveiled not only G-quadruplex targeting in non-human animal species but also the potential that G-quadruplexes can provide a ready tool for understanding the phenotypic effects of targeting certain important genes involved in differentiation and developmental processes in a living eukaryotic organism.

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (phi, ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (phi, ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (phi, ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.

Cyclic nucleotide binding and structural changes in the isolated GAF domain of Anabaena adenylyl cyclase, CyaB2

GAF domains are a large family of regulatory domains, and a subset are found associated with enzymes involved in cyclic nucleotide (cNMP) metabolism such as adenylyl cyclases and phosphodiesterases. CyaB2, an adenylyl cyclase from Anabaena, contains two GAF domains in tandem at the N-terminus and an adenylyl cyclase domain at the C-terminus. Cyclic AMP, but not cGMP, binding to the GAF domains of CyaB2 increases the activity of the cyclase domain leading to enhanced synthesis of cAMP. Here we show that the isolated GAFb domain of CyaB2 can bind both cAMP and cGMP, and enhanced specificity for cAMP is observed only when both the GAFa and the GAFb domains are present in tandem(GAFab domain). In silico docking and mutational analysis identified distinct residues important for interaction with either cAMP or cGMP in the GAFb domain. Structural changes associated with ligand binding to the GAF domains could not be detected by bioluminescence resonance energy transfer (BRET) experiments. However, amide hydrogen-deuterium exchange mass spectrometry (HDXMS) experiments provided insights into the structural basis for cAMP-induced allosteric regulation of the GAF domains, and differences in the changes induced by cAMP and cGMP binding to the GAF domain. Thus, our findings could allow the development of molecules that modulate the allosteric regulation by GAF domains present in pharmacologically relevant proteins.

SNP Discovery in European Anchovy (Engraulis encrasicolus, L) by High-Throughput Transcriptome and Genome Sequencing

12 p.

Altered expression of Alzheimer's disease-related genes in the cerebellum of autistic patients: a model for disrupted brain connectome and therapy

Autism and Alzheimer's disease (AD) are, respectively, neurodevelopmental and degenerative diseases with an increasing epidemiological burden. The AD-associated amyloid-beta precursor protein-alpha has been shown to be elevated in severe autism, leading to the 'anabolic hypothesis' of its etiology. Here we performed a focused microarray analysis of genes belonging to NOTCH and WNT signaling cascades, as well as genes related to AD and apoptosis pathways in cerebellar samples from autistic individuals, to provide further evidence for pathological relevance of these cascades for autism. By using the limma package from R and false discovery rate, we demonstrated that 31% (116 out of 374) of the genes belonging to these pathways displayed significant changes in expression (corrected P-values <0.05), with mitochondria- related genes being the most downregulated. We also found upregulation of GRIN1, the channel-forming subunit of NMDA glutamate receptors, and MAP3K1, known activator of the JNK and ERK pathways with anti-apoptotic effect. Expression of PSEN2 (presinilin 2) and APBB1 (or F65) were significantly lower when compared with control samples. Based on these results, we propose a model of NMDA glutamate receptor-mediated ERK activation of alpha-secretase activity and mitochondrial adaptation to apoptosis that may explain the early brain overgrowth and disruption of synaptic plasticity and connectome in autism. Finally, systems pharmacology analyses of the model that integrates all these genes together (NOWADA) highlighted magnesium (Mg2+) and rapamycin as most efficient drugs to target this network model in silico. Their potential therapeutic application, in the context of autism, is therefore discussed.

Estudo in vitro de derivados sintéticos da Isoniazida e da Pirazinamida sobre Leishmania (Viannia) braziliensis

As leishmanioses são um grupo de doenças causadas por protozoários do gênero Leishmania spp que afetam 98 países. No Brasil, no ano de 2013, foram relatados 3.253 casos de leishmaniose visceral e 18.226 casos de Leishmaniose Tegumentar Americana. O tratamento de primeira escolha continua sendo realizado com antimoniais pentavalentes, e em casos de insucessos os fármacos de segunda escolha são a pentamidina e a anfotericina B. Tais medicamentos causam intensos efeitos adversos e ultimamente têm surgido cepas resistentes aos mesmos. Em áreas endêmicas têm sido cada vez mais comum o surgimento da co-infecção Leishmania com Mycobacterium tuberculosis. O tratamento para a tuberculose com pirazinamida (PZA) e isoniazida (INZ), controla a leishmaniose. Esses dados sugerem atividade anti-leishmania da PZA e da INZ. O objetivo deste trabalho foi avaliar a atividade in vitro da INZ e da PZA e seus compostos derivados (série G e série R, respectivamente) sobre Leishmania (Viannia) braziliensis. As moléculas foram testadas em monocamadas de macrófagos peritoneais de camunongos infectados com L. (V) braziliensis durante 48h. Todas as moléculas testadas inibiram o índice de infecção de forma dose dependente em comparação aos controles. As moléculas da série R foram mais ativas do que a PZA, porém o resultado foi significativo somente para a R02 (p < 0,005). Apenas a molécula R05 (76,64M) foi relativamente tóxica para macrófagos. Os compostos mais ativos foram R02, G01 e G02, cujos índices de seletividade foram 14,31, 19 e 30, respectivamente. A dosagem de nitrito foi feita em sobrenadantes de monocamadas de macrófagos peritoniais infectados e tratados com as substâncias nas concentrações 10 e 100M. A G01 e a G02 estimularam a produção de NO2 nas duas concentrações, entretanto o resultado foi estatisticamente significativo para a G02 em 100M (p < 0,0001), a G05 só estimulou óxido nítrico na maior concentração. Todos os compostos da série R estimularam NO2, contudo, o resultado foi estatisticamente significativo para a R03 e R05 a 100M (p < 0,001). Adicionalmente, foi realizado uma análise preditiva in sílico de parâmetros farmacocinéticos das moléculas mais ativas in vitro, utilizando o software admetSAR. Os dados obtidos mostraram que de forma semelhante às suas moléculas originais a G01, G02 e R02 apresentaram alta capacidade de serem absorvidas pelo trato gastrointestinal, baixo potencial hepatotoxico e carcinogênico. Juntos, esses dados demonstram que essas moléculas são seletivamente tóxicas para o parasito com potencial para serem testadas pela via oral em estudos em modelo experimental de infecção.

The role of MC1R gene in buffalo coat color

Melanocortin-1 receptor (MC1R) plays a major role in pigmentation in many species. To investigate if the MC1R gene is associated with coat color in water buffalo, the coding region of MC1R gene of 216 buffalo samples was sequenced, which included 49 black river buffalo (Murrah and Nili-Ravi), 136 swamp buffalo (Dehong, Diandongnan, Dechang, Guizhou, and Xilin) with white and gray body, and 31 hybrid offspring of river buffalo Nili-Ravi (or Murrah) and swamp buffalo. Among the three variation sites found, SNP684 was synonymous, while SNP310 and SNP384 were nonsynonymous, leading to p.S104G and p.I128M changes, respectively. Only Individuals carrying homozygote E-BR/E-BR were black. The genotype and phenotype analysis of the hybrid offspring of black river buffalo and gray swamp buffalo further revealed that the river buffalo type allele E-BR or the allele carrying the amino acid p.104S was important for the full function of MC1R. The in silico functional analysis showed that the amino acid substitutions p.G104S and p.M128I had significant impact on the function of MC1R. Above results indicate that the allele E-BR or the allele carrying the amino acid p.104S was associated with the black coat color in buffalo.

Prevalent positive epistasis in Escherichia coli and Saccharomyces cerevisiae metabolic networks

Epistasis refers to the interaction between genes. Although high-throughput epistasis data from model organisms are being generated and used to construct genetic networks(1-3), the extent to which genetic epistasis reflects biologically meaningful interactions remains unclear(4-6). We have addressed this question through in silico mapping of positive and negative epistatic interactions amongst biochemical reactions within the metabolic networks of Escherichia coli and Saccharomyces cerevisiae using flux balance analysis. We found that negative epistasis occurs mainly between nonessential reactions with overlapping functions, whereas positive epistasis usually involves essential reactions, is highly abundant and, unexpectedly, often occurs between reactions without overlapping functions. We offer mechanistic explanations of these findings and experimentally validate them for 61 S. cerevisiae gene pairs.