Multidimensional hypoelliptic diffusions arise naturally in different fields, for example to model neuronal activity. Estimation in those models is complex because of the degenerate structure of the diffusion coefficient. In this paper we consider hypoelliptic diffusions, given as a solution of two-dimensional stochastic differential equations (SDEs), with the discrete time observations of both coordinates being available on an interval T = n∆ n , with ∆ n the time step between the observations. The estimation is studied in the asymptotic setting, with T → ∞ as ∆ n → 0. We build a consistent estimator of the drift and variance parameters with the help of a discretized log-likelihood of the continuous process. We discuss the difficulties generated by the hypoellipticity and provide a proof of the consistency and the asymptotic normality of the estimator. We test our approach numerically on the hypoelliptic FitzHugh-Nagumo model, which describes the firing mechanism of a neuron.