Nitrobenzene derivatives present an environmental concern to aquatic organisms, as well as a potential adverse effects on human health due to their toxicity. Therefore, monitoring nitrobenzene contamination levels is critical to protect both the ecosystem and the public. This study aims to quantitatively correlate the toxicity of nitrobenzene derivatives to their chemical structures using computational chemistry methods such as 3D‐QSPR, HQSPR, Molecular Docking and MD simulation. Additionally, to investigate the type of interactions and stability of nitrobenzene derivatives with the PP2 protein involved in the Phloem protein 2 (PP2). From each computational method a prediction model was extracted HQSPR (Atomic, Bonds, Donor/Acceptor, 3D‐QSPR (Steric, Electrostatic, Hydrophobic, Acceptor), and used for the design of X1 and X2 new nitrobenzene derivatives, which showed remarkable stability and better interactions. This study contributes to the understanding of the nitrobenzene derivatives’ toxicity and offers materials towards the understanding of the mechanism of action of these substances in protein–ligand interactions.