A molecular dynamics simulation study of the solvent isotope effect on copper plastocyanin

The effect of heavy water on the structure and dynamics of copper plastocyanin as well as on some aspects of the solvent dynamics at the protein-solvent interfacial region have been investigated by molecular dynamics simulation. The simulated system has been analyzed in terms of the atomic root mean square deviation and fluctuations, intraprotein H-bond pattern, dynamical cross-correlation map and the results have been compared with those previously obtained for plastocyanin in H2O (Ciocchetti et al. Biophys. Chem. 69 (1997), 185-198). The simulated plastocyanin structure in the two solvents, averaging 1 ns, is very similar along the beta-structure regions, while the most significant differences are registered, analogous to the turns and the regions likely involved in the electron transfer pathway. Moreover, plastocyanin in D2O shows an increase in the number of both the intraprotein H-bonds and the residues involved in correlated motions. An analysis of the protein-solvent coupling evidenced that D2O makes the H-bond formation more difficult with the solvent molecules for positively charged and polar residues, while an opposite trend is observed for negatively charged residues. On the other hand, the frequency of exchange of the solvent molecules involved in the protein-solvent H-bond formation is significantly depressed in D2O. The results are discussed also in connection with protein functionality and briefly with some experimental results connected with the thermostability of proteins in D2O.