Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

PIM2 Induced COX-2 and MMP-9 Expression in Macrophages Requires PI3K and Notch1 Signaling

Activation of inflammatory immune responses during granuloma formation by the host upon infection of mycobacteria is one of the crucial steps that is often associated with tissue remodeling and breakdown of the extracellular matrix. In these complex processes, cyclooxygenase-2 (COX-2) plays a major role in chronic inflammation and matrix metalloproteinase-9 (MMP-9) significantly in tissue remodeling. In this study, we investigated the molecular mechanisms underlying Phosphatidyl-myo-inositol dimannosides (PIM2), an integral component of the mycobacterial envelope, triggered COX-2 and MMP-9 expression in macrophages. PIM2 triggers the activation of Phosphoinositide-3 Kinase (PI3K) and Notch1 signaling leading to COX-2 and MMP-9 expression in a Toll-like receptor 2 (TLR2)-MyD88 dependent manner. Notch1 signaling perturbations data demonstrate the involvement of the cross-talk with members of PI3K and Mitogen activated protein kinase pathway. Enforced expression of the cleaved Notch1 in macrophages induces the expression of COX-2 and MMP-9. PIM2 triggered significant p65 nuclear factor-kappa B (NF-kappa B) nuclear translocation that was dependent on activation of PI3K or Notch1 signaling. Furthermore, COX-2 and MMP-9 expression requires Notch1 mediated recruitment of uppressor of Hairless (CSL) and NF-kappa B to respective promoters. Inhibition of PIM2 induced COX-2 resulted in marked reduction in MMP-9 expression clearly implicating the role of COX-2 dependent signaling events in driving the MMP-9 expression. Taken together, these data implicate PI3K and Notch1 signaling as obligatory early proximal signaling events during PIM2 induced COX-2 and MMP-9 expression in macrophages.

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.

Phenotypic and Functional Characterization of Human Mammary Stem/Progenitor Cells in Long Term Culture

Background: Cancer stem cells exhibit close resemblance to normal stem cells in phenotype as well as function. Hence, studying normal stem cell behavior is important in understanding cancer pathogenesis. It has recently been shown that human breast stem cells can be enriched in suspension cultures as mammospheres. However, little is known about the behavior of these cells in long-term cultures. Since extensive self-renewal potential is the hallmark of stem cells, we undertook a detailed phenotypic and functional characterization of human mammospheres over long-term passages. Methodology: Single cell suspensions derived from human breast `organoids' were seeded in ultra low attachment plates in serum free media. Resulting primary mammospheres after a week (termed T1 mammospheres) were subjected to passaging every 7th day leading to the generation of T2, T3, and T4 mammospheres. Principal Findings: We show that primary mammospheres contain a distinct side-population (SP) that displays a CD24(low)/CD44(low) phenotype, but fails to generate mammospheres. Instead, the mammosphere-initiating potential rests within the CD44(high)/CD24(low) cells, in keeping with the phenotype of breast cancer-initiating cells. In serial sphere formation assays we find that even though primary (T1) mammospheres show telomerase activity and fourth passage T4 spheres contain label-retaining cells, they fail to initiate new mammospheres beyond T5. With increasing passages, mammospheres showed an increase in smaller sized spheres, reduction in proliferation potential and sphere forming efficiency, and increased differentiation towards the myoepithelial lineage. Significantly, staining for senescence-associated beta-galactosidase activity revealed a dramatic increase in the number of senescent cells with passage, which might in part explain the inability to continuously generate mammospheres in culture. Conclusions: Thus, the self-renewal potential of human breast stem cells is exhausted within five in vitro passages of mammospheres, suggesting the need for further improvisation in culture conditions for their long-term maintenance.

Location of Pathogenic Bacteria during Persistent Infections: Insights from an Analysis Using Game Theory

Bacterial persistent infections are responsible for a significant amount of the human morbidity and mortality. Unlike acute bacterial infections, it is very difficult to treat persistent bacterial infections (e.g. tuberculosis). Knowledge about the location of pathogenic bacteria during persistent infection will help to treat such conditions by designing novel drugs which can reach such locations. In this study, events of bacterial persistent infections were analyzed using game theory. A game was defined where the pathogen and the host are the two players with a conflict of interest. Criteria for the establishment of Nash equilibrium were calculated for this game. This theoretical model, which is very simple and heuristic, predicts that during persistent infections pathogenic bacteria stay in both intracellular and extracellular compartments of the host. The result of this study implies that a bacterium should be able to survive in both intracellular and extracellular compartments of the host in order to cause persistent infections. This explains why persistent infections are more often caused by intracellular pathogens like Mycobacterium and Salmonella. Moreover, this prediction is in consistence with the results of previous experimental studies.

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

lac Repressor Is an Antivirulence Factor of Salmonella enterica: Its Role in the Evolution of Virulence in Salmonella

The genus Salmonella includes many pathogens of great medical and veterinary importance. Bacteria belonging to this genus are very closely related to those belonging to the genus Escherichia. lacZYA operon and lacI are present in Escherichia coli, but not in Salmonella enterica. It has been proposed that Salmonella has lost lacZYA operon and lacI during evolution. In this study, we have investigated the physiological and evolutionary significance of the absence of lacI in Salmonella enterica. Using murine model of typhoid fever, we show that the expression of Lacl causes a remarkable reduction in the virulence of Salmonella enterica. Lacl also suppresses the ability of Salmonella enterica to proliferate inside murine macrophages. Microarray analysis revealed that Lacl interferes with the expression of virulence genes of Salmonella pathogenicity island 2. This effect was confirmed by RT-PCR and Western blot analysis. Interestingly, we found that SBG0326 of Salmonella bongori is homologous to lacI of Escherichia coli. Salmonella bongori is the only other species of the genus Salmonella and it lacks the virulence genes of Salmonella pathogenicity island 2. Overall, our results demonstrate that Lacl is an antivirulence factor of Salmonella enterica and suggest that absence of lacI has facilitated the acquisition of virulence genes of Salmonella pathogenicity island 2 in Salmonella enterica making it a successful systemic pathogen.

A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Metabolome based reaction graphs of M.tuberculosis and M.leprae: A comparative network analysis

Background. Several types of networks, such as transcriptional, metabolic or protein-protein interaction networks of various organisms have been constructed, that have provided a variety of insights into metabolism and regulation. Here, we seek to exploit the reaction-based networks of three organisms for comparative genomics. We use concepts from spectral graph theory to systematically determine how differences in basic metabolism of organisms are reflected at the systems level and in the overall topological structures of their metabolic networks. Methodology/Principal Findings. Metabolome-based reaction networks of Mycobacterium tuberculosis, Mycobacterium leprae and Escherichia coli have been constructed based on the KEGG LIGAND database, followed by graph spectral analysis of the network to identify hubs as well as the sub-clustering of reactions. The shortest and alternate paths in the reaction networks have also been examined. Sub-cluster profiling demonstrates that reactions of the mycolic acid pathway in mycobacteria form a tightly connected sub-cluster. Identification of hubs reveals reactions involving glutamate to be central to mycobacterial metabolism, and pyruvate to be at the centre of the E. coli metabolome. The analysis of shortest paths between reactions has revealed several paths that are shorter than well established pathways. Conclusions. We conclude that severe downsizing of the leprae genome has not significantly altered the global structure of its reaction network but has reduced the total number of alternate paths between its reactions while keeping the shortest paths between them intact. The hubs in the mycobacterial networks that are absent in the human metabolome can be explored as potential drug targets. This work demonstrates the usefulness of constructing metabolome based networks of organisms and the feasibility of their analyses through graph spectral methods. The insights obtained from such studies provide a broad overview of the similarities and differences between organisms, taking comparative genomics studies to a higher dimension.

Tool for Identification of Duplicate Records Downloaded from Multiple Cd-Roms. A Case Study with Spirs Based Databases

As research becomes more and more interdisciplinary, literature search from CD-ROM databases is often carried out on more than one CD-ROM database. This results in retrieving duplicate records due to same literature being covered (indexed) in more than one database. The retrieval software does not identify such duplicate records. Three different programs have been written to accomplish the task of identifying the duplicate records. These programs are executed from a shell script to minimize manual intervention. The various fields that have been used (extracted) to identify the duplicate records include the article title, year, volume number, issue number and pagination. The shell script when executed prompts for input file that may contain duplicate records. The programs identify the duplicate records and write them to a new file.

Sesbania Mosaic Virus (SeMV) Infectious Clone: Possible Mechanism of 3 ` and 5 ` End Repair and Role of Polyprotein Processing in Viral Replication

Sesbania mosaic virus (SeMV) is a positive stranded RNA virus belonging to the genus Sobemovirus. Construction of an infectious clone is an essential step for deciphering the virus gene functions in vivo. Using Agrobacterium based transient expression system we show that SeMV icDNA is infectious on Sesbania grandiflora and Cyamopsis tetragonoloba plants. The efficiency of icDNA infection was found to be significantly high on Cyamopsis plants when compared to that on Sesbania grandiflora. The coat protein could be detected within 6 days post infiltration in the infiltrated leaves. Different species of viral RNA (double stranded and single stranded genomic and subgenomic RNA) could be detected upon northern analysis, suggesting that complete replication had taken place. Based on the analysis of the sequences at the genomic termini of progeny RNA from SeMV icDNA infiltrated leaves and those of its 3' and 5' terminal deletion mutants, we propose a possible mechanism for 3' and 5' end repair in vivo. Mutation of the cleavage sites in the polyproteins encoded by ORF 2 resulted in complete loss of infection by the icDNA, suggesting the importance of correct polyprotein processing at all the four cleavage sites for viral replication. Complementation analysis suggested that ORF 2 gene products can act in trans. However, the trans acting ability of ORF 2 gene products was abolished upon deletion of the N-terminal hydrophobic domain of polyprotein 2a and 2ab, suggesting that these products necessarily function at the replication site, where they are anchored to membranes.

Extent of Structural Asymmetry in Homodimeric Proteins: Prevalence and Relevance

Most homodimeric proteins have symmetric structure. Although symmetry is known to confer structural and functional advantage, asymmetric organization is also observed. Using a non-redundant dataset of 223 high-resolution crystal structures of biologically relevant homodimers, we address questions on the prevalence and significance of asymmetry. We used two measures to quantify global and interface asymmetry, and assess the correlation of several molecular and structural parameters with asymmetry. We have identified rare cases (11/223) of biologically relevant homodimers with pronounced global asymmetry. Asymmetry serves as a means to bring about 2:1 binding between the homodimer and another molecule; it also enables cellular signalling arising from asymmetric macromolecular ligands such as DNA. Analysis of these cases reveals two possible mechanisms by which possible infinite array formation is prevented. In case of homodimers associating via non-topologically equivalent surfaces in their tertiary structures, ligand-dependent mechanisms are used. For stable dimers binding via large surfaces, ligand-dependent structural change regulates polymerisation/depolymerisation; for unstable dimers binding via smaller surfaces that are not evolutionarily well conserved, dimerisation occurs only in the presence of the ligand. In case of homodimers associating via interaction surfaces with parts of the surfaces topologically equivalent in the tertiary structures, steric hindrance serves as the preventive mechanism of infinite array. We also find that homodimers exhibiting grossly symmetric organization rarely exhibit either perfect local symmetry or high local asymmetry. Binding of small ligands at the interface does not cause any significant variation in interface asymmetry. However, identification of biologically relevant interface asymmetry in grossly symmetric homodimers is confounded by the presence of similar small magnitude changes caused due to artefacts of crystallisation. Our study provides new insights regarding accommodation of asymmetry in homodimers.

Structural and Mutational Studies on Substrate Specificity and Catalysis of Salmonella typhimurium D-Cysteine Desulfhydrase

Salmonella typhimurium DCyD (StDCyD) is a fold type II pyridoxal 5' phosphate (PLP)-dependent enzyme that catalyzes the degradation of D-Cys to H2S and pyruvate. It also efficiently degrades beta-chloro-D-alanine (beta CDA). D-Ser is a poor substrate while the enzyme is inactive with respect to L-Ser and 1-amino-1-carboxy cyclopropane (ACC). Here, we report the X-ray crystal structures of StDCyD and of crystals obtained in the presence of D-Cys, beta CDA, ACC, D-Ser, L-Ser, D-cycloserine (DCS) and L-cycloserine (LCS) at resolutions ranging from 1.7 to 2.6 angstrom. The polypeptide fold of StDCyD consisting of a small domain (residues 48-161) and a large domain (residues 1-47 and 162-328) resembles other fold type II PLP dependent enzymes. The structures obtained in the presence of D-Cys and beta CDA show the product, pyruvate, bound at a site 4.0-6.0 angstrom away from the active site. ACC forms an external aldimine complex while D- and L-Ser bind non-covalently suggesting that the reaction with these ligands is arrested at C alpha proton abstraction and transimination steps, respectively. In the active site of StDCyD cocrystallized with DCS or LCS, electron density for a pyridoxamine phosphate (PMP) was observed. Crystals soaked in cocktail containing these ligands show density for PLP-cycloserine. Spectroscopic observations also suggest formation of PMP by the hydrolysis of cycloserines. Mutational studies suggest that Ser78 and Gln77 are key determinants of enzyme specificity and the phenolate of Tyr287 is responsible for C alpha proton abstraction from D-Cys. Based on these studies, a probable mechanism for the degradation of D-Cys by StDCyD is proposed.