This paper presents a novel method for predicting the failure probability of structural or mechanical systems subjected to random loads and material properties involving multiple design points. The method involves Multicut High Dimensional Model Representation (Multicut-HDMR) technique in conjunction with moving least squares to approximate the original implicit limit state/performance function with an explicit function. Depending on the order chosen sometimes truncated Cut-HDMR expansion is unable to approximate the original implicit limit state/performance function when multiple design points exist on the limit state/performance function or when the problem domain is large. Multicut-HDMR addresses this problem by using multiple reference points to improve accuracy of the approximate limit state/performance function. Numerical examples show the accuracy and efficiency of the proposed approach in estimating the failure probability.