Parameter estimation and uncertainty quantification in physiological modelling is a vital step towards personalised medicine. Current state-of-the-art in this research area performs parameter optimisation, with very limited uncertainty quantification. This paper demonstrates the advantage of novel sampling methods, applied on a complex biological PDE system of the pulmonary circulation. The aim is to find an efficient and accurate method for the inference and uncertainty quantification of unknown parameters, relevant for disease diagnosis (pulmonary hypertension) from limited and noisy blood pressure data. The data likelihood is expensive to evaluate as it requires solving numerically a system of PDEs. Therefore, having a model that best trades off accuracy and computational efficiency is of uppermost importance. In this study, we employ fast Bayesian methods, namely MCMC algorithms coupled with emulation using Gaussian Processes, to achieve a computational speed-up. We compare the Delayed Rejection Adaptive Metropolis algorithm in a History Matching framework to the delayed acceptance Adaptive Metropolis algorithm. The first algorithm draws samples from the approximate posterior distribution, while the latter is guaranteed to generate samples from the exact posterior distribution asymptotically. In this paper we propose and derive the n-steps ahead delayed acceptance Metropolis-Hastings algorithm, which is a generalisation of the classical 1-step ahead delayed acceptance Metropolis-Hastings. We show the superiority of the n-steps ahead algorithm over the 1-step ahead method.