Computational Compensatory Mutation Discovery Approach: Predicting a PARP1 Variant Rescue Mutation

The prediction of protein mutations that affect function may be exploited for multiple uses. In the context of disease variants, the prediction of compensatory mutations that reestablish functional phenotypes could aid in the development of genetic therapies. In this work, we present an integrated approach that combines coevolutionary analysis and molecular dynamics (MD) simulations to discover functional compensatory mutations. This approach is employed to investigate possible rescue mutations of a poly(ADP-ribose) polymerase 1 (PARP1) variant, PARP1 V762A, associated with lung cancer and follicular lymphoma. MD simulations show PARP1 V762A exhibits noticeable changes in structural and dynamical behavior compared with wild type PARP1. Our integrated approach predicts A755E as a possible compensatory mutation based on coevolutionary information, and molecular simulations indicate that the PARP1 A755E/V762A double mutant exhibits similar structural and dynamical behavior to WT PARP1. Our methodology can be broadly applied to a large number of systems where SNPs have been identified as connected to disease and can shed light on the biophysical effects of such changes as well as provide a way to discover potential mutants that could restore wild type-like functionality. This can in turn be further utilized in the design of molecular therapeutics that aim to mimic such compensatory effect.

Plasmid fitness costs are caused by specific genetic conflicts enabling resolution by compensatory mutation

Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid carriage are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or gene expression level. By combining the results of experimental evolution with genetics and transcriptomics, we show here that fitness costs of 2 divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid carriage. We identify one of these compensatory loci, the chromosomal gene PFLU4242, as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by up-regulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained up-regulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer (HGT) due to their propensity for amelioration by single compensatory mutations, helping to explain why plasmids are so common in bacterial genomes.

Respiratory Syncytial Virus Phosphoprotein Residue S156 Plays a Role in Regulating Genome Transcription and Replication

Respiratory syncytial virus (RSV) is a single-stranded, negative-sense, RNA virus in the family Pneumoviridae and genus Orthopneumoviridae that can cause severe disease in infants, immunocompromised adults, and the elderly. The RSV viral RNA-dependent RNA polymerase (vRdRp) complex is composed of the phosphoprotein (P) and the large polymerase protein (L). The P protein is constitutively phosphorylated by host kinases and has 41 serine (S) and threonine (T) residues as potential phosphorylation sites. To identify important phosphorylation residues in the P protein, we systematically and individually mutated all serine S and T residues to alanine (A) and first analyzed their effect on genome transcription and replication using a minigenome system. We found that the mutation of eight residues resulted in significantly reduced minigenome activity compared to wild-type P. We then incorporated these mutations (T210A, S203A, T151A, S156A, T160A, S23A, T188A, and T105A) into full-length genome cDNA to rescue recombinant RSV. We were able to recover four recombinant viruses (T151A, S156A, T160A, and S23A), suggesting RSV-P residues T210, S203, T188, and T105 are essential for viral RNA replication. Among the four rescued, rRSV-T160A caused a minor growth defect compared to its parental virus while rRSV-S156A had severely restricted replication due to decreased levels of genomic RNA. During infection, P-S156A phosphorylation was decreased, and when passaged, the S156A virus acquired a known compensatory mutation in L (L795I) that enhanced both WT-P and P-S156A minigenome activity and was able to partially rescue the S156A viral growth defect. This work demonstrates that residues T210, S203, T188, and T105 are critical for RSV replication, and S156 plays a critical role in viral RNA synthesis. Importance RSV-P is a heavily phosphorylated protein that is required for RSV replication. In this study, we identified several residues, including P-S156, as phosphorylation sites that play critical roles in efficient viral growth and genome replication. Future studies to identify the specific kinase(s) that phosphorylate these residues can lead to kinase inhibitors and anti-viral drugs for this important human pathogen.

An Absolutely Conserved Tryptophan in the Stem of the Envelope Protein E of Flaviviruses Is Essential for the Formation of Stable Particles

The major envelope protein E of flaviviruses contains an ectodomain that is connected to the transmembrane domain by the so-called “stem” region. In mature flavivirus particles, the stem is composed of two or three mostly amphipathic α-helices and a conserved sequence element (CS) with an undefined role in the viral life cycle. A tryptophan is the only residue within this region which is not only conserved in all vector-borne flaviviruses, but also in the group with no known vector. We investigated the importance of this residue in different stages of the viral life cycle by a mutagenesis-based approach using tick-borne encephalitis virus (TBEV). Replacing W421 by alanine or histidine strongly reduced the release of infectious virions and their thermostability, whereas fusion-related entry functions and virus maturation were still intact. Serial passaging of the mutants led to the emergence of a same-site compensatory mutation to leucine that largely restored these properties of the wildtype. The conserved tryptophan in CS (or another big hydrophobic amino acid at the same position) is thus essential for the assembly and infectivity of flaviviruses by being part of a network required for conferring stability to infectious particles.

Expression dysregulation as a mediator of fitness costs in antibiotic resistance

Antimicrobial resistance (AMR) poses a threat to global health and the economy. Rifampicin-resistant Mycobacterium tuberculosis accounts for a third of the global AMR burden. Gaining the upper hand on AMR requires a deeper understanding of the physiology of resistance. AMR often results in a fitness cost in absence of drug. Identifying the molecular mechanisms underpinning this cost could help strengthen future treatment regimens. Here, we used a collection of M. tuberculosis strains providing an evolutionary and phylogenetic snapshot of rifampicin resistance, and subjected them to genome-wide transcriptomic and proteomic profiling to identify key perturbations of normal physiology. We found that the clinically most common rifampicin resistance-conferring mutation RpoB Ser450Leu imparts considerable gene expression changes, many of which are mitigated by the compensatory mutation in RpoC Leu516Pro. However, our data also provide evidence for pervasive epistasis: the same resistance mutation imposed a different fitness cost and functionally distinct changes to gene expression in genetically unrelated clinical strains. Finally, we report a likely post-transcriptional modulation of gene expression that is shared in most of the tested strains carrying RpoB Ser450Leu, resulting in an increased abundance of proteins involved in central carbon metabolism. These changes contribute to a more general trend, in which the disruption of the composition of the proteome correlates with the fitness cost of the RpoB Ser450Leu mutation in different strains.

Oncolytic virotherapy induced CSDE1 neo-antigenesis restricts VSV replication but can be targeted by immunotherapy

In our clinical trials of oncolytic vesicular stomatitis virus expressing interferon beta (VSV-IFNβ), several patients achieved initial responses followed by aggressive relapse. We show here that VSV-IFNβ-escape tumors predictably express a point-mutated CSDE1P5S form of the RNA-binding Cold Shock Domain-containing E1 protein, which promotes escape as an inhibitor of VSV replication by disrupting viral transcription. Given time, VSV-IFNβ evolves a compensatory mutation in the P/M Inter-Genic Region which rescues replication in CSDE1P5S cells. These data show that CSDE1 is a major cellular co-factor for VSV replication. However, CSDE1P5S also generates a neo-epitope recognized by non-tolerized T cells. We exploit this predictable neo-antigenesis to drive, and trap, tumors into an escape phenotype, which can be ambushed by vaccination against CSDE1P5S, preventing tumor escape. Combining frontline therapy with escape-targeting immunotherapy will be applicable across multiple therapies which drive tumor mutation/evolution and simultaneously generate novel, targetable immunopeptidomes associated with acquired treatment resistance.

Septal Class A Penicillin-Binding Protein Activity and LD-Transpeptidases Mediate Selection of Colistin-Resistant Lipooligosaccharide-Deficient Acinetobacter baumannii

ABSTRACT Despite dogma suggesting that lipopolysaccharide/lipooligosaccharide (LOS) was essential for viability of Gram-negative bacteria, several Acinetobacter baumannii clinical isolates produced LOS− colonies after colistin selection. Inactivation of the conserved class A penicillin-binding protein, PBP1A, was a compensatory mutation that supported isolation of LOS− A. baumannii, but the impact of PBP1A mutation was not characterized. Here, we show that the absence of PBP1A causes septation defects and that these, together with ld-transpeptidase activity, support isolation of LOS− A. baumannii. PBP1A contributes to proper cell division in A. baumannii, and its absence induced cell chaining. Only isolates producing three or more septa supported selection of colistin-resistant LOS− A. baumannii. PBP1A was enriched at the midcell, where the divisome complex facilitates daughter cell formation, and its localization was dependent on glycosyltransferase activity. Transposon mutagenesis showed that genes encoding two putative ld-transpeptidases (LdtJ and LdtK) became essential in the PBP1A mutant. Both LdtJ and LdtK were required for selection of LOS− A. baumannii, but each had distinct enzymatic activities in the cell. Together, these findings demonstrate that defects in PBP1A glycosyltransferase activity and ld-transpeptidase activity remodel the cell envelope to support selection of colistin-resistant LOS− A. baumannii. IMPORTANCE The increasing prevalence of antibiotic treatment failure associated with Gram-negative bacterial infections highlights an urgent need to develop new alternative therapeutic strategies. The last-line antimicrobial colistin (polymyxin E) targets the ubiquitous outer membrane lipopolysaccharide (LPS)/LOS membrane anchor, lipid A, which is essential for viability of most diderms. However, several LOS− Acinetobacter baumannii clinical isolates were recovered after colistin selection, suggesting a conserved resistance mechanism. Here, we characterized a role for penicillin-binding protein 1A in A. baumannii septation and intrinsic β-lactam susceptibility. We also showed that defects in PBP1A glycosyltransferase activity and ld-transpeptidase activity support isolation of colistin-resistant LOS− A. baumannii.

Reversion is most likely under high mutation supply, when compensatory mutations don’t fully restore fitness costs

Adaptive mutations are often associated with a fitness cost. These costs can be compensated for through the acquisition of additional mutations, or the adaptations can be lost through reversion, in settings where they are no longer favored. While the dynamics of adaptation, reversion and compensation have been central features in several studies of microbial evolution, few studies have attempted to resolve the population genetics underlying how and when either compensation or reversion occur. Specifically, questions remain regarding how certain actors—the evolution of mutators and whether compensatory mutations alleviate costs fully or partially—may influence evolutionary dynamics of compensation and reversion. In this study, we attempt to explain findings from an experimental evolution study by utilizing computational and theoretical approaches towards a more refined understanding of how mutation rate and the fitness effects of compensatory mutation influence evolutionary dynamics. We find that high mutation rates increase the probability of reversion of deleterious adaptations when compensation is only partial. The existence of even a single fully compensatory mutation is associated with a dramatically decreased probability of reversion. Experimental results suggest that, in some contexts, compensatory mutations are not able to fully alleviate costs associated with adaption. Our findings emphasize the role of both mutation rate and the fitness effects of compensatory mutation in crafting evolutionary dynamics, and highlight the importance of population genetic theory for explaining findings from experimental evolution.

Structural characterization of HIV-1 matrix mutants implicated in envelope incorporation

ABSTRACTDuring the late phase of HIV-1 infection, viral Gag polyproteins are targeted to the plasma membrane (PM) for assembly. Gag localization at the PM is a prerequisite for the incorporation of the envelope protein (Env) into budding particles. Gag assembly and Env incorporation are mediated by the N-terminal myristoylated matrix (MA) domain of Gag. Nonconservative mutations in the trimer interface of MA (A45E, T70R, and L75G) were found to impair Env incorporation and infectivity, leading to the hypothesis that MA trimerization is an obligatory step for Env incorporation. Conversely, Env incorporation can be rescued by a compensatory mutation in the MA trimer interface (Q63R). The impact of these MA mutations on the structure and trimerization properties of MA is not known. In this study, we employed NMR spectroscopy, x-ray crystallography, and sedimentation techniques to characterize the structure and trimerization properties of HIV-1 MA A45E, Q63R, T70R, and L75G mutant proteins. NMR data revealed that these point mutations did not alter the overall structure and folding of MA but caused minor structural perturbations in the trimer interface. Analytical ultracentrifugation data indicated that mutations had a minimal effect on the MA monomer–trimer equilibrium. The high-resolution x-ray structure of the unmyristoylated MA Q63R protein revealed hydrogen bonding between the side chains of Arg-63 and Ser-67 located in the center of the trimer interface, providing the first structural evidence for a stabilization of the trimer form. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.

Radial Expansion Facilitates the Maintenance of Double Antibiotic Resistances

ABSTRACT Most microbes live in spatially confined subpopulations. Under spatial structure conditions, the efficacy of natural selection is often reduced (relative to homogeneous conditions) due to the increased importance of genetic drift and local competition. Additionally, under spatial structure conditions, the fittest genotype may not always be the one with better access to the heterogeneous distribution of nutrients. The effect of radial expansion may be particularly relevant for the elimination of antibiotic resistance mutations, as their dynamics within bacterial populations are strongly dependent on their growth rate. Here, we use Escherichia coli to systematically compare the allele frequency of streptomycin, rifampin, and fluoroquinolone single and double resistance mutants after 24 h of coexistence with a susceptible strain under radial expansion (local competition) and homogeneous (global competition) conditions. We show that there is a significant effect of structure on the maintenance of double resistances which is not observed for single resistances. Radial expansion also facilitates the persistence of double resistances when competing against their single counterparts. Importantly, we found that spatial structure reduces the rate of compensation of the double mutant RpsLK43T RpoBH526Y and that a strongly compensatory mutation in homogeneous conditions becomes deleterious under spatial structure conditions. Overall, our results unravel the importance of spatial structure for facilitating the maintenance and accumulation of multiple resistances over time and for determining the identity of compensatory mutations.