Atomistic description of the folding of a dimeric protein

Equilibrium molecular dynamics simulations are increasingly being used to describe the folding of individual proteins and protein domains at an atomic level of detail. Isolated protein domains, however, are rarely observed in vivo, where multidomain proteins and multimeric assemblies are far more common. It is clear that the folding of such proteins is often inextricably coupled with the process of dimerization; indeed, many protein monomers and protein domains are not stable in isolation, and fold to their native structures only when stabilized by interactions with other members of a protein complex. Here, we use equilibrium molecular dynamics simulations with an aggregate simulation length of 4 ms to elucidate key aspects of the folding mechanism, and of the associated free-energy surface, of the Top7-CFr dimer, a 114-amino-acid protein homodimer with a mixed α/β structure. In these simulations, we observed a number of folding and unfolding events. Each folding event was characterized by the assembly of two unfolded Top7-CFr monomers to form a stable, folded dimer. We found that the isolated monomer is unstable but that, early in the folding pathway, nascent native structure is stabilized by contacts between the two monomer subunits. These contacts are in some cases native, as in an induced-folding model, and in other cases non-native, as in a fly-casting mechanism. Occasionally, folding by conformational selection, in which both subunits form independently before dimerization, was also observed. Folding then proceeds through the sequential addition of strands to the protein β sheet. Although the long-time-scale relaxation of the folding process can be well described by a single exponential, these simulations reveal the presence of a number of kinetic traps, characterized by structures in which individual strands are added in an incorrect order.