Structural effects of charge destabilization and amino acid substitutions in amyloid fragments of CsgA

The most studied functional amyloid is the CsgA, major curli subunit protein, which is produced by numerous strains of Enterobacteriaceae. Although CsgA sequences are highly conserved, they exhibit species diversity, which reflects the specific evolutionary and functional adaptability of the major curli subunit. Herein, we performed bioinformatics analyses to uncover the differences in the amyloidogenic properties of the R4 fragments in Escherichia coli and Salmonella enterica and proposed four mutants for more detailed studies: M1, M2, M3, and M4. The mutated sequences were characterized by various experimental techniques, such as circular dichroism, ATR-FTIR, FT-Raman, thioflavin T, transmission electron microscopy and confocal microscopy. Additionally, molecular dynamics simulations were performed to determine the role of buffer ions in the aggregation process. Our results demonstrated that the aggregation kinetics, fibril morphology, and overall structure of the peptide were significantly affected by the positions of charged amino acids within the repeat sequences of CsgA. Notably, substituting glycine with lysine resulted in the formation of distinctive spherically packed globular aggregates. The differences in morphology observed are attributed to the influence of phosphate ions, which disrupt the local electrostatic interaction network of the polypeptide chains. This study provides knowledge on the preferential formation of amyloid fibrils based on charge states within the polypeptide chain.

Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR geneComputational analysis of bovine PKLR gene

Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins' stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle.Communicated by Ramaswamy H. Sarma.

Selective pressure mediated by influenza virus M158-66 epitope-specific CD8+T cells promotes accumulation of extra-epitopic amino acid substitutions associated with viral resistance to these T cells

Influenza viruses are notorious for their capacity to evade host immunity. Not only can they evade recognition by virus-neutralizing antibodies, there is also evidence that they accumulate mutations in epitopes recognized by virus-specific CD8+T cells. In addition, we have shown previously that human influenza A viruses were less well recognized than avian influenza viruses by CD8+T cells directed to the highly conserved, HLA-A*02:01 restricted M158-66 epitope located in the Matrix 1 (M1) protein. Amino acid differences at residues outside the epitope were responsible for the differential recognition, and it was hypothesized that this reflected immune adaptation of human influenza viruses to selective pressure exerted by M158-66-specific CD8+T cells in the human population. In the present study, we tested this hypothesis and investigated if selective pressure exerted by M158-66 epitope-specific CD8+T cells could drive mutations at the extra-epitopic residues in vitro. To this end, isogenic influenza A viruses with the M1 gene of a human or an avian influenza virus were serially passaged in human lung epithelial A549 cells that transgenically express the HLA-A*02:01 molecule or not, in the presence or absence of M158-66 epitope-specific CD8+T cells. Especially in the virus with the M1 gene of an avian influenza virus, variants emerged with mutations at the extra-epitopic residues associated with reduced recognition by M158-66-specific T cells as detected by Next Generation Sequencing. Although the emergence of these variants was observed in the absence of selective pressure exerted by M158-66 epitope-specific CD8+T cells, their proportion was much larger in the presence of this selective pressure.

Toxic to the touch: The makings of lethal mantles in pitohui birds and poison dart frogs

How do chemically defended animals resist their own toxins? This intriguing question on the concept of autotoxicity is at the heart of how species interactions evolve. In this issue of Molecular Ecology (Molecular Ecology, 2024, 33), Bodawatta and colleagues report on how Papua New Guinean birds coopted deadly neurotoxins to create lethal mantles that protect against predators and parasites. Combining chemical screening of the plumage of a diverse collection of passerine birds with genome sequencing, the researchers unlocked a deeper understanding of how some birds sequester deadly batrachotoxin (BTX) from their food without poisoning themselves. They identified that birds impervious to BTX bear amino acid substitutions in the toxin-binding site of the voltage-gated sodium channel Nav1.4, whose function is essential for proper contraction and relaxation of vertebrate muscles. Comparative genetic and molecular docking analyses show that several of the substitutions associated with insensitivity to BTX may have become prevalent among toxic birds through positive selection. Intriguingly, poison dart frogs that also co-opted BTX in their lethal mantles were found to harbour similar toxin insensitivity substitutions in their Nav1.4 channels. Taken together, this sets up a powerful model system for studying the mechanisms behind convergent molecular evolution and how it may drive biological diversity.

PrimPol Variant V102A with Altered Primase and Polymerase Activities

PrimPol is a human DNA primase-polymerase which restarts DNA synthesis beyond DNA lesions and non-B DNA structures blocking replication. Disfunction of PrimPol in cells leads to slowing of DNA replication rates in mitochondria and nucleus, accumulation of chromosome aberrations, cell cycle delay, and elevated sensitivity to DNA-damaging agents. A defective PrimPol has been suggested to be associated with the development of ophthalmic diseases, elevated mitochondrial toxicity of antiviral drugs and increased cell resistance to chemotherapy. Here, we describe a rare missense PrimPol variant V102A with altered biochemical properties identified in patients suffering from ovarian and cervical cancer. The Val102 to Ala substitution dramatically reduced both the primase and DNA polymerase activities of PrimPol as well as specifically decreased its ability to incorporate ribonucleotides. Structural analysis indicates that the V102A substitution can destabilize the hydrophobic pocket adjacent to the active site, affecting dNTP binding and catalysis.

Substitution of both histidines in the active site of yeast alcohol dehydrogenase 1 exposes underlying pH dependencies

Histidine residues 44 and 48 in yeast alcohol dehydrogenase (ADH) bind to the coenzymes NAD(H) and contribute to catalysis. The individual H44R and H48Q substitutions alter the kinetics and pH dependencies, and now the roles of other ionizable groups in the enzyme were studied in the doubly substituted H44R/H48Q ADH. The substitutions make the enzyme more resistant to inactivation by diethyl pyrocarbonate, modestly improve affinity for coenzymes, and substantially decrease catalytic efficiencies for ethanol oxidation and acetaldehyde reduction. The pH dependencies for several kinetic parameters are shifted from pK values for wild-type ADH of 7.3-8.1 to values for H44R/H48Q ADH of 8.0-9.6, and are assigned to the water or alcohol bound to the catalytic zinc. It appears that the rate of binding of NAD+ is electrostatically favored with zinc-hydroxide whereas binding of NADH is faster with neutral zinc-water. The pH dependencies of catalytic efficiencies (V/EtKm) for ethanol oxidation and acetaldehyde reduction are similarly controlled by deprotonation and protonation, respectively. The substitutions make an enzyme that resembles the homologous horse liver H51Q ADH, which has Arg-47 and Gln-51 and exhibits similar pK values. In the wild-type ADHs, it appears that His-48 (or His-51) in the proton relay systems linked to the catalytic zinc ligands modulate catalytic efficiencies.

A comprehensive analysis of SARS-CoV-2 missense mutations indicates that all possible amino acid replacements in the viral proteins occurred within the first two-and-a-half years of the pandemic

The surveillance of COVID-19 pandemic has led to the determination of millions of genome sequences of the SARS-CoV-2 virus, with the accumulation of a wealth of information never collected before for an infectious disease. Exploring the information retrieved from the GISAID database reporting at that time >13 million genome sequences, we classified the 141,639 unique missense mutations detected in the first two-and-a-half years (up to October 2022) of the pandemic. Notably, our analysis indicates that 98.2 % of all possible conservative amino acid replacements occurred. Even non-conservative mutations were highly represented (73.9 %). For a significant number of residues (3 %), all possible replacements with the other nineteen amino acids have been observed. These observations strongly indicate that, in this time interval, the virus explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain and that those that are not observed severely affect SARS-CoV-2 integrity. The implications of the present findings go well beyond the structural biology of SARS-CoV-2 as the huge amount of information here collected and classified may be valuable for the elucidation of the sequence-structure-function relationships in proteins.

Effect of S-region mutations on HBsAg in HBsAg-negative HBV-infected patients

BACKGROUND

Occult HBV infection (OBI) is a special form of hepatitis B virus (HBV) infection that may cause Liver cirrhosis and hepatocellular carcinoma, causing significant harm to patients. Given the insidious nature of OBI, it is usually not easy to be detected. Most of the samples currently studied are concentrated on blood donors, however, patients in this special state have not been fully studied. This project aimed to study the effect of HBV S region mutations on HBsAg in patients with clinical OBI.

METHODS

Collect 107 HBsAg-/HBV DNA + blood samples from Beijing Youan Hospital, Capital Medical University from August 2022 to April 2023. Next, the successfully extracted and amplified HBV DNA S regions were sequenced. Construct mutant plasmids to verify the cell function of the high-frequency mutation sites and explore the possible molecular mechanism.

RESULTS

Sixty-eight HBsAg-negative samples were sequenced, revealing high-frequency amino acid substitution sites in the HBV S protein, including immune escape mutations (i.e., sY100C、sK122R、sI126T、sT131P、and sS114T) and TMD (Transmembrane domain) region substitutions (i.e., sT5A、sG10D、sF20S、and sS3N). We constructed a portion of the mutant plasmids and found that sT5A, sF20S, sG10D, sS3N, sI68T, and sI126T single point mutations or combined mutations may decrease HBsAg expression or change the antigenicity of HBsAg leading to detection failure.

CONCLUSIONS

HBsAg-negative patients may show various mutations and amino acid replacement sites at high frequency in the HBV S-region, and these mutations may lead to undetectable Hepatitis B surface antigen (HBsAg), HBsAg antigenic changes or secretion inhibition.

Characterizing glucokinase variant mechanisms using a multiplexed abundance assay

BACKGROUND

Amino acid substitutions can perturb protein activity in multiple ways. Understanding their mechanistic basis may pinpoint how residues contribute to protein function. Here, we characterize the mechanisms underlying variant effects in human glucokinase (GCK) variants, building on our previous comprehensive study on GCK variant activity.

RESULTS

Using a yeast growth-based assay, we score the abundance of 95% of GCK missense and nonsense variants. When combining the abundance scores with our previously determined activity scores, we find that 43% of hypoactive variants also decrease cellular protein abundance. The low-abundance variants are enriched in the large domain, while residues in the small domain are tolerant to mutations with respect to abundance. Instead, many variants in the small domain perturb GCK conformational dynamics which are essential for appropriate activity.

CONCLUSIONS

In this study, we identify residues important for GCK metabolic stability and conformational dynamics. These residues could be targeted to modulate GCK activity, and thereby affect glucose homeostasis.

The Impact of SNP-Induced Amino Acid Substitutions L19P and G66R in the dRP-Lyase Domain of Human DNA Polymerase β on Enzyme Activities

Base excision repair (BER), which involves the sequential activity of DNA glycosylases, apurinic/apyrimidinic endonucleases, DNA polymerases, and DNA ligases, is one of the enzymatic systems that preserve the integrity of the genome. Normal BER is effective, but due to single-nucleotide polymorphisms (SNPs), the enzymes themselves-whose main function is to identify and eliminate damaged bases-can undergo amino acid changes. One of the enzymes in BER is DNA polymerase β (Polβ), whose function is to fill gaps in DNA. SNPs can significantly affect the catalytic activity of an enzyme by causing an amino acid substitution. In this work, pre-steady-state kinetic analyses and molecular dynamics simulations were used to examine the activity of naturally occurring variants of Polβ that have the substitutions L19P and G66R in the dRP-lyase domain. Despite the substantial distance between the dRP-lyase domain and the nucleotidyltransferase active site, it was found that the capacity to form a complex with DNA and with an incoming dNTP is significantly altered by these substitutions. Therefore, the lower activity of the tested polymorphic variants may be associated with a greater number of unrepaired DNA lesions.

Deep mutational scanning of the RNase III-like domain in Trypanosoma brucei RNA editing protein KREPB4

Kinetoplastid pathogens including Trypanosoma brucei, T. cruzi, and Leishmania species, are early diverged, eukaryotic, unicellular parasites. Functional understanding of many proteins from these pathogens has been hampered by limited sequence homology to proteins from other model organisms. Here we describe the development of a high-throughput deep mutational scanning approach in T. brucei that facilitates rapid and unbiased assessment of the impacts of many possible amino acid substitutions within a protein on cell fitness, as measured by relative cell growth. The approach leverages several molecular technologies: cells with conditional expression of a wild-type gene of interest and constitutive expression of a library of mutant variants, degron-controlled stabilization of I-SceI meganuclease to mediate highly efficient transfection of a mutant allele library, and a high-throughput sequencing readout for cell growth upon conditional knockdown of wild-type gene expression and exclusive expression of mutant variants. Using this method, we queried the effects of amino acid substitutions in the apparently non-catalytic RNase III-like domain of KREPB4 (B4), which is an essential component of the RNA Editing Catalytic Complexes (RECCs) that carry out mitochondrial RNA editing in T. brucei. We measured the impacts of thousands of B4 variants on bloodstream form cell growth and validated the most deleterious variants containing single amino acid substitutions. Crucially, there was no correlation between phenotypes and amino acid conservation, demonstrating the greater power of this method over traditional sequence homology searching to identify functional residues. The bloodstream form cell growth phenotypes were combined with structural modeling, RECC protein proximity data, and analysis of selected substitutions in procyclic form T. brucei. These analyses revealed that the B4 RNaseIII-like domain is essential for maintenance of RECC integrity and RECC protein abundances and is also involved in changes in RECCs that occur between bloodstream and procyclic form life cycle stages.

Fit for Purpose Approach To Evaluate Detection of Amino Acid Substitutions in Shotgun Proteomics

Amino acid substitutions (AASs) alter proteins from their genome-expected sequences. Accumulation of substitutions in proteins underlies numerous diseases and antibiotic mechanisms. Accurate global detection of AASs and their frequencies is crucial for understanding these mechanisms. Shotgun proteomics provides an untargeted method for measuring AASs but introduces biases when extrapolating from the genome to identify AASs. To characterize these biases, we created a "ground-truth" approach using the similarities betweenEscherichia coli and Salmonella typhimurium to model the complexity of AAS detection. Shotgun proteomics on mixed lysates generated libraries representing ∼100,000 peptide-spectra and 4161 peptide sequences with a single AAS and defined stoichiometry. Identifying S. typhimurium peptide-spectra with only the E. coli genome resulted in 64.1% correctly identified library peptides. Specific AASs exhibit variable identification efficiencies. There was no inherent bias from the stoichiometry of the substitutions. Short peptides and AASs localized near peptide termini had poor identification efficiency. We identify a new class of "scissor substitutions" that gain or lose protease cleavage sites. Scissor substitutions also had poor identification efficiency. This ground-truth AAS library reveals various sources of bias, which will guide the application of shotgun proteomics to validate AAS hypotheses.

Clinical isolates of Streptococcus mitis/oralis-related species with reduced carbapenem susceptibility, harboring amino acid substitutions in penicillin-binding proteins in Japan

Streptococcus mitis/oralis group isolates with reduced carbapenem susceptibility have been reported, but its isolation rate in Japan is unknown. We collected 356 clinical α-hemolytic streptococcal isolates and identified 142 of them as S. mitis/oralis using partial sodA sequencing. The rate of meropenem non-susceptibility was 17.6% (25/142). All 25 carbapenem-non-susceptible isolates harbored amino acid substitutions in/near the conserved motifs in PBP1A, PBP2B, and PBP2X. Carbapenem non-susceptibility is common among S. mitis/oralis group isolates in Japan.

Construction of an infectious clone for enterovirus A89 and mutagenesis analysis of viral infection and cell binding

Enterovirus A89 (EV-A89) is an unconventional strain belonging to the Enterovirus A species. Limited research has been conducted on EV-A89, leaving its biological and pathogenic properties unclear. Developing reverse genetic tools for EV-A89 would help to unravel its infection mechanisms and aid in the development of vaccines and anti-viral drugs. In this study, an infectious clone for EV-A89 was successfully constructed and recombinant enterovirus A89 (rEV-A89) was generated. The rEV-A89 exhibited similar characteristics such as growth curve, plaque morphology, and dsRNA expression with parental strain. Four amino acid substitutions were identified in the EV-A89 capsid, which were found to enhance viral infection. Mechanistic studies revealed that these substitutions increased the virus's cell-binding ability. Establishing reverse genetic tools for EV-A89 will significantly contribute to understanding viral infection and developing anti-viral strategies.IMPORTANCEEnterovirus A species contain many human pathogens and have been classified into conventional cluster and unconventional cluster. Most of the research focuses on various conventional members, while understanding of the life cycle and infection characteristics of unconventional viruses is still very limited. In our study, we constructed the infectious cDNA clone and single-round infectious particles for the unconventional EV-A89, allowing us to investigate the biological properties of recombinant viruses. Moreover, we identified key amino acids residues that facilitate EV-A89 infection and elucidate their roles in enhancing viral binding to host cells. The establishment of the reverse genetics system will greatly facilitate future study on the life cycle of EV-A89 and contribute to the development of prophylactic vaccines and anti-viral drugs.

Resistance of the Ginseng Gray Mold Pathogen, Botrytis cinerea, to Boscalid and Pyraclostrobin Fungicides in China

Gray mold caused by Botrytis cinerea severely threatens the yield of ginseng (Panax ginseng). Various categories of fungicides have been utilized to control gray mold on this crop. In this study, the resistance of 102 isolates of B. cinerea from 11 commercial ginseng-growing regions in China to fungicides was examined. A total of 32.4% were resistant to boscalid, with EC50 values that ranged from 12.26 to 235.87 μg/ml, and 94.1% were resistant to pyraclostrobin, with EC50 values that ranged from 5.88 to 487.72 μg/ml. Except for sdhA and sdhD, the amino acid substitutions of P225F, P225L, N230I, H272Y, and H272R in the sdhB subunit from 24 (4 sensitive [S] and 20 resistant [R]), 5 (1 S and 4 R), 1 (S), 1 (R), and 8 (4 S and 4 R) strains, respectively, and the concurrent amino acid substitutions of G85A + I93V + M158V + V168I in the sdhC subunit from 5 (4 S and 1 R) strains were identified. A G143A substitution in cytochrome b was identified in 96 isolates that were resistant to pyraclostrobin and three that were sensitive to it. The Bcbi-143/144 intron was identified in the other three isolates sensitive to pyraclostrobin, but it was absent in the isolates that harbored the G143A mutation. The results showed that the populations of B. cinerea on ginseng have developed strong resistance to pyraclostrobin. Therefore, it is not recommended to continue using this fungicide to control gray mold on P. ginseng. Boscalid is still effective against most isolates. However, to prevent fungicide resistance, it is recommended to use a mixture of boscalid with other categories of fungicides.

Recognition of the HLA-C*07:359 allele in Taiwanese individuals

One nucleotide substitution in codon 264 of HLA-C*07:02:01:01 results in a novel allele, HLA-C*07:359.

Cell Culture Adaptive Amino Acid Substitutions in FMDV Structural Proteins: A Key Mechanism for Altered Receptor Tropism

The foot-and-mouth disease virus is a highly contagious and economically devastating virus of cloven-hooved animals, including cattle, buffalo, sheep, and goats, causing reduced animal productivity and posing international trade restrictions. For decades, chemically inactivated vaccines have been serving as the most effective strategy for the management of foot-and-mouth disease. Inactivated vaccines are commercially produced in cell culture systems, which require successful propagation and adaptation of field isolates, demanding a high cost and laborious time. Cell culture adaptation is chiefly indebted to amino acid substitutions in surface-exposed capsid proteins, altering the necessity of RGD-dependent receptors to heparan sulfate macromolecules for virus binding. Several amino acid substations in VP1, VP2, and VP3 capsid proteins of FMDV, both at structural and functional levels, have been characterized previously. This literature review combines frequently reported amino acid substitutions in virus capsid proteins, their critical roles in virus adaptation, and functional characterization of the substitutions. Furthermore, this data can facilitate molecular virologists to develop new vaccine strains against the foot-and-mouth disease virus, revolutionizing vaccinology via reverse genetic engineering and synthetic biology.

Influence of Amino Acid Substitutions in Capsid Proteins of Coxsackievirus B5 on Free Chlorine and Thermal Inactivation

The sensitivity of enteroviruses to disinfectants varies among genetically similar variants and coincides with amino acid changes in capsid proteins, although the effect of individual substitutions remains unknown. Here, we employed reverse genetics to investigate how amino acid substitutions in coxsackievirus B5 (CVB5) capsid proteins affect the virus' sensitivity to free chlorine and heat treatment. Of ten amino acid changes observed in CVB5 variants with free chlorine resistance, none significantly reduced the chlorine sensitivity, indicating a minor role of the capsid composition in chlorine sensitivity of CVB5. Conversely, a subset of these amino acid changes located at the C-terminal region of viral protein 1 led to reduced heat sensitivity. Cryo-electron microscopy revealed that these changes affect the assembly of intermediate viral states (altered and empty particles), suggesting that the mechanism for reduced heat sensitivity could be related to improved molecular packing of CVB5, resulting in greater stability or altered dynamics of virus uncoating during infection.

Network of epistatic interactions in an enzyme active site revealed by large-scale deep mutational scanning

Cooperative interactions between amino acids are critical for protein function. A genetic reflection of cooperativity is epistasis, which is when a change in the amino acid at one position changes the sequence requirements at another position. To assess epistasis within an enzyme active site, we utilized CTX-M β-lactamase as a model system. CTX-M hydrolyzes β-lactam antibiotics to provide antibiotic resistance, allowing a simple functional selection for rapid sorting of modified enzymes. We created all pairwise mutations across 17 active site positions in the β-lactamase enzyme and quantitated the function of variants against two β-lactam antibiotics using next-generation sequencing. Context-dependent sequence requirements were determined by comparing the antibiotic resistance function of double mutations across the CTX-M active site to their predicted function based on the constituent single mutations, revealing both positive epistasis (synergistic interactions) and negative epistasis (antagonistic interactions) between amino acid substitutions. The resulting trends demonstrate that positive epistasis is present throughout the active site, that epistasis between residues is mediated through substrate interactions, and that residues more tolerant to substitutions serve as generic compensators which are responsible for many cases of positive epistasis. Additionally, we show that a key catalytic residue (Glu166) is amenable to compensatory mutations, and we characterize one such double mutant (E166Y/N170G) that acts by an altered catalytic mechanism. These findings shed light on the unique biochemical factors that drive epistasis within an enzyme active site and will inform enzyme engineering efforts by bridging the gap between amino acid sequence and catalytic function.

The Structural Role of N170 in Substrate-Assisted Deacylation in KPC-2 β-Lactamase

The amino acid substitutions in Klebsiella pneumoniae carbapenemase 2 (KPC-2) that have arisen in the clinic are observed to lead to the development of resistance to ceftazidime-avibactam, a preferred treatment for KPC bearing Gram-negative bacteria. Specific substitutions in the omega loop (R164-D179) result in changes in the structure and function of the enzyme, leading to alterations in substrate specificity, decreased stability, and more recently observed, increased resistance to ceftazidime/avibactam. Using accelerated rare-event sampling well-tempered metadynamics simulations, we explored in detail the structural role of R164 and D179 variants that are described to confer ceftazidime/avibactam resistance. The buried conformation of D179 substitutions produce a pronounced structural disorder in the omega loop - more than R164 mutants, where the crystallographic omega loop structure remains mostly intact. Our findings also reveal that the conformation of N170 plays an underappreciated role impacting drug binding and restricting deacylation. The results further support the hypothesis that KPC-2 D179 variants employ substrate-assisted catalysis for ceftazidime hydrolysis, involving the ring amine of the aminothiazole group to promote deacylation and catalytic turnover. Moreover, the shift in the WT conformation of N170 contributes to reduced deacylation and an altered spectrum of enzymatic activity.

Unravelling the link between SARS-CoV-2 mutation frequencies, patient comorbidities, and structural dynamics

Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.

Enhancing prime editor activity by directed protein evolution in yeast

Prime editing is a highly versatile genome editing technology that enables the introduction of base substitutions, insertions, and deletions. However, compared to traditional Cas9 nucleases prime editors (PEs) are less active. In this study we use OrthoRep, a yeast-based platform for directed protein evolution, to enhance the editing efficiency of PEs. After several rounds of evolution with increased selection pressure, we identify multiple mutations that have a positive effect on PE activity in yeast cells and in biochemical assays. Combining the two most effective mutations - the A259D amino acid substitution in nCas9 and the K445T substitution in M-MLV RT - results in the variant PE_Y18. Delivery of PE_Y18, encoded on DNA, mRNA or as a ribonucleoprotein complex into mammalian cell lines increases editing rates up to 3.5-fold compared to PEmax. In addition, PE_Y18 supports higher prime editing rates when delivered in vivo into the liver or brain. Our study demonstrates proof-of-concept for the application of OrthoRep to optimize genome editing tools in eukaryotic cells.

Structure-guided engineering of biased-agonism in the human niacin receptor via single amino acid substitution

The Hydroxycarboxylic acid receptor 2 (HCA2), also known as the niacin receptor or GPR109A, is a prototypical GPCR that plays a central role in the inhibition of lipolytic and atherogenic activities. Its activation also results in vasodilation that is linked to the side-effect of flushing associated with dyslipidemia drugs such as niacin. GPR109A continues to be a target for developing potential therapeutics in dyslipidemia with minimized flushing response. Here, we present cryo-EM structures of the GPR109A in complex with dyslipidemia drugs, niacin or acipimox, non-flushing agonists, MK6892 or GSK256073, and recently approved psoriasis drug, monomethyl fumarate (MMF). These structures elucidate the binding mechanism of agonists, molecular basis of receptor activation, and insights into biased signaling elicited by some of the agonists. The structural framework also allows us to engineer receptor mutants that exhibit G-protein signaling bias, and therefore, our study may help in structure-guided drug discovery efforts targeting this receptor.

Detection of CTX-M-15 ESBL in XDR Haemophilus parainfluenzae from a urethral swab

OBJECTIVES

Haemophilus parainfluenzae is an opportunistic pathogen causing respiratory tract infection and sexually transmitted diseases. The emergence of multidrug resistance in this species is particularly worrisome, especially since the recent description of CTX-M-15 ESBL-producing isolates in Spain. The aim of this study was to characterize a CTX-M-15-producing H. parainfluenzae clinical isolate, HP01, obtained from a urethral swab.

METHODS

MICs were determined with gradient strips for this isolate. Hydrolysis assays were performed with the β LACTA test. Genomic DNA from HP01 was subjected to Illumina and Oxford Nanopore sequencing to investigate the genetic environment of blaCTX-M-15. Phylogenetic analysis was performed with available H. parainfluenzae genomes from the NCBI database, including CTX-M-15 producers.

RESULTS

HP01, an XDR isolate, was resistant to penicillin, third-generation cephalosporins, fluoroquinolones, macrolides, cyclines and co-trimoxazole and susceptible only to carbapenems and rifampicin. HP01 carried blaTEM-1, blaCTX-M-15, tet(M), catS and mef(E)/mel and harboured amino acid substitutions in PBP3, PBP5, GyrA, ParC and FolA implicated in resistance. Genomic analysis revealed that blaCTX-M-15 was carried by a Tn3-like transposon inserted into a novel integrative and conjugative element (ICE), ICEHpaSLS, present on the chromosome and belonging to the ICEHin1056 family described in Haemophilus influenzae. The tet(M)-MEGA element was also detected on the chromosome. No plasmid was found. The phylogenetic analysis showed that four H. parainfluenzae producing CTX-M-15 clustered in the same clade.

CONCLUSIONS

Here we report the description of an XDR H. parainfluenzae producing blaCTX-M-15 isolated from a urethral swab. The blaCTX-M-15 gene was inserted into an ICE structure similar to those recently described in CTX-M-15 producers in Spain. The emergence of XDR H. parainfluenzae producing blaCTX-M-15 is a matter of great concern. Careful surveillance is required to prevent its spread.

Humanization and functional characterization of enhanced coagulation factor IX variants identified through ancestral sequence reconstruction

BACKGROUND

Laboratory resurrection of ancient coagulation factor (F) IX variants generated through ancestral sequence reconstruction led to the discovery of a FIX variant, designated An96, which possesses enhanced specific activity independent of and additive to that provided by human p.Arg384Lys, referred to as FIX-Padua.

OBJECTIVES

The goal of the current study was to identify the amino acid substitution(s) responsible for the enhanced activity of An96 and create a humanized An96 FIX transgene for gene therapy application.

METHODS

Reductionist screening approaches, including domain swapping and scanning residue substitution, were used and guided by one-stage FIX activity assays. In vitro characterization of top candidates included recombinant high-purity preparation, specific activity determination, and enzyme kinetic analysis. Final candidates were packaged into adeno-associated viral (AAV) vectors and delivered to hemophilia B mice.

RESULTS

Five of 42 total amino acid substitutions in An96 appear sufficient to retain the enhanced activity of An96 in an otherwise human FIX variant. Additional substitution of the Padua variant further increased the specific activity 5-fold. This candidate, designated ET9, demonstrated 51-fold greater specific activity than hFIX. AAV2/8-ET9 treated hemophilia B mice produced plasma FIX activities equivalent to those observed previously for AAV2/8-An96-Padua, which were 10-fold higher than AAV2/8-hFIX-Padua.

CONCLUSION

Starting from computationally inferred ancient FIX sequences, novel amino acid substitutions conferring activity enhancement were identified and translated into an AAV-FIX gene therapy cassette demonstrating high potency. This ancestral sequence reconstruction discovery and sequence mapping refinement approach represents a promising platform for broader protein drug and gene therapy candidate optimization.

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

Channel Gating in Kalium Channelrhodopsin Slow Mutants

Kalium channelrhodopsin 1 from Hyphochytrium catenoides (HcKCR1) is the first discovered natural light-gated ion channel that shows higher selectivity to K+ than to Na+ and therefore is used to silence neurons with light (optogenetics). Replacement of the conserved cysteine residue in the transmembrane helix 3 (Cys110) with alanine or threonine results in a >1,000-fold decrease in the channel closing rate. The phenotype of the corresponding mutants in channelrhodopsin 2 is attributed to breaking of a specific interhelical hydrogen bond (the "DC gate"). Unlike CrChR2 and other ChRs with long distance "DC gates", the HcKCR1 structure does not reveal any hydrogen bonding partners to Cys110, indicating that the mutant phenotype is likely caused by disruption of direct interaction between this residue and the chromophore. In HcKCR1_C110A, fast photochemical conversions corresponding to channel gating were followed by dramatically slower absorption changes. Full recovery of the unphotolyzed state in HcKCR1_C110A was extremely slow with two time constants 5.2 and 70 min. Analysis of the light-minus-dark difference spectra during these slow processes revealed accumulation of at least four spectrally distinct blue light-absorbing photocycle intermediates, L, M1 and M2, and a UV light-absorbing form, typical of bacteriorhodopsin-like channelrhodopsins from cryptophytes. Our results contribute to better understanding of the mechanistic links between the chromophore photochemistry and channel conductance, and provide the basis for using HcKCR1_C110A as an optogenetic tool.

Single-molecule force stability of the SARS-CoV-2-ACE2 interface in variants-of-concern

Mutations in SARS-CoV-2 have shown effective evasion of population immunity and increased affinity to the cellular receptor angiotensin-converting enzyme 2 (ACE2). However, in the dynamic environment of the respiratory tract, forces act on the binding partners, which raises the question of whether not only affinity but also force stability of the SARS-CoV-2-ACE2 interaction might be a selection factor for mutations. Using magnetic tweezers, we investigate the impact of amino acid substitutions in variants of concern (Alpha, Beta, Gamma and Delta) and on force-stability and bond kinetic of the receptor-binding domain-ACE2 interface at a single-molecule resolution. We find a higher affinity for all of the variants of concern (>fivefold) compared with the wild type. In contrast, Alpha is the only variant of concern that shows higher force stability (by 17%) compared with the wild type. Using molecular dynamics simulations, we rationalize the mechanistic molecular origins of this increase in force stability. Our study emphasizes the diversity of contributions to the transmissibility of variants and establishes force stability as one of the several factors for fitness. Understanding fitness advantages opens the possibility for the prediction of probable mutations, allowing a rapid adjustment of therapeutics, vaccines and intervention measures.

A New δ-Globin Gene Variant: Hb A2-Yulin [δ46(CD5)Gly→Arg,HBD: C.139G > A]

We report a new δ-chain hemoglobin (Hb) variant observed in a 5-year-old female living in Yulin, Guangxi, China. Capillary electrophoresis revealed splitting of the Hb A2 peak into two fractions (Hb A2 and Hb A2 variant), and the Hb A2 variant was also detected by high-performance liquid chromatography. However, it could not be detected using matrix-assisted laser desorption lonization-time of flight mass spectrometry. CD41-42 (-TCTT) heterozygosity was observed on the HBB gene by PCR and reverse dot-blot hybridization. Sanger sequencing showed a new transition (G > A) at codon 46 of the HBD gene, resulting in glycine changing to arginine. Based on the patient's place of residence, the new variant was named Hb A2-Yulin [δ46(CD5)Gly→Arg,HBD:c.139G > A].

A New α1-Globin Variant, Hb Ormylia [HBA1:c.63C > G; p.His21Gln]. Report of Eleven Cases in Northern Greece

The first identification of a novel α1-Globin variant, Hb Ormylia in 11 Greeks originating from a small village, Ormylia, Chalkidiki, Greece is reported. The new genetic variant leads to the production of a hemoglobin variant that can be identified and quantified by High-Performance Liquid Chromatography. Capillary and classic electrophoresis were not informative. Direct DNA sequencing revealed a new mutation C > G mutation at codon 21 of α1 gene (His > Gln). The new variant has been named Hb Ormylia and this is the first description of this genetic variant of α1 gene in the literature.

AAV-mediated hepatic expression of SLC30A10 and the Thr95Ile variant attenuates manganese excess and other phenotypes in Slc30a10-deficient mice

The manganese (Mn) export protein SLC30A10 is essential for Mn excretion via the liver and intestines. Patients with SLC30A10 deficiency develop Mn excess, dystonia, liver disease, and polycythemia. Recent genome-wide association studies revealed a link between the SLC30A10 variant T95I and markers of liver disease. The in vivo relevance of this variant has yet to be investigated. Using in vitro and in vivo models, we explore the impact of the T95I variant on SLC30A10 function. While SLC30A10 I95 expressed at lower levels than T95 in transfected cell lines, both T95 and I95 variants protected cells similarly from Mn-induced toxicity. Adeno-associated virus 8-mediated expression of T95 or I95 SLC30A10 using the liver-specific thyroxine binding globulin promoter normalized liver Mn levels in mice with hepatocyte Slc30a10 deficiency. Furthermore, Adeno-associated virus-mediated expression of T95 or I95 SLC30A10 normalized red blood cell parameters and body weights and attenuated Mn levels and differential gene expression in livers and brains of mice with whole body Slc30a10 deficiency. While our in vivo data do not indicate that the T95I variant significantly compromises SLC30A10 function, it does reinforce the notion that the liver is a key site of SLC30A10 function. It also supports the idea that restoration of hepatic SLC30A10 expression is sufficient to attenuate phenotypes in SLC30A10 deficiency.

Impact of amino acid substitutions in hepatitis C virus core region on the severe oxidative stress

Hepatitis C virus (HCV) infection is a major cause of chronic liver disease, leading to liver steatosis, fibrosis, and hepatocellular carcinoma (HCC). Despite the accumulation of clinical data showing the impact of amino acid substitutions at positions 70 (R70Q/H) and/or 91 (L91M) in the HCV core protein in progressive liver diseases, including HCC, the underlying mechanisms have not been elucidated. We analyzed 72 liver biopsy specimens from patients with chronic HCV genotype 1b (HCV-1b) infection prior to antiviral treatment. Levels of 8-hydroxy-2'-deoxyguanosine (8-OHdG) and nuclear factor erythroid 2-related factor 2 (NRF2) in the nucleus were quantified using liver tissue immunohistochemistry. The effects of amino acid substitutions in the HCV core region on hepatocellular oxidative stress were investigated using wild-type or double-mutant (R70Q/H+L91M) HCV-1b core transfection and stable expression in human hepatoma HuH-7 cells. Overall, 24, 19, 11, and 18 patients had the wild-type, R70Q/H, L91M, and R70Q/H+L91M genotypes, respectively, in the HCV core. A significantly higher accumulation of hepatocellular 8-OHdG and a lower NRF2/8-OHdG ratio were observed in patients with R70Q/H+L91M than in those with the wild-type disease. Increased levels of intracellular superoxide and hydrogen peroxide in the cytoplasm and mitochondria, mRNA expression of enzymes generating oxidative stress, and nuclear expression of nicotinamide adenine dinucleotide phosphate oxidase 4 were augmented in cells treated with R70Q+L91M. HCV core proteins harboring either or both substitutions of R70Q/H or L91M enhanced hepatocellular oxidative stress in vivo and in vitro. These amino acid substitutions may affect HCC development by enhancing hepatic oxidative stress in patients with chronic HCV-1b infection.

Estimating amino acid substitution models from genome datasets: a simulation study on the performance of estimated models

Estimating parameters of amino acid substitution models is a crucial task in bioinformatics. The maximum likelihood (ML) approach has been proposed to estimate amino acid substitution models from large datasets. The quality of newly estimated models is normally assessed by comparing with the existing models in building ML trees. Two important questions remained are the correlation of the estimated models with the true models and the required size of the training datasets to estimate reliable models. In this article, we performed a simulation study to answer these two questions based on simulated data. We simulated genome datasets with different numbers of genes/alignments based on predefined models (called true models) and predefined trees (called true trees). The simulated datasets were used to estimate amino acid substitution model using the ML estimation methods. Our experiments showed that models estimated by the ML methods from simulated datasets with more than 100 genes have high correlations with the true models. The estimated models performed well in building ML trees in comparison with the true models. The results suggest that amino acid substitution models estimated by the ML methods from large genome datasets are a reliable tool for analyzing amino acid sequences.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

Substitution Models of Protein Evolution with Selection on Enzymatic Activity

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

Prevalence of persistent SARS-CoV-2 in a large community surveillance study

Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as 'persistent infections' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

Twenty natural amino acid substitution screening at the last residue 121 of influenza A virus NS2 protein reveals the critical role of NS2 in promoting virus genome replication by coordinating with viral polymerase

Both influenza A virus genome transcription (vRNA→mRNA) and replication (vRNA→cRNA→vRNA), catalyzed by the influenza RNA polymerase (FluPol), are dynamically regulated across the virus life cycle. It has been reported that the last amino acid I121 of the viral NS2 protein plays a critical role in promoting viral genome replication in influenza mini-replicon systems. Here, we performed a 20 natural amino acid substitution screening at residue NS2-I121 in the context of virus infection. We found that the hydrophobicity of the residue 121 is essential for virus survival. Interestingly, through serial passage of the rescued mutant viruses, we further identified adaptive mutations PA-K19E and PB1-S713N on FluPol which could effectively compensate for the replication-promoting defect caused by NS2-I121 mutation in the both mini-replicon and virus infection systems. Structural analysis of different functional states of FluPol indicates that PA-K19E and PB1-S713N could stabilize the replicase conformation of FluPol. By using a cell-based NanoBiT complementary reporter assay, we further demonstrate that both wild-type NS2 and PA-K19E/PB1-S713N could enhance FluPol dimerization, which is necessary for genome replication. These results reveal the critical role NS2 plays in promoting viral genome replication by coordinating with FluPol.IMPORTANCEThe intrinsic mechanisms of influenza RNA polymerase (FluPol) in catalyzing viral genome transcription and replication have been largely resolved. However, the mechanisms of how transcription and replication are dynamically regulated remain elusive. We recently reported that the last amino acid of the viral NS2 protein plays a critical role in promoting viral genome replication in an influenza mini-replicon system. Here, we conducted a 20 amino acid substitution screening at the last residue 121 in virus rescue and serial passage. Our results demonstrate that the replication-promoting function of NS2 is important for virus survival and efficient multiplication. We further show evidence that NS2 and NS2-I121 adaptive mutations PA-K19E/PB1-S713N regulate virus genome replication by promoting FluPol dimerization. This work highlights the coordination between NS2 and FluPol in fulfilling efficient genome replication. It further advances our understanding of the regulation of viral RNA synthesis for influenza A virus.

A naturally occurring HA-stabilizing amino acid (HA1-Y17) in an A(H9N2) low-pathogenic influenza virus contributes to airborne transmission

IMPORTANCE

Despite the accumulation of evidence showing that airborne transmissible influenza A virus (IAV) typically has a lower pH threshold for hemagglutinin (HA) fusion activation, the underlying mechanism for such a link remains unclear. In our study, by using a pair of isogenic recombinant A(H9N2) viruses with a phenotypical difference in virus airborne transmission in a ferret model due to an acid-destabilizing mutation (HA1-Y17H) in the HA, we demonstrate that an acid-stable A(H9N2) virus possesses a multitude of advantages over its less stable counterpart, including better fitness in the ferret respiratory tract, more effective aerosol emission from infected animals, and improved host susceptibility. Our study provides supporting evidence for the requirement of acid stability in efficient airborne transmission of IAV and sheds light on fundamental mechanisms for virus airborne transmission.

The impact of PA/I38 substitutions and PA polymorphisms on the susceptibility of zoonotic influenza A viruses to baloxavir

Genetic reassortment of avian, swine, and human influenza A viruses (IAVs) poses potential pandemic risks. Surveillance is important for influenza pandemic preparedness, but the susceptibility of zoonotic IAVs to the cap-dependent endonuclease inhibitor baloxavir acid (BXA) has not been thoroughly researched. Although an amino acid substitution at position 38 in the polymerase acidic protein (PA/I38) in seasonal IAVs reduces BXA susceptibility, PA polymorphisms at position 38 are rarely seen in zoonotic IAVs. Here, we examined the impact of PA/I38 substitutions on the BXA susceptibility of recombinant A(H5N1) viruses. PA mutants that harbored I38T, F, and M were 48.2-, 24.0-, and 15.5-fold less susceptible, respectively, to BXA than wild-type A(H5N1) but were susceptible to the neuraminidase inhibitor oseltamivir acid and the RNA polymerase inhibitor favipiravir. PA mutants exhibited significantly impaired replicative fitness in Madin-Darby canine kidney cells at 24 h postinfection. In addition, in order to investigate new genetic markers for BXA susceptibility, we screened geographically and temporally distinct IAVs isolated worldwide from birds and pigs. The results showed that BXA exhibited antiviral activity against avian and swine viruses with similar levels to seasonal isolates. All viruses tested in the study lacked the PA/I38 substitution and were susceptible to BXA. Isolates harboring amino acid polymorphisms at positions 20, 24, and 37, which have been implicated in the binding of BXA to the PA endonuclease domain, were also susceptible to BXA. These results suggest that monitoring of the PA/I38 substitution in animal-derived influenza viruses is important for preparedness against zoonotic influenza virus outbreaks.

Isolation, characterization, and phylogenetic analysis of two new porcine parvovirus 1 isolates from Northern China

Porcine parvovirus (PPV) is a pathogen of infectious reproductive disease, which can cause stillbirth, mummification, embryo death, and infertility (SMEDI) syndrome in pigs. The objective of this study was to gain new insights into the evolution and phylogeny of the PPV1 genome. In this study, we isolated two new PPV1 (HLJ202108-Y and SDLC202109) from northern China and sequenced their whole genomes. The new isolates were found to have three amino acid substitutions (K195R, K562R, and S578P) in nonstructural protein 1. The VP2 amino acid site contained nine nonsynonymous substitutions, including six substitutions of the Kresse strain corresponding to the NADL-2 strain and three substitutions of A414S, S436T, and N555K. Genetic evolution analysis was conducted on 107 reference sequences available in the GenBank database, and 4-5 PPV1 taxa were defined. The new isolates were in the same phylogenetic cluster as strain 27a. The changes in the cluster, specifically marker amino acids, and their potential role in enhancing pathogenicity are discussed in this study. Furthermore, the evolutionary tree map results showed that the strains in China were evolving in two directions: one was becoming increasingly similar to early NADL-2 strains, while the other was evolving toward 27a-like strains. We also compared the proliferation ability of the isolated strains in susceptible cells by analyzing the multistep growth curves. The results showed that the virulence titer of the mutant strain was high. In summary, this study introduced the latest changes in PPV and discussed the virus characteristics that were considered to affect virulence.

Hb Hebei [α20 (B1) His→Leu; HBA2:C.62A > T]: A Novel Hemoglobin Variant Found during Measurement of Glycated Hemoglobin

We report a novel hemoglobin (Hb) variant found in a 34-year-old Chinese male during a routine measurement of glycated hemoglobin. The variant resulted in a P3 peak of 27.5% of the total Hb on high performance liquid chromatography (HPLC) with a glycated hemoglobin mode. However, no abnormal Hb peaks were observed in capillary electrophoresis (CE) with 3.1% Hb A2 and 96.9% Hb A. The amino acid substitution was determined by Sanger sequencing as α20 (B1) His→Leu; the corresponding DNA mutation was identified as CAC > CTC at the first position of codon 20 of the α-chain. This is the first description of the mutation, and we have named it Hb Hebei for the region of origin of the proband.

The phosphorylation of the movement protein TGBp1 regulates the accumulation of the Bamboo mosaic virus

Phosphorylation and dephosphorylation of viral movement proteins plays a crucial role in regulating virus movement. Our study focused on investigating the movement protein TGBp1 of Bamboo mosaic virus (BaMV), which is a single-stranded positive-sense RNA virus. Specifically, we examined four potential phosphorylation sites (S15, S18, T58, and S247) within the TGBp1 protein. To study the impact of phosphorylation, we introduced amino acid substitutions at the selected sites. Alanine substitutions were used to prevent phosphorylation, while aspartate substitutions were employed to mimic phosphorylation. Our findings suggest that mimicking phosphorylation at S15, S18 and T58 of TGBp1 might be linked to silencing suppressor activities. The phosphorylated form at these sites exhibits a loss of silencing suppressor activity, leading to reduced viral accumulation in the inoculated leaves. Furthermore, mimicking phosphorylation at residues S15 and S18 could diminish viral accumulation at the single-cell level, while doing so at residue T58 could influence virus movement. However, mimicking phosphorylation at residue S247 does not appear to be relevant to both functions of TGBp1. Overall, our study provides insights into the functional significance of specific phosphorylation sites in BaMV TGBp1, illuminating the regulatory mechanisms involved in virus movement and silencing suppression.

Identification of a buried β-strand as a novel disease-related motif in the human polysialyltransferases

The polysialyltransferases ST8SIA2 and ST8SIA4 and their product, polysialic acid (polySia), are known to be related to cancers and mental disorders. ST8SIA2 and ST8SIA4 have conserved amino acid (AA) sequence motifs essential for the synthesis of the polySia structures on the neural cell adhesion molecule. To search for a new motif in the polysialyltransferases, we adopted the in silico Individual Meta Random Forest program that can predict disease-related AA substitutions. The Individual Meta Random Forest program predicted a new eight-amino-acids sequence motif consisting of highly pathogenic AA residues, thus designated as the pathogenic (P) motif. A series of alanine point mutation experiments in the pathogenic motif (P motif) showed that most P motif mutants lost the polysialylation activity without changing the proper enzyme expression levels or localization in the Golgi. In addition, we evaluated the enzyme stability of the P motif mutants using newly established calculations of mutation energy, demonstrating that the subtle change of the conformational energy regulates the activity. In the AlphaFold2 model, we found that the P motif was a buried β-strand underneath the known surface motifs unique to ST8SIA2 and ST8SIA4. Taken together, the P motif is a novel buried β-strand that regulates the full activity of polysialyltransferases from the inside of the molecule.

Structural Evolution of Primate Glutamate Dehydrogenase 2 as Revealed by In Silico Predictions and Experimentally Determined Structures

Glutamate dehydrogenase (GDH) interconverts glutamate to a-ketoglutarate and ammonia, interconnecting amino acid and carbohydrate metabolism. In humans, two functional GDH genes, GLUD1 and GLUD2, encode for hGDH1 and hGDH2, respectively. GLUD2 evolved from retrotransposition of the GLUD1 gene in the common ancestor of modern apes. These two isoenzymes are involved in the pathophysiology of human metabolic, neoplastic, and neurodegenerative disorders. The 3D structures of hGDH1 and hGDH2 have been experimentally determined; however, no information is available about the path of GDH2 structure changes during primate evolution. Here, we compare the structures predicted by the AlphaFold Colab method for the GDH2 enzyme of modern apes and their extinct primate ancestors. Also, we analyze the individual effect of amino acid substitutions emerging during primate evolution. Our most important finding is that the predicted structure of GDH2 in the common ancestor of apes was the steppingstone for the structural evolution of primate GDH2s. Two changes with a strong functional impact occurring at the first evolutionary step, Arg443Ser and Gly456Ala, had a destabilizing and stabilizing effect, respectively, making this step the most important one. Subsequently, GDH2 underwent additional modifications that fine-tuned its enzymatic properties to adapt to the functional needs of modern-day primate tissues.

A Single Amino Acid Substitution in the Transmembrane Domain of Glycoprotein H Functionally Compensates for the Absence of gL in Pseudorabies Virus

Herpesvirus entry requires the coordinated action of at least four viral glycoproteins. Virus-specific binding to a cellular receptor triggers a membrane fusion cascade involving the conserved gH/gL complex and gB. Although gB is the genuine herpesvirus fusogen, it requires gH/gL for fusion, but how activation occurs is still unclear. To study the underlying mechanism, we used a gL-deleted pseudorabies virus (PrV) mutant characterized by its limited capability to directly infect neighboring cells that was exploited for several independent serial passages in cell culture. Unlike previous revertants that acquired mutations in the gL-binding N-terminus of gH, we obtained a variant, PrV-ΔgLPassV99, that unexpectedly contained two amino acid substitutions in the gH transmembrane domain (TMD). One of these mutations, I662S, was sufficient to compensate for gL function in virus entry and in in vitro cell-cell fusion assays in presence of wild type gB, but barely for cell-to-cell spread. Additional expression of receptor-binding PrV gD, which is dispensable for cell-cell fusion mediated by native gB, gH and gL, resulted in hyperfusion in combination with gH V99. Overall, our results uncover a yet-underestimated role of the gH TMD in fusion regulation, further shedding light on the complexity of herpesvirus fusion involving all structural domains of the conserved entry glycoproteins.

Amino acid substitution L232F in non-structural protein 6 identified as a possible human-adaptive mutation in clade B MERS coronaviruses

IMPORTANCE

Viral host adaptation plays an important role in inter-species transmission of coronaviruses and influenza viruses. Multiple human-adaptive mutations have been identified in influenza viruses but not so far in MERS-CoV that circulates widely in dromedary camels in the Arabian Peninsula leading to zoonotic transmission. Here, we analyzed clade B MERS-CoV sequences and identified an amino acid substitution L232F in nsp6 that repeatedly occurs in human MERS-CoV. Using a loss-of-function reverse genetics approach, we found the nsp6 L232F conferred increased viral replication competence in vitro, in cultures of the upper human respiratory tract ex vivo, and in lungs of mice infected in vivo. Our results showed that nsp6 L232F may be an adaptive mutation associated with zoonotic transmission of MERS-CoV. This study highlighted the capacity of MERS-CoV to adapt to transmission to humans and also the need for continued surveillance of MERS-CoV in camels.

Amino acid substitutions in norovirus VP1 dictate host dissemination via variations in cellular attachment

IMPORTANCE

All viruses initiate infection by utilizing receptors to attach to target host cells. These virus-receptor interactions can therefore dictate viral replication and pathogenesis. Understanding the nature of virus-receptor interactions could also be important for the development of novel therapies. Noroviruses are non-enveloped icosahedral viruses of medical importance. They are a common cause of acute gastroenteritis with no approved vaccine or therapy and are a tractable model for studying fundamental virus biology. In this study, we utilized the murine norovirus model system to show that variation in a single amino acid of the major capsid protein alone can affect viral infectivity through improved attachment to suspension cells. Modulating plasma membrane mobility reduced infectivity, suggesting an importance of membrane mobility for receptor recruitment and/or receptor conformation. Furthermore, different substitutions at this site altered viral tissue distribution in a murine model, illustrating how in-host capsid evolution could influence viral infectivity and/or immune evasion.

Probing of the reactive center loop region of alpha-1-antitrypsin by mutagenesis predicts new type-2 dysfunctional variants

Lung disease in alpha-1-antitrypsin deficiency (AATD) mainly results from insufficient control of the serine proteases neutrophil elastase (NE) and proteinase-3 due to reduced plasma levels of alpha-1-antitrypsin (AAT) variants. Mutations in the specificity-determining reactive center loop (RCL) of AAT would be predicted to minimally affect protein folding and secretion by hepatocytes but can impair anti-protease activity or alter the target protease. These properly secreted but dysfunctional 'type-2' variants would not be identified by common diagnostic protocols that are predicated on a reduction in circulating AAT. This has potential clinical relevance: in addition to the dysfunctional Pittsburgh and Iners variants reported previously, several uncharacterized RCL variants are present in genome variation databases. To prospectively evaluate the impact of RCL variations on secretion and anti-protease activity, here we performed a systematic screening of amino acid substitutions occurring at the AAT-NE interface. Twenty-three AAT variants that can result from single nucleotide polymorphisms in this region, including 11 present in sequence variation databases, were expressed in a mammalian cell model. All demonstrated unaltered protein folding and secretion. However, when their ability to form stable complexes with NE was evaluated by western blot, enzymatic assays, and a novel ELISA developed to quantify AAT-NE complexes, substrate-like and NE-binding deficient dysfunctional variants were identified. This emphasizes the ability of the RCL to accommodate inactivating substitutions without impacting the integrity of the native molecule and demonstrates that this class of molecule violates a generally accepted paradigm that equates circulating levels with functional protection of lung tissue.

The atypical sphingolipid SPB 18:1(14Z);O2 is a biomarker for DEGS1 related hypomyelinating leukodystrophy

Sphingolipids (SL) represent a structurally diverse class of lipids that are central to cellular physiology and neuronal development and function. Defects in the sphingolipid metabolism are typically associated with nervous system disorders. The C4-dihydroceramide desaturase (DEGS1) catalyzes the conversion of dihydroceramide to ceramide, the final step in the SL de-novo synthesis. Loss of function mutations in DEGS1 cause a hypomyelinating leukodystrophy, which is associated with increased plasma dihydrosphingolipids (dhSL) and with the formation of an atypical SPB 18:1(14Z);O2 metabolite. Here, we characterize two novel DEGS1 variants of unknown significance (VUS), provide a structural model with a predicted substrate binding site, and propose a regulatory link between DEGS1 and fatty acid desaturase 3 (FADS3). Both VUS involve single amino acid substitutions near the C-terminus within conserved regions of the enzyme. Patient 1 (p.R311K variant) shows severe progressive tetraspasticity, intellectual disability, and epilepsy in combination with brain magnetic resonance imaging (MRI) findings, typical for DEGS1-related leukodystrophy. Patient 2 (p.G270E variant) presents with delayed psychomotor development, oculomotor apraxia, and a normal brain MRI. Plasma from the p.R311K carrier showed a significantly elevated dhSL species and the presence of SPB 18:1(14Z);O2, while the plasma SL profile for the p.G270E variant was not altered. This suggests the p.R331K variant is pathogenic, while the p.G270E appears benign. As an increase in dihydroSL species is also seen in other pathological disorders of the SL metabolism, the SPB 18:1(14Z);O2 seems to be a more specific biomarker to discriminate between pathogenic and benign DEGS1 variants.