Frequency of Fusion Inhibitor Resistance Mutations Among Therapy-Naïve HIV Patients

Glycoprotein 41 (gp41) of the human immunodeficiency virus type 1 (HIV-1) protein plays a critical role in membrane fusion. Gp41 binds to proteins in the plasma membrane of CD4+ T cells, particularly the T-cell antigen receptor (TCR). These findings indicate that gp41 is involved in the assembly of HIV-1 at the plasma membrane of T cells and affects the stimulation of the TCR. To control HIV-1, new inhibitors were introduced to target the gp41 protein. However, mutations in this region might reduce their efficacy. The Gp41 region was amplified from the sera of 30 patients using nested polymerase chain reaction. The sequences were analyzed by bioinformatics tools to identify mutations and gp41 structural features. Subtyping and the interaction between fusion inhibitors and gp41 proteins were also examined. As the first report from Iran, docking analysis between fusion inhibitors and Iranian gp41 proteins showed that mutations in gp41 could not reduce the efficacy of the fusion inhibitors. Most of the patients were infected with CRF35-AD. Several post-modification positions, including glycosylation and phosphorylation sites, were identified in the gp41 protein. Our findings revealed no known multinational drug resistance to gp41 inhibitors; thus, fusion inhibitors can effectively inhibit HIV in Iranian patients. In addition, the present study introduced a new gp41 region (36-44 aa), which considerably influences the interactions between gp41 inhibitors and the gp41 protein. This region may play a pivotal role in suppressing gp41 inhibitors in CFR35-AD. Furthermore, gp41 can be considered a good target for subtyping analysis via the phylogenetic method.

Interaction Between Vpx and SAMHD1, Vital for SAMHD1 Inhibition

It was confirmed that the sterile alpha motif and HD domain 1 (SAMHD1) limits human immunodeficiency virus type 1 (HIV-1) replication. In contrast, viral protein x (Vpx) in HIV-2 and some simian immunodeficiency viruses can counteract this effect. The possible interaction between SAMHD1 and Vpx was suggested by previous studies; however, there are no data to confirm this interaction. Therefore, this study aimed to study the interaction between two proteins and the properties of Vpx protein for the first time using bioinformatic tools. Vpx and SAMHD1 sequences were obtained from the National Center for Biotechnology Information GenBank. Several software were used to define Vpx properties and the interaction between Vpx and different SAMHD1 isoforms. Our findings indicated the difference in interaction sites among different Vpx. However, in all Vpx proteins, this region is from amino acids 4 to 90. In addition, two regions (26-31 and 134-139) and two amino acids 425 and 429 in SAMHD1 are vital in the possible interaction. In addition, our analysis determined the physicochemical and immunological properties of the Vpx. Considering all factors, this study could confirm that Vpx interacts with SAMHD1, which could inhibit SAMHD1. Moreover, our findings can pave the way for future studies to express and purify Vpx in the laboratory and study this protein in vitro.

Transmitted Drug Resistance Against Integrase Strand Transfer Inhibitors in Iranian HIV-Infected Naïve Patients

Background

Human Immunodeficiency Virus (HIV) has claimed the lives of millions of people during the past decades. While several antiretroviral drugs like Integrase Strand Transfer Inhibitors (INSTIs) have been introduced to control HIV, Transmitted Drug Resistance (TDR) in HIV genome caused failure in treatment. This study aimed to investigate TDR and natural occurring mutations (NOPs) in HIV integrase gene in Iranian HIV patients.

Methods

In this cross-sectional study, blood samples of 30 HIV-positive patients who had never taken integrase inhibitors were considered for CD4 T cell count, RT real-time PCR, and, Nested PCR. The sequencing results were analyzed by CLC sequence viewer software and Stanford University HIV Drug Resistance Database.

Results

In all samples, nine NOPs with a high prevalence were found; however, we did not find any drug resistance mutations, except for a mutation in one sample, which showed a low resistance level. Subtype A1 was dominant in all samples.

Conclusion

Based on the findings and compared to our previous study, all patients were sustainable to main integrase inhibitors, including bictegravir, raltegravir, bictegravir, elvitegravir and dolutegravir. It seems the resistant mutation pattern attributed to integrase inhibitors was not diffent among studied patients; hence, the prescription of such inhibitors helps physicians to control HIV infection in Iranian HIV-infected patients.

First report of computational protein-ligand docking to evaluate susceptibility to HIV integrase inhibitors in HIV-infected Iranian patients

Background

Iran has recently included integrase (INT) inhibitors (INTIs) in the first-line treatment regimen in human immunodeficiency virus (HIV)-infected patients. However, there is no bioinformatics data to elaborate the impact of resistance-associated mutations (RAMs) and naturally occurring polymorphisms (NOPs) on INTIs treatment outcome in Iranian patients.

Method

In this cross-sectional survey, 850 HIV-1-infected patients enrolled; of them, 78 samples had successful sequencing results for INT gene. Several analyses were performed including docking screening, genotypic resistance, secondary/tertiary structures, post-translational modification (PTM), immune epitopes, etc.

Result

The average docking energy (E value) of different samples with elvitegravir (EVG) and raltegravir (RAL) was more than other INTIs. Phylogenetic tree analysis and Stanford HIV Subtyping program revealed HIV-1 CRF35-AD was the predominant subtype (94.9%) in our cases; in any event, online subtyping tools confirmed A1 as the most frequent subtype. For the first time, CRF-01B and BF were identified as new subtypes in Iran. Decreased CD4 count was associated with several factors: poor or unstable adherence, naïve treatment, and drug user status.

Conclusion

As the first bioinformatic report on HIV-integrase from Iran, this study indicates that EVG and RAL are the optimal INTIs in first-line antiretroviral therapy (ART) in Iranian patients. Some conserved motifs and specific amino acids in INT-protein binding sites have characterized that mutation(s) in them may disrupt INT-drugs interaction and cause a significant loss in susceptibility to INTIs. Good adherence, treatment of naïve patients, and monitoring injection drug users are fundamental factors to control HIV infection in Iran effectively.

First report on molecular docking analysis and drug resistance substitutions to approved HCV NS5A and NS5B inhibitors amongst Iranian patients

Background

NS5A and NS5B proteins of hepatitis C virus (HCV) are the main targets of compounds that directly inhibit HCV infections. However, the emergence of resistance-associated substitutions (RASs) may cause substantial reductions in susceptibility to inhibitors.

Methods

Viral load and genotyping were determined in eighty-seven naïve HCV-infected patients, and the amplified NS5A and NS5B regions were sequenced by Sanger sequencing. In addition, physicochemical properties, structural features, immune epitopes, and inhibitors-protein interactions of sequences were analyzed using several bioinformatics tools.

Results

Several amino acid residue changes were found in NS5A and NS5B proteins; however, we did not find any mutations related to resistance to the treatment in NS5B. Different phosphorylation and few glycosylation sites were assessed. Disulfide bonds were identified in both proteins that had a significant effect on the function and structure of HCV proteins. Applying reliable software to predict B-cell epitopes, 3 and 5 regions were found for NS5A and NS5B, respectively, representing a considerable potential to induce the humoral immune system. Docking analysis determined amino acids involved in the interaction of inhibitors and mentioned proteins may not decrease the drug efficiency.

Conclusions

Strong interactions between inhibitors, NS5A and NS5B proteins and the lack of efficient drug resistance mutations in the analyzed sequences may confirm the remarkable ability of NS5A and NS5B inhibitors to control HCV infection amongst Iranian patients. The results of bioinformatics analysis could unveil all features of both proteins, which can be beneficial for further investigations on HCV drug resistance and designing novel vaccines.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12876-021-01988-y.

Assessment of New E2 Protein Domain Interaction with PKR Protein to Control IFN Signaling

Introduction:

Controversy exists regarding the impact of Phosphorylation Homology Domain (PePHD) of Hepatitis C Virus (HCV) E2 protein on the interruption of the antiviral signaling pathway. A mechanism by which the virus evades the antiviral effect of interferon (IFN) alpha involves protein kinase (PKR) eukaryotic Initiation Factor 2 alpha (eIF2a) PePHD. By binding to PKR, PePHD inhibits its activity and, therefore, causes the virus to evade the antiviral activity of IFN. This study aimed to clarify the inconsistency of different conclusions reached in previous studies using reliable bioinformatics tools.

Methods:

Fifty-eight Iranian patients infected with HCV genotypes 1a and 3a and 58 healthy control individuals were examined. Plasma viral RNA was used to amplify and sequence the HCV E2 gene; also, HCV viral load, genotyping, IL-28B genotyping, alanine Aminotransferase (ALT), and aspartate aminotransferase (AST) test were determined. Bioinformatics tools determined the physicochemical properties, B-cell epitopes, post-modification changes, and secondary/tertiary structures, and also evaluated the interactions between E2 and PePHD regions.

Results:

The results showed a new domain that is responsible for binding to PKR protein and is important to the intrigued antiviral response. Physicochemical features and post-modifications were defined, showing that E2 is highly phosphorylated and there were numerous possible disulfide bonds. Secondary and tertiary structures for E2 protein were constructed. No significant relationship between ALT and AST and treatment failure was detected.

Conclusions:

Docking analysis showed a new domain in E2 protein that can be involved in the interaction between PKR and E2 protein. This finding may justify our results revealing the non-significant relationship between mutations in PePHD and treatment failure.

Subtype Classification by Polymerase and Gag Genes of HIV-1 Iranian Sequences Registered in the NCBI GenBank

Background:

The rate of Human Immunodeficiency Virus type 1 (HIV-1) infection in Iran has increased dramatically in the last few years.

Objective:

The aim of this study was to investigate the HIV subtype amongst all Iranian HIV sequences, using 8 websites.

Methods:

In this study, 637 sequences of polymerase, and gag genes of HIV-1 were obtained from NCBI. HIV-1 subtyping was done, using 8 reliable software.

Results:

The final results of the 8 online tools indicated that the majority of sequences were HIV-1 subtype CRF35 AD. However, it appeared that in some genes, a few programs could not determine specific subtypes and in some cases they described different subtypes.

Conclusion:

Considering the CRF35 AD diagram, it was clear that integrase was not an appropriate region to define this subtype. Also the full length of gag gene should be used for subtyping. For CRF1, AE envelop gene is a reliable region to define this subtype. Stanford software was used to determine the drug resistance prevalence and in 5.7% of the sequences, drug resistance mutations were found.

First report of HHV-8 viral load and seroprevalence of major blood-borne viruses in Iranian patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is characterized by autoimmune manifestations, and viral infections may have a key role in the development and progression of it. This study aimed to investigate the seroprevalence of major blood-borne viruses and HHV-8 viral load in Iranian SSc patients.

METHODS

In this cross-sectional study, 90 patients with a confirmed history of SSc and 90 healthy blood donors were enrolled. The frequency of HHV-8, CMV, EBV, HIV, HBV, and HCV antibodies and HHV-8 viral load were evaluated by enzyme-linked immunosorbent assay and real-time PCR assay, respectively.

RESULTS

HHV-8 IgG antibody was diagnosed in 61 (67.8%) patients and 3 (3.3%) healthy individuals (p<0.0001), but its genomic DNA was not detected in the patients or healthy blood donors. CMV and EBV antibodies were detected in 100% and 88.9% of SSc patients without any significant difference with healthy population (p>0.05). None of the patients or healthy population was positive for HBsAg and HIVAb; however, HCVAb was detected in two patients.

CONCLUSION

According to the results, HHV-8 antibody was uniquely increased in SSc population while its frequency in healthy population was very low. Since none of the SSc patients were positive for HHV-8 genomic DNA, the high prevalence of HHV-8 antibody in this group was not related to the real history of infection. Therefore, antibody-mediated epitope mimicry can play a role to get the high rate of seropositivity and lead to pathogeneses of SSc. Besides, CMV and EBV viral load monitoring in SSc patients can help the physician to prescribe the viral drugs to suppress the viral replication and avoid the crucial effect of reactivation.

Association of Mutations in the NS5A-PKRBD Region and IFNL4 Genotypes with Hepatitis C Interferon Responsiveness and its Functional and Structural Analysis

Background:

The cellular antiviral responses induced by interferons require some cellular protein kinase for its activation. Evidence indicated that a number of Hepatitis C Virus (HCV) proteins can repress double-stranded (ds) RNA-dependent Protein Kinase (PKR) function and help HCV to escape. However, the reports are controversial, some researchers have suggested that a region in Nonstructural 5A (NS5A) gene called Protein Kinase R-Binding Domain (PKR-BD) is associated with HCV sensitivity to the antiviral effects of Interferon (IFN). In addition, the other factor that might be associated with response to PEGylated-IFNα (Peg-IFNα) and Ribavirin (RBV) combination therapy, is IFNL4 genotypes.

Objective:

The aim of this study was to investigate the association between amino acid (aa) substitutions in the NS5A region and the IFNL4 genotypes in two Single Nucleotide Polymorphism (SNP) (rs8099917. rs12979860) in patients with HCV genotypes 1a and 3a. We also examined their response to combination therapy and the effect of these mutations on the function and structure of PKR-BD.

Methods:

Eighty-six patients with hepatitis C were recruited and follow-up for 6 months. Several tests, including alanine aminotransferase (ALT), aspartate aminotransferase (AST), viral load, IFNL4 genotyping, and PKR-BD sequencing were performed. Using several well-known and trustworthy bioinformatics tools, sequences were analyzed to define physio-chemical properties, structural features, immune epitopes and protein-protein interaction.

Results:

Of the 86 patients, 65.1% had high viral load at baseline, 64% had CT genotype for rs12979860 and 57% had GT genotype for rs8099917. Several aa residues changes were found in the PKR-BD region. We could not find any link between mutations in the PKR-BD region and different genotypes of IFNL4 in response to antiviral therapy. Regardless of pI, PKR-BD 1a and 3a showed similar physio-chemical properties, and 2 phosphorylation sites and one glycosylation site were estimated for both PKR-BD 1a and 3a. Trustworthy software were employed in order to predict B-cell epitopes, 3 regions (6-17, 26-32, 34-41) were found for both proteins, indicating a huge potential of PKR-BD protein to induce humoral immune system. Docking analysis determined non-responder sequences in both 1a and 3a genotypes to have higher energy value and are more compatible with PKR.

Conclusion:

To sum up, our results could not determine any significant relationship between mutations of PKR-BD and genotypes of IFNL4 with other factors; ALT, AST, viral load. However, docking results showed strengthened interaction between PKR-BD and PKR in non-responders that could have a momentous impact on the illness severity.

Evaluation of miR-122 Serum Level and IFN-λ3 Genotypes in Patients with Chronic HCV and HCV-Infected Liver Transplant Candidate

BACKGROUND

Alanine aminotransferase (ALT) and aspartate aminotransferase (AST) are the most common markers of liver damage, but serum level interpretation can be complicated. In hepatocytes, microRNA-122 (miR-122) is the most abundant miRs and its high expression in the serum is a characteristic of liver disease.

OBJECTIVE

We aimed to compare the circulatory level of miR-122 in patients with Chronic Hepatitis C (CHC), Hepatitis C Virus (HCV) infected Liver Transplant Candidates (LTC) and healthy controls to determine if miR-122 can be considered as an indicator of chronic and advanced stage of liver disease.

METHODS

MiR-122 serum level was measured in 170 Interferon-naïve (IFN-naïve) CHC patients, 62 LTC patients, and 132 healthy individuals via TaqMan real-time PCR. Serum levels of miR-122 were normalized to the serum level of Let-7a and miR-221. Also, the ALT and AST levels were measured.

RESULTS

ALT and AST activities and the expression of circulatory miR-122 were similar in the CHC and LTC groups, but it had significantly increased compared to healthy individuals (P<0.001 and P<0.001, respectively). Up-regulation of miR-122 in the sample of patients with normal ALT and AST activities was also observed, indicating that miR-122 is a good marker with high sensitivity and specificity for diagnosing liver damage.

CONCLUSION

miR-122 seemed to be more specific for liver diseases in comparison with the routine ALT and AST liver enzymes. Since the lower levels of circulating miR-122 were observed in the LTC group compared to the CHC group, advanced liver damages might reduce the release of miR-122 from the hepatocytes, as a sign of liver function deficiency.

Incidence of Hepatotoxicity in Iranian Patients With HIV on Antiretroviral Therapies and Its Correlation with Virologic Response to HIV Treatment

OBJECTIVE

To investigate hepatotoxicity in Iranian patients with HIV to assess the association between virologic response to HIV treatment and serum alanine aminotransferase (ALT).

METHODS

This study was conducted with 200 control patients, 75 patients with HIV naïve to antiretroviral therapy (ART), and 443 patients who received ARTs with virologic response (≤1000 copies/mL) or virologic treatment failure (>1000 copies/mL). Serum ALT level and HIV viral load were determined in all patients.

RESULTS

Patient ALT levels were significantly higher than those of control patients (45.1 ± 44.4 IU/L vs 23.8 ± 5.4 IU/L). Compared to patients who were ART-naïve, patients with ART experience had significantly higher ALT levels (38.2 ± 26.2 IU/L vs 46.3 ± 46.7 IU/L), and severe hepatotoxicity was only detected in those with ART experience (8 patients, 1.8%). Mean ALT had no significant difference between virologic response/failure groups. The ALT activity and HIV load had a negative correlation coefficient, but it was not significant.

CONCLUSION

Periodic monitoring for the possibility of hepatotoxicity is highly recommended in all patients with HIV, especially in those receiving ART treatment.

The Predictive Factors for Favorable Outcomes of Peginterferon and Ribavirin Combination Therapy in HCV-Infected Patients

We aimed to investigate the association of pretreatment host and/or viral related factors with sustained virological response (SVR) rate in chronic hepatitis C (CHC) infected patients. This cohort study was performed on 200 IFN-naïve Iranian CHC patients who were treated with pegylated interferon-α (PEG-IFN-α) plus ribavirin (RBV). Pretreatment levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), fasting blood sugar (FBS), HCV load and genotype were determined, and the pattern of changes was monitored throughout the course of treatment. The baseline FBS value in the non-responder group was significantly higher than that of the SVR group. The SVR group showed a rapid and continuous decline of ALT/AST activity from the beginning of the treatment, while the ALT level was fluctuating in non-responder and relapse groups. Persistent normalization of transaminases during combination antiviral therapy was significantly associated with SVR rate. Besides, age and FBS levels had the greatest impact on SVR. Minocycline seems to be a safe and effective adjuvant in the management of patients with schizophrenia.

Impact of IL28 Genotypes and Modeling the Interactions of HCV Core Protein on Treatment of Hepatitis C

BACKGROUND

Mutations in the core CVR region of hepatitis C virus (HCV) and polymorphisms of interleukin 28B (IL28B) are associated with progression toward liver disease and in response to therapy. In addition, interactions of the core protein with some cell interactors can be related to HCV liver damage.

AIM

This study aimed to evaluate the effect of core mutations as well as IL28B polymorphism on clinical features, sustained virological response (SVR) in 1a and 3a HCV genotypes amongst Iranian HCV infected patients, and the impact of mutations on core protein properties, antigenic properties, and interactions with HCV inhibitors, using several bioinformatics tools.

METHODS

Seventy-nine Iranian patients infected with HCV genotypes 1a and 3a and diagnosed with chronic active hepatitis were examined. Plasma viral RNA was used to amplify and sequence the HCV Core gene; also, HCV viral load, molecular genotyping, and the liver enzymes were determined for all samples. The sequencing results were analyzed by several reliable bioinformatics tools to determine the physicochemical properties, B cell epitopes, post-modification changes, and secondary/tertiary structures; and evaluate the interactions with 4 drugs by docking method.

RESULT

There were some substitutions in core CVR related to ALT and AST enzymes that can lead to HCV advanced liver disease. The most prevalent mutation for 3a genotypes was a substitution in aa 162 (I to V) while we did not find any mutation in 1a responder group. Polymorphism of the rs8099917 showed that the majority of patients had TG heterozygous and carried CT genotype at the rs12979860. Analysis indicated several phosphorylation sits for core protein as well as two important disulfide bonds. Immunogenic prediction showed that core protein can strongly induce the immune system. Interaction analysis, using the docking method revealed two potential interactors (Vitronectin and SETD2).

CONCLUSION

Generally, mutations in all core CVR regions in all patients showed a relationship between such substitutions and higher liver enzymes that can result in advanced liver disease progression in HCV infected patients. Furthermore, immunoinformatics analysis determined the possible immunodominant regions to be considered in HCV vaccine designs. Furthermore, no association between SVR and IL28B polymorphism was shown. In silico analysis determined modification sites, structures, B-cell epitopes of core protein and interactions with several interactors can lead to persistent HCV infection in the cell and the progress of liver diseases.

Interaction of Human Herpesvirus 8 Viral Interleukin-6 with Human Interleukin-6 Receptor Using In Silico Approach: The Potential Role in HHV-8 Pathogenesis

Introduction:

Human Herpesvirus 8 (HHV-8) causes classical, endemic (African), and Acquired Immunodeficiency Syndrome (AIDS)-related Kaposi’s Sarcoma (KS), Body Cavity-Based Primary Effusion Lymphomas (BCBL), HHV-8-associated peritoneal Primary Effusion Lymphoma (PEL), and Multicentric Castleman’s Disease (MCD). HHV8 genome encodes several structural and non-structural proteins, among which vIL6 is a functional homologue of Interleukin-6 (IL-6). It has been established that vIL6 plays a vital role in HHV8 infections; also, it has been suggested that its function was mediated through gp130, rather than the gp80 (IL-6 receptor [IL-6R]). This study aimed to investigate the physicochemical and structural properties as well as the immunological features, and finally the interaction between vIL6 and IL6 receptor (IL6R) by using several bioinformatics tools which could provide both valuable insight into vIL6 protein and advantageous data for further studies on HHV8 inhibitors and new vaccines.

Material and Methods:

vIL6, human IL6 (hIL6), and IL6R were obtained from NCBI GenBank and Uniport, which were aligned by The CLC Genomics Workbench. "Signal-BLAST" and “predisi" were employed to define signal peptide; also, “Expasy’sProtParam” was used to predict physicochemical properties as well as "DiANNA", and "SCRATCH" predicted the disulfide bonds. “NetPhosK”, “DISPHOS”, “NetPhos”, ”NetNGlyc”, and ”GlycoEP” were involved to determine post-modification sites. To define immunoinformatics analysis, “BcePred”, “ABCpred”, “Bepipred”, “AlgPred”, and "VaxiJen" were used. “SOPMA”, “I-TASSER”, “GalaxyRefine”, and “3D-Refine” predicted and refined the secondary and tertiary structures. TM-align server was used to align 3D structures. In addition, docking analysis was done by “Hex 5.0.”, and finally the results were illustrated by “Discovery Studio”.

Results:

A signal peptide (1-22) was defined in the vIL6 sequences and analysis has shown that vIL6 is an acidic protein which is significantly stable in all organisms. Three Disulfide bonds were predicted and immunoinformatics analysis showed 5 distinct B-cell epitopes. vIL6 is predicted as a non-allergen protein and the majority of its structure consists of Alpha helix. TM-align pointed the significant similarity between vIL6 and hIL6 in protein folding. The high energy value between vIL6 protein and IL6R was calculated and further analysis illustrated 5 conserved regions as well as 4 conserved amino acids which had a significant role in vIL6 and IL6R interaction.

Discussion:

An in silico study by numerous software determined the possible interaction between vIL6 and IL6R and the possible role of this interaction in HHV8 pathogenesis and the progress of infection. These have been overlooked by previous studies and will be beneficial to gain a more comprehensive understanding of vIL6 function during HHV8 lifecycle and infections. Structural analysis showed the significant similarity between vIL6 and hIL6 folding which can describe the similarity of the functions or interactions of both proteins. Furthermore, several conserved regions in the interaction site which interestingly were highly conserved among all vIL6 sequences can be used as new target for vIL6 inhibitors. Moreover, our results could predict immunological properties of vIL6 which suggested the ability of this protein in induction of the humoral immune response. Such a protein may be used for further studies on therapeutic vaccine fields.

The possible regions to design Human Papilloma Viruses vaccine in Iranian L1 protein

Human Papilloma Virus (HPV) genome encodes several proteins, as L1is major capsid protein and L2 is minor capsid protein. Among all HPV types HPV-16 and HPV-18 are the most common high-risk HPV (HR-HPV) types globally and the majority of cases are infected with these types. HPV entry and the initial interaction with the host cell are mainly related to the L1 protein which is the main component of HPV vaccines. The aim of this research was comparison analysis among all Iranian L1 protein sequences submitted in NCBI GenBank to find the major substitutions as well as structural and immune properties of this protein. All sequences HPV L1 protein from Iranian isolates from 2014 to 2016 were selected and obtained from NCBI data bank. "CLC Genomics Workbench" was used to translate alignment. To predict B cell epitopes, we employed several programs. Modification sites such as phosphorylation, glycosylation, and disulfide bonds were determined. Secondary and tertiary structures of all sequences were analyzed. Several mutations were found and major mutations were in amino acid residues 102, 202, 207, 292, 379, and 502. The mentioned mutations showed the minor effect on B cell and physicochemical properties of the L1 protein. Six disulfide bonds were determined in L1 protein and also in several N-link glycosylation and phosphorylation sites. Five L1 loops were determined, which had great potential to be B cell epitopes with high antigenic properties. All in all, this research as the first report from Iran described the tremendous potential of two L1 loops (BC and FG) to induce immune system which can be used as the descent candidate to design a new vaccine against HPV in the Iranian population. In addition, some differences between the reference sequence and Iranian patients' sequences were determined. It is essential to consider these differences to monitor the effectiveness and efficacy of the vaccine for the Iranian population. Our results provide a vast understanding of L1 protein that can be useful for further studies on HPV infections and new vaccine generations.

Bioinformatics Analysis of Domain 1 of HCV-Core Protein: Iran

Hepatitis C virus (HCV) infection is a serious global health problem and a cause of chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma (HCC). Bioinformatics software has been an effective tool to study the HCV genome as well as core domains. Our research was based on employing several bioinformatics software applications to find important mutations in domain 1 of core protein in Iranian HCV infected samples from 2006 to 2017, and an investigation of general properties, B-cell and T-cell epitopes, modification sites, and structure of domain 1. Domain 1 sequences of 188 HCV samples isolated from 2006 to 2017, Iran, were retrieved from NCBI gene bank. Using several tools, all sequences were analyzed for determination of mutations, physicochemical analysis, B-cell epitopes prediction, T-cell and CTL epitopes prediction, post modification, secondary and tertiary structure prediction. Our analysis determined several mutations in some special positions (70, 90, 91, and 110) that are associated with HCC and hepatocarcinogenesis, efficacy of triple therapy and sustained virological response, and interaction between core and CCR6. Several B-cell, T-cell, and CTL epitopes were recognized. Secondary and tertiary structures were mapped fordomain1 and core proteins. Our study, as a first report, offered inclusive data about frequent mutation in HCV-core gene domain 1 in Iranian sequences that can provide helpful analysis on structure and function of domain 1 of the core gene.

Functional and structural characterization of Ebola virus glycoprotein (1976–2015) — Anin silicostudy

Ebola virus (EBOV) is the causative agent of a severe hemorrhagic fever disease associated with high mortality rates in humans. This virus has five strains of which Zaire Ebola virus (ZEBOV) is the first and most important strain. It can be transmitted through contact with contaminated surfaces and objects. The genome of EBOV codes one non-structural and seven structural proteins consisting of two forms of glycoprotein (GP): soluble glycoprotein (sGP) and GP (spike). In this paper, we attempted to characterize and predict physicochemical properties, B-cell epitopes, mutation sites, phosphorylation sites, glycosylation sites, and different protein structures of EBOV GP to provide comprehensive data about changes of this GP during a 40-years course (1976–2015). GP sequences were obtained from NCBI gene bank from 1976–2015. Expasy’sProtParam, PROTSCALE, immuneepitope, Bepipred, BcePred, ABCpred, VaxiJen, DISPHOS, NetPhos, and numerous programs were used to predict and analyze all sequences. More variety of mutations were found in 2015 sequences and mutations were related to huge changes in B-cell epitopes, phosphorylation and glycosylation sites. In addition, our results determined different sites of disulfide bonds and an important mutation related to IgE epitope as well as four potent B-cell epitopes (380–387, 318–338, 405–438 and 434–475). In this study, we suggested the effect of mutations on GP properties, six positions for disulfide bonds and four phosphorylation sites for protein kinase C enzyme. Selected sequences were shown different sites for O-link and N-link glycosylation. A mutation that changed GP to an allergen protein and has a significant role in vaccine designing as well as four potent B-cell epitopes in GP protein were found.