We present a multiscale coupling method to address contact problems. The components of the model are a molecular dynamics engine, a finite element program and a coupling scheme. We validate the approach on an Hertzian contact. Various measures are provided to highlight limitations and new opportunities in conducting large scale simulations of contact brought by the proposed multiscale approach.