Karine Leitão Lima Thiers1 & João Hermínio Martins da Silva2 & Geraldo Rodrigues Sartori2 & Clesivan Pereira dos Santos1 & Kátia Daniella da Cruz Saraiva1,3 &André Luiz Maia Roque1 & Birgit Arnholdt-Schmitt1,4,5 & José Hélio Costa1 

Plant plastoquinol oxidase (PTOX) is a chloroplast oxidoreductase involved in carotenoid biosynthesis, chlororespiration, and response to environmental stresses. The present study aimed to gain insight of the potential role of nucleotide/amino acid changes linked to environmental adaptation in PTOX gene/protein from Arabidopsis thaliana accessions. SNPs in the single-copy PTOX gene were identified in 1190 accessions of Arabidopsis using the Columbia-0 PTOX as a reference. The identified SNPs were correlated with geographical distribution of the accessions according to altitude, climate, and rainfall. Among the 32 identified SNPs in the coding region of the PTOX gene, 16 of these were characterized as non-synonymous SNPs (in which an AA is altered). A higher incidence of AA changes occurred in the mature protein at positions 78 (31%), 81 (31.4%), and 323 (49.9%). Three-dimensional structure prediction indicated that the AA change at position 323 (D323N) leads to a PTOX structure with the most favorable interaction with the substrate plastoquinol, when compared with the reference PTOX structure (Columbia-0). Molecular docking analysis suggested that the most favorable D323N PTOX-plastoquinol interaction is due to a better enzyme- substrate binding affinity. The molecular dynamics revealed that plastoquinol should be more stable in complex with D323N PTOX, likely due a restraint mechanism in this structure that stabilize plastoquinol inside of the reaction center. The integrated analysis made from accession geographical distribution and PTOX SNPs indicated that AA changes in PTOX are related to altitude and rainfall, potentially due to an adaptive positive environmental selection.