Coordinated evolution of the Hepatitis B Virus Polymerase

The detection of compensatory mutations that abrogate negative fitness effects of drug-resistance and vaccine-escape mutations indicates the important role of epistatic connectivity in evolution of viruses, especially under the strong selection pressures. Mapping of epistatic connectivity in the form of coordinated substitutions should help to characterize molecular mechanisms shaping viral evolution and provides a tool for the development of novel anti-viral drugs and vaccines. We analyzed coordinated variation among amino acid sites in 370 the hepatitis B virus (HBV) polymerase sequences using Bayesian networks. Among the HBV polymerase domains the spacer domain separating terminal protein from the reverse-transcriptase domain, showed the highest network centrality. Coordinated substitutions preserve the hydrophobicity and charge of Spacer. Maximum likelihood estimates of codon selection showed that Spacer contains the highest number of positively selected sites. Identification of 67% of the domain lacking an ordered structure suggests that Spacer belongs to the class of intrinsically disordered domains and proteins whose crucial functional role in the regulation of transcription, translation and cellular signal transduction has only recently been recognized. Spacer plays a central role in the epistatic network associating substitutions across the HBV genome, including those conferring viral virulence, drug resistance and vaccine escape. The data suggest that Spacer is extensively involved in coordination of HBV evolution.

[1]  Sergei L. Kosakovsky Pond,et al.  Spidermonkey: rapid detection of co-evolving sites using Bayesian graphical models , 2008, Bioinform..

[2]  Sergei L. Kosakovsky Pond,et al.  Datamonkey: rapid detection of selective pressure on individual sites of codon alignments , 2005, Bioinform..

[3]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[4]  J. Villeneuve,et al.  Genotypic succession of mutations of the hepatitis B virus polymerase associated with lamivudine resistance. , 2000, Journal of hepatology.

[5]  Ching-Lung Lai,et al.  Long-term safety of lamivudine treatment in patients with chronic hepatitis B. , 2003, Gastroenterology.

[6]  W. Atchley,et al.  Solving the protein sequence metric problem. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[7]  W. Atchley,et al.  Correlations among amino acid sites in bHLH protein domains: an information theoretic analysis. , 2000, Molecular biology and evolution.

[8]  A. Makhov,et al.  Prediction of terminal protein and ribonuclease H domains in the gene P product of hepadnaviruses , 1989, FEBS letters.

[9]  J. Felsenstein.  Phylogenies and the Comparative Method , 1985, The American Naturalist.

[10]  Christopher Bystroff,et al.  Fully automated ab initio protein structure prediction using I-STES, HMMSTR and ROSETTA , 2002, ISMB.

[11]  Y. Khudyakov.  Coevolution and HBV drug resistance. , 2010, Antiviral therapy.

[12]  N. Hayashi,et al.  Supportive role played by precore and preS2 genomic changes in the establishment of lamivudine-resistant hepatitis B virus. , 2008, Journal of Infectious Diseases.

[13]  P. Tompa,et al.  The pairwise energy content estimated from amino acid composition discriminates between folded and intrinsically unstructured proteins. , 2005, Journal of molecular biology.

[14]  V. Uversky.  Intrinsically Disordered Proteins , 2014 .

[15]  A. Horovitz,et al.  Detection and reduction of evolutionary noise in correlated mutation analysis. , 2005, Protein engineering, design & selection : PEDS.

[16]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[17]  G. Chelvanayagam,et al.  An analysis of simultaneous variation in protein structures. , 1997, Protein engineering.

[18]  D. Richman,et al.  Nomenclature for antiviral‐resistant human hepatitis B virus mutations in the polymerase region , 2001, Hepatology.

[19]  Vladimir Batagelj,et al.  Pajek - Analysis and Visualization of Large Networks , 2001, Graph Drawing Software.

[20]  J. Lara,et al.  Coordinated evolution of the hepatitis C virus , 2008, Proceedings of the National Academy of Sciences.

[21]  Sergei L. Kosakovsky Pond,et al.  Not so different after all: a comparison of methods for detecting amino acid sites under selection. , 2005, Molecular biology and evolution.

[22]  Y. Khudyakov,et al.  Evaluation of Intra-Host Variants of the Entire Hepatitis B Virus Genome , 2011, PloS one.

[23]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[24]  Nikolay A. Kolchanov,et al.  Detection of conserved physico-chemical characteristics of proteins by analyzing clusters of positions with co-ordinated substitutions , 2001, Bioinform..

[25]  R. Aurora,et al.  Genome-wide hepatitis C virus amino acid covariance networks can predict response to antiviral therapy in humans. , 2008, The Journal of clinical investigation.

[26]  S F Altschul,et al.  Weights for data related by a tree. , 1989, Journal of molecular biology.

[27]  Sergei L. Kosakovsky Pond,et al.  An Evolutionary-Network Model Reveals Stratified Interactions in the V3 Loop of the HIV-1 Envelope , 2007, PLoS Comput. Biol..

[28]  D. Holtzman,et al.  Genotypic distribution of hepatitis B virus (HBV) among acute cases of HBV infection, selected United States counties, 1999-2005. , 2011, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[29]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[30]  T. Harrison,et al.  Hepatitis B Virus: Molecular Virology and Common Mutants , 2006, Seminars in liver disease.

[31]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[32]  De-min Yu,et al.  Early changes of hepatitis B virus quasispecies during lamivudine treatment and the correlation with antiviral efficacy. , 2009, Journal of hepatology.

[33]  G. Radziwill,et al.  Mutational analysis of the hepatitis B virus P gene product: domain structure and RNase H activity , 1990, Journal of virology.

[34]  K. Dill,et al.  An iterative method for extracting energy-like quantities from protein structures. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[35]  D. Baker,et al.  Prediction of local structure in proteins using a library of sequence-structure motifs. , 1998, Journal of molecular biology.

[36]  Nikolay A. Kolchanov,et al.  CRASP: a program for analysis of coordinated substitutions in multiple alignments of protein sequences , 2004, Nucleic Acids Res..

[37]  Hiroyuki Ogata,et al.  AAindex: Amino Acid Index Database , 1999, Nucleic Acids Res..

[38]  F. Zoulim,et al.  Hepatitis B virus resistance to nucleos(t)ide analogues. , 2009, Gastroenterology.

[39]  Ziheng Yang,et al.  Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes. , 2002, Molecular biology and evolution.