Nowadays, in the Scientific Machine Learning (SML) research field, the traditional machine learning (ML) tools and scientific computing approaches are fruitfully intersected for solving problems modelled by Partial Differential Equations (PDEs) in science and engineering applications. Challenging SML methodologies are the new computational paradigms named Physics-Informed Neural Networks (PINNs). PINN has revolutionized the classical adoption of ML in scientific computing, representing a novel class of promising algorithms where the learning process is constrained to satisfy known physical laws described by differential equations. In this paper, we propose a PINN-based computational study to deal with a non-linear partial differential equations system. In particular, using this approach, we solve the Gray-Scott model, a reaction–diffusion system that involves an irreversible chemical reaction between two reactants. In the unstable region of the model, we consider some a priori information related to dynamical behaviors, i. e. a supervised approach that relies on a finite difference method (FDM). Finally, simulation results show that PINNs can successfully provide an approximated Grey-Scott system solution, reproducing the characteristic Turing patterns for different parameter configurations.