Analyzing Machupo virus-receptor binding by molecular dynamics simulations

In many biological applications, we would like to be able to computationally predict mutational effects on affinity in protein-protein interactions. However, many commonly used methods to predict these effects perform poorly in important test cases. In particular, the effects of multiple mutations, non-alanine substitutions, and flexible loops are difficult to predict with available tools and protocols. We present here an existing method applied in a novel way to a new test case; we interrogate affinity differences resulting from mutations in a host-virus protein-protein interface. We use steered molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then approximate affinity using the maximum applied force of separation and the area under the force-versus-distance curve. We find, even without the rigor and planning required for free energy calculations, that these quantities can provide novel biophysical insight into the GP1/hTfR1 interaction. First, with no prior knowledge of the system we can differentiate among wild type and mutant complexes. Moreover, we show that this simple SMD scheme correlates well with relative free energy differences computed via free energy perturbation. Second, although the static co-crystal structure shows two large hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate that one of them may not be important for tight binding. Third, one viral site known to be critical for infection may mark an important evolutionary suppressor site for infection-resistant hTfR1 mutants. Finally, our approach provides a framework to compare the effects of multiple mutations, individually and jointly, on protein-protein interactions.Comment: 33 pages, 8 figures, 5 table

Identification of a Genomic Reservoir for New TRIM Genes in Primate Genomes

Tripartite Motif (TRIM) ubiquitin ligases act in the innate immune response against viruses. One of the best characterized members of this family, TRIM5α, serves as a potent retroviral restriction factor with activity against HIV. Here, we characterize what are likely to be the youngest TRIM genes in the human genome. For instance, we have identified 11 TRIM genes that are specific to humans and African apes (chimpanzees, bonobos, and gorillas) and another 7 that are human-specific. Many of these young genes have never been described, and their identification brings the total number of known human TRIM genes to approximately 100. These genes were acquired through segmental duplications, most of which originated from a single locus on chromosome 11. Another polymorphic duplication of this locus has resulted in these genes being copy number variable within the human population, with a Han Chinese woman identified as having 12 additional copies of these TRIM genes compared to other individuals screened in this study. Recently, this locus was annotated as one of 34 “hotspot” regions that are also copy number variable in the genomes of chimpanzees and rhesus macaques. Most of the young TRIM genes originating from this locus are expressed, spliced, and contain signatures of positive natural selection in regions known to determine virus recognition in TRIM5α. However, we find that they do not restrict the same retroviruses as TRIM5α, consistent with the high degree of divergence observed in the regions that control target specificity. We propose that this recombinationally volatile locus serves as a reservoir from which new TRIM genes arise through segmental duplication, allowing primates to continually acquire new antiviral genes that can be selected to target new and evolving pathogens

Special Issue: 2020 Research Data Access and Preservation Summit

The Journal of eScience Librarianship has partnered with the Research Data Access & Preservation (RDAP) Association for a third year to publish selected conference proceedings. This issue highlights the research presented at the RDAP 2020 Summit and the community it has fostered

Discordant Evolution of the Adjacent Antiretroviral Genes TRIM22 and TRIM5 in Mammals

TRIM5α provides a cytoplasmic block to retroviral infection, and orthologs encoded by some primates are active against HIV. Here, we present an evolutionary comparison of the TRIM5 gene to its closest human paralogs: TRIM22, TRIM34, and TRIM6. We show that TRIM5 and TRIM22 have a dynamic history of gene expansion and loss during the evolution of mammals. The cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. In contrast, TRIM6 and TRIM34 have been strictly preserved as single gene orthologs in human, dog, and cow. A more focused analysis of primates reveals that, while TRIM6 and TRIM34 have evolved under purifying selection, TRIM22 has evolved under positive selection as was previously observed for TRIM5. Based on TRIM22 sequences obtained from 27 primate genomes, we find that the positive selection of TRIM22 has occurred episodically for approximately 23 million years, perhaps reflecting the changing pathogenic landscape. However, we find that the evolutionary episodes of positive selection that have acted on TRIM5 and TRIM22 are mutually exclusive, with generally only one of these genes being positively selected in any given primate lineage. We interpret this to mean that the positive selection of one gene has constrained the adaptive flexibility of its neighbor, probably due to genetic linkage. Finally, we find a striking congruence in the positions of amino acid residues found to be under positive selection in both TRIM5α and TRIM22, which in both proteins fall predominantly in the β2-β3 surface loop of the B30.2 domain. Astonishingly, this same loop is under positive selection in the multiple cow TRIM5 genes as well, indicating that this small structural loop may be a viral recognition motif spanning a hundred million years of mammalian evolution

An expanded clade of rodent Trim5 genes

AbstractTrim5α from primates (including humans), cows, and rabbits has been shown to be an active antiviral host gene that acts against a range of retroviruses. Although this suggests that Trim5α may be a common antiviral restriction factor among mammals, the status of Trim5 genes in rodents has been unclear. Using genomic and phylogenetic analyses, we describe an expanded paralogous cluster of at least eight Trim5-like genes in mice (including the previously described Trim12 and Trim30 genes), and three Trim5-like genes in rats. Our characterization of the rodent Trim5 locus, and comparison to the Trim5 locus in humans, cows, and rabbits, indicates that Trim5 has undergone independent evolutionary expansions within species. Evolutionary analysis shows that rodent Trim5 genes have evolved under positive selection, suggesting evolutionary conflicts consistent with important antiviral function. Sampling six rodent Trim5 genes failed to reveal antiviral activities against a set of eight retroviral challenges, although we predict that such activities exist

Rapid Evolution of BRCA1 and BRCA2 in Humans and Other Primates

The maintenance of chromosomal integrity is an essential task of every living organism and cellular repair mechanisms exist to guard against insults to DNA. Given the importance of this process, it is expected that DNA repair proteins would be evolutionarily conserved, exhibiting very minimal sequence change over time. However, BRCA1, an essential gene involved in DNA repair, has been reported to be evolving rapidly despite the fact that many protein-altering mutations within this gene convey a significantly elevated risk for breast and ovarian cancers. Results: To obtain a deeper understanding of the evolutionary trajectory of BRCA1, we analyzed complete BRCA1 gene sequences from 23 primate species. We show that specific amino acid sites have experienced repeated selection for amino acid replacement over primate evolution. This selection has been focused specifically on humans and our closest living relatives, chimpanzees (Pan troglodytes) and bonobos (Pan paniscus). After examining BRCA1 polymorphisms in 7 bonobo, 44 chimpanzee, and 44 rhesus macaque (Macaca mulatta) individuals, we find considerable variation within each of these species and evidence for recent selection in chimpanzee populations. Finally, we also sequenced and analyzed BRCA2 from 24 primate species and find that this gene has also evolved under positive selection. Conclusions: While mutations leading to truncated forms of BRCA1 are clearly linked to cancer phenotypes in humans, there is also an underlying selective pressure in favor of amino acid-altering substitutions in this gene. A hypothesis where viruses are the drivers of this natural selection is discussed.National Institutes of Health R01-GM-093086, 8U42OD011197-13National Science Foundation BCS-07115972Burroughs Wellcome FundMolecular Bioscience

Widespread Sequence Variations in VAMP1 across Vertebrates Suggest a Potential Selective Pressure from Botulinum Neurotoxins

Botulinum neurotoxins (BoNT/A-G), the most potent toxins known, act by cleaving three SNARE proteins required for synaptic vesicle exocytosis. Previous studies on BoNTs have generally utilized the major SNARE homologues expressed in brain (VAMP2, syntaxin 1, and SNAP-25). However, BoNTs target peripheral motor neurons and cause death by paralyzing respiratory muscles such as the diaphragm. Here we report that VAMP1, but not VAMP2, is the SNARE homologue predominantly expressed in adult rodent diaphragm motor nerve terminals and in differentiated human motor neurons. In contrast to the highly conserved VAMP2, BoNT-resistant variations in VAMP1 are widespread across vertebrates. In particular, we identified a polymorphism at position 48 of VAMP1 in rats, which renders VAMP1 either resistant (I48) or sensitive (M48) to BoNT/D. Taking advantage of this finding, we showed that rat diaphragms with I48 in VAMP1 are insensitive to BoNT/D compared to rat diaphragms with M48 in VAMP1. This unique intra-species comparison establishes VAMP1 as a physiological toxin target in diaphragm motor nerve terminals, and demonstrates that the resistance of VAMP1 to BoNTs can underlie the insensitivity of a species to members of BoNTs. Consistently, human VAMP1 contains I48, which may explain why humans are insensitive to BoNT/D. Finally, we report that residue 48 of VAMP1 varies frequently between M and I across seventeen closely related primate species, suggesting a potential selective pressure from members of BoNTs for resistance in vertebrates

Dengue viruses cleave STING in humans but not in nonhuman primates, their presumed natural reservoir.

Human dengue viruses emerged from primate reservoirs, yet paradoxically dengue does not reach high titers in primate models. This presents a unique opportunity to examine the genetics of spillover versus reservoir hosts. The dengue virus 2 (DENV2) - encoded protease cleaves human STING, reducing type I interferon production and boosting viral titers in humans. We find that both human and sylvatic (reservoir) dengue viruses universally cleave human STING, but not the STING of primates implicated as reservoir species. The special ability of dengue to cleave STING is thus specific to humans and a few closely related ape species. Conversion of residues 78/79 to the human-encoded 'RG' renders all primate (and mouse) STINGs sensitive to viral cleavage. Dengue viruses may have evolved to increase viral titers in the dense and vast human population, while maintaining decreased titers and pathogenicity in the more rare animals that serve as their sustaining reservoir in nature

TRIM5alpha Restricts Flavivirus Replication by Targeting the Viral Protease for Proteasomal Degradation

Tripartite motif-containing protein 5alpha (TRIM5alpha) is a cellular antiviral restriction factor that prevents early events in retrovirus replication. The activity of TRIM5alpha is thought to be limited to retroviruses as a result of highly specific interactions with capsid lattices. In contrast to this current understanding, we show that both human and rhesus macaque TRIM5alpha suppress replication of specific flaviviruses. Multiple viruses in the tick-borne encephalitis complex are sensitive to TRIM5alpha-dependent restriction, but mosquito-borne flaviviruses, including yellow fever, dengue, and Zika viruses, are resistant. TRIM5alpha suppresses replication by binding to the viral protease NS2B/3 to promote its K48-linked ubiquitination and proteasomal degradation. Importantly, TRIM5alpha contributes to the antiviral function of IFN-I against sensitive flaviviruses in human cells. Thus, TRIM5alpha possesses remarkable plasticity in the recognition of diverse virus families, with the potential to influence human susceptibility to emerging flaviviruses of global concern

Fine-mapping identifies multiple prostate cancer risk loci at 5p15, one of which associates with TERT expression

Associations between single nucleotide polymorphisms (SNPs) at 5p15 and multiple cancer types have been reported. We have previously shown evidence for a strong association between prostate cancer (PrCa) risk and rs2242652 at 5p15, intronic in the telomerase reverse transcriptase (TERT) gene that encodes TERT. To comprehensively evaluate the association between genetic variation across this region and PrCa, we performed a fine-mapping analysis by genotyping 134 SNPs using a custom Illumina iSelect array or Sequenom MassArray iPlex, followed by imputation of 1094 SNPs in 22 301 PrCa cases and 22 320 controls in The PRACTICAL consortium. Multiple stepwise logistic regression analysis identified four signals in the promoter or intronic regions of TERT that independently associated with PrCa risk. Gene expression analysis of normal prostate tissue showed evidence that SNPs within one of these regions also associated with TERT expression, providing a potential mechanism for predisposition to disease