In this study, molecular dynamics simulation has been applied to investigate Heulandite family zeolites (Heulandite and Clinoptilolite differing in Si/Al content) adsorption capacity of Favipiravir acting as both a potential treatment candidate for COVID-19 and the pharmaceutical environmental pollutant using Materials Studio package. The Heulandite zeolite structure was taken from the package database and its substitutions were carried out according to Lowenstein’s and Dempsey’s rules to create Clinoptilolite unit cell. The simulations are performed in the constant-volume constant-temperature (NVT) ensemble. Respectively, adsorption energy values of -83.99 and -68.33 kcal/mol were obtained for the Heulandite and Clinoptilolite-based structures following the adsorption of the initial drug. These values indicate that obtained compounds are suitable adsorbents for removing unwanted Favipiravir drug emitted in the environment. To illustrate the diffusion of the drug within the framework structure, mean-squared displacement (MSD) calculations were utilized, revealing that Favipiravir had a greater propensity to diffuse through the Heulandite pore network in comparison to Clinoptilolite.