a b s t r a c tThis article presents an overview of the simulation studies of the behaviour of multihaem cytochromes using theoretical/computational methodologies, with an emphasis on cytochrome c 3 . It starts with the first studies using rigid molecules and continuum electrostatic models, where protonation and redox events were treated as independent. The gradual addition of physical details is then described, from the inclusion of proton isomerism, to the proper treatment of the thermodynamics of electron-proton coupling, to the explicit inclusion of the solvent and protein structural reorganization into the models, culminating with the method for molecular dynamics simulations at constant pH and reduction potential, where the solvation, conformational, protonation and redox features are all simulated in a fully integrated and coupled way. We end with a discussion of the strategies used to study the interaction between multihaem cytochromes, taking into account the further coupling effect introduced by the molecular association.Ó 2011 Federation of European Biochemical Societies. Published by Elsevier B.V.
ForewordAntónio V. Xavier dedicated most of his scientific life to the tetrahaem cytochrome c 3 , and he was the reason why we started to study this system, back in the middle of the nineties. The thermodynamic characterization of the redox behaviour of this small, but feature rich cytochrome by himself and co-workers (see [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20] for some examples) opened thought provoking questions concerning multi-electron transfer and the coupling between electron and proton transfer. Using this thermodynamic data and the available X-ray and NMR structures, we started to implement computational models in order to understand these phenomena at the molecular level. This wealth of experimental data was of outmost importance for these modelling efforts, since it allowed us to design models and methodologies with our feet firmly stepped on the ground, not necessarily to predict exact numbers, but rather to explain phenomena. This small paper is a summary of this journey through this scientific research that still goes on, and a tribute to the man that prompt all this, one spring morning in his office in Oeiras, and many more mornings in the years that followed.Cytochrome c 3 is a paradigmatic multihaem protein, having four c-type haems with bishistidinyl-axial coordination packed together in a small protein core, constituted by little more than 100 residues (Fig. 1). Its biological function is to serve as an electron exchanger protein in sulphate reducing bacteria, being the physiological partner of periplasmatic hydrogenases [21] and other redox systems in this compartment (see [22,23] for reviews). Multiple X-ray and NMR structures exist for this protein, from different species and different redox and pH states. This makes it the prototypal multihaem protein, and one of the best characterised in its class. Several factors contribute it to be a good model for s...