We address the virtual network embedding problem (VNE) which, given a physical (substrate) network and a collection of virtual networks (VNs), calls for an embedding of the most profitable subset of VNs onto the physical substrate, subject to capacity constraints. In practical applications, node and link demands of the different VNs are, typically, uncertain and difficult to know a priori. To face this issue, we first model VNE as a chance-constrained Mixed-Integer Linear Program (MILP) where the uncertain demands are assumed to be random variables. We then propose a $$\varGamma $$ Γ -robust optimization approach to approximate the original chance-constrained formulation, capable of yielding solutions with a large profit that are feasible for almost all the possible realizations of the uncertain demands. To solve larger scale instances, for which the exact approach is computationally too demanding, we propose two MILP-based heuristics: a parametric one, which relies on a parameter setting chosen a priori, and an adaptive one, which does not. We conclude by reporting on extensive computational experiments where the different methods and approaches are compared.