Estimating Frequencies of Minority Nevirapine-Resistant Strains in Chronically HIV-1-Infected Individuals Naive to Nevirapine by Using Stochastic Simulations and a Mathematical Model

Nevirapine forms the mainstay of our efforts to curtail the pediatric AIDS epidemic through prevention of mother-to-child transmission of HIV-1. A key limitation, however, is the rapid selection of HIV-1 strains resistant to nevirapine following the administration of a single dose. This rapid selection of resistance suggests that nevirapine-resistant strains preexist in HIV-1 patients and may adversely affect outcomes of treatment. The frequencies of nevirapine-resistant strains in vivo, however, remain poorly estimated, possibly because they exist as a minority below current assay detection limits. Here, we employ stochastic simulations and a mathematical model to estimate the frequencies of strains carrying different combinations of the common nevirapine resistance mutations K103N, V106A, Y181C, Y188C, and G190A in chronically infected HIV-1 patients naive to nevirapine. We estimate the relative fitness of mutant strains from an independent analysis of previous competitive growth assays. We predict that single mutants are likely to preexist in patients at frequencies (similar to 0.01% to 0.001%) near or below current assay detection limits (>0.01%), emphasizing the need for more-sensitive assays. The existence of double mutants is subject to large stochastic variations. Triple and higher mutants are predicted not to exist. Our estimates are robust to variations in the recombination rate, cellular superinfection frequency, and the effective population size. Thus, with 10(7) to 10(8) infected cells in HIV-1 patients, even when undetected, nevirapine-resistant genomes may exist in substantial numbers and compromise efforts to prevent mother-to-child transmission of HIV-1, accelerate the failure of subsequent antiretroviral treatments, and facilitate the transmission of drug resistance.

Systems biology of tuberculosis

Systems biology seeks to study biological systems as a whole, by adopting an integrated approach to study and understand the function of biological systems, particularly, the response of such systems to various perturbations. In this article, we focus on the Indian efforts towards systems-level studies of Mycobacterium tuberculosis and its interaction with the host. Availability of a variety of genome-scale experimental data, providing first level `omics' descriptions of the pathogen, render it feasible to study it at a systems level. Various aspects of the pathogen, from metabolic pathways to protein-protein interaction networks have been modelled and simulated, while host-pathogen interactions have been studied experimentally using siRNA-based techniques. These studies have been useful in obtaining a global perspective of the pathogen and its interactions with the host in many ways. For example, significant insights have been gained about different aspects such as proteins essential for bacterial survival, proteins that are highly influential in the network, pathways that are highly connected, host factors responsible for maintaining the TB infection and key factors involved in autophagy and pathogenesis. A rational pipeline developed for drug target identification incorporating analyses of the interactome, reactome, genome, pocketome and the transcriptome is discussed. Finally, exploring host factors as drug targets and insights about the emergence of drug resistance are also discussed. (C) 2011 Elsevier Ltd. All rights reserved.

RAD51 as a potential biomarker and therapeutic target for pancreatic cancer

Chemotherapy is a very important therapeutic strategy for cancer treatment. The failure of conventional and molecularly targeted chemotherapeutic regimes for the treatment of pancreatic cancer highlights a desperate need for novel therapeutic interventions. Chemotherapy often fails to eliminate all tumor cells because of intrinsic or acquired drug resistance, which is the most common cause of tumor recurrence. Overexpression of RAD51 protein, a key player in DNA repair/recombination has been observed in many cancer cells and its hyperexpression is implicated in drug resistance. Recent studies suggest that RAD51 overexpression contributes to the development, progression and drug resistance of pancreatic cancer cells. Here we provide a brief overview of the available pieces of evidence in support of the role of RAD51 in pancreatic tumorigenesis and drug resistance, and hypothesize that RAD51 could serve as a potential biomarker for diagnosis of pancreatic cancer. We discuss the possible involvement of RAD51 in the drug resistance associated with epithelial to mesenchymal transition and with cancer stem cells. Finally, we speculate that targeting RAD51 in pancreatic cancer cells may be a novel approach for the treatment of pancreatic cancer. (C) 2011 Elsevier B.V. All rights reserved.

An Antitubulin Agent BCFMT Inhibits Proliferation of Cancer Cells and Induces Cell Death by Inhibiting Microtubule Dynamics

Using cell based screening assay, we identified a novel anti-tubulin agent (Z)-5-((5-(4-bromo-3-chlorophenyl)furan-2-yl)methylene)-2-thioxothiazoli din-4-one (BCFMT) that inhibited proliferation of human cervical carcinoma (HeLa) (IC50, 7.2 +/- 1.8 mu M), human breast adenocarcinoma (MCF-7) (IC50, 10.0 +/- 0.5 mu M), highly metastatic breast adenocarcinoma (MDA-MB-231) (IC50, 6.0 +/- 1 mu M), cisplatin-resistant human ovarian carcinoma (A2780-cis) (IC50, 5.8 +/- 0.3 mu M) and multi-drug resistant mouse mammary tumor (EMT6/AR1) (IC50, 6.5 +/- 1 mu M) cells. Using several complimentary strategies, BCFMT was found to inhibit cancer cell proliferation at G2/M phase of the cell cycle apparently by targeting microtubules. In addition, BCFMT strongly suppressed the dynamics of individual microtubules in live MCF-7 cells. At its half maximal proliferation inhibitory concentration (10 mu M), BCFMT reduced the rates of growing and shortening phases of microtubules in MCF-7 cells by 37 and 40%, respectively. Further, it increased the time microtubules spent in the pause (neither growing nor shortening detectably) state by 135% and reduced the dynamicity (dimer exchange per unit time) of microtubules by 70%. In vitro, BCFMT bound to tubulin with a dissociation constant of 8.3 +/- 1.8 mu M, inhibited tubulin assembly and suppressed GTPase activity of microtubules. BCFMT competitively inhibited the binding of BODIPY FL-vinblastine to tubulin with an inhibitory concentration (K-i) of 5.2 +/- 1.5 mu M suggesting that it binds to tubulin at the vinblastine site. In cultured cells, BCFMT-treatment depolymerized interphase microtubules, perturbed the spindle organization and accumulated checkpoint proteins (BubR1 and Mad2) at the kinetochores. BCFMT-treated MCF-7 cells showed enhanced nuclear accumulation of p53 and its downstream p21, which consequently activated apoptosis in these cells. The results suggested that BCFMT inhibits proliferation of several types of cancer cells including drug resistance cells by suppressing microtubule dynamics and indicated that the compound may have chemotherapeutic potential.

Stochastic Simulations Suggest that HIV-1 Survives Close to Its Error Threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, mu(c), however, is not known. Application of the quasispecies theory to determine mu(c) poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and mu(c). We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated mu(c) to be 7 x 10(-5) -1 x 10(-4) substitutions/site/replication, similar to 2-6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, mu(c) increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of mu(c) may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Multimodal therapy for complete regression of malignant melanoma using constrained nonlinear optimal dynamic inversion

Using a realistic nonlinear mathematical model for melanoma dynamics and the technique of optimal dynamic inversion (exact feedback linearization with static optimization), a multimodal automatic drug dosage strategy is proposed in this paper for complete regression of melanoma cancer in humans. The proposed strategy computes different drug dosages and gives a nonlinear state feedback solution for driving the number of cancer cells to zero. However, it is observed that when tumor is regressed to certain value, then there is no need of external drug dosages as immune system and other therapeutic states are able to regress tumor at a sufficiently fast rate which is more than exponential rate. As model has three different drug dosages, after applying dynamic inversion philosophy, drug dosages can be selected in optimized manner without crossing their toxicity limits. The combination of drug dosages is decided by appropriately selecting the control design parameter values based on physical constraints. The process is automated for all possible combinations of the chemotherapy and immunotherapy drug dosages with preferential emphasis of having maximum possible variety of drug inputs at any given point of time. Simulation study with a standard patient model shows that tumor cells are regressed from 2 x 107 to order of 105 cells because of external drug dosages in 36.93 days. After this no external drug dosages are required as immune system and other therapeutic states are able to regress tumor at greater than exponential rate and hence, tumor goes to zero (less than 0.01) in 48.77 days and healthy immune system of the patient is restored. Study with different chemotherapy drug resistance value is also carried out. (C) 2014 Elsevier Ltd. All rights reserved.

Bmi1 regulates self-renewal and epithelial to mesenchymal transition in breast cancer cells through Nanog

Background: The Bmi1 polycomb ring finger oncogene, a transcriptional repressor belonging to the Polycomb group of proteins plays an important role in the regulation of stem cell self-renewal and is elevated in several cancers. In the current study, we have explored the role of Bmi1 in regulating the stemness and drug resistance of breast cancer cells. Methods: Using real time PCR and immunohistochemistry primary breast tissues were analyzed. Retro-and lentiviruses were utilized to overexpress and knockdown Bmi1, RT-PCR and Western blot was performed to evaluate mRNA and protein expression. Stemness properties were analyzed by flow cytometry and sphere-formation and tumor formation was determined by mouse xenograft experiments. Dual luciferase assay was employed to assess promoter activity and MTT assay was used to analyze drug response. Results: We found Bmi1 overexpression in 64% of grade III invasive ductal breast adenocarcinomas compared to normal breast tissues. Bmi1 overexpression in immortalized and transformed breast epithelial cells increased their sphere-forming efficiency, induced epithelial to mesenchymal transition ( EMT) with an increase in the expression of stemness-related genes. Knockdown of Bmi1 in tumorigenic breast cells induced epithelial morphology, reduced expression of stemness-related genes, decreased the IC50 values of doxorubicin and abrogated tumor-formation. Bmi1-high tumors showed elevated Nanog expression whereas the tumors with lower Bmi1 showed reduced Nanog levels. Overexpression of Bmi1 increased Nanog levels whereas knockdown of Bmi1 reduced its expression. Dual luciferase promoter-reporter assay revealed Bmi1 positively regulated the Nanog and NF kappa B promoter activity. RT-PCR analysis showed that Bmi1 overexpression activated the NF kappa B pathway whereas Bmi1 knockdown reduced the expression of NF kappa B target genes, suggesting that Bmi1 might regulate Nanog expression through the NF kappa B pathway. Conclusions: Our study showed that Bmi1 is overexpressed in several high-grade, invasive ductal breast adenocarcinomas, thus supporting its role as a prognostic marker. While Bmi1 overexpression increased self-renewal and promoted EMT, its knockdown reversed EMT, reduced stemness, and rendered cells drug sensitive, thus highlighting a crucial role for Bmi1 in regulating the stemness and drug response of breast cancer cells. Bmi1 may control self-renewal through the regulation of Nanog expression via the NF kappa B pathway.

Exploring the potential of adjunct therapy in tuberculosis

A critical unmet need for treatment of drug-resistant tuberculosis (TB) is to find novel therapies that are efficacious, safe, and shorten the duration of treatment. Drug discovery approaches for TB primarily target essential genes of the pathogen Mycobacterium tuberculosis (Mtb) but novel strategies such as host-directed therapies and nonmicrobicidal targets are necessary to bring about a paradigm shift in treatment. Drugs targeting the host pathways and nonmicrobicidal proteins can be used only in conjunction with existing drugs as adjunct therapies. Significantly, host-directed adjunct therapies have the potential to decrease duration of treatment, as they are less prone to drug resistance, target the immune responses, and act via novel mechanism of action. Recent advances in targeting host-pathogen interactions have implicated pathways such as eicosanoid regulation and angiogenesis. Furthermore, several approved drugs such as metformin and verapamil have been identified that appear suitable for repurposing for the treatment of TB. These findings and the challenges in the area of host- and/or pathogen-directed adjunct therapies and their implications for TB therapy are discussed.

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.

Constraint effects on multiple void interaction in pressure sensitive plastic solids

In this work, effects of pressure sensitive yielding and plastic dilatancy on void growth and void interaction mechanisms in fracture specimens displaying high and low constraint levels are investigated. To this end, large deformation finite element simulations are carried out with discrete voids ahead of the notch. It is observed that multiple void interaction mechanism which is favored by high initial porosity is further accelerated by pressure sensitive yielding, but is retarded by loss of constraint. The resistance curves predicted based on a simple void coalescence criterion show enhancement in fracture resistance when constraint level is low and when pressure sensitivity is suppressed.

Influence of surface chemistry of mesoporous alumina with wide pore distribution on controlled drug release

The crucial role of the drug carrier surface chemical moeities on the uptake and in vitro release of drug is discussed here in a systematic manner. Mesoporous alumina with a wide pore size distribution (2-7 nm) functionalized with various hydrophilic and hydrophobic surface chemical groups was employed as the carrier for delivery of the model drug ibuprofen. Surface functionalization with hydrophobic groups resulted in low degree of drug loading (approximately 20%) and fast rate of release (85% over a period of 5 h) whereas hydrophilic groups resulted in a significantly higher drug payloads (21%-45%) and slower rate of release (12%-40% over a period of 5 h). Depending on the chemical moiety, the diffusion controlled (proportional to time(-0.5)) drug release was additionally observed to be dependent on the mode of arrangement of the functional groups on the alumina surface as well as on the pore characteristics of the matrix. For all mesoporous alumina systems the drug dosages were far lower than the maximum recommended therapeutic dosages (MRTD) for oral delivery. We envisage that the present study would aid in the design of delivery systems capable of sustained release of multiple drugs.

Involvement of Cytochrome P450 in Conferring Chloroquine Resistance to the Malarial Parasite, Plasmodium falciparum

The higher levels of cytochrone P-450 dependent enzyme activities reported earlier are traced to higher levels of cytochrome P-450 (CYPIIB1/B2 like) messenger RNA in the chloroquine resistant than the sensitive strains. The messenger RNA is also induced by phenobarbitone in the sensitive strain. Pretreatment with phenobarbitone affords partial protection to chloroquine toxicity in the sensitive strain and this is not due to a differential accumulation of the drug.

Open source drug discovery- A new paradigm of collaborative research in tuberculosis drug development

It is being realized that the traditional closed-door and market driven approaches for drug discovery may not be the best suited model for the diseases of the developing world such as tuberculosis and malaria, because most patients suffering from these diseases have poor paying capacity. To ensure that new drugs are created for patients suffering from these diseases, it is necessary to formulate an alternate paradigm of drug discovery process. The current model constrained by limitations for collaboration and for sharing of resources with confidentiality hampers the opportunities for bringing expertise from diverse fields. These limitations hinder the possibilities of lowering the cost of drug discovery. The Open Source Drug Discovery project initiated by Council of Scientific and Industrial Research, India has adopted an open source model to power wide participation across geographical borders. Open Source Drug Discovery emphasizes integrative science through collaboration, open-sharing, taking up multi-faceted approaches and accruing benefits from advances on different fronts of new drug discovery. Because the open source model is based on community participation, it has the potential to self-sustain continuous development by generating a storehouse of alternatives towards continued pursuit for new drug discovery. Since the inventions are community generated, the new chemical entities developed by Open Source Drug Discovery will be taken up for clinical trial in a non-exclusive manner by participation of multiple companies with majority funding from Open Source Drug Discovery. This will ensure availability of drugs through a lower cost community driven drug discovery process for diseases afflicting people with poor paying capacity. Hopefully what LINUX the World Wide Web have done for the information technology, Open Source Drug Discovery will do for drug discovery. (C) 2011 Elsevier Ltd. All rights reserved.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

Ductile fracture by multiple void growth and interaction ahead of a notch tip in polycrystalline plastic solids

The objectives of this paper are to study the effects of plastic anisotropy and evolution in crystallographic texture with deformation on the ductile fracture behaviour of polycrystalline solids. To this end, numerical simulations of multiple void growth and interaction ahead of a notch tip are performed under mode I, plane strain, small scale yielding conditions using two approaches. The first approach is based on the Hill yield theory, while the second employs crystal plasticity constitutive equations and a Taylor-type homogenization in order to represent the ductile polycrystalline solid. The initial textures pertaining to continuous cast Al-Mg AA5754 sheets in recrystallized and cold rolled conditions are considered. The former is nearly-isotropic, while the latter displays pronounced anisotropy. The results indicate distinct changes in texture in the ligaments bridging the voids ahead of the notch tip with increase in load level which gives rise to retardation in porosity evolution and increase in tearing resistance for both materials.