Most algorithms that are used to predict the effects of variants rely on evolutionary conservation. However, a majority of such techniques compute evolutionary conservation by solely using the alignment of multiple sequences while overlooking the evolutionary context of substitution events. We had introduced PHACT, a scoring-based pathogenicity predictor for missense mutations that can leverage phylogenetic trees, in our previous study. By building on this foundation, we now propose PHACTboost, a gradient boosting tree-based classifier that combines PHACT scores with information from multiple sequence alignments, phylogenetic trees, and ancestral reconstruction. By learning from data PHACTboost outperforms PHACT. Furthermore, the results of comprehensive experiments on carefully constructed sets of variants demonstrated that PHACTboost can outperform 40 prevalent pathogenicity predictors reported in the dbNSFP, including conventional tools, meta-predictors, and deep learning-based approaches as well as more recent tools such as, AlphaMissense, EVE, and CPT-1. The superiority of PHACTboost over these methods was particularly evident in case of hard variants for which different pathogenicity predictors offered conflicting results. We provide predictions of 215 million amino acid alterations over 20,191 proteins. PHACTboost is available at https://github.com/CompGenomeLab/PHACTboost. PHACTboost can improve our understanding of genetic diseases and facilitate more accurate diagnoses.

Monitoring the genetic variability of human respiratory syncytial virus (hRSV) is of paramount importance, especially for the potential implication of key antigenic mutations on the emergence of immune escape variants. Thus, to describe the genetic diversity and evolutionary dynamics of hRSV circulating in Sicily (Italy), a total of 153 hRSV whole-genome sequences collected from 770 hRSV-positive subjects between 2017 and 2023, before the introduction of expanded immunization programs into the population, were investigated. The phylogenetic analyses indicated that the genotypes GA.2.3.5 (ON1) for hRSV-A and GB.5.0.5a (BA9) for hRSV-B co-circulated in our region. Amino acid (AA) substitutions in the surface and internal proteins were evaluated, including the F protein antigenic sites, as the major targets of immunoprophylactic monoclonal antibodies and vaccines. Overall, the proportion of AA changes ranged between 1.5% and 22.6% among hRSV-A, whereas hRSV-B varied in the range 0.8-16.9%; the latter was more polymorphic than hRSV-A within the key antigenic sites. No AA substitutions were found at site III of both subgroups. Although several non-synonymous mutations were found, none of the polymorphisms known to potentially affect the efficacy of current preventive measures were documented. These findings provide new insights into the global hRSV molecular epidemiology and highlight the importance of defining a baseline genomic picture to monitor for future changes that might be induced by the selective pressures of immunological preventive measures, which will soon become widely available.

When designing live-attenuated respiratory syncytial virus (RSV) vaccine candidates, attenuating mutations can be developed through biologic selection or reverse-genetic manipulation and may include point mutations, codon and gene deletions, and genome rearrangements. Attenuation typically involves the reduction in virus replication, due to direct effects on viral structural and replicative machinery or viral factors that antagonize host defense or cause disease. However, attenuation must balance reduced replication and immunogenic antigen expression. In the present study, we explored a new approach in order to discover attenuating mutations. Specifically, we used protein structure modeling and computational methods to identify amino acid substitutions in the RSV nonstructural protein 1 (NS1) predicted to cause various levels of structural perturbation. Twelve different mutations predicted to alter the NS1 protein structure were introduced into infectious virus and analyzed in cell culture for effects on viral mRNA and protein expression, interferon and cytokine expression, and caspase activation. We found the use of structure-based machine learning to predict amino acid substitutions that reduce the thermodynamic stability of NS1 resulted in various levels of loss of NS1 function, exemplified by effects including reduced multi-cycle viral replication in cells competent for type I interferon, reduced expression of viral mRNAs and proteins, and increased interferon and apoptosis responses.

In our prior investigations, we elucidated the role of the tryptophan-to-tyrosine substitution at the 61st position in the nonstructural protein NSsW61Y in diminishing the interaction between nonstructural proteins (NSs) and nucleoprotein (NP), impeding viral replication. In this study, we focused on the involvement of NSs in replication via the modulation of autophagosomes. Initially, we examined the impact of NP expression levels, a marker for replication, upon the infection of HeLa cells with severe fever thrombocytopenia syndrome virus (SFTSV), with or without the inhibition of NP binding. Western blot analysis revealed a reduction in NP levels in NSsW61Y-expressing conditions. Furthermore, the expression levels of the canonical autophagosome markers p62 and LC3 decreased in HeLa cells expressing NSsW61Y, revealing the involvement of individual viral proteins on autophagy. Subsequent experiments confirmed that NSsW61Y perturbs autophagy flux, as evidenced by reduced levels of LC3B and p62 upon treatment with chloroquine, an inhibitor of autophagosome-lysosome fusion. LysoTracker staining demonstrated a decrease in lysosomes in cells expressing the NS mutant compared to those expressing wild-type NS. We further explored the mTOR-associated regulatory pathway, a key regulator affected by NS mutant expression. The observed inhibition of replication could be linked to conformational changes in the NSs, impairing their binding to NP and altering mTOR regulation, a crucial upstream signaling component in autophagy. These findings illuminate the intricate interplay between NSsW61Y and the suppression of host autophagy machinery, which is crucial for the generation of autophagosomes to facilitate viral replication.

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

It has been shown previously that a set of three modifications-termed S1, Crystal Kappa, and elbow-act synergistically to improve the crystallizability of an antigen-binding fragment (Fab) framework. Here, we prepared a phage-displayed library and performed crystallization screenings to identify additional substitutions-located near the heavy-chain elbow region-which cooperate with the S1, Crystal Kappa, and elbow modifications to increase expression and improve crystallizability of the Fab framework even further. One substitution (K141Q) supports the signature Crystal Kappa-mediated Fab:Fab crystal lattice packing interaction. Another substitution (E172G) improves the compatibility of the elbow modification with the Fab framework by alleviating some of the strain incurred by the shortened and bulkier elbow linker region. A third substitution (F170W) generates a split-Fab conformation, resulting in a powerful crystal lattice packing interaction comprising the biological interaction interface between the variable heavy and light chain domains. In sum, we have used K141Q, E172G, and F170W substitutions-which complement the S1, Crystal Kappa, and elbow modifications-to generate a set of highly crystallizable Fab frameworks that can be used as chaperones to enable facile elucidation of Fab:antigen complex structures by x-ray crystallography.

Conserved tryptophan residues are critical for the structure and the stability of β/γ-crystallin in the lenses of vertebrates. During aging, in which the lenses are continuously exposed to ultraviolet irradiation and other environmental stresses, oxidation of tryptophan residues in β/γ-crystallin is triggered and impacts the lens proteins to varying degrees. Kynurenine derivatives, formed by oxidation of tryptophan, accumulate, resulting in destabilization and insolubilization of β/γ-crystallin, which correlates with age-related cataract formation. To understand the contribution of tryptophan modification on the structure and stability of human βB2-crystallin, five tryptophan residues were mutated to phenylalanine considering its similarity in structure and hydrophilicity to kynurenine. Among all mutants, W59F and W151F altered the stability and homo-oligomerization of βB2-crystallin-W59F promoted tetramerization whereas W151F blocked oligomerization. Most W59F dimers transformed into tetramer in a month, and the separated dimer and tetramer of W59F demonstrated different structures and hydrophobicity, implying that the biochemical properties of βB2-crystallin vary over time. By using SAXS, we found that the dimer of βB2-crystallin in solution resembled the lattice βB1-crystallin dimer (face-en-face), whereas the tetramer of βB2-crystallin in solution resembled its lattice tetramer (domain-swapped). Our results suggest that homo-oligomerization of βB2-crystallin includes potential inter-subunit reactions, such as dissociation, unfolding, and re-formation of the dimers into a tetramer in solution. The W>F mutants are useful in studying different folding states of βB2-crystallin in lens.

BACKGROUND: HIV-1 produces Tat, a crucial protein for transcription, viral replication, and CNS neurotoxicity. Tat interacts with TAR, enhancing HIV reverse transcription. Subtype C Tat variants (C31S, R57S, Q63E) are associated with reduced transactivation and neurovirulence compared to subtype B. However, their precise impact on Tat-TAR binding is unclear. This study investigates how these substitutions affect Tat-TAR interaction.
METHODS: We utilized molecular modelling techniques, including MODELLER, to produce precise three-dimensional structures of HIV-1 Tat protein variants. We utilized Tat subtype B as the reference or wild type, and generated Tat variants to mirror those amino acid variants found in Tat subtype C. Subtype C-specific amino acid substitutions were selected based on their role in the neuropathogenesis of HIV-1. Subsequently, we conducted molecular docking of each Tat protein variant to TAR using HDOCK, followed by molecular dynamic simulations.
RESULTS: Molecular docking results indicated that Tat subtype B (TatWt) showed the highest affinity for the TAR element (-262.07), followed by TatC31S (-261.61), TatQ63E (-256.43), TatC31S/R57S/Q63E (-238.92), and TatR57S (-222.24). However, binding free energy analysis showed higher affinities for single variants TatQ63E (-349.2 ± 10.4 kcal/mol) and TatR57S (-290.0 ± 9.6 kcal/mol) compared to TatWt (-247.9 ± 27.7 kcal/mol), while TatC31S and TatC31S/R57SQ/63E showed lower values. Interactions over the protein trajectory were also higher for TatQ63E and TatR57S compared to TatWt, TatC31S, and TatC31S/R57SQ/63E, suggesting that modifying amino acids within the Arginine/Glutamine-rich region notably affects TAR interaction. Single amino acid mutations TatR57S and TatQ63E had a significant impact, while TatC31S had minimal effect. Introducing single amino acid variants from TatWt to a more representative Tat subtype C (TatC31S/R57SQ/63E) resulted in lower predicted binding affinity, consistent with previous findings.
CONCLUSIONS: These identified amino acid positions likely contribute significantly to Tat-TAR interaction and the differential pathogenesis and neuropathogenesis observed between subtype B and subtype C. Additional experimental investigations should prioritize exploring the influence of these amino acid signatures on TAR binding to gain a comprehensive understanding of their impact on viral transactivation, potentially identifying them as therapeutic targets.

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK\'s Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a \"stability rheostat.\" The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

We report a novel HLA-DRB3*03 allele, now named DRB3*03:65, identified by next-generation sequencing.