Structural implications of BK polyomavirus sequence variations in the major viral capsid protein Vp1 and large T-antigen: a computational study

• Published in: mSphere 2024-04, Vol.9 (4), p.e0079923-e0079923
• Main Authors: Durairaj, Janani, Follonier, Océane M, Leuzinger, Karoline, Alexander, Leila T, Wilhelm, Maud, Pereira, Joana, Hillenbrand, Caroline A, Weissbach, Fabian H, Schwede, Torsten, Hirsch, Hans H
• Format: Article
• Language: English
• Subjects: Amino acids; Antigen T (large); Antigens; Antiviral agents; BK virus; Capsid protein; Computational Biology; Computer applications; Cytotoxicity; Deoxyribonucleic acid; DNA; DNA biosynthesis; DNA viruses; Epitopes; Genetic variability; Genomes; genotype; Genotype & phenotype; Genotypes; Hemorrhagic cystitis; Histocompatibility antigen HLA; Immune evasion; Lymphocytes T; Monoclonal antibodies; mutant; Mutation; Nephropathy; polyomavirus; Proteins; Research Article; serotype; Taxonomy; Transplants & implants; Urothelial cancer; Vaccines; variant; VP1 protein
• ISSN: 2379-5042; 2379-5042
• DOI: 10.1128/msphere.00799-23

BK polyomavirus (BKPyV) is a double-stranded DNA virus causing nephropathy, hemorrhagic cystitis, and urothelial cancer in transplant patients. The BKPyV-encoded capsid protein Vp1 and large T-antigen (LTag) are key targets of neutralizing antibodies and cytotoxic T-cells, respectively. Our single-center data suggested that variability in Vp1 and LTag may contribute to failing BKPyV-specific immune control and impact vaccine design. We, therefore, analyzed all available entries in GenBank (1516 ; 742 ) and explored potential structural effects using computational approaches. BKPyV-genotype (gt)1 was found in 71.18% of entries, followed by BKPyV-gt4 (19.26%), BKPyV-gt2 (8.11%), and BKPyV-gt3 (1.45%), but rates differed according to country and specimen type. Vp1-mutations matched a serotype different than the assigned one or were serotype-independent in 43%, 18% affected more than one amino acid. Notable Vp1-mutations altered antibody-binding domains, interactions with sialic acid receptors, or were predicted to change conformation. LTag-sequences were more conserved, with only 16 mutations detectable in more than one entry and without significant effects on LTag-structure or interaction domains. However, LTag changes were predicted to affect HLA-class I presentation of immunodominant 9mers to cytotoxic T-cells. These global data strengthen single center observations and specifically our earlier findings revealing mutant 9mer epitopes conferring immune escape from HLA-I cytotoxic T cells. We conclude that variability of BKPyV-Vp1 and LTag may have important implications for diagnostic assays assessing BKPyV-specific immune control and for vaccine design. Type and rate of amino acid variations in BKPyV may provide important insights into BKPyV diversity in human populations and an important step toward defining determinants of BKPyV-specific immunity needed to protect vulnerable patients from BKPyV diseases. Our analysis of BKPyV sequences obtained from human specimens reveals an unexpectedly high genetic variability for this double-stranded DNA virus that strongly relies on host cell DNA replication machinery with its proof reading and error correction mechanisms. BKPyV variability and immune escape should be taken into account when designing further approaches to antivirals, monoclonal antibodies, and vaccines for patients at risk of BKPyV diseases.