Parallel surrogate-based optimization using Mesh Adaptive Direct Search