Mixed-effect modeling of codon usage

Logistic mixed effects models are used to determine whether optimal codons associate with two specific properties of the expressed protein: solvent accessibility, aggregation propensity, or evolutionary conservation. Both random components and fixed structures in the models are decided by following certain selection procedures. More models are also developed by considering different factor combinations using the same selection procedure. The results show that evolutionary conservation is the most important factor for predicting for the optimal codon usage for most amino acids; aggregation propensity is also an important factor, and solvent accessibility is the least important factor for most amino acids.The results of this analysis are consistent with the previous literature, provide more straightforward way to study the research question and also more information for the insight relationships.