The in silico studies were performed in order to assess the binding affinity of selected organophosphorus compounds towards the acetylcholinesterase enzyme (AChE). Quantum mechanical calculations, molecular docking and molecular dynamics (MD) with molecular mechanics Generalized-Born surface area (MM/GBSA) were applied to assess quantitatively differences between the binding energies of acetylcholine (ACh; the natural agonist of AChE) and neurotoxic, synthetic correlatives (so-called “Novichoks”, and selected compounds from the G- and V-series). Several additional quantitative descriptors like root mean square fluctuation (RMSF) and the solvent accessible surface area (SASA) were briefly discussed to give – to the best of our knowledge – the first quantitative in silico description of AChE – Novichok non-covalent binding process and thus facilitate the search for an efficient and effective treatment for Novichok intoxication and in a broader sense - intoxication with other warfare nerve agents as well.