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.

The Mycobacterial MsDps2 Protein Is a Nucleoid-Forming DNA Binding Protein Regulated by Sigma Factors sigma(A) and sigma(B)

The Dps (DNA-binding protein from starved cells) proteins from Mycobacterium smegmatis MsDps1 and MsDps2 are both DNA-binding proteins with some differences. While MsDps1 has two oligomeric states, with one of them responsible for DNA binding, MsDps2 has only one DNA-binding oligomeric state. Both the proteins however, show iron-binding activity. The MsDps1 protein has been shown previously to be induced under conditions of starvation and osmotic stress and is regulated by the extra cellular sigma factors sigma(H) and sigma(F). We show here, that the second Dps homologue in M. smegmatis, namely MsDps2, is purified in a DNA-bound form and exhibits nucleoid-like structures under the atomic force microscope. It appears that the N-terminal sequence of Dps2 plays a role in nucleoid formation. MsDps2, unlike MsDps1, does not show elevated expression in nutritionally starved or stationary phase conditions; rather its promoter is recognized by RNA polymerase containing sigma(A) or sigma(B), under in vitro conditions. We propose that due to the nucleoid-condensing ability, the expression of MsDps2 is tightly regulated inside the cells.

Toward the Discovery of Vaccine Adjuvants: Coupling In Silico Screening and In Vitro Analysis of Antagonist Binding to Human and Mouse CCR4 Receptors

Background: Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses. Methodology: Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cellmigration, including that of CCR4(+) Tregs. Significance: Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.

The Role of UPF0157 in the Folding of M-tuberculosis Dephosphocoenzyme A Kinase and the Regulation of the Latter by CTP

Background: Targeting the biosynthetic pathway of Coenzyme A (CoA) for drug development will compromise multiple cellular functions of the tubercular pathogen simultaneously. Structural divergence in the organization of the penultimate and final enzymes of CoA biosynthesis in the host and pathogen and the differences in their regulation mark out the final enzyme, dephosphocoenzyme A kinase (CoaE) as a potential drug target. Methodology/Principal Findings: We report here a complete biochemical and biophysical characterization of the M. tuberculosis CoaE, an enzyme essential for the pathogen's survival, elucidating for the first time the interactions of a dephosphocoenzyme A kinase with its substrates, dephosphocoenzyme A and ATP; its product, CoA and an intrinsic yet novel inhibitor, CTP, which helps modulate the enzyme's kinetic capabilities providing interesting insights into the regulation of CoaE activity. We show that the mycobacterial enzyme is almost 21 times more catalytically proficient than its counterparts in other prokaryotes. ITC measurements illustrate that the enzyme follows an ordered mechanism of substrate addition with DCoA as the leading substrate and ATP following in tow. Kinetic and ITC experiments demonstrate that though CTP binds strongly to the enzyme, it is unable to participate in DCoA phosphorylation. We report that CTP actually inhibits the enzyme by decreasing its Vmax. Not surprisingly, a structural homology search for the modeled mycobacterial CoaE picks up cytidylmonophosphate kinases, deoxycytidine kinases, and cytidylate kinases as close homologs. Docking of DCoA and CTP to CoaE shows that both ligands bind at the same site, their interactions being stabilized by 26 and 28 hydrogen bonds respectively. We have also assigned a role for the universal Unknown Protein Family 0157 (UPF0157) domain in the mycobacterial CoaE in the proper folding of the full length enzyme. Conclusions/Significance: In view of the evidence presented, it is imperative to assign a greater role to the last enzyme of Coenzyme A biosynthesis in metabolite flow regulation through this critical biosynthetic pathway.

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.

Characterization of the Entamoeba histolytica Ornithine Decarboxylase-Like Enzyme

Background: The polyamines putrescine, spermidine, and spermine are organic cations that are required for cell growth and differentiation. Ornithine decarboxylase (ODC), the first and rate-limiting enzyme in the polyamine biosynthetic pathway, is a highly regulated enzyme. Methodology and Results: To use this enzyme as a potential drug target, the gene encoding putative ornithine decarboxylase (ODC)-like sequence was cloned from Entamoeba histolytica, a protozoan parasite causing amoebiasis. DNA sequence analysis revealed an open reading frame (ORF) of similar to 1,242 bp encoding a putative protein of 413 amino acids with a calculated molecular mass of 46 kDa and a predicted isoelectric point of 5.61. The E. histolytica putative ODC-like sequence has 33% sequence identity with human ODC and 36% identity with the Datura stramonium ODC. The ORF is a single-copy gene located on a 1.9-Mb chromosome. The recombinant putative ODC protein (48 kDa) from E. histolytica was heterologously expressed in Escherichia coli. Antiserum against recombinant putative ODC protein detected a band of anticipated size similar to 46 kDa in E. histolytica whole-cell lysate. Difluoromethylornithine (DFMO), an enzyme-activated irreversible inhibitor of ODC, had no effect on the recombinant putative ODC from E. histolytica. Comparative modeling of the three-dimensional structure of E. histolytica putative ODC shows that the putative binding site for DFMO is disrupted by the substitution of three amino acids-aspartate-332, aspartate-361, and tyrosine-323-by histidine-296, phenylalanine-305, and asparagine-334, through which this inhibitor interacts with the protein. Amino acid changes in the pocket of the E. histolytica enzyme resulted in low substrate specificity for ornithine. It is possible that the enzyme has evolved a novel substrate specificity. Conclusion: To our knowledge this is the first report on the molecular characterization of putative ODC-like sequence from E. histolytica. Computer modeling revealed that three of the critical residues required for binding of DFMO to the ODC enzyme are substituted in E. histolytica, resulting in the likely loss of interactions between the enzyme and DFMO.

Transcription of Human Resistin Gene Involves an Interaction of Sp1 with Peroxisome Proliferator-Activating Receptor Gamma (PPAR gamma)

Background: Resistin is a cysteine rich protein, mainly expressed and secreted by circulating human mononuclear cells. While several factors responsible for transcription of mouse resistin gene have been identified, not much is known about the factors responsible for the differential expression of human resistin.Methodology/Principal Finding: We show that the minimal promoter of human resistin lies within similar to 80 bp sequence upstream of the transcriptional start site (-240) whereas binding sites for cRel, CCAAT enhancer binding protein alpha (C/EBP-alpha), activating transcription factor 2 (ATF-2) and activator protein 1 (AP-1) transcription factors, important for induced expression, are present within sequences up to -619. Specificity Protein 1(Sp1) binding site (-276 to -295) is also present and an interaction of Sp1 with peroxisome proliferator activating receptor gamma (PPAR gamma) is necessary for constitutive expression in U937 cells. Indeed co-immunoprecipitation assay demonstrated a direct physical interaction of Sp1 with PPAR gamma in whole cell extracts of U937 cells. Phorbol myristate acetate (PMA) upregulated the expression of resistin mRNA in U937 cells by increasing the recruitment of Sp1, ATF-2 and PPAR gamma on the resistin gene promoter. Furthermore, PMA stimulation of U937 cells resulted in the disruption of Sp1 and PPAR gamma interaction. Chromatin immunoprecipitation (ChIP) assay confirmed the recruitment of transcription factors phospho ATF-2, Sp1, Sp3, PPAR gamma, chromatin modifier histone deacetylase 1 (HDAC1) and the acetylated form of histone H3 but not cRel, C/EBP-alpha and phospho c-Jun during resistingene transcription.Conclusion: Our findings suggest a complex interplay of Sp1 and PPAR gamma along with other transcription factors that drives the expression of resistin in human monocytic U937 cells.