  This paper reports on the determination of $\alpha_S$ from lattice
simulations with 2+1+1 twisted-mass dynamical flavours {\it via} the
computation of the ghost-gluon coupling renormalized in the MOM Taylor scheme.
A high-statistics sample of gauge configurations, used to evaluate the coupling
from ghost and gluon propagators, allows for the appropriate update of previous
results, now performing an improved analysis of data with reduced statistical
errors and the systematical uncertainties under a better control.
