Molecular dynamics simulation of the effect of heat on the conformation of bovine β-lactoglobulin A: A comparison of conventional and accelerated methods

Three molecular dynamics (MD) simulation methods are used to follow the thermal unfolding of bovine beta-lactoglobulin (beta-lac). The methods used, classical MD simulation at different temperatures in the range 300-500 K, essential dynamics and replica exchange MD, were chosen to give a range of conventional and accelerated methods. At 350 K, just above the experimentally determined denaturation temperature of beta-lac only small changes to the tertiary fold and secondary structure were seen during a 110 ns simulation. At 400 K and 450 K more unfolding was observed, but it was not until the temperature was increased to 500 K that substantial disruption to the protein structure was seen. For the time that the heated simulated beta-lac molecules occupy the same conformation space, they appear to follow similar unfolding pathways. ED simulations biased along the first 5 eigenvectors (determined at 350 K) give rise to conformations that have a more elongated tertiary fold compared to heated simulations, and do not follow similar unfolding pathways to the heated simulations. The REMD simulation, which is the equivalent of an approximately 0.7 mu s simulation at 350 K, shows a small degree of tertiary structure unfolding, but very little secondary structure change. The results are discussed in terms of the structural changes that have been observed to occur experimentally in b-lac. (C) 2012 Elsevier Ltd. All rights reserved.