Fibroblast growth factor 2 (FGF2) is an attractive biomaterial for pharmaceuticals and functional cosmetics. To improve the thermo-stability of FGF2, we designed two mutants harboring four-point mutations: FGF2-M1 (D28E/C78L/C96I/S137P) and FGF2-M2 (D28E/C78I/C96I/S137P) through bioinformatics, molecular thermodynamics, and molecular modeling. The D28E mutation reduced fragmentation of the FGF2 wild type during preparation, and the substitution of a whale-specific amino acid, S137P, enhanced the thermal stability of FGF2. Surface-exposed cysteines that participate in oligomerization through intermolecular disulfide bond formation were substituted with hydrophobic residues (C78L/C78I and C96I) using the in silico method. High-resolution crystal structures revealed at the atomic level that the introduction of mutations stabilizes each local region by forming more favorable interactions with neighboring residues. In particular, P137 forms CH-π interactions with the side chain indole ring of W123, which seems to stabilize a β-hairpin structure, containing a heparin-binding site of FGF2. Compared to the wild type, both FGF2-M1 and FGF2-M2 maintained greater solubility after a week at 45 °C, with their Tm values rising by ~ 5 °C. Furthermore, the duration for FGF2-M1 and FGF2-M2 to reach 50% residual activity at 45 °C extended to 8.8- and 8.2-fold longer, respectively, than that of the wild type. Interestingly, the hydrophobic substitution of surface-exposed cysteine in both FGF2 mutants makes them more resistant to proteolytic cleavage by trypsin, subtilisin, proteinase K, and actinase than the wild type and the Cys → Ser substitution. The hydrophobic replacements can influence protease resistance as well as oligomerization and thermal stability. It is notable that hydrophobic substitutions of surface-exposed cysteines, as well as D28E and S137P of the FGF2 mutants, were designed through various approaches with structural implications. Therefore, the engineering strategies and structural insights adopted in this study could be applied to improve the stability of other proteins.

The ubiquitously expressed small GTPase Ras-related protein 1B (RAP1B) acts as a molecular switch that regulates cell signaling, cytoskeletal remodeling, and cell trafficking and activates integrins in platelets and lymphocytes. The residue G12 in the P-loop is required for the RAP1B-GTPase conformational switch. Heterozygous germline RAP1B variants have been described in patients with syndromic thrombocytopenia. However, the causality and pathophysiological impact remained unexplored. We report a boy with neonatal thrombocytopenia, combined immunodeficiency, neutropenia, and monocytopenia caused by a heterozygous de novo single nucleotide substitution, c.35G>A (p.G12E) in RAP1B. We demonstrate that G12E and the previously described G12V and G60R were gain-of-function variants that increased RAP1B activation, talin recruitment, and integrin activation, thereby modifying late responses such as platelet activation, T cell proliferation, and migration. We show that in our patient, G12E was a somatic variant whose allele frequency decreased over time in the peripheral immune compartment, but remained stable in bone marrow cells, suggesting a differential effect in distinct cell populations. Allogeneic hematopoietic stem cell transplantation fully restored the patient\'s hemato-immunological phenotype. Our findings define monoallelic RAP1B gain-of-function variants as a cause for constitutive immunodeficiency and thrombocytopenia. The phenotypic spectrum ranged from isolated hematological manifestations in our patient with somatic mosaicism to complex syndromic features in patients with reported germline RAP1B variants.

The epidemic and outbreaks of influenza B Victoria lineage (Bv) during 2019-2022 led to an analysis of genetic, epitopes, charged amino acids and Bv outbreaks. Based on the National Influenza Surveillance Network (NISN), the Bv 72 strains isolated during 2019-2022 were selected by spatio-temporal sampling, then were sequenced. Using the Compare Means, Correlate and Cluster, the outbreak data were analyzed, including the single nucleotide variant (SNV), amino acid (AA), epitope, evolutionary rate (ER), Shannon entropy value (SV), charged amino acid and outbreak. With the emergence of COVID-19, the non-pharmaceutical interventions (NPIs) made Less distant transmission and only Bv outbreak. The 2021-2022 strains in the HA genes were located in the same subset, but were distinct from the 2019-2020 strains (P < 0.001). The codon G → A transition in nucleotide was in the highest ratio but the transversion of C → A and T → A made the most significant contribution to the outbreaks, while the increase in amino acid mutations characterized by polar, acidic and basic signatures played a key role in the Bv epidemic in 2021-2022. Both ER and SV were positively correlated in HA genes (R = 0.690) and NA genes (R = 0.711), respectively, however, the number of mutations in the HA genes was 1.59 times higher than that of the NA gene (2.15/1.36) from the beginning of 2020 to 2022. The positively selective sites 174, 199, 214 and 563 in HA genes and the sites 73 and 384 in NA genes were evolutionarily selected in the 2021-2022 influenza outbreaks. Overall, the prevalent factors related to 2021-2022 influenza outbreaks included epidemic timing, Tv, Ts, Tv/Ts, P137 (B → P), P148 (B → P), P199 (P → A), P212 (P → A), P214 (H → P) and P563 (B → P). The preference of amino acid mutations for charge/pH could influence the epidemic/outbreak trends of infectious diseases. Here was a good model of the evolution of infectious disease pathogens. This study, on account of further exploration of virology, genetics, bioinformatics and outbreak information, might facilitate further understanding of their deep interaction mechanisms in the spread of infectious diseases.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new emerging coronavirus that caused coronavirus disease 2019 (COVID-19). Whole-genome tracking of SARS-CoV-2 enhanced our understanding of the mechanism of the disease, control, and prevention of COVID-19.
METHODS: we analyzed 3368 SARS-CoV-2 protein sequences from Iran and compared them with 15.6 million global sequences in the GISAID database, using the Wuhan-Hu-1 strain as a reference.
RESULTS: Our investigation revealed that NSP12-P323L, ORF9c-G50N, NSP14-I42V, membrane-A63T, Q19E, and NSP3-G489S were found to be the most frequent mutations among Iranian SARS-CoV-2 sequences. Furthermore, it was observed that more than 94% of the SARS-CoV-2 genome, including NSP7, NSP8, NSP9, NSP10, NSP11, and ORF8, had no mutations when compared to the Wuhan-Hu-1 strain. Finally, our data indicated that the ORF3a-T24I, NSP3-G489S, NSP5-P132H, NSP14-I42V, envelope-T9I, nucleocapsid-D3L, membrane-Q19E, and membrane-A63T mutations might be responsible factors for the surge in the SARS-CoV-2 Omicron variant wave in Iran.
CONCLUSIONS: real-time genomic surveillance is crucial for detecting new SARS-CoV-2 variants, updating diagnostic tools, designing vaccines, and understanding adaptation to new environments.

The West Nile virus (WNV) subtype Kunjin virus (WNVKUN) is endemic to Australia. Here, we characterized the classical WNVKUN strain, OR393. The original OR393 strain contained two types of viruses: small plaque-forming virus (SP) and large plaque-forming virus (LP). The amino acid residues at positions 156 and 332 in the E protein (E156 and E332) of SP were Ser and Lys (E156S/332K), respectively, whereas those in LP were Phe and Thr (E156F/332T). SP grew slightly faster than LP in vitro. The E protein of SP was N-glycosylated, whereas that of LP was not. Analysis using two recombinant single-mutant LP viruses, rKUNV-LP-EF156S and rKUNV-LP-ET332K, indicated that E156S enlarged plaques formed by LP, but E332K potently reduced them, regardless of the amino acid at E156. rKUNV-LP-EF156S showed significantly higher neuroinvasive ability than LP, SP, and rKUNV-LP-ET332K. Our results indicate that the low-pathogenic classical WNVKUN can easily change its pathogenicity through only a few amino acid substitutions in the E protein. It was also found that Phe at E156 of the rKUNV-LP-ET332K was easily changed to Ser during replication in vitro and in vivo, suggesting that E156S is advantageous for the propagation of WNVKUN in mammalian cells.

Eotaxin-3 is a key chemokine with a relevant role in eosinophilic esophagitis, a rare chronic immune/antigen-mediated inflammatory disorder. Eotaxin-3 is a potent activator of eosinophil emergence and migration, which may lead to allergic airway inflammation. We investigated, using bioinformatics tools, the protein structure and the possible effects of the known variations reported in public databases. Following a procedure already established, we created a 3D model of the whole protein and modeled the structure of 105 protein variants due to known point mutations. The effects of the amino acid substitution at the level of impact on protein structure, stability, and possibly function were detected by the bioinformatics procedure and described in detail. A web application was implemented to browse the results of the analysis and visualize the 3D models, with the opportunity of downloading the models and analyzing them using their own software. Among 105 amino acid substitutions investigated, the study evidenced in 44 cases at least one change in any of the investigated structural parameters. Other six variations are also relevant, although a structural effect was not detected by our analysis, because they affected amino acids highly conserved, which suggests a possible function role. All these variations should be the object of particular attention, as they may induce a loss of functionality in the protein.

Recombinant mutant holotoxin BoNTs (rBoNTs) are being evaluated as possible vaccines against botulism. Previously, several rBoNTs containing 2-3 amino acid mutations in the light chain (LC) showed significant decreases in toxicity (2.5-million-fold-12.5-million-fold) versus wild-type BoNT/A1, leading to their current exclusion from the Federal Select Agent list. In this study, we added four additional mutations in the receptor-binding domain, translocation domain, and enzymatic cleft to further decrease toxicity, creating 7M rBoNT/A1. Due to poor expression in E. coli, 7M rBoNT/A1 was produced in an endogenous C. botulinum expression system. This protein had higher residual toxicity (LD50: 280 ng/mouse) than previously reported for the catalytically inactive rBoNT/A1 containing only three of the mutations (>10 µg/mouse). To investigate this discrepancy, several additional rBoNT/A1 constructs containing individual sets of amino acid substitutions from 7M rBoNT/A1 and related mutations were also endogenously produced. Similarly to endogenously produced 7M rBoNT/A1, all of the endogenously produced mutants had ~100-1000-fold greater toxicity than what was reported for their original heterologous host counterparts. A combination of mutations in multiple functional domains resulted in a greater but not multiplicative reduction in toxicity. This report demonstrates the impact of production systems on residual toxicity of genetically inactivated rBoNTs.

Human terminal deoxynucleotidyl transferase (TdT) can catalyze template-independent DNA synthesis during the V(D)J recombination and DNA repair through nonhomologous end joining. The capacity for template-independent random addition of nucleotides to single-stranded DNA makes this polymerase useful in various molecular biological applications involving sequential stepwise synthesis of oligonucleotides using modified dNTP. Nonetheless, a serious limitation to the applications of this enzyme is strong selectivity of human TdT toward dNTPs in the order dGTP > dTTP ≈ dATP > dCTP. This study involved molecular dynamics to simulate a potential impact of amino acid substitutions on the enzyme\'s selectivity toward dNTPs. It was found that the formation of stable hydrogen bonds between a nitrogenous base and amino acid residues at positions 395 and 456 is crucial for the preferences for dNTPs. A set of single-substitution and double-substitution mutants at these positions was analyzed by molecular dynamics simulations. The data revealed two TdT mutants-containing either substitution D395N or substitutions D395N+E456N-that possess substantially equalized selectivity toward various dNTPs as compared to the wild-type enzyme. These results will enable rational design of TdT-like enzymes with equalized dNTP selectivity for biotechnological applications.

The omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified in 2021 as a variant with heavy amino acid mutations in the spike protein, which is targeted by most vaccines, compared to previous variants. Amino acid substitutions in the spike proteins may alter their affinity for host viral receptors and the host interactome. Here, we found that the receptor-binding domain (RBD) of the omicron variant of SARS-CoV-2 exhibited an increased affinity for human angiotensin-converting enzyme 2, a viral cell receptor, compared to the prototype RBD. Moreover, we identified β- and γ-actin as omicron-specific binding partners of RBD. Protein complex predictions revealed that many omicron-specific amino acid substitutions affected the affinity between RBD of the omicron variant and actin. Our findings indicate that proteins localized to different cellular compartments exhibit strong binding to the omicron RBD.

Mosquito-borne Zika virus (ZIKV) from sub-Saharan Africa has recently gained attention due to its epidemic potential and its capacity to be highly teratogenic. To improve our knowledge on currently circulating strains of African ZIKV, we conducted protein sequence alignment and identified contemporary West Africa NS1 (NS1CWA) protein as a highly conserved viral protein. Comparison of NS1CWA with the NS1 of the historical African ZIKV strain MR766 (NS1MR766), revealed seven amino acid substitutions. The effects of NS1 mutations on protein expression, virus replication, and innate immune activation were assessed in human cells using recombinant NS1 proteins and a chimeric viral clone MR766 with NS1CWA replacing NS1MR766. Our data indicated higher secretion efficiency of NS1CWA compared to NS1MR766 associated with a change in subcellular distribution. A chimeric MR766 virus with NS1CWA instead of authentic protein displayed a greater viral replication efficiency, leading to more pronounced cell death compared to parental virus. Enhanced viral growth was associated with reduced activation of innate immunity. Our data raise questions of the importance of NS1 protein in the pathogenicity of contemporary ZIKV from sub-Saharan Africa and point to differences within viral strains of African lineage.