This article proposes a new method for the resolution of parameterized coupled systems of equations, such as typical in multidisciplinary analysis (MDA).It is based on the use of disciplinary surrogate models within a multi-query context. The main idea is to replace the costly disciplinary solvers of the MDA by Proper Orthogonal Decomposition and Interpolation models. The main challenge we address is the high-dimensional coupling variables whose ranges are unknown. To overcome this issue, a training strategy is developed by uncoupling the disciplinary solvers from the MDA context. This new surrogate MDA called Disciplinary Proper Orthogonal Decomposition and Interpolation is iteratively enriched to solve the analysis with an estimation of the error made by the surrogate models. The disciplinary surrogate model which has the most influence on this error is determined by a sensitivity analysis and thus enriched. This approach allows to uncouple the disciplinary solvers during the training and enrichment phases. The approach is applied to various aeroelastic problems of an aircraft wing and allows to reduce by a factor of 5 the mean number of disciplinary solver calls needed for the resolution of the MDA.