RESUMO
In the present work, the distribution of the electronic charge density of the natural sex pheromone, the (Z)-13-hexadecen-11-ynyl acetate, in the female processionary moth, Thaumetopoea pytiocampa, and its nine analogue derivatives was studied within the framework of the Density Functional Theory and the Atoms in Molecules (AIM) Theory at B3LYP/6-31G *//B3LYP/6-31++G * * level. Additionally, molecular electrostatic potential (MEP) maps of the previously mentioned compounds were computed and compared. Furthermore, the substitution of hydrogen atoms from the methyl group in the acetate group by electron withdrawing substituents (i.e., halogen atoms) as well as the replacement effect of hydrogen by electron donor substituents (+I effect) as methyl group, were explored. The key feature of the topological distribution of the charge density in analogue compounds, such as the variations of the topological properties encountered in the region formed by neighbouring atoms from the substitution site were presented and discussed. Using topological parameters, such as electronic charge density, Laplacian, kinetic energy density, and potential energy density evaluated at bond critical points (BCP), we provide here a detailed analysis of the nature of the chemical bonding of these molecules. In addition, the atomic properties (population, charge, energy, volume, and dipole moment) were determined on selected atoms. These properties were analyzed at the substitution site (with respect to the natural sex pheromone) and related to the biological activity and to the possible binding site with the pheromone binding protein, (PBP). Moreover, the Laplacian function of the electronic density was used to locate electrophilic regions susceptible to be attacked (by deficient electron atoms or donor hydrogen). Our results indicate that the change in the atomic properties, such as electronic population and atomic volume, are sensitive indicators of the loss of the biological activity in the analogues studied here. The crucial interaction between the acetate group of the natural sex pheromone and the PBP is most likely to be a hydrogen bonding and the substitution of hydrogen atoms by electronegative atoms in the pheromone molecule reduces the hydrogen acceptor capacity. This situation is mirrored by the diminish of the electronic population on carbon and oxygen atoms at the carbonylic group in the halo-acetate group. Additionally, the modified acetate group (with electronegative atoms) shows new charge concentration critical points or regions of concentration of charge density in which an electrophilic attack can also occur. Finally, the use of the topological analysis based in the charge density distribution and its Laplacian function, in conjunction with MEP maps provides valuable information about the steric volume and electronic requirement of the sex pheromone for binding to the PBP.