Full genome ultra-deep pyrosequencing associates G-to-A hypermutation of the hepatitis B virus genome with the natural progression of hepatitis B.

Human APOBEC3 (A3) cytosine deaminases are antiviral restriction factors capable of editing the genome the hepatitis B virus (HBV). Despite the importance of the human A3 protein family for the innate immune response little is known about the clinical relevance for hepatitis B. The aim of this study was to utilize ultra-deep pyrosequencing (UDPS) data to analyse the phenomenon of G-to-A hypermutation of the complete HBV genome and to relate it to fundamental characteristics of patients with chronic hepatitis B. By analysing the viral population of 80 treatment naïve patients (47 HBeAg-positive and 33 HBeAg-negative), we identified an unequal distribution of G-to-A hypermutations across the genome. Our data indicate that G-to-A hypermutation occurs predominantly in a region between nucleotide positions 600 and 1800 a region which is usually single stranded in matured HBV particles. This implies that A3 likely edits HBV in the virion. Hypermutation rates for HBeAg-negative patients were more than 10-fold higher than those of HBeAg-positive patients. For HBeAg-negative patients higher hypermutation rates were significantly associated with the degree of fibrosis. Additionally, we found that for HBeAg-positive chronic hepatitis G-to-A hypermutation rates were significantly associated with the relative prevalence of the G1764A mutation, which is related to HBeAg seroconversion. In total, our data imply an important association of hypermutation mediated by A3 deaminases with the natural progression of chronic hepatitis B infections both in terms of HBeAg seroconversion and disease progression towards cirrhosis.

[Importance of bioinformatics in personalised medicine]. / Stellenwert der Bioinformatik für die personalisierte Medizin.

Medicine is experiencing a period of change: Extensive molecular biological data on the patient are increasingly included in diagnosis and treatment. This trend is based on the development of targeted drugs and accompanying diagnostics, which serve the purpose of providing advance evidence that the medication promises therapy success for the patient. According to this concept drugs are often given in combination. The sizes of patient groups for which a given therapy out of many possible alternatives can be expected to be successful are quite limited. The relationship between the molecular data pertaining to a patient and their disease phenotype are complex and cannot be determined manually. Thus, computer-based bioinformatics methods play a central role in interpreting the molecular data and as an instrument for providing recommendations for the practicing physician. Bioinformatics is an essential component in basic research, in the development of new concepts for diagnosis and therapy as well as in clinical practice, in which these concepts are applied to treating patients. This article discusses the role of bioinformatics in both basic research and clinical practice. We present the example of treatment of HIV patients, for which bioinformatics-assisted therapy selection has already entered clinical practice. Such a therapy concept is also predestined for other diseases (e.g., cancer). The article concludes with remarks on the prerequisites for society as a whole for ensuring success of this concept of personalised medicine as a factor of medical progress.

Prediction of response to antiretroviral therapy by human experts and by the EuResist data-driven expert system (the EVE study).

OBJECTIVES: The EuResist expert system is a novel data-driven online system for computing the probability of 8-week success for any given pair of HIV-1 genotype and combination antiretroviral therapy regimen plus optional patient information. The objective of this study was to compare the EuResist system vs. human experts (EVE) for the ability to predict response to treatment. METHODS: The EuResist system was compared with 10 HIV-1 drug resistance experts for the ability to predict 8-week response to 25 treatment cases derived from the EuResist database validation data set. All current and past patient data were made available to simulate clinical practice. The experts were asked to provide a qualitative and quantitative estimate of the probability of treatment success. RESULTS: There were 15 treatment successes and 10 treatment failures. In the classification task, the number of mislabelled cases was six for EuResist and 6-13 for the human experts [mean±standard deviation (SD) 9.1±1.9]. The accuracy of EuResist was higher than the average for the experts (0.76 vs. 0.64, respectively). The quantitative estimates computed by EuResist were significantly correlated (Pearson r=0.695, P<0.0001) with the mean quantitative estimates provided by the experts. However, the agreement among experts was only moderate (for the classification task, inter-rater κ=0.355; for the quantitative estimation, mean±SD coefficient of variation=55.9±22.4%). CONCLUSIONS: With this limited data set, the EuResist engine performed comparably to or better than human experts. The system warrants further investigation as a treatment-decision support tool in clinical practice.

Effect of ribavirin on the frequency of RNase L cleavage sites within the hepatitis C viral genome.

The mechanisms of synergy in antiviral activity of interferon-alpha and ribavirin in treating chronic hepatitis C virus (HCV) infection are still unknown. Interferon-alpha indirectly induces cleavage of viral RNA by RNase L at UU/UA dinucleotides. There is evidence that HCV genomes with a higher number of UU/UA dinucleotides are more sensitive to interferon-alpha. As a guanosine analogue, ribavirin exerts a mutagenic effect promoting G-to-A and C-to-U transitions. This study investigates whether ribavirin-induced mutagenesis causes a higher frequency of UU/UA dinucleotides in the viral progeny sequences. Increased mutational frequencies in favour of G-to-A and C-to-U transitions during ribavirin treatment was reported by Hofmann et al. (Gastroenterology 2007;132:921-930). Overall, 937 nucleotide sequences from that publication were reanalysed for RNase L cleavage sites. These included HCV NS3 quasispecies from three patients with ribavirin monotherapy and NS5B quasispecies from patients who received ribavirin alone (n = 7) or in combination with interferon-alpha (n = 7) at baseline and during treatment; NS5B quasispecies from a subgenomic HCV replicon system after 24, 48 and 72 h of cultivation with or without ribavirin or with levovirin. For NS3 quasispecies during ribavirin monotherapy and NS5B quasispecies from patients who received ribavirin alone or in combination with interferon-alpha, analysis of RNase L cleavage sites did not reveal changes during treatment or differences between treatment regimes. Similarly, RNaseL cleavage sites from NS5B quasispecies of the HCV replicon did not differ significantly between time points or treatments. In conclusion, Ribavirin-induced mutagenesis did not increase RNase L cleavage sites (UU/UA dinucleotides) within the HCV NS3 or NS5B encoding regions.

Genotypic coreceptor analysis.

HIV infects target cells by binding of its envelope gp120 protein to CD4 and a coreceptor on the cell surface. In vivo, the different HIV-strains use either CCR5 or CXCR4 as coreceptor. CCR5-using strains are named R5 viruses, while CXCR4-using strains are named X4. X4 viruses usually occur in the later stages. Coreceptor usage is a marker for disease progression. Additionally interest on coreceptors continually raises as a consequence of the development of a new class of antiretroviral drugs, namely the coreceptor antagonists or blockers. These specific drugs block the CCR5 or the CXCR4 coreceptors. So far, the CXCR4 blockers are not allowed to be used in the clinical practice due to their severe side effects. On the other hand, CCR5 blockers are currently in clinical practice, although they can only be administered after a baseline determination of the coreceptor usage of the predominant viral strain. Most of the coreceptor analyses in clinical cohorts have been performed with commercially available phenotypic assays. As for resistance testing of NRTIs, NNRTIs and PIs, efforts have also been made to predict the coreceptor usage from the genotype of the viruses. Different rules have been published based on the amino acid sequence of the Env-V3 region of HIV-gp120, which is known to be the major determinant of coreceptor usage. Among these, the most widely used is the 11/25 rule. Recently, bioinformatics driven prediction systems have been developed. Three of the interpretation systems are freely available via internet: WetCat, WebPSSM, geno2pheno[coreceptor]. All three systems focus on the Env-V3 region and take the amino acid sequence only into account. They learn from phenotypic and corresponding genotypic data. So far, two cohorts have been analyzed with such a genotypic approach and provided frequencies of R5 virus strains that are within the range of those reported with phenotypic assays. For one of the systems, geno2pheno[coreceptor], additional clinical data (e.g. CD4+T-cell counts) or structural information can be used to improve the prediction. Such genotypic systems provide the possibility for rapid screening of patients who may be administered with CCR5 blockers like the recently licensed Maraviroc.

Correlation of amino acid variations within nonstructural 4B protein with initial viral kinetics during interferon-alpha-based therapy in HCV-1b-infected patients.

Chronic hepatitis C is a major cause of liver cirrhosis leading to chronic liver failure and hepatocellular carcinoma. Different hepatitis C virus (HCV) proteins have been associated with resistance to interferon-alpha-based therapy. However, the exact mechanisms of virus-mediated interferon resistance are not completely understood. The importance of amino acid (aa) variations within the HCV nonstructural (NS)4B protein for replication efficiency and viral decline during the therapy is unknown. We investigated pretreatment sera from 42 patients with known outcome to interferon-based therapy. The complete NS4B gene was amplified and sequenced. Mutational analyses of predicted conformational, functional, structural and phylogenetic properties of the deduced aa sequences were performed. The complete NS4B protein was highly conserved with a median frequency of 0.015 +/- 0.009 aa exchanges (median +/- SD, 4.00 +/- 2.31). Especially within the predicted transmembranous domains of the NS4B protein, the mean number of aa variations was low (median frequency, 0.013 +/- 0.013). Neither the number of aa variations nor specific aa exchanges were correlated with HCV RNA serum concentration at baseline. A rapid initial HCV RNA decline of >/=1.5 log(10) IU/mL at week 2 of interferon-based therapy was associated with a higher frequency of nonconservative aa exchanges within the complete NS4B protein in comparison with patients with a nonrapid HCV RNA decline (median frequency, 0.011 +/- 0.005 vs 0.004 +/- 0.003, P = 0.006). Overall, the aa sequence of the NS4B protein was highly conserved, indicating an important role for replication in vivo. Amino acid variations with relevant changes of physicochemical properties may influence replication efficiency, associated with a rapid early virological response.

Predicting protein structure classes from function predictions.

MOTIVATION: We introduce a new approach to using the information contained in sequence-to-function prediction data in order to recognize protein template classes, a critical step in predicting protein structure. The data on which our method is based comprise probabilities of functional categories; for given query sequences these probabilities are obtained by a neural net that has previously been trained on a variety of functionally important features. On a training set of sequences we assess the relevance of individual functional categories for identifying a given structural family. Using a combination of the most relevant categories, the likelihood of a query sequence to belong to a specific family can be estimated. RESULTS: The performance of the method is evaluated using cross-validation. For a fixed structural family and for every sequence, a score is calculated that measures the evidence for family membership. Even for structural families of small size, family members receive significantly higher scores. For some examples, we show that the relevant functional features identified by this method are biologically meaningful. The proposed approach can be used to improve existing sequence-to-structure prediction methods. AVAILABILITY: Matlab code is available on request from the authors. The data are available at http://www.mpisb.mpg.de/~sommer/Fun2Struc/

The Helmholtz Network for Bioinformatics: an integrative web portal for bioinformatics resources.

SUMMARY: The Helmholtz Network for Bioinformatics (HNB) is a joint venture of eleven German bioinformatics research groups that offers convenient access to numerous bioinformatics resources through a single web portal. The 'Guided Solution Finder' which is available through the HNB portal helps users to locate the appropriate resources to answer their queries by employing a detailed, tree-like questionnaire. Furthermore, automated complex tool cascades ('tasks'), involving resources located on different servers, have been implemented, allowing users to perform comprehensive data analyses without the requirement of further manual intervention for data transfer and re-formatting. Currently, automated cascades for the analysis of regulatory DNA segments as well as for the prediction of protein functional properties are provided. AVAILABILITY: The HNB portal is available at http://www.hnbioinfo.de

Centralization: a new method for the normalization of gene expression data.

Microarrays measure values that are approximately proportional to the numbers of copies of different mRNA molecules in samples. Due to technical difficulties, the constant of proportionality between the measured intensities and the numbers of mRNA copies per cell is unknown and may vary for different arrays. Usually, the data are normalized (i.e., array-wise multiplied by appropriate factors) in order to compensate for this effect and to enable informative comparisons between different experiments. Centralization is a new two-step method for the computation of such normalization factors that is both biologically better motivated and more robust than standard approaches. First, for each pair of arrays the quotient of the constants of proportionality is estimated. Second, from the resulting matrix of pairwise quotients an optimally consistent scaling of the samples is computed.

Derivation of a scoring function for crystal structure prediction.

The ever increasing number of experimentally resolved crystal structures supports the possibility of fully empirical crystal structure prediction for small organic molecules. Empirical methods promise to be significantly more efficient than methods that attempt to solve the same problem from first principles. However, the transformation from data to empirical knowledge and further to functional algorithms is not trivial and the usefulness of the result depends strongly on the quantity and the quality of the data. In this work, a simple scoring function is parameterized to discriminate between the correct structure and a set of decoys for a large number of different molecular systems. The method is fully automatic and has the advantage that the complete scoring function is parametrized at once, leading to a self-consistent set of parameters. The obtained scoring function is tested on an independent set of crystal structures taken from the P1 and P1; space groups. With the trained scoring function and FlexCryst, a program for small-molecule crystal structure prediction, it is shown that approximately 73% of the 239 tested molecules in space group P1 are predicted correctly. For the more complex space group P1;, the success rate is 26%. Comparison with force-field potentials indicates the physical content of the obtained scoring function, a result of direct importance for protein threading where such database-based potentials are being applied.