AGNpy is an open-source python package designed to model the radiative processes at work
in the jets of active galactic nuclei (AGN) and blazars, whose high energy emission can arise
due to both leptonic and hadronic processes. The former assumes the emission to be dominated
by electrons and positrons, while the latter considers a significant contribution by protons. A
first version of the package allowed a fully leptonic interpretation of the sources spectral energy
distribution (SED). However, the evidence for multi-messenger photon and neutrino emission
from the blazar TXS 0506+056, in 2017, suggested blazars as potential neutrino – and thus
cosmic-ray – sources giving relevance to the role of protons in blazar jets. In this contribution we
present the implementation of the first hadronic processes in AGNpy: synchrotron emission from
accelerated protons and radiation of secondaries from photo-meson interactions, in particular.
The user is now able to model the sources SED also by means of hadronic and lepto-hadronic
processes and to estimate the spectra of the produced muons and neutrinos. We show the software
validation through comparison with other codes and few examples of application. These additions
make agnpy the first open-source software allowing also a lepto-hadronic interpretation of AGNs
emission. The development of hadronic processes in an open-source package opens the modelling
effort to a larger community of astrophysicists and ensures the reproducibility and validation of
results, becoming essential with the amount of open access data that will be provided by the next
generation high-energy observatories.