Two broad paradigms exist for inferring dN=dS; the ratio of nonsynonymous to synonymous substitution rates, from coding sequences: (i) a one-rate approach, where dN=dS is represented with a single parameter, or (ii) a two-rate approach, where d N and d S are estimated separately. The performances of these two approaches have been well studied in the specific context of proper model specification, i.e., when the inference model matches the simulation model. By contrast, the relative performances of one-rate vs. two-rate parameterizations when applied to data generated according to a different mechanism remain unclear. Here, we compare the relative merits of one-rate and two-rate approaches in the specific context of model misspecification by simulating alignments with mutation-selection models rather than with dN=dS-based models. We find that one-rate frameworks generally infer more accurate dN=dS point estimates, even when d S varies among sites. In other words, modeling d S variation may substantially reduce accuracy of dN=dS point estimates. These results appear to depend on the selective constraint operating at a given site. For sites under strong purifying selection (dN=dS ≲ 0:3), one-rate and two-rate models show comparable performances. However, one-rate models significantly outperform two-rate models for sites under moderate-to-weak purifying selection. We attribute this distinction to the fact that, for these more quickly evolving sites, a given substitution is more likely to be nonsynonymous than synonymous. The data will therefore be relatively enriched for nonsynonymous changes, and modeling d S contributes excessive noise to dN=dS estimates. We additionally find that high levels of divergence among sequences, rather than the number of sequences in the alignment, are more critical for obtaining precise point estimates.KEYWORDS dN/dS; mutation-selection models; evolutionary rate; sequence simulation; molecular evolution A variety of computational approaches have been developed to infer selection pressure from protein-coding sequences in a phylogenetically aware context. Among the most commonly used methods are those that compute the evolutionary rate ratio dN=dS; which represents the ratio of nonsynonymous to synonymous substitution rates. Beginning in the mid-1990s, this value has been calculated using maximum-likelihood (ML) approaches (Goldman and Yang 1994;Muse and Gaut 1994), and since then, a wide variety of inference frameworks have been developed to infer dN=dS at individual sites in protein-coding sequences (Nielsen and Yang 1998;Yang et al. 2000;Yang and Nielsen 2002;Yang and Swanson 2002;Lemey et al. 2012;Murrell et al. 2012bMurrell et al. , 2013.Typically, the performance of evolutionary models is examined using simulation-based approaches wherein sequences are simulated according to the model being examined, and inferences are subsequently performed on the simulated data. This approach ensures that simulated and inferred parameters correspond. Although useful, this strategy fundam...