Molecular dynamics simulations of metalloproteins: A folding study of rubredoxin from Pyrococcus furiosus

The constant increase of computational power has made feasible to investigate the folding mechanism of small proteins using molecular dynamics (MD). Metal-binding proteins (metalloproteins) are usually complicated to model, largely due to the presence of the metal cofactor. Thus, the study of metal-coupled folding is still challenging. In this work, we addressed the folding process of Pyrococcus furiosus rubredoxin (PfRd), a 53-residue protein binding a single iron ion, using different MD methods. Starting from an extended conformation of the polypeptide chain where we preserved the coordination of the metal ion, a classical MD simulation and an extensive accelerated MD run were performed to reconstruct the folding process of the metal-bound protein. For comparison, we simulated also the dynamics of folded PfRd devoid of the metal cofactor (apo-form), starting from the folded structure. For these MD trajectories, we computed various structural and biochemical properties. In addition, we took advantage of available experimental data to quantify the degree to which our simulations sampled conformations close to the native folded state. We observed that the compaction of the hydrophobic core is the main feature driving the folding of the structure. However, we could not reach a fully folded conformation within our trajectories, because of the incomplete removal of the solvent from the core. Altogether, the various MD simulations, including that of the folded apo-form of the protein, suggest that an improvement in the accuracy of the protein force-field is still needed.