Positive Selection or Free to Vary? Assessing the Functional Significance of Sequence Change Using Molecular Dynamics

PLoS One. 2016 Feb 12;11(2):e0147619. doi: 10.1371/journal.pone.0147619. eCollection 2016.

Abstract

Evolutionary arms races between pathogens and their hosts may be manifested as selection for rapid evolutionary change of key genes, and are sometimes detectable through sequence-level analyses. In the case of protein-coding genes, such analyses frequently predict that specific codons are under positive selection. However, detecting positive selection can be non-trivial, and false positive predictions are a common concern in such analyses. It is therefore helpful to place such predictions within a structural and functional context. Here, we focus on the p19 protein from tombusviruses. P19 is a homodimer that sequesters siRNAs, thereby preventing the host RNAi machinery from shutting down viral infection. Sequence analysis of the p19 gene is complicated by the fact that it is constrained at the sequence level by overprinting of a viral movement protein gene. Using homology modeling, in silico mutation and molecular dynamics simulations, we assess how non-synonymous changes to two residues involved in forming the dimer interface-one invariant, and one predicted to be under positive selection-impact molecular function. Interestingly, we find that both observed variation and potential variation (where a non-synonymous change to p19 would be synonymous for the overprinted movement protein) does not significantly impact protein structure or RNA binding. Consequently, while several methods identify residues at the dimer interface as being under positive selection, MD results suggest they are functionally indistinguishable from a site that is free to vary. Our analyses serve as a caveat to using sequence-level analyses in isolation to detect and assess positive selection, and emphasize the importance of also accounting for how non-synonymous changes impact structure and function.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Amino Acid Sequence
  • Binding Sites
  • Crystallography, X-Ray
  • Gene Expression
  • Immune Evasion / genetics
  • Molecular Dynamics Simulation*
  • Molecular Sequence Data
  • Mutation*
  • Phylogeny
  • Plant Diseases / immunology
  • Plant Diseases / virology
  • Protein Binding
  • Protein Interaction Domains and Motifs
  • Protein Multimerization
  • Protein Structure, Secondary
  • RNA / chemistry*
  • Selection, Genetic
  • Sequence Alignment
  • Sequence Homology, Amino Acid
  • Solanum lycopersicum / immunology
  • Solanum lycopersicum / virology
  • Tombusvirus / classification
  • Tombusvirus / genetics*
  • Tombusvirus / immunology
  • Viral Proteins / chemistry*
  • Viral Proteins / genetics

Substances

  • P19 protein, tomato bushy stunt virus
  • Viral Proteins
  • RNA

Grants and funding

The authors received funding from the following sources: Royal Society of New Zealand Marsden Fund Fast Start award (13-MAU-039) to JRA; Royal Society of New Zealand Rutherford Discovery Fellowship (RDF-11-UOC-013) to AMP; Biomolecular Interaction Centre, University of Canterbury (to AMP, JRA); Massey University Research Fund Early Career Grant (to JRA); and Swedish Research Council (to AMP). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.