An in vitro and computational validation of a novel loss-of-functional mutation in PAX9 associated with non-syndromic tooth agenesis

Congenital tooth agenesis (CTA) is one of the most common craniofacial anomalies. Its frequency varies among different population depending upon the genetic heterogeneity. CTA could be of familial or sporadic and syndromic or non-syndromic. Five major genes are found to be associated with non-syndromic CTA namely, PAX9, MSX1, EDA1, AXIN2 and WNT10A. In this study, an India family with CTA was investigated and a novel c.336C>G variation was identified in the exon 3 of PAX9, leading to substitution of evolutionary conserved Cys with Trp at 112 amino acid position located at the functionally significant DNA binding paired domain region. Functional analysis revealed that p.Cys112Trp mutation did not prevent the nuclear localization although mutant protein had higher cytoplasmic retention. EMSA using e5 probe revealed that mutant protein was unable to bind with the paired-domain binding site. Subsequently, GST pull-down assay revealed lower binding activity of the mutant protein with its known interactor MSX1. Further RNA-sequencing of PAX9 over-expressed HEK293, identified two potential novel targets, WNT4 and WNT7b those are up-regulated by wild-type PAX9 but not by mutant. These in vitro results were consistent with the computational results. The in vitro and computational observations altogether suggest that c.336C>G (p.Cys112Trp) variation leads to loss-of-function of PAX9 leading to CTA in this family.

RSK1 SUMOylation is required for KSHV lytic replication

RSK1, a downstream kinase of the MAPK pathway, has been shown to regulate multiple cellular processes and is essential for lytic replication of a variety of viruses, including Kaposi’s sarcoma-associated herpesvirus (KSHV). Besides phosphorylation, it is not known whether other post-translational modifications play an important role in regulating RSK1 function. We demonstrate that RSK1 undergoes robust SUMOylation during KSHV lytic replication at lysine residues K110, K335, and K421. SUMO modification does not alter RSK1 activation and kinase activity upon KSHV ORF45 co-expression, but affects RSK1 downstream substrate phosphorylation. Compared to wild-type RSK1, the overall phosphorylation level of RxRxxS*/T* motif is significantly declined in RSK1K110/335/421R expressing cells. Specifically, SUMOylation deficient RSK1 cannot efficiently phosphorylate eIF4B. Sequence analysis showed that eIF4B has one SUMO-interacting motif (SIM) between the amino acid position 166 and 170 (166IRVDV170), which mediates the association between eIF4B and RSK1 through SUMO-SIM interaction. These results indicate that SUMOylation regulates the phosphorylation of RSK1 downstream substrates, which is required for efficient KSHV lytic replication.

Prediction of Protein–Protein Interaction Sites Based on Stratified Attentional Mechanisms

Proteins are the basic substances that undertake human life activities, and they often perform their biological functions through interactions with other biological macromolecules, such as cell transmission and signal transduction. Predicting the interaction sites between proteins can deepen the understanding of the principle of protein interactions, but traditional experimental methods are time-consuming and labor-intensive. In this study, a new hierarchical attention network structure, named HANPPIS, by adding six effective features of protein sequence, position-specific scoring matrix (PSSM), secondary structure, pre-training vector, hydrophilic, and amino acid position, is proposed to predict protein–protein interaction (PPI) sites. The experiment proved that our model has obtained very effective results, which was better than the existing advanced calculation methods. More importantly, we used the double-layer attention mechanism to improve the interpretability of the model and to a certain extent solved the problem of the “black box” of deep neural networks, which can be used as a reference for location positioning on the biological level.

SARS-CoV-2 Variants, RBD Mutations, Binding Affinity, and Antibody Escape

Since 2020, the receptor-binding domain (RBD) of the spike protein of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been constantly mutating, producing most of the notable missense mutations in the context of “variants of concern”, probably in response to the vaccine-driven alteration of immune profiles of the human population. The Delta variant, in particular, has become the most prevalent variant of the epidemic, and it is spreading in countries with the highest vaccination rates, causing the world to face the risk of a new wave of the contagion. Understanding the physical mechanism responsible for the mutation-induced changes in the RBD’s binding affinity, its transmissibility, and its capacity to escape vaccine-induced immunity is the “urgent challenge” in the development of preventive measures, vaccines, and therapeutic antibodies against the coronavirus disease 2019 (COVID-19) pandemic. In this study, entropy–enthalpy compensation and the Gibbs free energy change were used to analyze the impact of the RBD mutations on the binding affinity of SARS-CoV-2 variants with the receptor angiotensin converting enzyme 2 (ACE2) and existing antibodies. Through the analysis, we found that the existing mutations have already covered almost all possible detrimental mutations that could result in an increase of transmissibility, and that a possible mutation in amino-acid position 498 of the RBD can potentially enhance its binding affinity. A new calculation method for the binding energies of protein–protein complexes is proposed based on the entropy–enthalpy compensation rule. All known structures of RBD–antibody complexes and the RBD–ACE2 complex comply with the entropy–enthalpy compensation rule in providing the driving force behind the spontaneous protein–protein docking. The variant-induced risk of breakthrough infections in vaccinated people is attributed to the L452R mutation’s reduction of the binding affinity of many antibodies. Mutations reversing the hydrophobic or hydrophilic performance of residues in the spike RBD potentially cause breakthrough infections of coronaviruses due to the changes in geometric complementarity in the entropy–enthalpy compensations between antibodies and the virus at the binding sites.

In silico molecular and morphological analysis of rice blast resistant gene Pi-ta in Sri Lankan rice germplasm

Background Pi-ta is a major blast resistant gene, introgressed from indica rice varieties. In this study, diversity of the Pi-ta gene of 47 Sri Lankan rice accessions was studied by bioinformatics, and the results were validated with molecular and disease reaction assays. Sequences of rice accessions at the locus Os12g0281300 were retrieved from Rice SNP-Seek Database, and the coding sequence of reference Pi-ta gene of cultivar Tetep (accession no. GQ918486.1) was obtained from GenBank. Comparisons were made at nucleotide, amino acid, and protein structure level, and the 3D models predicted using Phyre2 software were superimposed using TM-align software. Results In silico analysis revealed that 10 accessions possessed resistant allele of the Pi-ta gene. The remaining accessions recorded high polymorphism in the leucine-rich domain resulting in 9 allele types, leading to single–amino acid substitutions at 27 different positions including a functional mutation of alanine to serine at the 918th amino acid position. None of the genotypes led to truncations in the amino acid sequence. The in silico analysis results were validated on 23 accessions comprising resistant and susceptible genotypes and another 25 cultivars from Northern Sri Lanka, by molecular assay using YL183/YL87 and YL155/YL87 resistant and susceptible allele-specific markers. Resistance of Pi-ta gene for the causal fungus, Magnaporthe oryzae, was further validated through pathogenicity assay. Conclusion The Pi-ta gene, especially the LRD region, revealed significant variations within Sri Lankan rice cultivars leading to high levels of resistance against blast. This information would be highly useful in breeding programmes for resistance against rice blast.

A Cross-Reactive Monoclonal Antibody Against Neuraminidases of Both H9N2 and H3N2 Influenza Viruses Shows Protection in Mice Challenging Models

Neuraminidases (NAs) of H9N2 avian influenza virus (AIV) and H3N2 human seasonal influenza virus (HSIV) share similar antigenic structures. However, there are few reports on epitopes shared by these two NAs. We previously reported a monoclonal antibody (mAb) 1G8 against the NA of H9N2 AIV with neuraminidase inhibition (NI) ability. In this study, 1G8 was shown to cross-react with and inhibit the NA of H3N2 HSIV. In a passive transfer experiment, 1G8 provided protection to mice challenged with rescued H1N2 viruses carrying H9N2 NA or H3N2 NA. Mutation at amino acid position 199 was also selected and proved to be crucial for H3N2 HSIV to escape from mAb 1G8. Moreover, we found that residue 199 contributed to inducing broad protective antibodies without the influence of the N-linked glycosylation at amino acid position 200 in NAs. Residues as residue 199, which are not shielded by glycosylation modification, would form ideal epitopes for developing universal vaccine and protective antibodies.

Mutations at a split codon in the GTPase-encoding domain of OPA1 cause dominant optic atrophy through different molecular mechanisms

Exonic (i.e. coding) variants in genes associated with disease can exert pathogenic effects both at the protein and mRNA level, either by altering the amino acid sequence or by affecting pre-mRNA splicing. The latter is often neglected due to the lack of RNA analyses in genetic diagnostic testing. In this study we considered both pathomechanisms and performed a comprehensive analysis of nine exonic nucleotide changes in OPA1, which is the major gene underlying autosomal dominant optic atrophy (DOA) and is characterized by pronounced allelic heterogeneity. We focused on the GTPase-encoding domain of OPA1, which harbors most of the missense variants associated with DOA. Given that the consensus splice sites extend into the exons, we chose a split codon, namely codon 438, for our analyses. Variants at this codon are the second most common cause of disease in our large cohort of DOA patients harboring disease-causing variants in OPA1. In silico splice predictions, heterologous splice assays, analysis of patient’s RNA when available, and protein modeling revealed different molecular outcomes for variants at codon 438. The wildtype aspartate residue at amino acid position 438 is directly involved in the dimerization of OPA1 monomers. We found that six amino acid substitutions at codon 438 (i.e. all substitutions of the first and second nucleotide of the codon) destabilized dimerization while only substitutions of the first nucleotide of the codon caused exon skipping. Our study highlights the value of combining RNA analysis and protein modeling approaches to accurately assign patients to future precision therapies.

Spike Protein NTD mutation G142D in SARS-CoV-2 Delta VOC lineages is associated with frequent back mutations, increased viral loads, and immune evasion

The significantly greater infectivity of the SARS-CoV-2 Delta variants of concern (VOC) is hypothesized to be driven by key mutations that result in increased transmissibility, viral load and/or evasion of host immune response. We surveyed the mutational profiles of Delta VOC genomes between September 2020 and mid-August 2021 and identified a previously unreported mutation pattern at amino acid position 142 in the N-terminal domain (NTD) of the spike protein which demonstrated multiple rounds of mutation from G142 to D142 and back. This pattern of frequent back mutations was observed at multiple time points and across Delta VOC sub-lineages. The etiology for these recurrent mutations is unclear but raises the possibility of host-directed editing of the SARS-CoV-2 genome. Within Delta VOC this mutation is associated with higher viral load, further enhanced in the presence of another NTD mutation (T95I) which was also frequently observed in these cases. Protein modeling of both mutations predicts alterations of the surface topography of the NTD by G142D, specifically disturbance of the ‘super site’ epitope that binds NTD-directed neutralizing antibodies (NAbs). The appearance of frequent and repeated G142D followed by D142G back mutations is previously unreported in SARS-CoV-2 and may represent viral adaptation to evolving host immunity characterized by increasing frequency of spike NAbs, from both prior infection and vaccine-based immunity. The emergence of alterations of the NTD in and around the main NAb epitope is a concerning development in the ongoing evolution of SARS-CoV-2 which may contribute to increased infectivity, immune evasion and ‘breakthrough infections’ characteristic of Delta VOC. Future vaccine and therapy development may benefit by recognizing the emergence of these novel spike NTD mutations and considering their impact on antibody recognition, viral neutralization, infectivity, replication, and viral load.

MLST/MVLST Analysis and Antibiotic Resistance of Vibrio cholerae in Shandong Province of China

Background: Vibrio cholerae is an important bacterium causing profuse watery diarrhea. Cholera had swept the whole Shandong province from 1975 to 2013. Methods: From epidemiological data and pulsed-field gel electrophoresis data, we selected 86 V. cholerae isolates appearing in Shandong Province in China from 1975 to 2013 and characterized them by multilocus sequence typing (MLST)/multi-virulence locus sequence typing (MVLST), antibiogram and analysis of genes related to antibiotic resistance. Results: Combined MLST/MVLST data revealed 33 sequence types and a major group. Within the group, 3 subgroups (ST1, ST24 and ST29) were revealed, prevalent in the strains isolated during the 1980s, 1990s and 21st century, respectively. All the O1 isolates after 1990 were found to be El Tor variants harboring the classical ctxB gene. The tcpA gene of O139 strains had a mutation at amino acid position 62 (N→D). Antibiotic resistance of V. cholerae increased over time. Most El Tor variants between 1998 and 1999 were resistant to trimethoprim/sulfamethoxazole. The O139 strain, since its appearance in 1997, had significantly broader spectrum of antibiotic resistance than O1 variants. The presence of the SXT element corresponds to the trend of growing drug resistance. Conclusion: The analysis of genotypic polymorphism and enhanced resistance of V. cholerae indicated continuous variation and evolution of this pathogenic agent in Shandong Province.

Delta plus variant with neutral mutation is less virulent compared to wild type SARS-CoV2

BackgroundWith the start of the Coronavirus disease19 (COVID19) pandemic, the Coronavirus has mutated constantly. Recently a new variant called Delta plus has been reported in few countries, including South Africa, Brazil and India. The Delta plus variant contains an additional mutation called K417N on the coronavirus spike. The present study aims to determine the virulence and transmissibility of the Delta plus variant and to check the efficiency of different antibodies on its neutralization.Materials and MethodsDifferent computational tools such as PROVEAN, an online tool, HOPE server, simulation using CABS Flex, Clus pro, an online docking tool, were used to predict the structure and function of Delta plus variant by performing a comparative study with wild type protein. Also, to find an effective antibody against Delta plus variant, antigen-antibody docking studies were conducted through Clus pro server. Furthermore, we performed a 2D interaction diagram analysis to find the amino acid residue's interaction against antibodies.Results PROVEAN and HOPE showed the mutation (K417N) in the S-glycoprotein of Delta plus as NEUTRAL mutation. This mutation causes the loss of cysteine bonds leading to the destabilization of the 3D structure of spike protein. Furthermore, the RMSF plot emphasizing the 17th amino acid position of wild and Delta plus mutant revealed the high fluctuation of mutant protein structure compared to the wild protein structure. Further, a comparative docking study against hACE2 shows higher binding energy of wild-type RBD (-751.7 kcal/mol) than mutant RBD (-750.1 kcal/mol). Moreover, antigen-antibody docking study revealed higher affinity of BD-23 Fab antibodies with greater interaction energy ( -997 kcal/mol) compared to other antibodies and thus may prove to be a promising therapeutic against Delta plus variant.ConclusionDelta plus variant is less stable, has a lower binding affinity to hACE2 and has less virulence than wild type. However, the BD-23 Fab antibody has shown a more significant association for this variant and can be used in its treatment.