We propose a Bayesian nonparametrics method, including algorithm for posterior computation, for the sparse regression problem. Our method applies in a general setting, when there are direct or indirect noisy observations of the signal. We try to make a wide focus on smoothness properties and sparsity of the approximate. As an example, we consider the ill-posed inverse problem of Quantum Homodyne Tomography.