A systems perspective of host-pathogen interactions: predicting disease outcome in tuberculosis

The complex web of interactions between the host immune system and the pathogen determines the outcome of any infection. A computational model of this interaction network, which encodes complex interplay among host and bacterial components, forms a useful basis for improving the understanding of pathogenesis, in filling knowledge gaps and consequently to identify strategies to counter the disease. We have built an extensive model of the Mycobacterium tuberculosis host-pathogen interactome, consisting of 75 nodes corresponding to host and pathogen molecules, cells, cellular states or processes. Vaccination effects, clearance efficiencies due to drugs and growth rates have also been encoded in the model. The system is modelled as a Boolean network. Virtual deletion experiments, multiple parameter scans and analysis of the system's response to perturbations, indicate that disabling processes such as phagocytosis and phagolysosome fusion or cytokines such as TNF-alpha and IFN-gamma, greatly impaired bacterial clearance, while removing cytokines such as IL-10 alongside bacterial defence proteins such as SapM greatly favour clearance. Simulations indicate a high propensity of the pathogen to persist under different conditions.

A systems perspective of host-pathogen interactions: predicting disease outcome in tuberculosis

The complex web of interactions between the host immune system and the pathogen determines the outcome of any infection. A computational model of this interaction network, which encodes complex interplay among host and bacterial components, forms a useful basis for improving the understanding of pathogenesis, in filling knowledge gaps and consequently to identify strategies to counter the disease. We have built an extensive model of the Mycobacterium tuberculosis host-pathogen interactome, consisting of 75 nodes corresponding to host and pathogen molecules, cells, cellular states or processes. Vaccination effects, clearance efficiencies due to drugs and growth rates have also been encoded in the model. The system is modelled as a Boolean network. Virtual deletion experiments, multiple parameter scans and analysis of the system's response to perturbations, indicate that disabling processes such as phagocytosis and phagolysosome fusion or cytokines such as TNF-alpha and IFN-gamma, greatly impaired bacterial clearance, while removing cytokines such as IL-10 alongside bacterial defence proteins such as SapM greatly favour clearance. Simulations indicate a high propensity of the pathogen to persist under different conditions.

Analysis of a Genomic DNA Expression Library of Mycobacterium tuberculosis Using Tuberculosis Patient Sera: Evidence for Modulation of Host Immune Response

DNA obtained from a human sputum isolate of Mycobacterium tuberculosis, NTI-64719, which showed extensive dissemination in the guinea pig model resulting in a high score for virulence was used to construct an expression library in the lambda ZAP vector. The size of DNA inserts in the library ranged from 1 to 3 kb, and recombinants represented 60% of the total plaques obtained. When probed with pooled serum from chronically infected tuberculosis patients, the library yielded 176 recombinants with a range of signal intensities. Among these, 93 recombinants were classified into 12 groups on the basis of DNA hybridization experiments, The polypeptides synthesized by the recombinants were predominantly LacZ fusion proteins, Serum obtained from patients who were clinically diagnosed to be in the early phase of M. tuberculosis infection was used to probe the 176 recombinants obtained. interestingly, some recombinants that gave very strong signals in the original screen did not react with early-phase serum; conversely, others whose signals were extremely weak in the original screen gave very intense signals with serum from recently infected patients, This indicates the differential nature of either the expression of these antigens or the immune response elicited by them as a function of disease progression.

Genetic characterization of Indian type O FMD virus 3A region in context with host cell preference

The 3A region of foot-and-mouth disease virus has been implicated in host range and virulence. For example, amino acid deletions in the porcinophilic strain (O/TAW/97) at 93-102 aa of the 153 codons long 3A protein have been recognized as the determinant of species specificity. In the present study, 18 type 0 FMDV isolates from India were adapted in different cell culture systems and the 3A sequence was analyzed. These isolates had complete 3A coding sequence (153 aa) and did not exhibit growth restriction in cells based on species of origin. The 3A region was found to be highly conserved at N-terminal half (1-75 aa) but exhibited variability or substitutions towards C-terminal region (80-153). Moreover the amino acid substitutions were more frequent in recent Indian buffalo isolates but none of the Indian isolates showed deletion in 3A protein, which may be the reason for the absence of host specificity in vitro. Further inclusive analysis of 3A region will reveal interesting facts about the variability of FMD virus 3A region in an endemic environment. (C) 2010 Elsevier B.V. All rights reserved.

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Prediction of protein-protein interactions between human host and a pathogen and its application to three pathogenic bacteria

Molecular understanding of disease processes can be accelerated if all interactions between the host and pathogen are known. The unavailability of experimental methods for large-scale detection of interactions across host and pathogen organisms hinders this process. Here we apply a simple method to predict protein-protein interactions across a host and pathogen organisms. We use homology detection approaches against the protein-protein interaction databases. DIP and iPfam in order to predict interacting proteins in a host-pathogen pair. In the present work, we first applied this approach to the test cases involving the pairs phage T4 - Escherichia coli and phage lambda - E. coli and show that previously known interactions could be recognized using our approach. We further apply this approach to predict interactions between human and three pathogens E. coli, Salmonella enterica typhimurium and Yersinia pestis. We identified several novel interactions involving proteins of host or pathogen that could be thought of as highly relevant to the disease process. Serendipitously, many interactions involve hypothetical proteins of yet unknown function. Hypothetical proteins are predicted from computational analysis of genome sequences with no laboratory analysis on their functions yet available. The predicted interactions involving such proteins could provide hints to their functions. (C) 2011 Elsevier B.V. All rights reserved.

Progress in tuberculosis vaccine development and host-directed therapies-a state of the art review

Tuberculosis continues to kill 1.4 million people annually. During the past 5 years, an alarming increase in the number of patients with multidrug-resistant tuberculosis and extensively drug-resistant tuberculosis has been noted, particularly in eastern Europe, Asia, and southern Africa. Treatment outcomes with available treatment regimens for drug-resistant tuberculosis are poor. Although substantial progress in drug development for tuberculosis has been made, scientific progress towards development of interventions for prevention and improvement of drug treatment outcomes have lagged behind. Innovative interventions are therefore needed to combat the growing pandemic of multidrug-resistant and extensively drug-resistant tuberculosis. Novel adjunct treatments are needed to accomplish improved cure rates for multidrug-resistant and extensively drug-resistant tuberculosis. A novel, safe, widely applicable, and more effective vaccine against tuberculosis is also desperately sought to achieve disease control. The quest to develop a universally protective vaccine for tuberculosis continues. So far, research and development of tuberculosis vaccines has resulted in almost 20 candidates at different stages of the clinical trial pipeline. Host-directed therapies are now being developed to refocus the anti-Mycobacterium tuberculosis-directed immune responses towards the host; a strategy that could be especially beneficial for patients with multidrug-resistant tuberculosis or extensively drug-resistant tuberculosis. As we are running short of canonical tuberculosis drugs, more attention should be given to host-directed preventive and therapeutic intervention measures.

Differences in host species relationships and biogeographic influences produce contrasting patterns of prevalence, community composition and genetic structure in two genera of avian malaria parasites in southern Melanesia

1. Host-parasite interactions have the potential to influence broadscale ecological and evolutionary processes, levels of endemism, divergence patterns and distributions in host populations. Understanding the mechanisms involved requires identification of the factors that shape parasite distribution and prevalence. 2. A lack of comparative information on community-level host-parasite associations limits our understanding of the role of parasites in host population divergence processes. Avian malaria (haemosporidian) parasites in bird communities offer a tractable model system to examine the potential for pathogens to influence evolutionary processes in natural host populations. 3. Using cytochrome b variation, we characterized phylogenetic diversity and prevalence of two genera of avian haemosporidian parasites, Plasmodium and Haemoproteus, and analysed biogeographic patterns of lineages across islands and avian hosts, in southern Melanesian bird communities to identify factors that explain patterns of infection. 4. Plasmodium spp. displayed isolation-by-distance effects, a significant amount of genetic variation distributed among islands but insignificant amounts among host species and families, and strong local island effects with respect to prevalence. Haemoproteus spp. did not display isolation-by-distance patterns, showed marked structuring of genetic variation among avian host species and families, and significant host species prevalence patterns. 5. These differences suggest that Plasmodium spp. infection patterns were shaped by geography and the abiotic environment, whereas Haemoproteus spp. infection patterns were shaped predominantly by host associations. Heterogeneity in the complement and prevalence of parasite lineages infecting local bird communities likely exposes host species to a mosaic of spatially divergent disease selection pressures across their naturally fragmented distributions in southern Melanesia. Host associations for Haemoproteus spp. indicate a capacity for the formation of locally co-adapted host-parasite relationships, a feature that may limit intraspecific gene flow or range expansions of closely related host species.

Heme biosynthesis by the malarial parasite - Import of delta-aminolevulinate dehydrase from the host red cell

The mouse and human malarial parasites, Plasmodium berghei and Plasmodium falciparum, respectively, synthesize heme de novo following the standard pathway observed in animals despite the availability of large amounts of heme, derived from red cell hemoglobin, which is stored as hemozoin pigment, The enzymes, delta-aminolevulinate dehydrase (ALAD), coproporphyrinogen oxidase, and ferrochelatase are present at strikingly high levels in the P, berghei infected mouse red cell in vivo, The isolated parasite has low levels of ALAD and the data clearly indicate it to be of red cell origin. The purified enzyme preparations from the uninfected red cell and the parasite are identical in kinetic properties, subunit molecular weight, cross-reaction with antibodies to the human enzyme, and N-terminal amino acid sequence. Immunogold electron microscopy of the infected culture indicates that the enzyme is present inside the parasite and, therefore, is not a contaminant, The parasite derives functional ALAD from the host and the enzyme binds specifically to isolated parasite membrane in vitro, suggestive of the involvement of a receptor in its translocation into the parasite, While, ALAD, coproporphyrinogen oxidase, and ferrochelatase from the parasite and the uninfected red cell supernatant have identical subunit molecular weights on SDS-polyacrylamide gel electrophoresis and show immunological cross-reaction with antibodies to the human enzymes, as revealed by Western analysis, the first enzyme of the pathway, namely, delta-aminolevulinate synthase (ALAS) in the parasite, unlike that of the red cell host, does not cross-react with antibodies to the human enzyme, However, ALAS enzyme activity in the parasite is higher than that of the infected red cell supernatant. We therefore conclude that the parasite, while making its own ALAS, imports ALAD and perhaps most of the other enzymes of the pathway from the host to synthesize heme de novo, and this would enable it to segregate this heme from the heme derived from red cell hemoglobin degradation, ALAS of the parasite and the receptor(s) involved in the translocation of the host enzymes into the parasite would be unique drug targets.

How does contrasting dependence of impurity-atom diffusivity on the density of host disordered medium arise?

We report results of molecular dynamics investigations into neutral impurity diffusing within an amorphous solid as a function of the size of the diffusant and density of the host amorphous matrix. We find that self diffusivity exhibits an anomalous maximum as a function of the size of the impurity species. An analysis of properties of the impurity atom with maximum diffusivity shows that it is associated with lower mean square force, reduced backscattering of velocity autocorrelation function, near-exponential decay of the intermediate scattering function (as compared to stretched-exponential decay for other sizes of the impurity species) and lower activation energy. These results demonstrate the existence of size-dependent diffusivity maximum in disordered solids. Further, we show that the diffusivity maximum is observed at lower impurity diameters with increase in density. This is explained in terms of the Levitation parameter and the void structure of the amorphous solid. We demonstrate that these results imply contrasting dependence of self diffusivity (D) on the density of the amorphous matrix, p. D increases with p for small sizes of the impurity but shows an increase followed by a decrease for intermediate sizes of the impurity atom. For large sizes of the impurity atom, D decreases with increase in p. These contrasting dependence arises naturally from the existence of Levitation Effect.

Genetic analysis of foot-and-mouth disease virus serotype A of Indian origin and detection of positive selection and recombination in leader protease- and capsid-coding regions

The leader protease (L-pro) and capsid-coding sequences (P1) constitute approximately 3 kb of the foot-and-mouth disease virus (FMDV). We studied the phylogenetic relationship of 46 FMDV serotype A isolates of Indian origin collected during the period 1968-2005 and also eight vaccine strains using the neighbour-joining tree and Bayesian tree methods. The viruses were categorized under three major groups - Asian, Euro-South American and European. The Indian isolates formed a distinct genetic group among the Asian isolates. The Indian isolates were further classified into different genetic subgroups (<5% divergence). Post-1995 isolates were divided into two subgroups while a few isolates which originated in the year 2005 from Andhra Pradesh formed a separate group. These isolates were closely related to the isolates of the 1970s. The FMDV isolates seem to undergo reverse mutation or onvergent evolution wherein sequences identical to the ancestors are present in the isolates in circulation. The eight vaccine strains included in the study were not related to each other and belonged to different genetic groups. Recombination was detected in the L-pro region in one isolate (A IND 20/82) and in the VP1 coding 1D region in another isolate (A RAJ 21/96). Positive selection was identified at aa positions 23 in the L-pro (P<0.05; 0.046*) and at aa 171 in the capsid protein VP1 (P<0.01; 0.003**).