Rational in silico design of aptamers for organophosphates based on the example of paraoxon

Poisoning by organophosphates (OPs) takes one of the leading places in the total number of exotoxicoses. Detoxication of OPs at the first stage of the poison entering the body could be achieved with the help of DNA- or RNA-aptamers, which are able to bind poisons in the bloodstream. The aim of the research was to develop an approach to rational in silico design of aptamers for OPs based on the example of paraoxon. From the published sequence of an aptamer binding organophosphorus pesticides, its three-dimensional model has been constructed. The most probable binding site for paraoxon was determined by molecular docking and molecular dynamics (MD) methods. Then the nucleotides of the binding site were mutated consequently and the values of free binding energy have been calculated using MD trajectories and MM-PBSA approach. On the basis of the energy values, two sequences that bind paraoxon most efficiently have been selected. The value of free binding energy of paraoxon with peripheral anionic site of acetylcholinesterase (AChE) has been calculated as well. It has been revealed that the aptamers found bind paraoxon more effectively than AChE. The peculiarities of paraoxon interaction with the aptamers nucleotides have been analyzed. The possibility of improving in silico approach for aptamer selection is discussed.