Precise genome editing underlines the distinct contributions of mutations in ERG11, ERG3, MRR1, and TAC1 genes to antifungal resistance in Candida parapsilosis

Candida parapsilosis has recently emerged as a major threat due to the worldwide emergence of fluconazole-resistant strains causing clonal outbreaks in hospitals and poses a therapeutic challenge due to the limited antifungal armamentarium. Here, we used precise genome editing using CRISPR-Cas9 to gain further insights into the contribution of mutations in ERG11, ERG3, MRR1, and TAC1 genes and the influence of allelic dosage to antifungal resistance in C. parapsilosis. Seven of the most common amino acid substitutions previously reported in fluconazole-resistant clinical isolates (including Y132F in ERG11) were engineered in two fluconazole-susceptible C. parapsilosis lineages (ATCC 22019 and STZ5). Each mutant was then challenged in vitro against a large array of antifungals, with a focus on azoles. Any possible change in virulence was also assessed in a Galleria mellonella model. We successfully generated a total of 19 different mutants, using CRISPR-Cas9. Except for R398I (ERG11), all remaining amino acid substitutions conferred reduced susceptibility to fluconazole. However, the impact on fluconazole in vitro susceptibility varied greatly according to the engineered mutation, the stronger impact being noted for G583R acting as a gain-of-function mutation in MRR1. Cross-resistance with newer azoles, non-medical azoles, but also non-azole antifungals such as flucytosine, was occasionally noted. Posaconazole and isavuconazole remained the most active in vitro. Except for G583R, no fitness cost was associated with the acquisition of fluconazole resistance. We highlight the distinct contributions of amino acid substitutions in ERG11, ERG3, MRR1, and TAC1 genes to antifungal resistance in C. parapsilosis.

Elucidating the Impact of Deleterious Mutations on IGHG1 and Their Association with Huntington's Disease

Huntington's disease (HD) is a chronic, inherited neurodegenerative condition marked by chorea, dementia, and changes in personality. The primary cause of HD is a mutation characterized by the expansion of a triplet repeat (CAG) within the huntingtin gene located on chromosome 4. Despite substantial progress in elucidating the molecular and cellular mechanisms of HD, an effective treatment for this disorder is not available so far. In recent years, researchers have been interested in studying cerebrospinal fluid (CSF) as a source of biomarkers that could aid in the diagnosis and therapeutic development of this disorder. Immunoglobulin heavy constant gamma 1 (IGHG1) is one of the CSF proteins found to increase significantly in HD. Considering this, it is reasonable to study the potential involvement of deleterious mutations in IGHG1 in the pathogenesis of this disorder. In this study, we explored the potential impact of deleterious mutations on IGHG1 and their subsequent association with HD. We evaluated 126 single-point amino acid substitutions for their impact on the structure and functionality of the IGHG1 protein while exploiting multiple computational resources such as SIFT, PolyPhen-2, FATHMM, SNPs&Go mCSM, DynaMut2, MAESTROweb, PremPS, MutPred2, and PhD-SNP. The sequence- and structure-based tools highlighted 10 amino acid substitutions that were deleterious and destabilizing. Subsequently, out of these 10 mutations, eight variants (Y32C, Y32D, P34S, V39E, C83R, C83Y, V85M, and H87Q) were identified as pathogenic by disease phenotype predictors. Finally, two pathogenic variants (Y32C and P34S) were found to reduce the solubility of the protein, suggesting their propensity to form protein aggregates. These variants also exhibited higher residual frustration within the protein structure. Considering these findings, the study hypothesized that the identified variants of IGHG1 may compromise its function and potentially contribute to HD pathogenesis.

Predicting the Impact of OTOF Gene Missense Variants on Auditory Neuropathy Spectrum Disorder

Auditory neuropathy spectrum disorder (ANSD) associated with mutations of the OTOF gene is one of the common types of sensorineural hearing loss of a hereditary nature. Due to its high genetic heterogeneity, ANSD is considered one of the most difficult hearing disorders to diagnose. The dataset from 270 known annotated single amino acid substitutions (SAV) related to ANSD was created. It was used to estimate the accuracy of pathogenicity prediction using the known (from dbNSFP4.4) method and a new one. The new method (ConStruct) for the creation of the protein-centric classification model is based on the use of Random Forest for the analysis of missense variants in exons of the OTOF gene. A system of predictor variables was developed based on the modern understanding of the structure and function of the otoferlin protein and reflecting the location of changes in the tertiary structure of the protein due to mutations in the OTOF gene. The conservation values of nucleotide substitutions in genomes of 100 vertebrates and 30 primates were also used as variables. The average prediction of balanced accuracy and the AUC value calculated by the 5-fold cross-validation procedure were 0.866 and 0.903, respectively. The model shows good results for interpreting data from the targeted sequencing of the OTOF gene and can be implemented as an auxiliary tool for the diagnosis of ANSD in the early stages of ontogenesis. The created model, together with the results of the pathogenicity prediction of SAVs via other known accurate methods, were used for the evaluation of a manually created set of 1302 VUS related to ANSD. Based on the analysis of predicted results, 16 SAVs were selected as the new most probable pathogenic variants.

Mutation Analysis of SARS-CoV-2 Variants Isolated from Symptomatic Cases from Andhra Pradesh, India

There has been a continuous evolution in the SARS-CoV-2 genome; therefore, it is necessary to monitor the shifts in the SARS-CoV-2 variants. This study aimed to detect various SARS-CoV-2 variants circulating in the state of Andhra Pradesh, India. The study attempted to sequence the complete S-gene of SARS-CoV-2 of 104 clinical samples using Sanger's method to analyze and compare the mutations with the global prevalence. The method standardized in this study was able to amplify the complete length of the S-gene (3822 bp). The resulting nucleotide and amino acid mutations were analyzed and compared with the local and global SARS-CoV-2 databases using Nextclade and GISAID tools. The Delta variant was the most common variant reported in the present study, followed by the Omicron variant. A variant name was not assigned to thirteen samples using the Nextclade tool. There were sixty-nine types of amino acid substitutions reported (excluding private mutations) throughout the spike gene. The T95I mutation was observed predominantly in Delta variants (15/38), followed by Kappa (3/8) and Omicron (1/31). Nearly all Alpha and Omicron lineages had the N501Y substitution; Q493R was observed only in the Omicron lineage; and other mutations (L445, F486, and S494) were not observed in the present study. Most of these mutations found in the Omicron variant are located near the furin cleavage site, which may play a role in the virulence, pathogenicity, and transmission of the virus. Phylogenetic analysis showed that the 104 complete CDS of SARS-CoV-2 belonged to different phylogenetic clades like 20A, 20B, 20I (Alpha), 21A (Delta), 21B (Kappa), 21I (Delta), 21J (Delta), and 21L (Omicron).

Amino Acid Substitutions at P1 Position Change the Inhibitory Activity and Specificity of Protease Inhibitors BmSPI38 and BmSPI39 from Bombyx mori

It was found that silkworm serine protease inhibitors BmSPI38 and BmSPI39 were very different from typical TIL-type protease inhibitors in sequence, structure, and activity. BmSPI38 and BmSPI39 with unique structure and activity may be good models for studying the relationship between the structure and function of small-molecule TIL-type protease inhibitors. In this study, site-directed saturation mutagenesis at the P1 position was conducted to investigate the effect of P1 sites on the inhibitory activity and specificity of BmSPI38 and BmSPI39. In-gel activity staining and protease inhibition experiments confirmed that BmSPI38 and BmSPI39 could strongly inhibit elastase activity. Almost all mutant proteins of BmSPI38 and BmSPI39 retained the inhibitory activities against subtilisin and elastase, but the replacement of P1 residues greatly affected their intrinsic inhibitory activities. Overall, the substitution of Gly54 in BmSPI38 and Ala56 in BmSPI39 with Gln, Ser, or Thr was able to significantly enhance their inhibitory activities against subtilisin and elastase. However, replacing P1 residues in BmSPI38 and BmSPI39 with Ile, Trp, Pro, or Val could seriously weaken their inhibitory activity against subtilisin and elastase. The replacement of P1 residues with Arg or Lys not only reduced the intrinsic activities of BmSPI38 and BmSPI39, but also resulted in the acquisition of stronger trypsin inhibitory activities and weaker chymotrypsin inhibitory activities. The activity staining results showed that BmSPI38(G54K), BmSPI39(A56R), and BmSPI39(A56K) had extremely high acid-base and thermal stability. In conclusion, this study not only confirmed that BmSPI38 and BmSPI39 had strong elastase inhibitory activity, but also confirmed that P1 residue replacement could change their activity and inhibitory specificity. This not only provides a new perspective and idea for the exploitation and utilization of BmSPI38 and BmSPI39 in biomedicine and pest control, but also provides a basis or reference for the activity and specificity modification of TIL-type protease inhibitors.

[Variability of genes encoding nonstructural proteins of rotavirus А (Reoviridae: Rotavirus: Rotavirus A) genotype G9P[8] during the period of dominance in the territory of Nizhny Novgorod (central part of Russia) (2011-2020)]

INTRODUCTION

In Russia, rotavirus A is the main cause of severe viral gastroenteritis in young children. The molecular features that allow a rotavirus of a particular genotype to gain an evolutionary advantage remain unclear, therefore, the study of the genetic diversity of rotaviruses based on genes encoding nonstructural proteins (NSPs) responsible for the reproduction of the virus in the cell is an urgent task.

OBJECTIVE

To study the genetic diversity of rotaviruses of genotype G9P[8], which dominated Nizhny Novgorod in 20112020, based on genes encoding nonstructural proteins.

MATERIALS AND METHODS

Rotavirus-positive samples were subjected to PCR-genotyping and sequencing of NSP1 NSP5 genes. Phylogenetic analysis was carried out in the MEGA X program.

RESULTS

In the period 20112020, G9P[8] rotaviruses with four variants of the NSP2 gene were co-circulating in Nizhny Novgorod. New alleles were noted in 2012 (N1-a-III), 2016 (N1-a-IV) and in 2019 (N1-a-II). The appearance of new variants of other genes occurred in 2014 (E1-3, NSP4), 2018 (T1-a3-III, NSP3) and in 2019 (A1-b-II, NSP1). NSP2 gene had the most variable amino acid sequence (16 substitutions), 2 to 7 substitutions were observed in NSP1, NSP3 and NSP4, NSP5 was conservative.

DISCUSSION

The results obtained are consistent with the literature data and indicate the participation of NSP genes in maintaining the heterogeneity of the rotavirus population.

CONCLUSION

Until 2018, the genetic diversity of rotaviruses in Nizhny Novgorod was determined by the circulation of strains carrying several alleles of the NSP2 gene and conservative genes NSP1, NSP3NSP5. By the end of the study period, new variants of the genotype G9P[8] were formed in the population, carrying previously unknown combinations of alleles of nonstructural genes.

Peptides That Block RAS-p21 Protein-Induced Cell Transformation

This is a review of approaches to the design of peptides and small molecules that selectively block the oncogenic RAS-p21 protein in ras-induced cancers. Single amino acid substitutions in this protein, at critical positions such as at Gly 12 and Gln 61, cause the protein to become oncogenic. These mutant proteins cause over 90 percent of pancreatic cancers, 40-50 percent of colon cancers and about one third of non-small cell cancers of the lung (NSCCL). RAS-p21 is a G-protein that becomes activated when it exchanges GDP for GTP. Several promising approaches have been developed that target mutant (oncogenic) RAS-p21 proteins in these different cancers. These approaches comprise: molecular simulations of mutant and wild-type proteins to identify effector domains, for which peptides can be made that selectively inhibit the oncogenic protein that include PNC-1 (ras residues 115-126), PNC-2 (ras residues 96-110) and PNC7 (ras residues 35-47); the use of contiguous RAS-p21 peptide sequences that can block ras signaling; cyclic peptides from large peptide libraries and small molecule libraries that can be identified in high throughput assays that can selectively stabilize inactive forms of RAS-p21; informatic approaches to discover peptides and small molecules that dock to specific domains of RAS-p21 that can block mitogenic signal transduction by oncogenic RAS-p21; and the use of cell-penetrating peptides (CPPs) that are attached to the variable domains of the anti-RAS-p21 inactivating monoclonal antibody, Y13 259, that selectively enters oncogenic RAS-p21-containing cancer cells, causing these cells to undergo apoptosis. Several new anti-oncogenic RAS-p21 agents, i.e., Amgen's AMG510 and Mirati Therapeutics' MRTX849, polycyclic aromatic compounds, have recently been FDA-approved and are already being used clinically to treat RAS-p21-induced NSCCL and colorectal carcinomas. These new drugs target the inactive form of RAS-p21 bound to GDP with G12C substitution at the critical Gly 12 residue by binding to a groove bordered by specific domains in this mutant protein into which these compounds insert, resulting in the stabilization of the inactive GDP-bound form of RAS-p21. Other peptides and small molecules have been discovered that block the G12D-RAS-p21 oncogenic protein. These agents can treat specific mutant protein-induced cancers and are excellent examples of personalized medicine. However, many oncogenic RAS-p21-induced tumors are caused by other mutations at positions 12, 13 and 61, requiring other, more general anti-oncogenic agents that are being provided using alternate methods.

Functional Expression of Recombinant Candida auris Proteins in Saccharomyces cerevisiae Enables Azole Susceptibility Evaluation and Drug Discovery

Candida auris infections are difficult to treat due to acquired drug resistance against one or multiple antifungal drug classes. The most prominent resistance mechanisms in C. auris are overexpression and point mutations in Erg11, and the overexpression of efflux pump genes CDR1 and MDR1. We report the establishment of a novel platform for molecular analysis and drug screening based on acquired azole-resistance mechanisms found in C. auris. Constitutive functional overexpression of wild-type C. auris Erg11, Erg11 with amino acid substitutions Y132F or K143R and the recombinant efflux pumps Cdr1 and Mdr1 has been achieved in Saccharomyces cerevisiae. Phenotypes were evaluated for standard azoles and the tetrazole VT-1161. Overexpression of CauErg11 Y132F, CauErg11 K143R, and CauMdr1 conferred resistance exclusively to the short-tailed azoles Fluconazole and Voriconazole. Strains overexpressing the Cdr1 protein were pan-azole resistant. While CauErg11 Y132F increased VT-1161 resistance, K143R had no impact. Type II binding spectra showed tight azole binding to the affinity-purified recombinant CauErg11 protein. The Nile Red assay confirmed the efflux functions of CauMdr1 and CauCdr1, which were specifically inhibited by MCC1189 and Beauvericin, respectively. CauCdr1 exhibited ATPase activity that was inhibited by Oligomycin. The S. cerevisiae overexpression platform enables evaluation of the interaction of existing and novel azole drugs with their primary target CauErg11 and their susceptibility to drug efflux.

The influence of viral protein R amino acid substitutions on clinical outcomes in people living with HIV: A systematic review

BACKGROUND

The HIV viral protein R (Vpr) is a multifunction protein involved in the pathophysiology of HIV-1. Recent evidence has suggested that Vpr amino acid substitutions influence the pathophysiology of HIV-1 and clinical outcomes in people living with HIV (PLWH). Several studies have linked Vpr amino acid substitutions to clinical outcomes in PLWH; however, there is no clear consensus as to which amino acids or amino acid substitutions are most important in the pathophysiology and clinical outcomes in PLWH. We, therefore, conducted a systematic review of studies investigating Vpr amino acid substitutions and clinical outcomes in PLWH.

METHODS

PubMed, Scopus and Web of Science databases were searched according to PRISMA guidelines using a search protocol designed specifically for this study.

RESULTS

A total of 22 studies were included for data extraction, comprising 14 cross-sectional and 8 longitudinal studies. Results indicated that Vpr amino acid substitutions were associated with specific clinical outcomes, including disease progressions, neurological outcomes and treatment status. Studies consistently showed that the Vpr substitution 63T was associated with slower disease progression, whereas 77H and 85P were associated with no significant contribution to disease progression.

CONCLUSIONS

Vpr-specific amino acid substitutions may be contributors to clinical outcomes in PLWH, and future studies should consider investigating the Vpr amino acid substitutions highlighted in this review.

Mouse adaptation of H6 avian influenza viruses and their molecular characteristics

H6 avian influenza viruses (AIVs) not only continue to circulate in both domestic poultry and wild waterfowl, but also have occasionally caused spillovers infections in pigs and humans, posing a potential threat to public health. However, the molecular mechanism of H6 AIV adaptation to mammals remains largely unknown. In this study, two mouse-adapted (MA) H6 AIV strains, named as MA E-Teal/417 and MA GWF-Goose/740, were generated through blind passages in BALB/c mice. The two MA H6 strains replicated more efficiently and showed higher virulence than the corresponding wild type (WT) H6 strains in mice. Genome sequencing revealed that MA E-Teal/417 and MA GWF-Goose/740 carried six amino acid mutations (PB2-T224A/E627K, HA-G124R, NA-F167L/Y356H and M1-M92R), and four amino acid mutations (PB1-K577E, PA-T97I/D514E and HA-T276K), respectively, when compared to the corresponding WT virus. Receptor binding assay showed MA E-Teal/417 had stronger binding activity to α-2,3 SA than WT E-Teal/417. Moreover, the polymerase activity analysis found the RNP polymerase activity of both MA H6 viruses was significantly higher than that of the corresponding WT virus in 293T cells. All these demonstrate that H6 AIV can acquire limit amino acid substitutions to adapt to mammals and increase virulence, highlighting the significance of monitoring such mutations of H6 AIV in the field for alarming the potential of its cross-transmission and pathogenesis in mammals.

Comparative Analysis of SARS-CoV-2 Variants of Concern, Including Omicron, Highlights Their Common and Distinctive Amino Acid Substitution Patterns, Especially at the Spike ORF

In order to gain a deeper understanding of the recently emerged and highly divergent Omicron variant of concern (VoC), a study of amino acid substitution (AAS) patterns was performed and compared with those of the other four successful variants of concern (Alpha, Beta, Gamma, Delta) and one closely related variant of interest (VoI-Lambda). The Spike ORF consistently emerges as an AAS hotspot in all six lineages, but in Omicron this enrichment is significantly higher. The progenitors of each of these VoC/VoI lineages underwent positive selection in the Spike ORF. However, once they were established, their Spike ORFs have been undergoing purifying selection, despite the application of global vaccination schemes from 2021 onwards. Our analyses reject the hypothesis that the heavily mutated receptor binding domain (RBD) of the Omicron Spike was introduced via recombination from another closely related Sarbecovirus. Thus, successive point mutations appear as the most parsimonious scenario. Intriguingly, in each of the six lineages, we observed a significant number of AAS wherein the new residue is not present at any homologous site among the other known Sarbecoviruses. Such AAS should be further investigated as potential adaptations to the human host. By studying the phylogenetic distribution of AAS shared between the six lineages, we observed that the Omicron (BA.1) lineage had the highest number (8/10) of recurrent mutations.

Evolutionary analysis of rotavirus G1P[8] strains from Chennai, South India

Rotaviruses by virtue of its segmented genome generate numerous genotypes. G1P[8] is the most common genotype reported globally. We intend to identify the evolutionary differences among G1P[8] strains from the study with vaccine strains. Stool samples collected from children <5 years were screened for rotavirus antigen by enzyme linked immunosorbent assay. The samples that tested positive for rotavirus were subjected to VP7 and VP4 semi-nested RT-PCR. Sanger sequencing was performed in randomly chosen VP7 and VP4 rotavirus strains. Phylogenetic analysis showed less homology between study strains and vaccine strains and they were placed in different lineages. The VP7 and VP4 proteins of rotavirus were analyzed by two different platforms to identify the amino acid substitutions in the epitope regions. Nine amino acid substitutions with respect to Rotarix®, RotaTeq® and Rotasiil®-V66A, A/T68S, Q72R, N94S, D100E, T113I, S123N, M217T, and I281T were observed in VP7. There were five amino acid substitutions-S145G, N/D195G, N113D, N/I78T, E150D in VP4 (VP8 portion) with respect to Rotarix® and RotaTeq® vaccine strains. M217T substitution in VP7 (epitope 7-2) and N113D, D195G substitution in VP4 (epitope 8-3, 8-1) confer changes in polarity/electrical charge with respect to vaccine strains, thus indicating the need for continued surveillance.

Improving the Activity of Antimicrobial Peptides Against Aquatic Pathogen Bacteria by Amino Acid Substitutions and Changing the Ratio of Hydrophobic Residues

With the increasing number of drug-resistant bacteria, there is an urgent need for new antimicrobial agents, and antimicrobial peptides (AMPs), which exist in the human non-specific immune system, are one of the most promising candidates. It is an effective optimization strategy to modify antimicrobial peptides (AMPs) according to the distribution of amino acids and hydrophobic characteristics. The addition of bacterial pheromones to the N short peptide can increase the ability to recognize bacteria. In this study, we designed and synthesized AMP1–6 by amino acid substitution of mBjAMP1. Additionally, P-6, S-6, and L-6 were designed and synthesized by adding bacterial pheromones based on 1–6. Functional tests showed that the four AMPs had the ability to kill Gram-negative Vibrio anguillarum, Pseudomonas mendocina, and Vibrio parahaemolyticus, and Gram-positive Micrococcus luteus and Listeria monocytogenes. Additionally, all four AMPs induced permeabilization and depolarization of bacterial cell membranes and increased intracellular reactive oxygen species (ROS) levels. Importantly, they had little or no mammalian cytotoxicity. At the same time, 1–6 and L-6 protected the stability of intestinal flora in Sebastes schlegelii and increased the relative abundance of Lactobacillaceae. In summary, our results indicate that the designed AMPs have broad application prospects as a new type of polypeptide antimicrobial agent.

Impact of single amino acid substitution on the structure and function of TANK-binding kinase-1

Tank-binding kinase 1 (TBK1) is a serine/threonine protein kinase involved in various signaling pathways and subsequently regulates cell proliferation, apoptosis, autophagy, antiviral and antitumor immunity. Dysfunction of TBK1 can cause many complex diseases, including autoimmunity, neurodegeneration, and cancer. This dysfunction of TBK1 may result from single amino acid substitutions and subsequent structural alterations. This study analyzed the effect of substituting amino acids on TBK1 structure, function, and subsequent disease using advanced computational methods and various tools. In the initial assessment, a total of 467 mutations were found to be deleterious. After that, in detailed structural and sequential analyses, 13 mutations were found to be pathogenic. Finally, based on the functional importance, two variants (K38D and S172A) of the TBK1 kinase domain were selected and studied in detail by utilizing all-atom molecular dynamics (MD) simulation for 200 ns. MD simulation, including correlation matrix and principal component analysis, helps to get deeper insights into the TBK1 structure at the atomic level. We observed a substantial change in variants' conformation, which may be possible for structural alteration and subsequent TBK1 dysfunction. However, substitution S172A shows a significant conformational change in TBK1 structure as compared to K38D. Thus, this study provides a structural basis to understand the effect of mutations on the kinase domain of TBK1 and its function associated with disease progression.

Molecular Epidemiological Survey of Canine Parvovirus Circulating in China from 2014 to 2019

The global distribution of canine parvovirus (CPV-2) derived from a closely related carnivore parvovirus poses a considerable threat to the dog population. The virus is continuously undergoing genetic evolution, giving rise to several variants. To investigate the prevalence of Chinese CPV-2 strains in recent years, a total of 30 CPV-2 strains were collected from 2018 to 2021 and the VP2 gene was sequenced and analyzed. Two variants, new CPV-2a (297Ala, 426Asn) and CPV-2c (426Glu), were identified. In contrast to previous reports, the CPV-2c variant has gained an epidemiological advantage over the new CPV-2a variant in China. To compensate for the relatively small sample size, 683 Chinese CPV-2 strains identified between 2014 and 2019 were retrieved from the GenBank database and previous publications, and analyses of these strains further supported our findings, which should be considered since the CPV-2c variant has been frequently associated with immune failure in adult dogs. VP2 protein sequence analysis revealed several amino acid substitutions, including Ala5Gly, Pro13Ser, Phe267Tyr, Tyr324Ile, Gln370Arg, Thr440Ala, and Lys570Arg. Phylogenetic analysis of full-length VP2 gene indicated a close relationship between Chinese strains and other Asian strains, suggesting mutual transmission between Asian countries. Furthermore, intercontinental transmission is a cause for concern. Surprisingly, two feline panleukopenia virus (FPV) strains with the Ile101Thr mutation in the VP2 protein were identified in canine fecal samples; FPV has been considered incapable of infecting dogs. This study clarified the epidemic characteristics of Chinese CPV-2 strains detected between 2014 and 2019, offering a reference for epidemic control. In addition, the detection of FPV in canine samples may provide information for future studies on the evolution of carnivore parvoviruses.

Geographical and Tick-Dependent Distribution of Flavi-Like Alongshan and Yanggou Tick Viruses in Russia

The genus Flavivirus includes related, unclassified segmented flavi-like viruses, two segments of which have homology with flavivirus RNA-dependent RNA polymerase NS5 and RNA helicase-protease NS3. This group includes such viruses as Jingmen tick virus, Alongshan virus, Yanggou tick virus and others. We detected the Yanggou tick virus in Dermacentor nuttalli and Dermacentor marginatus ticks in two neighbouring regions of Russia. The virus prevalence ranged from 0.5% to 8.0%. We detected RNA of the Alongshan virus in 44 individuals or pools of various tick species in eight regions of Russia. The virus prevalence ranged from 0.6% to 7.8%. We demonstrated the successful replication of the Yanggou tick virus and Alongshan virus in IRE/CTVM19 and HAE/CTVM8 tick cell lines without a cytopathic effect. According to the phylogenetic analysis, we divided the Alongshan virus into two groups: an Ixodes persulcatus group and an Ixodes ricinus group. In addition, the I. persulcatus group can be divided into European and Asian subgroups. We found amino acid signatures specific to the I. ricinus and I. persulcatus groups and also distinguished between the European and Asian subgroups of the I. persulcatus group.

The Effect of I155T, K156Q, K156E and N186K Mutations in Hemagglutinin on the Virulence and Reproduction of Influenza A/H5N1 Viruses

The continued circulation of influenza A virus subtype H5 may cause the emergence of new potential pandemic virus variants, which can be transmitted from person to person. The occurrence of such variants is mainly related to mutations in hemagglutinin (HA). Previously we discovered mutations in H5N1 influenza virus hemagglutinin, which contributes to virus immune evasion. The purpose of this work was to study the role of these mutations in changing other, non-antigenic properties of the virus and the possibility of their maintenance in the viral population. Mutations were introduced into the HA gene of a recombinant H5N1 influenza A virus (VNH5N1-PR8/CDC-RG) using site-specific mutagenesis. The "variant" viruses were investigated and compared with respect to replication kinetics in chicken embryos, thermostability, reproductive activity at different temperatures (33, 37 and 40°C), and virulence for mice. Amino acid substitutions I155T, K156Q, K156E+V138A, N186K led to a decrease in thermal stability, replication activity of the mutant viruses in chicken embryos, and virulence for mice, although these effects differed between the variants. The K156Q and N186K mutations reduced viral reproduction at elevated temperature (40°C). The analysis of the frequency of these mutations in natural isolates of H5N1 influenza viruses indicated that the K156E/Q and N186K mutations have little chance to gain a foothold during evolution, in contrast to the I155T mutation, which is the most responsible for antigenic drift. The A138V and N186K mutations seem to be adaptive in mammalian viruses.

[The Effect of I155T, K156Q, K156E and N186K Mutations in Hemagglutinin on the Virulence and Reproduction of Influenza A/H5N1 Viruses]

The continued circulation of influenza A virus subtype H5 may cause the emergence of new potential pandemic virus variants, which can be transmitted from person to person. The occurrence of such variants is mainly related to mutations in hemagglutinin (HA). Previously we discovered mutations in H5N1 influenza virus hemagglutinin, which contributes to virus immune evasion. The purpose of this work was to study the role of these mutations in changing other, non-antigenic properties of the virus and the possibility of their maintenance in the viral population. Mutations were introduced into the HA gene of a recombinant H5N1 influenza A virus (VNH5N1-PR8/CDC-RG) using site-specific mutagenesis. The "variant" viruses were investigated and compared with respect to replication kinetics in chicken embryos, thermostability, reproductive activity at different temperatures (33, 37 and 40°C), and virulence for mice. Amino acid substitutions I155T, K156Q, K156E+V138A, N186K led to a decrease in thermal stability, replication activity of the mutant viruses in chicken embryos, and virulence for mice, although these effects differed between the variants. The K156Q and N186K mutations reduced viral reproduction at elevated temperature (40°C). The analysis of the frequency of these mutations in natural isolates of H5N1 influenza viruses indicated that the K156E/Q and N186K mutations have little chance to gain a foothold during evolution, in contrast to the I155T mutation, which is the most responsible for antigenic drift. The A138V and N186K mutations seem to be adaptive in mammalian viruses.

Experimental Infection Using Mouse-Adapted Influenza B Virus in a Mouse Model

Every year, influenza B viruses (IBVs) contribute to annual illness, and infection can lead to serious respiratory disease among humans. More attention is needed in several areas, such as increasing virulence or pathogenicity of circulating B viruses and developing vaccines against current influenza. Since preclinical trials of anti-influenza drugs are mainly conducted in mice, we developed an appropriate infection model, using an antigenically-relevant IBV strain, for furtherance of anti-influenza drug testing and influenza vaccine protective efficacy analysis. A Victoria lineage (clade 1A) IBV was serially passaged 17 times in BALB/c mice, and adaptive amino acid substitutions were found in hemagglutinin (HA) (T214I) and neuraminidase (NA) (D432N). By electron microscopy, spherical and elliptical IBV forms were noted. Light microscopy showed that mouse-adapted IBVs caused influenza pneumonia on day 6 post inoculation. We evaluated the illness pathogenicity, viral load, and histopathological features of mouse-adapted IBVs and estimated anti-influenza drugs and vaccine efficiency in vitro and in vivo. Assessment of an investigational anti-influenza drug (oseltamivir ethoxysuccinate) and an influenza vaccine (Ultrix^®, SPBNIIVS, Saint Petersburg, Russia) showed effectiveness against the mouse-adapted influenza B virus.

Elucidating Binding Sites and Affinities of ERα Agonists and Antagonists to Human Alpha-Fetoprotein by In Silico Modeling and Point Mutagenesis

Alpha-fetoprotein (AFP) is a major embryo- and tumor-associated protein capable of binding and transporting a variety of hydrophobic ligands, including estrogens. AFP has been shown to inhibit estrogen receptor (ER)-positive tumor growth, which can be attributed to its estrogen-binding ability. Despite AFP having long been investigated, its three-dimensional (3D) structure has not been experimentally resolved and molecular mechanisms underlying AFP–ligand interaction remains obscure. In our study, we constructed a homology-based 3D model of human AFP (HAFP) with the purpose of molecular docking of ERα ligands, three agonists (17β-estradiol, estrone and diethylstilbestrol), and three antagonists (tamoxifen, afimoxifene and endoxifen) into the obtained structure. Based on the ligand-docked scoring functions, we identified three putative estrogen- and antiestrogen-binding sites with different ligand binding affinities. Two high-affinity binding sites were located (i) in a tunnel formed within HAFP subdomains IB and IIA and (ii) on the opposite side of the molecule in a groove originating from a cavity formed between domains I and III, while (iii) the third low-affinity binding site was found at the bottom of the cavity. Here, 100 ns molecular dynamics (MD) simulation allowed us to study their geometries and showed that HAFP–estrogen interactions were caused by van der Waals forces, while both hydrophobic and electrostatic interactions were almost equally involved in HAFP–antiestrogen binding. Molecular mechanics/Generalized Born surface area (MM/GBSA) rescoring method exploited for estimation of binding free energies (ΔG_bind) showed that antiestrogens have higher affinities to HAFP as compared to estrogens. We performed in silico point substitutions of amino acid residues to confirm their roles in HAFP–ligand interactions and showed that Thr132, Leu138, His170, Phe172, Ser217, Gln221, His266, His316, Lys453, and Asp478 residues, along with two disulfide bonds (Cys224–Cys270 and Cys269–Cys277), have key roles in both HAFP–estrogen and HAFP–antiestrogen binding. Data obtained in our study contribute to understanding mechanisms underlying protein–ligand interactions and anticancer therapy strategies based on ERα-binding ligands.

[Change of phenotypic properties of escape mutants and readaptants of influenza virus A (H1N1)pdm09 under the influence of selected mutations in the molecule of hemagglutinin.]

INTRODUCTION

After the emergence and spread of pandemic H1N1 viruses in 2009, antigenic epitopes recognized by neutralizing antibodies against the hemagglutinin of influenza A/Moscow/01/09(H1N1)pdm09 viruses were studied.

PURPOSE

The purpose of the study was to obtain readapted variants of the virus from a low-virulent escapemutant that has an increased affinity of the avian and the human types cellular receptors compared to the wild type and the comparative study of their antigenic and receptor specificity.

MATERIAL AND METHODS

Viruses were accumulated in 10-day-old chicken embryos. The MAB panel against HA of influenza virus strain A/IIV-Moscow/01/09(H1N1)sw1 was used in the form of ascites fluids from mice. Immunization of mice, HI testing, elution of viruses from chicken erythrocytes, PCR and sequencing of readapted variants were performed by standard methods.

RESULTS

The amino acid substitution A198E acquired in the process of readaptation leads to changes in the antigenic specificity. A correlation was found between a decrease in virulence of a low-virulent escape mutant associated with the substitution D190N in the hemagglutinin molecule and an increase in the hemagglutinating titer to inhibitors in normal mouse serum. Viruses with low affinity of cellular receptor analogs and carrying amino acid substitutions have an increased ability to elute from chicken erythrocytes.

DISCUSSION

The results discuss the effect of mutations in the HA molecule of the influenza A(H1N1) pdm09 virus to the change in antigen specificity; virulence for mice, adsorption-elution at cellular receptors.

CONCLUSION

A comparative study of the antigenic specificity and receptor-binding activity of the escape mutants was conducted for the hemagglutinin of the influenza virus A/Moscow/01/2009 (H1N1)swl, and the readapted variants obtained for one of the escape mutants with reduced virulence for mouse. Monitoring the pleiotropic effect of mutations in the hemagglutinin H1 molecule is necessary to predict variants of the virus with pandemic potential.

Amino Acid Substitutions and Differential Gene Expression of Outer Membrane Proteins in Adherent-Invasive Escherichia coli

Variations in the sequence and/or the expression of outer membrane proteins (OMPs) may modulate bacterial virulence. OmpA and OmpC have been involved in the interaction of adherent-invasive Escherichia coli (AIEC) strain LF82 with intestinal epithelial cells (IECs). Scarce data exist about OMPs sequence variants in a collection of AIEC strains, and no study of OMPs expression during infection exists. We aimed to determine whether particular mutations or differential expression of OMPs are associated with AIEC virulence. The ompA, ompC, and ompF genes in 14 AIEC and 30 non-AIEC strains were sequenced by Sanger method, and the protein expression profile was analyzed by urea-SDS-PAGE. Gene expression was determined during in vitro bacterial infection of intestine-407 cells by RT-qPCR. The distribution of amino acid substitutions in OmpA-A200V, OmpC-S89N, V220I, and W231D associated with pathotype and specific changes (OmpA-A200V, OmpC-V220I, D232A, OmpF-E51V, and M60K) correlated with adhesion and/or invasion indices but no particular variants were found specific of AIEC. OMPs protein levels did not differ according to pathotype when growing in Mueller-Hinton broth. Interestingly, higher OMPs gene expression levels were reported in non-AIEC growing in association with cells compared with those non-AIEC strains growing in the supernatants of infected cultures (p < 0.028), whereas in AIEC strains ompA expression was the only increased when growing in association with cells (p = 0.032), but they did not significantly alter ompC and ompF expression under this condition (p > 0.146). Despite no particular OMPs sequence variants have been found as a common and distinctive trait in AIEC, some mutations could facilitate a better interaction with the host. Moreover, the different behavior between pathotypes regarding OMPs gene expression at different stages of infection could be related with the virulence of the strains.

New Insights into the Cyp51 Contribution to Azole Resistance in Aspergillus Section Nigri

Invasive aspergillosis (IA) is a severe condition mainly caused by Aspergillus fumigatus, although other species of the genus, such as section Nigri members, can also be involved. Voriconazole (VRC) is the recommended treatment for IA; however, the prevalence of azole-resistant Aspergillus isolates has alarmingly increased in recent years, and the underlying resistance mechanisms in non-fumigatus species remain unclear. We have determined the in vitro susceptibility of 36 strains from section Nigri to VRC, posaconazole (POS), and itraconazole (ITC), and we have explored the role of Cyp51A and Cyp51B, both targets of azoles, in azole resistance. The three drugs were highly active; POS displayed the best in vitro activity, while ITC and VRC showed MICs above the established epidemiological cutoff values in 9 and 16% of the strains, respectively. Furthermore, expression studies of cyp51A and cyp51B in control condition and after VRC exposure were performed in 14 strains with different VRC susceptibility. We found higher transcription of cyp51A, which was upregulated upon VRC exposure, but no correlation between MICs and cyp51 transcription levels was observed. In addition, cyp51A sequence analyses revealed nonsynonymous mutations present in both, wild-type and non-wild-type strains of A. niger and A. tubingensis Nevertheless, a few mutations were exclusively present in non-wild-type A. tubingensis strains. Altogether, our results suggest that azole resistance in section Nigri is not clearly explained by Cyp51A protein alteration or by cyp51 gene upregulation, which indicates that other mechanisms might be involved.

[Three Mutations in the Stalk Region of Hemagglutinin Affect the pH of Fusion and Pathogenicity of H5N1 Influenza Virus]

Previously, an attenuated variant Ku/at was obtained from the highly pathogenic avian influenza virus A/chicken/Kurgan/3/2005 (H5N1) by a reverse selection method aimed at increasing the virus resistance to a proteolytic cleavage and acidic pH values. In the Ku/at, 10 mutations in proteins PB2, PB1, HA, NA, and NS1 occurred. In comparison with the parental strain, the pH of the conformational transition of the viral glycoprotein hemagglutinin (HA) and virulence for mice and chickens have decreased in an attenuated variant. The purpose of this work is to clarify the role of three mutations in the stalk region of HA: Asp54Asn in HA1 and Val48Ile and Lys131Thr in HA2 (H3 HA numbering). To attain these ends, analogous substitutions were introduced into HA with a deleted polybasic cleavage site (important for pathogenicity) of the recombinant A/Vietnam/1203/04-PR8/CDC-RG (H5N1) virus, and so we created the VN3x-PR variant. Viruses VN3x-PR and Ku/at with the same three mutations, but different proteolytic cleavage sites in HA, as well as the corresponding initial viruses, were tested for pathogenicity in mice and in the erythrocyte hemolysis test. Compared with the parental strains, the virulence of their mutant variants in the case of intranasal infection of BALB/c mice decreased by 4-5 orders of magnitude, and the pH of the conformational transition of HA decreased from 5.70-5.80 to 5.25-5.30, which is typical for low pathogenic natural isolates. Thus, as a result of the study, the attenuating role of these three mutations in HA has been proved, a correlation was established between the pH value of the HA conformational transition and the virulence of H5N1 influenza viruses, and it was shown that the polybasic cleavage site of the H5 HA does not always determine high pathogenicity of the virus.

Nuclear and mitochondrial RNA editing systems have opposite effects on protein diversity

RNA editing can yield protein products that differ from those directly encoded by genomic DNA. This process is pervasive in the mitochondria of many eukaryotes, where it predominantly results in the restoration of ancestral protein sequences. Nuclear mRNAs in metazoans also undergo editing (adenosine-to-inosine or 'A-to-I' substitutions), and most of these edits appear to be nonadaptive 'misfirings' of adenosine deaminases. However, recent analysis of cephalopod transcriptomes found that many editing sites are shared by anciently divergent lineages within this group, suggesting they play some adaptive role. Recent discoveries have also revealed that some fungi have an independently evolved A-to-I editing mechanism, resulting in extensive recoding of their nuclear mRNAs. Here, phylogenetic comparisons were used to determine whether RNA editing generally restores ancestral protein sequences or creates derived variants. Unlike in mitochondrial systems, RNA editing in metazoan and fungal nuclear transcripts overwhelmingly leads to novel sequences not found in inferred ancestral proteins. Even for the subset of RNA editing sites shared by deeply divergent cephalopod lineages, the primary effect of nuclear editing is an increase-not a decrease-in protein divergence. These findings suggest fundamental differences in the forces responsible for the evolution of RNA editing in nuclear versus mitochondrial systems.

Amino Acid Substitutions Associated with Avian H5N6 Influenza A Virus Adaptation to Mice

At least 15 cases of human beings infected with H5N6 have been reported since 2014, of which at least nine were fatal. The highly pathogenic avian H5N6 influenza virus may pose a serious threat to both public health and the poultry industry. However, the molecular features promoting the adaptation of avian H5N6 influenza viruses to mammalian hosts is not well understood. Here, we sequentially passaged an avian H5N6 influenza A virus (A/Northern Shoveler/Ningxia/488-53/2015) 10 times in mice to identify the adaptive amino acid substitutions that confer enhanced virulence to H5N6 in mammals. The 1st and 10th passages of the mouse-adapted H5N6 viruses were named P1 and P10, respectively. P1 and P10 displayed higher pathogenicity in mice than their parent strain. P10 showed significantly higher replication capability in vivo and could be detected in the brains of mice, whereas P1 displayed higher replication efficiency in their lungs but was not detectable in the brain. Similar to its parent strain, P10 remained no transmissible between guinea pigs. Using genome sequencing and alignment, multiple amino acid substitutions, including PB2 E627K, PB2 T23I, PA T97I, and HA R239H, were found in the adaptation of H5N6 to mice. In summary, we identified amino acid changes that are associated with H5N6 adaptation to mice.

Mutations in Cancer Cause Gain of Cysteine, Histidine, and Tryptophan at the Expense of a Net Loss of Arginine on the Proteome Level

Accumulation of somatic mutations is critical for the transition of a normal cell to become cancerous. Mutations cause amino acid substitutions that change properties of proteins. However, it has not been studied as to what extent the composition and accordingly chemical properties of the cell proteome is altered as a result of the increased mutation load in cancer. Here, we analyzed data on amino acid substitutions caused by mutations in about 2000 protein coding genes from the Cancer Cell Line Encyclopedia that contains information on nucleotide and amino acid alterations in 782 cancer cell lines, and validated the analysis with information on amino acid substitutions for the same set of proteins in the Catalogue of Somatic Mutations in Cancer (COSMIC; v78) in circa 18,000 tumor samples. We found that nonsynonymous single nucleotide substitutions in the analyzed proteome subset ultimately result in a net gain of cysteine, histidine, and tryptophan at the expense of a net loss of arginine. The extraordinary loss of arginine may be attributed to some extent to composition of its codons as well as to the importance of arginine in the functioning of prominent tumor suppressor proteins like p53.

[Changes in the phenotypic properties of highly pathogenic influenza A virus of H5N1 subtype induced by N186I and N186T point mutations in hemagglutinin]

The change in the phenotypic properties resulting from amino acid substitutions in the hemagglutinin (HA) molecule is an important link in the evolutionary process of influenza viruses. It is believed to be one of the mechanisms of the emergence of highly pathogenic strains of influenza A viruses, including subtype H5N1. Using the site-directed mutagenesis, we introduced mutations in the HA gene of the H5N1 subtype of influenza A virus. The obtained virus variants were analyzed and compared using the following parameters: optimal pH of conformational transition (according to the results of the hemolysis test), specificity of receptor binding (using a set of synthetic analogues of cell surface sialooligosaccharides), thermoresistance (heat-dependent reduction of hemagglutinin activity), virulence in mice, and the kinetics of replication in chicken embryos, and reproductive activity at different temperatures (RCT-based). N186I and N186T mutations in the HA protein increased the virulence of the original virus in mice. These mutations accelerated virus replication in the early stages of infection in chicken embryos and increased the level of replication at late stages. In addition, compared to the original virus, the mutant variants replicated more efficiently at lower temperatures. The obtained data clearly prove the effect of amino acid substitutions at the 186 position of HA on phenotypic properties of the H5N1 subtype of influenza A.

Loss and Gain of Human Acidic Mammalian Chitinase Activity by Nonsynonymous SNPs

Acidic mammalian chitinase (AMCase) is implicated in asthma, allergic inflammation, and food processing. Little is known about genetic and evolutional regulation of chitinolytic activity of AMCase. Here, we relate human AMCase polymorphisms to the mouse AMCase, and show that the highly active variants encoded by nonsynonymous single-nucleotide polymorphisms (nsSNPs) are consistent with the mouse AMCase sequence. The chitinolytic activity of the recombinant human AMCase was significantly lower than that of the mouse counterpart. By creating mouse-human chimeric AMCase protein we found that the presence of the N-terminal region of human AMCase containing conserved active site residues reduced the enzymatic activity of the molecule. We were able to significantly increase the activity of human AMCase by amino acid substitutions encoded by nsSNPs (N45, D47, and R61) with those conserved in the mouse homologue (D45, N47, and M61). For abolition of the mouse AMCase activity, introduction of M61R mutation was sufficient. M61 is conserved in most of primates other than human and orangutan as well as in other mammals. Orangutan has I61 substitution, which also markedly reduced the activity of the mouse AMCase, indicating that the M61 is a crucial residue for the chitinolytic activity. Altogether, our data suggest that human AMCase has lost its chitinolytic activity by integration of nsSNPs during evolution and that the enzyme can be reactivated by introducing amino acids conserved in the mouse counterpart.

Virological, epidemiological, clinic, and molecular genetic features of the influenza epidemic in 2015-2016: prevailing of the influenza A(H1N1)09 pdm virus in Russia and countries of the Northern hemisphere

This work describes the specific features of the influenza virus circulating in the period from October 2015 to March 2016 in 10 cities of Russia, the basic laboratories of CEEI at the D.I. Ivanovsky Institute of Virology "Federal Research Centre of Epidemilogy and Microbiology named after the honorary academician N.F. Gamaleya" of the Ministry of Health of the Russian Federation. The increase in the morbidity caused by influenza viruses was detected in January-February 2016. The duration of the morbidity peak was 4-5 weeks. The most vulnerable group included children at the age from 3 to 6; a high rate of hospitalization was also detected among people at the age of 15-64 (65%). In clinic symptoms there were middle and severe forms with high frequency of hospitalization as compared with the season of 2009-2010, but much higher in comparison with the season of 2014-2015. Some of the hospitalized patients had virus pneumonias, half of which were bilateral. Among these patients, 10% were children; 30%, adults. The mortality in the intensive care unit of the hospital was 46%. Almost all lethal cases were among unvaccinated patients in the case of late hospitalization and without early antiviral therapy. The predominance of the influenza A(H1N1)09pdm virus both in the Russian Federation and the major part of the countries in the Northern hemisphere was noted. The results of the study of the antigenic properties of influenza strains of A(H1N1)pdm09 virus did not reveal any differences with respect to the vaccine virus. The sequencing data showed the amino acid substitutions in hemagglutinin (receptor binding and Sa sites) and in genes encoding internal proteins (PA, NP, M1, NS1). Strains were sensitive to oseltamivir and zanamivir and maintained resistance to rimantadine. The participation of non-influenza ARI viruses was comparable to that in preliminary epidemic seasons.

Classification of Amino Acid Substitutions in Mismatch Repair Proteins Using PON-MMR2

Variations in mismatch repair (MMR) system genes are causative of Lynch syndrome and other cancers. Thousands of variants have been identified in MMR genes, but the clinical relevance is known for only a small proportion. Recently, the InSiGHT group classified 2,360 MMR variants into five classes. One-third of variants, majority of which is nonsynonymous variants, remain to be of uncertain clinical relevance. Computational tools can be used to prioritize variants for disease relevance investigations. Previously, we classified 248 MMR variants as likely pathogenic and likely benign using PON-MMR. We have developed a novel tool, PON-MMR2, which is trained on a larger and more reliable dataset. In performance comparison, PON-MMR2 outperforms both generic tolerance prediction methods as well as methods optimized for MMR variants. It achieves accuracy and MCC of 0.89 and 0.78, respectively, in cross-validation and 0.86 and 0.69, respectively, on an independent test dataset. We classified 354 class 3 variants in InSiGHT database as well as all possible amino acid substitutions in four MMR proteins. Likely harmful variants mainly appear in the protein core, whereas likely benign variants are on the surface. PON-MMR2 is a highly reliable tool to prioritize variants for functional analysis. It is freely available at http://structure.bmc.lu.se/PON-MMR2/.

Evolutionary patterns in coiled-coils

Models of protein evolution are used to describe evolutionary processes, for phylogenetic analyses and homology detection. Widely used general models of protein evolution are biased toward globular domains and lack resolution to describe evolutionary processes for other protein types. As three-dimensional structure is a major constraint to protein evolution, specific models have been proposed for other types of proteins. Here, we consider evolutionary patterns in coiled-coil forming proteins. Coiled-coils are widespread structural domains, formed by a repeated motif of seven amino acids (heptad repeat). Coiled-coil forming proteins are frequently rods and spacers, structuring both the intracellular and the extracellular spaces that often form protein interaction interfaces. We tested the hypothesis that due to their specific structure the associated evolutionary constraints differ from those of globular proteins. We showed that substitution patterns in coiled-coil regions are different than those observed in globular regions, beyond the simple heptad repeat. Based on these substitution patterns we developed a coiled-coil specific (CC) model that in the context of phylogenetic reconstruction outperforms general models in tree likelihood, often leading to different topologies. For multidomain proteins containing both a coiled-coil region and a globular domain, we showed that a combination of the CC model and a general one gives higher likelihoods than a single model. Finally, we showed that the model can be used for homology detection to increase search sensitivity for coiled-coil proteins. The CC model, software, and other supplementary materials are available at http://www.evocell.org/cgl/resources (last accessed January 29, 2015).

Sequential analysis of amino acid substitutions with hepatitis B virus in association with nucleoside/nucleotide analog treatment detected by deep sequencing

Taking nucleoside/nucleotide analogs is a major antiviral therapy for chronic hepatitis B infection. The problem with this treatment is the selection for drug-resistant mutants. Currently, identification of genotypic drug resistance is conducted by molecular cloning sequenced by the Sanger method. However, this methodology is complicated and time-consuming. These limitations can be overcome by deep sequencing technology. Therefore, we performed sequential analysis of the frequency of drug resistance in one individual, who was treated with lamivudine on-and-off therapy for 2 years, by deep sequencing. The lamivudine-resistant mutations at rtL180M and rtM204V and the entecavir-resistant mutation at rtT184L were detected in the first subject. The lamivudine- and entecavir-resistant strain was still detected in the last subject. However, in the deep sequencing analysis, rt180 of the first subject showed a mixture in 76.9% of the methionine and in 23.1% of the leucine, and rt204 also showed a mixture in 69.0% of the valine and 29.8% of the isoleucine. During the treatment, the ratio of resistant mutations increased. At rt184, the resistant variants were detectable in 58.7% of the sequence, with the replacement of leucine by the wild-type threonine in the first subject. Gradually, entecavir-resistant variants increased in 82.3% of the leucine in the last subject. In conclusion, we demonstrated the amino acid substitutions of the serial nucleoside/nucleotide analog resistants. We revealed that drug-resistant mutants appear unchanged at first glance, but actually there are low-abundant mutations that may develop drug resistance against nucleoside/nucleotide analogs through the selection of dominant mutations.

Complex dynamics of hepatitis B virus resistance to adefovir

In patients with hepatitis B e antigen-negative chronic hepatitis B, adefovir dipivoxil administration selects variants bearing reverse transcriptase rtN236T and/or rtA181V/T substitutions in 29% of cases after 5 years. The aim of this study was to characterize the dynamics of adefovir-resistant variant populations during adefovir monotherapy in order to better understand the molecular mechanisms underlying hepatitis B virus resistance to this class of nucleotide analogues. Patients included in a 240-week clinical trial of adefovir monotherapy who developed adefovir resistance-associated substitutions were studied. The dynamics of hepatitis B virus populations were analyzed over time, after generating nearly 4,000 full-length reverse transcriptase sequences, and compared with the replication kinetics of the virus during therapy. Whatever the viral kinetics pattern, adefovir resistance was characterized by exclusive detection of a dominant wild-type, adefovir-sensitive variant population at baseline and late and gradual selection by adefovir of several coexisting resistant viral populations, defined by the presence of amino acid substitutions at position rt236, position rt181, or both. The gain in fitness of one or the other of these resistant populations during adefovir administration was never associated with the selection of additional amino acid substitutions in the reverse transcriptase.

CONCLUSION

Our results suggest that adefovir administration selects poorly fit preexisting or emerging viral populations with low-level adefovir resistance, which subsequently compete to fill the replication space. Viral kinetics depends on the initial virological response to adefovir. Lamivudine add-on restores some antiviral efficacy, but adefovir-resistant variants remain predominant. Whether these adefovir resistance-associated substitutions may confer cross-resistance to tenofovir in vivo will need to be determined.

Structural features of peptide analogs of human histocompatibility leukocyte antigen class I epitopes that are more potent and immunogenic than wild-type peptide

Certain peptide analogs that carry substitutions at residues other than the main major histocompatibility complex anchors and are surprisingly much more antigenic than wild-type peptide (heteroclitic analogs). To date, it was unknown how frequently wild-type epitopes could be modified to obtain heteroclitic activity. In this study, we analyzed a large panel of analogs of two different human histocompatibility leukocyte antigen (HLA)-A2.1-restricted epitopes and found that heteroclitic analogs were associated with higher magnitude responses and increased (up to 10(7)-fold) sensitivity to antigen, and corresponded to conservative or semiconservative substitutions at odd-numbered positions in the middle of the peptide (positions 3, 5, or 7). These findings were validated by performing additional immunogenicity studies in murine and human systems with four additional epitopes. The biological relevance of heteroclitic analogs was underlined when predicted analogs of the p53.261 epitope was shown to induce cytotoxic T lymphocytes (CTLs) that recognize low concentrations of peptide (high avidity) in vivo and demonstrate in vitro antitumor recognition. Finally, in vitro immunization of human peripheral blood mononuclear cells with two heteroclitic analogs resulted in recruitment of more numerous CTLs which were associated with increased antigen sensitivity. In conclusion, heteroclitic analogs were identified in each of the six cases studied and structural features were defined which allow identification of such analogs. The strong CTL immunity elicited by heteroclitic epitopes suggest that they could be of significant value in vaccination against tolerant or weakly immunogenic tumor-associated and viral antigens.

Heat capacity changes upon burial of polar and nonpolar groups in proteins

In this paper we address the question of whether the burial of polar and nonpolar groups in the protein locale is indeed accompanied by the heat capacity changes, DeltaC(p), that have an opposite sign, negative for nonpolar groups and positive for polar groups. To accomplish this, we introduced amino acid substitutions at four fully buried positions of the ubiquitin molecule (Val5, Val17, Leu67, and Gln41). We substituted Val at positions 5 and 17 and Leu at position 67 with a polar residue, Asn. As a control, Ala was introduced at the same three positions. We also replaced the buried polar Gln41 with Val and Leu, nonpolar residues that have similar size and shape as Gln. As a control, Asn was introduced at Gln41 as well. The effects of these amino acid substitutions on the stability, and in particular, on the heat capacity change upon unfolding were measured using differential scanning calorimetry. The effect of the amino acid substitutions on the structure was also evaluated by comparing the (1)H-(15)N HSQC spectra of the ubiquitin variants. It was found that the Ala substitutions did not have a considerable effect on the heat capacity change upon unfolding. However, the substitutions of aliphatic side chains (Val or Leu) with a polar residue (Asn) lead to a significant (> 30%) decrease in the heat capacity change upon unfolding. The decrease in heat capacity changes does not appear to be the result of significant structural perturbations as seen from the HSQC spectra of the variants. The substitution of a buried polar residue (Gln41) to a nonpolar residue (Leu or Val) leads to a significant (> 25%) increase in heat capacity change upon unfolding. These results indicate that indeed the heat capacity change of burial of polar and nonpolar groups has an opposite sign. However, the observed changes in DeltaC(p) are several times larger than those predicted, based on the changes in water accessible surface area upon substitution.