The purpose of the present paper is to develop $C^1$ Virtual Elements in three dimensions for linear elliptic fourth order problems, motivated by the difficulties that standard conforming Finite Elements encounter in this framework. We focus the presentation on the lowest order case, the generalization to higher orders being briefly provided in the Appendix. The degrees of freedom of the proposed scheme are only 4 per mesh vertex, representing function values and gradient values. Interpolation error estimates for the proposed space are provided, together with a set of numerical tests to validate the method at the practical level.