In Silico Prediction of the Mechanobiological Response of Arterial Tissue: Application to Angioplasty and Stenting

One way to restore physiological blood flow to occluded arteries involves the deformation of plaque using an intravascular balloon and preventing elastic recoil using a stent. Angioplasty and stent implantation cause unphysiological loading of the arterial tissue, which may lead to tissue in-growth and reblockage; termed “restenosis.” In this paper, a computational methodology for predicting the time-course of restenosis is presented. Stress-induced damage, computed using a remaining life approach, stimulates inflammation (production of matrix degrading factors and growth stimuli). This, in turn, induces a change in smooth muscle cell phenotype from contractile (as exists in the quiescent tissue) to synthetic (as exists in the growing tissue). In this paper, smooth muscle cell activity (migration, proliferation, and differentiation) is simulated in a lattice using a stochastic approach to model individual cell activity. The inflammation equations are examined under simplified loading cases. The mechanobiological parameters of the model were estimated by calibrating the model response to the results of a balloon angioplasty study in humans. The simulation method was then used to simulate restenosis in a two dimensional model of a stented artery. Cell activity predictions were similar to those observed during neointimal hyperplasia, culminating in the growth of restenosis. Similar to experiment, the amount of neointima produced increased with the degree of expansion of the stent, and this relationship was found to be highly dependant on the prescribed inflammatory response. It was found that the duration of inflammation affected the amount of restenosis produced, and that this effect was most pronounced with large stent expansions. In conclusion, the paper shows that the arterial tissue response to mechanical stimulation can be predicted using a stochastic cell modeling approach, and that the simulation captures features of restenosis development observed with real stents. The modeling approach is proposed for application in three dimensional models of cardiovascular stenting procedures.

In silico prediction of the effects of mutations in the human UDP-galactose 4'-epimerase gene:Towards a predictive framework for type III galactosemia

The enzyme UDP-galactose 4'-epimerase (GALE) catalyses the reversible epimerisation of both UDP-galactose and UDP-N-acetyl-galactosamine. Deficiency of the human enzyme (hGALE) is associated with type III galactosemia. The majority of known mutations in hGALE are missense and private thus making clinical guidance difficult. In this study a bioinformatics approach was employed to analyse the structural effects due to each mutation using both the UDP-glucose and UDP-N-acetylglucosamine bound structures of the wild-type protein. Changes to the enzyme's overall stability, substrate/cofactor binding and propensity to aggregate were also predicted. These predictions were found to be in good agreement with previous in vitro and in vivo studies when data was available and allowed for the differentiation of those mutants that severely impair the enzyme's activity against UDP-galactose. Next this combination of techniques were applied to another twenty-six reported variants from the NCBI dbSNP database that have yet to be studied to predict their effects. This identified p.I14T, p.R184H and p.G302R as likely severely impairing mutations. Although severely impaired mutants were predicted to decrease the protein's stability, overall predicted stability changes only weakly correlated with residual activity against UDP-galactose. This suggests other protein functions such as changes in cofactor and substrate binding may also contribute to the mechanism of impairment. Finally this investigation shows that this combination of different in silico approaches is useful in predicting the effects of mutations and that it could be the basis of an initial prediction of likely clinical severity when new hGALE mutants are discovered.

Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

Proteomics and in silico approaches to extend understanding of the glutathione transferase superfamily of the tropical liver fluke Fasciola gigantica

Fasciolosis is an important foodborne, zoonotic disease of livestock and humans, with global annual health and economic losses estimated at several billion US$. Fasciola hepatica is the major species in temperate regions, while F. gigantica dominates in the tropics. In the absence of commercially available vaccines to control fasciolosis, increasing reports of resistance to current chemotherapeutic strategies and the spread of fasciolosis into new areas, new functional genomics approaches are being used to identify potential new drug targets and vaccine candidates. The glutathione transferase (GST) superfamily is both a candidate drug and vaccine target. This study reports the identification of a putatively novel Sigma class GST, present in a water-soluble cytosol extract from the tropical liver fluke F. gigantica. The GST was cloned and expressed as an enzymically active recombinant protein. This GST shares a greater identity with the human schistosomiasis GST vaccine currently at Phase II clinical trials than previously discovered F. gigantica GSTs, stimulating interest in its immuno-protective properties. In addition, in silico analysis of the GST superfamily of both F. gigantica and F. hepatica has revealed an additional Mu class GST, Omega class GSTs, and for the first time, a Zeta class member.

In silico analysis and DHPLC screening strategy identifies novel apoptotic gene targets of aberrant promoter hypermethylation in prostate cancer

BACKGROUND: Aberrant DNA methylation has been implicated as a key survival mechanism in cancer, whereby promoter hypermethylation silences genes essential for many cellular processes including apoptosis. Limited data is available on the methylation profile of apoptotic genes in prostate cancer (CaP). The aim of this study was to profile methylation of apoptotic-related genes in CaP using denaturing high performance liquid chromatography (DHPLC).

METHODS: Based on an in silico selection process, 13 genes were screened for methylation in CaP cell lines using DHPLC. Quantitative methylation specific PCR was employed to determine methylation levels in prostate tissue specimens (n = 135), representing tumor, histologically benign prostate, high-grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. Gene expression was measured by QRT-PCR in cell lines and tissue specimens.

RESULTS: The promoters of BIK, BNIP3, cFLIP, TMS1, DCR1, DCR2, and CDKN2A appeared fully or partially methylated in a number of malignant cell lines. This is the first report of aberrant methylation of BIK, BNIP3, and cFLIP in CaP. Quantitative methylation analysis in prostate tissues identified 5 genes (BNIP3, CDKN2A, DCR1, DCR2 and TMS1) which were frequently methylated in tumors but were unmethylated in 100% of benign tissues. Furthermore, 69% of tumors were methylated in at least one of the five-gene panel. In the case of all genes, except BNIP3, promoter hypermethylation was associated with concurrent downregulation of gene expression.

CONCLUSION: Future examination of this "CaP apoptotic methylation signature" in a larger cohort of patients is justified to further evaluate its value as a diagnostic and prognostic marker.

In silico mining identifies IGFBP3 as a novel target of methylation in prostate cancer

Promoter hypermethylation is central in deregulating gene expression in cancer. Identification of novel methylation targets in specific cancers provides a basis for their use as biomarkers of disease occurrence and progression. We developed an in silico strategy to globally identify potential targets of promoter hypermethylation in prostate cancer by screening for 5' CpG islands in 631 genes that were reported as downregulated in prostate cancer. A virtual archive of 338 potential targets of methylation was produced. One candidate, IGFBP3, was selected for investigation, along with glutathione-S-transferase pi (GSTP1), a well-known methylation target in prostate cancer. Methylation of IGFBP3 was detected by quantitative methylation-specific PCR in 49/79 primary prostate adenocarcinoma and 7/14 adjacent preinvasive high-grade prostatic intraepithelial neoplasia, but in only 5/37 benign prostatic hyperplasia (P < 0.0001) and in 0/39 histologically normal adjacent prostate tissue, which implies that methylation of IGFBP3 may be involved in the early stages of prostate cancer development. Hypermethylation of IGFBP3 was only detected in samples that also demonstrated methylation of GSTP1 and was also correlated with Gleason score > or =7 (P=0.01), indicating that it has potential as a prognostic marker. In addition, pharmacological demethylation induced strong expression of IGFBP3 in LNCaP prostate cancer cells. Our concept of a methylation candidate gene bank was successful in identifying a novel target of frequent hypermethylation in early-stage prostate cancer. Evaluation of further relevant genes could contribute towards a methylation signature of this disease.

Genetic regulation of the ramA locus and its expression in clinical isolates of Klebsiella pneumoniae

Tigecycline resistance has been attributed to ramA overexpression and subsequent acrA upregulation. The ramA locus, originally identified in Klebsiella pneumoniae, has homologues in Enterobacter and Salmonella spp. In this study, we identify in silico that the ramR binding site is also present in Citrobacter spp. and that Enterobacter, Citrobacter and Klebsiella spp. share key regulatory elements in the control of the romA–ramA locus. RACE (rapid amplification of cDNA ends) mapping indicated that there are two promoters from which romA–ramA expression can be regulated in K. pneumoniae. Correspondingly, electrophoretic binding studies clearly showed that purified RamA and RamR proteins bind to both of these promoters. Hence, there appear to be two RamR binding sites within the Klebsiella romA–ramA locus. Like MarA, RamA binds the promoter region, implying that it might be subject to autoregulation. We have identified changes within ramR in geographically distinct clinical isolates of K. pneumoniae. Intriguingly, levels of romA and ramA expression were not uniformly affected by changes within the ramR gene, thereby supporting the dual promoter finding. Furthermore, a subset of strains sustained no changes within the ramR gene but which still overexpressed the romA–ramA genes, strongly suggesting that a secondary regulator may control ramA expression.

Establishing Clonal Cell Lines with Endothelial-Like Potential from CD9(hi), SSEA-1(-) Cells in Embryonic Stem Cell-Derived Embryoid Bodies

Background. Differentiation of embryonic stem cells (ESCs) into specific cell types with minimal risk of teratoma formation could be efficiently directed by first reducing the differentiation potential of ESCs through the generation of clonal, self-renewing lineage-restricted stem cell lines. Efforts to isolate these stem cells are, however, mired in an impasse where the lack of purified lineage-restricted stem cells has hindered the identification of defining markers for these rare stem cells and, in turn, their isolation. Methodology/Principal Findings. We describe here a method for the isolation of clonal lineage-restricted cell lines with endothelial potential from ESCs through a combination of empirical and rational evidence-based methods. Using an empirical protocol that we have previously developed to generate embryo-derived RoSH lines with endothelial potential, we first generated E-RoSH lines from mouse ESC-derived embryoid bodies (EBs). Despite originating from different mouse strains, RoSH and E-RoSH lines have similar gene expression profiles (r(2) = 0.93) while that between E-RoSH and ESCs was 0.83. In silico gene expression analysis predicted that like RoSH cells, E-RoSH cells have an increased propensity to differentiate into vasculature. Unlike their parental ESCs, E-RoSH cells did not form teratomas and differentiate efficiently into endothelial-like cells in vivo and in vitro. Gene expression and FACS analysis revealed that RoSH and E-RoSH cells are CD9(hi), SSEA-1(-) while ESCs are CD9(lo), SSEA-1(+). Isolation of CD9(hi), SSEA-1(-) cells that constituted 1%-10% of EB-derived cultures generated an E-RoSH-like culture with an identical E-RoSH-like gene expression profile (r(2) = 0.95) and a propensity to differentiate into endothelial-like cells. Conclusions. By combining empirical and rational evidence-based methods, we identified definitive selectable surface antigens for the isolation and propagation of lineage-restricted stem cells with endothelial-like potential from mouse ESCs.

Towards accurate estimation of the proportion of true null hypotheses in multiple testing

Background

Biomedical researchers are now often faced with situations where it is necessary to test a large number of hypotheses simultaneously, eg, in comparative gene expression studies using high-throughput microarray technology. To properly control false positive errors the FDR (false discovery rate) approach has become widely used in multiple testing. The accurate estimation of FDR requires the proportion of true null hypotheses being accurately estimated. To date many methods for estimating this quantity have been proposed. Typically when a new method is introduced, some simulations are carried out to show the improved accuracy of the new method. However, the simulations are often very limited to covering only a few points in the parameter space.

Results

Here I have carried out extensive in silico experiments to compare some commonly used methods for estimating the proportion of true null hypotheses. The coverage of these simulations is unprecedented thorough over the parameter space compared to typical simulation studies in the literature. Thus this work enables us to draw conclusions globally as to the performance of these different methods. It was found that a very simple method gives the most accurate estimation in a dominantly large area of the parameter space. Given its simplicity and its overall superior accuracy I recommend its use as the first choice for estimating the proportion of true null hypotheses in multiple testing.

BRCA1 transcriptionally regulates genes associated with the basal-like phenotype in breast cancer

Expression profiling of BRCA1-deficient tumours has identified a pattern of gene expression similar to basal-like breast tumours. In this study, we examine whether a BRCA1-dependent transcriptional mechanism may underpin the link between BRCA1 and basal-like phenotype. In methods section, the mRNA and protein were harvested from a number of BRCA1 mutant and wild-type breast cancer cell lines and from matched isogenic controls. Microarray-based expression profiling was used to identify potential BRCA1-regulated transcripts. These gene targets were then validated (by in silico analysis of tumour samples) by real-time PCR and Western blot analysis. Chromatin immunoprecipitation (ChIP) assays were used to confirm recruitment of BRCA1 to specific promoters. In results, we demonstrate that functional BRCA1 represses the expression of cytokeratins 5(KRT5) and 17(KRT17) and p-Cadherin (CDH3) in HCC1937 and T47D breast cancer cell lines at both mRNA and protein level. ChIP assays demonstrate that BRCA1 is recruited to the promoters of KRT5, KRT17 and CDH3, and re-ChIP assays confirm that BRCA1 is recruited independently to form c-Myc and Sp1 complexes on the CDH3 promoter. We show that siRNA-mediated inhibition of endogenous c-Myc (and not Sp1) results in a marked increase in CDH3 expression analogous to that observed following the inhibition of endogenous BRCA1. The data provided suggest a model whereby BRCA1 and c-Myc form a repressor complex on the promoters of specific basal genes and represent a potential mechanism to explain the observed overexpression of key basal markers in BRCA1-deficient tumours.

The WaaL O-antigen lipopolysaccharide ligase has features in common with metal ion-independent inverting glycosyltransferases(*)

WaaL is a membrane enzyme that catalyzes a key step in lipopolysaccharide (LPS) synthesis: the glycosidic bonding of a sugar at the proximal end of the undecaprenyl-diphosphate (Und-PP) O-antigen with a terminal sugar of the lipid A-core oligosaccharide (OS). Utilizing an in vitro assay, we demonstrate here that ligation with purified Escherichia coli WaaL occurs without adenosine-5'-triphosphate (ATP) and magnesium ions. Furthermore, E. coli and Pseudomonas aeruginosa WaaL proteins cannot catalyze ATP hydrolysis in vitro. We also show that a lysine substitution of the arginine (Arg)-215 residue renders an active protein, whereas WaaL mutants with alanine replacements in the periplasmic-exposed residues Arg-215, Arg-288 and histidine (His)-338 and also the membrane-embedded aspartic acid-389 are nonfunctional. An in silico approach, combining predicted topological information with the analysis of sequence conservation, confirms the importance of a positive charge at the small periplasmic loop of WaaL, since an Arg corresponding to Arg-215 was found at a similar position in all the WaaL homologs. Also, a universally conserved H[NSQ]X(9)GXX[GTY] motif spanning the C-terminal end of the predicted large periplasmic loop and the membrane boundary of the transmembrane helix was identified. The His residue in this motif corresponds to His-338. A survey of LPS structures in which the linkage between O-antigen and lipid A-core OS was elucidated reveals that it is always in the beta-configuration, whereas the sugars bound to Und-PP are in the alpha-configuration. Together, our biochemical and in silico data argue that WaaL proteins use a common reaction mechanism and share features of metal ion-independent inverting glycosyltransferases.