Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MUTATIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

Plug & play directed evolution of proteins with gradient-based discrete MCMC

Abstract A long-standing goal of machine-learning-based protein engineering is to accelerate the discovery of novel mutations that improve the function of a known protein. We introduce a sampling framework for evolving proteins in silico that supports mixing and matching a variety of unsupervised models, such as protein language models, and supervised models that predict protein function from sequence. By composing these models, we aim to improve our ability to evaluate unseen mutations and constrain search to regions of sequence space likely to contain functional proteins. Our framework achieves this without any model fine-tuning or re-training by constructing a product of experts distribution directly in discrete protein space. Instead of resorting to brute force search or random sampling, which is typical of classic directed evolution, we introduce a fast Markov chain Monte Carlo sampler that uses gradients to propose promising mutations. We conduct in silico directed evolution experiments on wide fitness landscapes and across a range of different pre-trained unsupervised models, including a 650 M parameter protein language model. Our results demonstrate an ability to efficiently discover variants with high evolutionary likelihood as well as estimated activity multiple mutations away from a wild type protein, suggesting our sampler provides a practical and effective new paradigm for machine-learning-based protein engineering.

59 BASIC BIOLOGICAL SCIENCES↗

Orchestration of ErbB3 signaling through heterointeractions and homointeractions

Members of the ErbB family of receptor tyrosine kinases are capable of both homointeractions and heterointeractions. Because each receptor has a unique set of binding sites for downstream signaling partners and differential catalytic activity, subtle shifts in their combinatorial interplay may have a large effect on signaling outcomes. The overexpression and mutation of ErbB family members are common in numerous human cancers and shift the balance of activation within the signaling network. Here we report the development of a spatial stochastic model that addresses the dynamics of ErbB3 homodimerization and heterodimerization with ErbB2. The model is based on experimental measures for diffusion, dimer off-rates, kinase activity, and dephosphorylation. We also report computational analysis of ErbB3 mutations, generating the prediction that activating mutations in the intracellular and extracellular domains may be subdivided into classes with distinct underlying mechanisms. We show experimental evidence for an ErbB3 gain-of-function point mutation located in the C-lobe asymmetric dimerization interface, which shows enhanced phosphorylation at low ligand dose associated with increased kinase activity.

59 BASIC BIOLOGICAL SCIENCES↗

Promoter deletion in the soybean Compact mutant leads to overexpression of a gene with homology to the C 20 -gibberellin 2-oxidase family

Height is a critical component of plant architecture, significantly affecting crop yield. The genetic basis of this trait in soybean remains unclear. In this study, we report the characterization of the Compact mutant of soybean, which has short internodes. The candidate gene was mapped to chromosome 17, and the interval containing the causative mutation was further delineated using biparental mapping. Whole-genome sequencing of the mutant revealed an 8.7 kb deletion in the promoter of the Glyma.17g145200 gene, which encodes a member of the class III gibberellin (GA) 2-oxidases. The mutation has a dominant effect, likely via increased expression of the GA 2-oxidase transcript observed in green tissue, as a result of the deletion in the promoter of Glyma.17g145200. We further demonstrate that levels of GA precursors are altered in the Compact mutant, supporting a role in GA metabolism, and that the mutant phenotype can be rescued with exogenous GA3. We also determined that overexpression of Glyma.17g145200 in Arabidopsis results in dwarfed plants. Thus, gain of promoter activity in the Compact mutant leads to a short internode phenotype in soybean through altered metabolism of gibberellin precursors. Furthermore, these results provide an example of how structural variation can control an important crop trait and a role for Glyma.17g145200 in soybean architecture, with potential implications for increasing crop yield.

59 BASIC BIOLOGICAL SCIENCES↗

Structural insights into the substrate specificity of the endonuclease activity of the influenza virus cap-snatching mechanism

Abstract The endonuclease activity within the influenza virus cap-snatching process is a proven therapeutic target. The anti-influenza drug baloxavir is highly effective, but is associated with resistance mutations that threaten its clinical efficacy. The endonuclease resides within the N-terminal domain of the PA subunit (PAN) of the influenza RNA dependent RNA polymerase, and we report here complexes of PAN with RNA and DNA oligonucleotides to understand its specificity and the structural basis of baloxavir resistance mutations. The RNA and DNA oligonucleotides bind within the substrate binding groove of PAN in a similar fashion, explaining the ability of the enzyme to cleave both substrates. The individual nucleotides occupy adjacent conserved pockets that flank the two-metal active site. However, the 2′ OH of the RNA ribose moieties engage in additional interactions that appear to optimize the binding and cleavage efficiency for the natural substrate. The major baloxavir resistance mutation at position 38 is at the core of the substrate binding site, but structural studies and modeling suggest that it maintains the necessary virus fitness via compensating interactions with RNA. These studies will facilitate the development of new influenza therapeutics that spatially match the substrate and are less likely to elicit resistance mutations.

59 BASIC BIOLOGICAL SCIENCES↗

Structure and functional determinants of Rad6–Bre1 subunits in the histone H2B ubiquitin-conjugating complex

The conserved complex of the Rad6 E2 ubiquitin-conjugating enzyme and the Bre1 E3 ubiquitin ligase catalyzes histone H2B monoubiquitination (H2Bub1), which regulates chromatin dynamics during transcription and other nuclear processes. Here, we report a crystal structure of Rad6 and the non-RING domain N-terminal region of Bre1, which shows an asymmetric homodimer of Bre1 contacting a conserved loop on the Rad6 ‘backside’. This contact is distant from the Rad6 catalytic site and is the location of mutations that impair telomeric silencing in yeast. Mutational analyses validated the importance of this contact for the Rad6–Bre1 interaction, chromatin-binding dynamics, H2Bub1 formation and gene expression. Moreover, the non-RING N-terminal region of Bre1 is sufficient to confer nucleosome binding ability to Rad6 in vitro. Interestingly, Rad6 P43L protein, an interaction interface mutant and equivalent to a cancer mutation in the human homolog, bound Bre1 5-fold more tightly than native Rad6 in vitro, but showed reduced chromatin association of Bre1 and reduced levels of H2Bub1 in vivo. These surprising observations imply conformational transitions of the Rad6–Bre1 complex during its chromatin-associated functional cycle, and reveal the differential effects of specific disease-relevant mutations on the chromatin-bound and unbound states. Overall, our study provides structural insights into Rad6–Bre1 interaction through a novel interface that is important for their biochemical and biological responses.

59 BASIC BIOLOGICAL SCIENCES↗

New Delhi Metallo‐Beta‐Lactamase Variants NDM‐4 and NDM‐12 from E. coli Clinical Isolates Exhibit Increased Activity and Stability

Metallo‐Beta‐Lactamases (MBLs) confer resistance to carbapenems, cephalosporins, and penicillins in several clinically relevant Gram‐negative bacteria including Acinetobacter, Pseudomonas aeruginosa, and many Enterobacteriaceae. New Delhi Metallo‐Beta‐Lactamases (NDMs) are amongst some of the most worrisome and prevalent MBLs. This family of MBLs is capable of hydrolyzing all generations of bicyclic beta‐lactams. Two new NDM variants, NDM‐4 and ‐12, have recently been identified in E. coli clinical isolates. NDM‐4 contains an M154L substitution while NDM‐12 bears both M154L and G222D mutations. Considering the clinical significance of Gram‐negative species possessing these enzymes, there is a need to characterize the structures and activities of these variants in order to provide a better basis for drug discovery efforts. We have found that NDM‐4 and NDM‐12 possess higher hydrolytic activity compared to NDM‐1. Our differential scanning fluorimetry studies suggest that the M154L mutation confers increased thermal stability and improved affinity for active site catalytic Zn metal ions. The additional G222D mutation of NDM‐12 further improves stability and Zn affinity. M154L, found in the secondary coordination sphere, accounts for changes in stability by acting as a buttress to the loop containing His 122 , a direct Zn coordinating residue. G222D, on the other hand, is a mutation located in a loop just above the active site approximately 10 Å from the active site zinc ions. Together, our crystal structures of these NDM variants, including the first known structure of NDM‐12, biochemical analyses, and biophysical analyses provide a basis for understanding the roles of key residues surrounding the active site within NDM variants. A better understanding of clinically relevant NDM variants is essential in the fight against MBL‐mediated antibiotic resistance. Our structural, biochemical, and biophysical studies will help to inform drug design efforts to develop new MBL inhibitors. Support or Funding Information The authors acknowledge financial support from the US National Institutes of Health (Award No. R01 GM111926 to RCP, DLT, MWC, and WF) and institutional support from Miami University through the Robert H. and Nancy J. Blayney Professorship (to RCP). The Advanced Light Source is supported by the US Department of Energy under contract number DE‐AC03‐76SF00098 at Lawrence Berkeley National Laboratory.

Williams, Cameron↗

Genomic analysis of Klebsiella aerogenes circulating in New Mexico

Klebsiella aerogenes is an opportunistic pathogen and a growing cause of healthcare-associated infections, characterized by multidrug resistance and the emergence of global high-risk clones. However, regional genomic surveillance data remain limited. Here, we sought to characterize the population structure, transmission dynamics and resistance mechanisms of clinical K. aerogenes in Albuquerque, New Mexico. We sequenced 177 clinical isolates collected between 2021 and 2023. We also developed a novel, species-specific PopPUNK database to facilitate rapid, high-resolution typing. The New Mexico K. aerogenes population was diverse but dominated by two global pandemic lineages, ST93 (47.5%) and ST4 (7.9%), which were significantly enriched for the virulence factors yersiniabactin and colibactin. Genomic evidence for recent local transmission was rare, with only four putative transmission pairs identified. The resistome was characterized by intrinsic and adaptive mutations. Nearly all isolates possessed gyrA mutations associated with decreased fluoroquinolone susceptibility. Mutations in the AmpC regulator AmpD and the outer membrane porin Omp36 were common, particularly within the dominant ST93 lineage. These mutations have been associated with increased AmpC-mediated carbapenem resistance. Our findings underscore the critical importance of genomic surveillance to monitor the transmission and evolution of adaptive resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Improved exome prioritization of disease genes through cross-species phenotype comparison

Numerous new disease-gene associations have been identified by whole-exome sequencing studies in the last few years. However, many cases remain unsolved due to the sheer number of candidate variants remaining after common filtering strategies such as removing low quality and common variants and those deemed unlikely to be pathogenic (non-coding, not affecting splicing, synonymous or missense mutations annotated as non-pathogenic by prediction algorithms). The observation that each of our genomes contains about 100 genuine loss of function variants with ~20 genes completely inactivated makes identification of the causative mutation problematic when using these strategies alone. In some cases it may be possible to use multiple affected individuals, linkage data, identity-by-descent inference, de novo heterozygous mutations from trio analysis, or prior knowledge of affected pathways to narrow down to the causative variant. In cases where this is not possible or has proven unsuccessful we propose using the wealth of genotype to phenotype data that already exists from model organism studies to assess the potential impact of these exome variants. Here, we introduce PHenotypic Interpretation of Variants in Exomes (PHIVE), an algorithm that integrates the calculation of phenotype similarity between human diseases and genetically modified mouse models with evaluation of the variants according to allele frequency, pathogenicity and mode of inheritance approaches in our Exomiser tool. Large-scale validation of PHIVE analysis using 100,000 exomes containing known mutations demonstrated a substantial improvement (up to 54.1 fold) over purely variant-based (frequency and pathogenicity) methods with the correct gene recalled as the top hit in up to 83% of samples, corresponding to an area under the ROC curve of over 95%. We conclude that incorporation of phenotype data can play a vital role in translational bioinformatics and propose that exome sequencing projects should systematically capture clinical phenotypes to take advantage of the strategy presented here.

59 BASIC BIOLOGICAL SCIENCES↗

Scaling Out a Combinatorial Algorithm for Discovering Carcinogenic Gene Combinations to Thousands of GPUs

Cancer is a leading cause of death in the US, second only to heart disease. It is primarily a result of a combination of an estimated two-nine genetic mutations (multi-hit combinations). Although a body of research has identified hundreds of cancer-causing genetic mutations, we don’t know the specific combination of mutations responsible for specific instances of cancer for most cancer types. An approximate algorithm for solving the weighted set cover problem was previously adapted to identify combinations of genes with mutations that may be responsible for individual instances of cancer. However, the algorithm’s computational requirement scales exponentially with the number of genes, making it impractical for identifying more than three-hit combinations, even after the algorithm was parallelized and scaled up to a V100 GPU. Since most cancers have been estimated to require more than three hits, we scaled out the algorithm to identify combinations of four or more hits using 1000 nodes (6000 V100 GPUs with ≈48×106 processing cores) on the Summit supercomputer at Oak Ridge National Laboratory. Efficiently scaling out the algorithm required a series of algorithmic innovations and optimizations for balancing an exponentially divergent workload across processors and for minimizing memory latency and inter-node communication. We achieved an average strong scaling efficiency of 90.14% (80.96%–97.96% for 200 to 1000 nodes), compared to a 100 node run, with 84.18% scaling efficiency for 1000 nodes. With experimental validation, the multi-hit combinations identified here could provide further insight into the etiology of different cancer subtypes and provide a rational basis for targeted combination therapy.

Dash, Sajal↗

Plasmidic Expression of nemA and yafC * Increased Resistance of Ethanologenic Escherichia coli LY180 to Nonvolatile Side Products from Dilute Acid Treatment of Sugarcane Bagasse and Artificial Hydrolysate

Hydrolysate-resistant Escherichia coli SL100 was previously isolated from ethanologenic LY180 after sequential transfers in AM1 medium containing a dilute acid hydrolysate of sugarcane bagasse and was used as a source of resistance genes. Many genes that affect tolerance to furfural, the most abundant inhibitor, have been described previously. To identify genes associated with inhibitors other than furfural, plasmid clones were selected in an artificial hydrolysate that had been treated with a vacuum to remove furfural. Two new resistance genes were discovered from Sau3A1 libraries of SL100 genomic DNA: nemA (N-ethylmaleimide reductase) and a putative regulatory gene containing a mutation in the coding region, yafC *. The presence of these mutations in SL100 was confirmed by sequencing. A single mutation was found in the upstream regulatory region of nemR ( nemRA operon) in SL100. This mutation increased nemA activity 20-fold over that of the parent organism (LY180) in AM1 medium without hydrolysate and increased nemA mRNA levels >200-fold. Addition of hydrolysates induced nemA expression (mRNA and activity), in agreement with transcriptional control. NemA activity was stable in cell extracts (9 h, 37°C), eliminating a role for proteinase in regulation. LY180 with a plasmid expressing nemA or yafC * was more resistant to a vacuum-treated sugarcane bagasse hydrolysate and to a vacuum-treated artificial hydrolysate than LY180 with an empty-vector control. Neither gene affected furfural tolerance. The vacuum-treated hydrolysates inhibited the reduction of N -ethylmaleimide by NemA while also serving as substrates. In conclusion, expression of the nemA or yafC * plasmid in LY180 doubled the rate of ethanol production from the vacuum-treated sugarcane bagasse hydrolysate.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic, Genomic, and Transcriptomic Studies of Pyruvate Metabolism in Methanosarcina barkeri Fusaro

Pyruvate, a central intermediate in the carbon fixation pathway of methanogenic archaea, is rarely used as an energy source by these organisms. The sole exception to this rule is a genetically uncharacterized Methanosarcina barkeri mutant capable of using pyruvate as a sole energy and carbon source (the Pyr+ phenotype). Here, we provide evidence that suggests that the Pyr+ mutant is able to metabolize pyruvate by overexpressing pyruvate ferredoxin oxidoreductase (por) and mutating genes involved in central carbon metabolism. Genomic analysis showed that the Pyr+ strain has two mutations localized to Mbar_A1588, the biotin protein ligase subunit of the pyruvate carboxylase (pyc) operon, and Mbar_A2165, a putative transcriptional regulator. Mutants expressing the Mbar_A1588 mutation showed no growth defect compared to the wild type (WT), yet the strains lacked pyc activity. Recreation of the Mbar_A2165 mutation resulted in a 2-fold increase of Por activity and gene expression, suggesting a role in por transcriptional regulation. Further transcriptomic analysis revealed that Pyr+ strains also overexpress the gene encoding phosphoenolpyruvate carboxylase, indicating the presence of a previously uncharacterized route for synthesizing oxaloacetate in M. barkeri and explaining the unimpaired growth in the absence of Pyc. Surprisingly, stringent repression of the por operon was lethal, even when the media were supplemented with pyruvate and/or Casamino Acids, suggesting that por plays an unidentified essential function in M. barkeri.

59 BASIC BIOLOGICAL SCIENCES↗

Experimental Evolution of Extreme Resistance to Ionizing Radiation in Escherichia coli after 50 Cycles of Selection

In prior work (D. R. Harris et al., J Bacteriol 191:5240–5252, 2009, https://doi.org/10.1128/JB.00502-09; B. T. Byrne et al., Elife 3:e01322, 2014, https://doi.org/10.7554/eLife.01322), we demonstrated that Escherichia coli could acquire substantial levels of resistance to ionizing radiation (IR) via directed evolution. Major phenotypic contributions involved adaptation of organic systems for DNA repair. We have currently undertaken an extended effort to generate E. coli populations that are as resistant to IR as Deinococcus radiodurans . After an initial 50 cycles of selection using high-energy electron beam IR, four replicate populations exhibit major increases in IR resistance but have not yet reached IR resistance equivalent to D. radiodurans . Regular deep sequencing reveals complex evolutionary patterns with abundant clonal interference. Prominent IR resistance mechanisms involve novel adaptations to DNA repair systems and alterations in RNA polymerase. Adaptation is highly specialized to resist IR exposure, since isolates from the evolved populations exhibit highly variable patterns of resistance to other forms of DNA damage. Sequenced isolates from the populations possess between 184 and 280 mutations. IR resistance in one isolate, IR9-50-1, is derived largely from four novel mutations affecting DNA and RNA metabolism: RecD A90E, RecN K429Q, and RpoB S72N/RpoC K1172I. Additional mechanisms of IR resistance are evident. IMPORTANCE Some bacterial species exhibit astonishing resistance to ionizing radiation, with Deinococcus radiodurans being the archetype. As natural IR sources rarely exceed mGy levels, the capacity of Deinococcus to survive 5,000 Gy has been attributed to desiccation resistance. To understand the molecular basis of true extreme IR resistance, we are using experimental evolution to generate strains of Escherichia coli with IR resistance levels comparable to Deinococcus . Experimental evolution has previously generated moderate radioresistance for multiple bacterial species. Yet, these efforts could not take advantage of modern genomic sequencing technologies. In this report, we examine four replicate bacterial populations after 50 selection cycles. Genomic sequencing allows us to follow the genesis of mutations in populations throughout selection. Novel mutations affecting genes encoding DNA repair proteins and RNA polymerase enhance radioresistance. However, more contributors are apparent.

59 BASIC BIOLOGICAL SCIENCES↗

Growth Trade-Offs Accompany the Emergence of Glycolytic Metabolism in Shewanella oneidensis MR-1

ABSTRACT Bacteria increase their metabolic capacity via the acquisition of genetic material or by the mutation of genes already present in the genome. Here, we explore the mechanisms and trade-offs involved when Shewanella oneidensis , a bacterium that typically consumes small organic and amino acids, rapidly evolves to expand its metabolic capacity to catabolize glucose after a short period of adaptation to a glucose-rich environment. Using whole-genome sequencing and genetic approaches, we discovered that deletions in a region including the transcriptional repressor ( nagR ) that regulates the expression of genes associated with catabolism of N -acetylglucosamine are the common basis for evolved glucose metabolism across populations. The loss of nagR results in the constitutive expression of genes for an N -acetylglucosamine permease ( nagP ) and kinase ( nagK ). We demonstrate that promiscuous activities of both NagP and NagK toward glucose allow for the transport and phosphorylation of glucose to glucose-6-phosphate, the initial events of glycolysis otherwise thought to be absent in S. oneidensis . 13 C-based metabolic flux analysis uncovered that subsequent utilization was mediated by the Entner-Doudoroff pathway. This is an example whereby gene loss and preexisting enzymatic promiscuity, and not gain-of-function mutations, were the drivers of increased metabolic capacity. However, we observed a significant decrease in the growth rate on lactate after adaptation to glucose catabolism, suggesting that trade-offs may explain why glycolytic function may not be readily observed in S. oneidensis in natural environments despite it being readily accessible through just a single mutational event. IMPORTANCE Gains in metabolic capacity are frequently associated with the acquisition of novel genetic material via natural or engineered horizontal gene transfer events. Here, we explored how a bacterium that typically consumes small organic acids and amino acids expands its metabolic capacity to include glucose via a loss of genetic material, a process frequently associated with a deterioration of metabolic function. Our findings highlight how the natural promiscuity of transporters and enzymes can be a key driver in expanding metabolic diversity and that many bacteria may possess a latent metabolic capacity accessible through one or a few mutations that remove regulatory functions. Our discovery of trade-offs between growth on lactate and on glucose suggests why this easily gained trait is not observed in nature.

59 BASIC BIOLOGICAL SCIENCES↗

Structural Basis for Escape of Human Astrovirus from Antibody Neutralization: Broad Implications for Rational Vaccine Design

ABSTRACT Human astroviruses are recognized as a leading cause of viral diarrhea worldwide in children, immunocompromised patients, and the elderly. There are currently no vaccines available to prevent astrovirus infection; however, antibodies developed by healthy individuals during previous infection correlate with protection from reinfection, suggesting that an effective vaccine could be developed. In this study, we investigated the molecular mechanism by which several strains of human astrovirus serotype 2 (HAstV-2) are resistant to the potent HAstV-2-neutralizing monoclonal antibody PL-2 (MAb PL-2). Sequencing of the HAstV-2 capsid genes reveals mutations in the PL-2 epitope within the capsid's spike domain. To understand the molecular basis for resistance from MAb PL-2 neutralization, we determined the 1.35-Å-resolution crystal structure of the capsid spike from one of these HAstV-2 strains. Our structure reveals a dramatic conformational change in a loop within the PL-2 epitope due to a serine-to-proline mutation, locking the loop in a conformation that sterically blocks binding and neutralization by MAb PL-2. We show that mutation to serine permits loop flexibility and recovers MAb PL-2 binding. Importantly, we find that HAstV-2 capsid spike containing a serine in this loop is immunogenic and elicits antibodies that neutralize all HAstV-2 strains. Taken together, our results have broad implications for rational selection of vaccine strains that do not contain prolines in antigenic loops, so as to elicit antibodies against diverse loop conformations. IMPORTANCE Human astroviruses (HAstVs) infect nearly every person in the world during childhood and cause diarrhea, vomiting, and fever. In this study, we investigated how several strains of HAstV are resistant to a virus-neutralizing monoclonal antibody. We determined the crystal structure of the capsid protein spike domain from one of these HAstV strains and found that a single amino acid mutation induces a structural change in a loop that is responsible for antibody binding. Our findings reveal how viruses can escape antibody neutralization and provide insight for the rational design of vaccines to elicit diverse antibodies that provide broader protection from infection.

Bogdanoff, Walter A. (ORCID:0000000165321006)↗

N-Glycolylneuraminic Acid Binding of Avian and Equine H7 Influenza A Viruses

Influenza A viruses (IAV) initiate infection by binding to glycans with terminal sialic acids on the cell surface. Hosts of IAV variably express two major forms of sialic acid, N-acetylneuraminic acid (NeuAc) and N-glycolylneuraminic acid (NeuGc). NeuGc is produced in most mammals, including horses and pigs, but is absent in humans, ferrets, and birds. The only known naturally occurring IAV that exclusively bind NeuGc are extinct highly pathogenic equine H7N7 viruses. We determined the crystal structure of a representative equine H7 hemagglutinin (HA) in complex with NeuGc and observed high similarity in the receptor-binding domain with an avian H7 HA. To determine the molecular basis for NeuAc and NeuGc specificity, we performed systematic mutational analyses, based on the structural insights, on two distant avian H7 HAs and an H15 HA. We found that the A135E mutation is key for binding α2,3-linked NeuGc but does not abolish NeuAc binding. The additional mutations S128T, I130V, T189A, and K193R converted the specificity from NeuAc to NeuGc. We investigated the residues at positions 128, 130, 135, 189, and 193 in a phylogenetic analysis of avian and equine H7 HAs. This analysis revealed a clear distinction between equine and avian residues. The highest variability was observed at key position 135, of which only the equine glutamic acid led to NeuGc binding. These results demonstrate that genetically distinct H7 and H15 HAs can be switched from NeuAc to NeuGc binding and vice versa after the introduction of several mutations, providing insights into the adaptation of H7 viruses to NeuGc receptors.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of the I38T PA Substitution as a Resistance Marker for Next-Generation Influenza Virus Endonuclease Inhibitors

ABSTRACT The clinical severity and annual occurrence of influenza virus epidemics, combined with the availability of just a single class of antivirals to treat infections, underscores the urgent need to develop new anti-influenza drugs. The endonuclease activity within the viral acidic polymerase (PA) protein is an attractive target for drug discovery due to the critical role it plays in viral gene transcription. RO-7 is a next-generation PA endonuclease inhibitor of influenza A and B viruses, but its drug resistance potential is unknown. Through serial passage of influenza A(H1N1) viruses in MDCK cells under selective pressure of RO-7, we identified an I38T substitution within the PA endonuclease domain that conferred in vitro resistance to RO-7 (up to a 287-fold change in 50% effective concentration [EC 50 ]). I38T emerged between 5 and 10 passages, and when introduced into recombinant influenza A(H1N1) viruses, alone conferred RO-7 resistance (up to an 81-fold change in EC 50 ). Cocrystal structures of mutant and wild-type endonuclease domains with RO-7 provided the structural basis of resistance, where a key hydrophobic interaction between RO-7 and the Ile38 side chain is compromised when mutated to the polar threonine. While Ile38 does not have a crucial role in coordinating the endonuclease active site, the switch to threonine does affect the polymerase activity of some viruses and influences RO-7 affinity for the PA N target (i.e., the ≈200-residue N-terminal domain of PA). However, the change does not lead to a complete loss of replication activity in vitro . Our results predict that RO-7-resistant influenza viruses carrying the I38T substitution may emerge under treatment. This should be taken into consideration for clinical surveillance and in refinement of these drugs. IMPORTANCE The effectiveness of antiviral drugs can be severely compromised by the emergence of resistant viruses. Therefore, determination of the mechanisms by which viruses become resistant is critical for drug development and clinical use. RO-7 is a compound that potently inhibits influenza virus replication and belongs to a new class of drugs in late-stage clinical trials for treatment of influenza virus infection. Here we demonstrate that a single amino acid change acquired under prolonged virus exposure to RO-7 renders influenza viruses significantly less susceptible to its inhibitory effects. We have discovered how the mutation can simultaneously interfere with drug activity and still maintain efficient virus replication. These findings have important implications for the development of more effective derivatives of RO-7-like drugs and provide guidance for how to monitor the emergence of resistance.

Jones, Jeremy C.↗

The Frequency of Sex: Population Genomics Reveals Differences in Recombination and Population Structure of the Aflatoxin-Producing Fungus Aspergillus flavus

The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.

59 BASIC BIOLOGICAL SCIENCES↗