A better prediction of conformational changes of proteins using minimally connected network models

Elastic network models have recently been used for studying low-frequency collective motions of proteins. These models simplify the complexity that arises from normal mode analysis by considering a simplified potential involving a few parameters. Two common parameters in most of the elastic network models are cutoff radius and force constant. Although the latter has been studied extensively and even elaborate new models were introduced, for the former usually an ad-hoc cutoff radius is considered. Moreover, the quality of the network models is usually assessed by evaluating their prediction against experimental B-factors. In this work, we consider various common elastic network models with different cutoff radii and assess them by their ability to predict conformational changes of proteins in complexes from unbound to bound state. This prediction is performed by perturbing the unbound structure using a number of low-frequency normal modes of its network model to optimally fit the bound structure. We evaluated a dataset of 30 proteins with distinct unbound and bound structures using this criterion. The results showed that, opposed to the common calibration process based on B-factors, a meaningful relationship exists between the quality of the prediction and model parameters. It was shown that the cutoff radius has a major role in this prediction and minimally connected network models, which use the shortest cutoff radius for which the network is stable, give the best results. It was also shown that by considering the first ten normal modes, the conformational changes can be predicted by about 25 percent. Hence, the evaluation process was extended to the case of considering the contribution of all normal modes in the prediction. The results indicated that minimally connected network models are superior, despite their simplicity, when any number of modes are considered in the prediction.