Binding of S. cerevisiae iso-1 cytochrome c and its surface lysine-to-alanine variants to cardiolipin: charge effects and the role of the lipid to protein ratio.

The interaction of cytochrome c with cardiolipin (CL) is a critical step in the initial stages of apoptosis and is mediated by a positively charged region on the protein surface comprising several lysine residues (site A). Here, the interaction of wt S. cerevisiae cytochrome c (ycc) and its K72A/K73A, K72A/K79A, K73A/K79A and K72A/K73A/K79A variants with CL was studied through UV-Vis and MCD spectroscopies at pH 7 and molecular dynamics (MD) simulations, to clarify the role of the mutated lysines. Moreover, the influence of the lipid to protein ratio on the interaction mechanism was investigated using low (0.5-10) and high (5-60) CL/ycc molar ratios, obtained with small and gradual or large and abrupt CL additions, respectively. Although all proteins bind to CL, switching from the native low-spin His/Met-ligated form to a low-spin bis-His conformer and to a high-spin species at larger CL concentrations, the two schemes of CL addition show relevant differences in the CL/ycc molar ratios at which the various conformers appear, due to differences in the interaction mechanism. Extended lipid anchorage and peripheral binding appear to prevail at low and high CL/ycc molar ratios, respectively. Simultaneous deletion of two or three surface positive charges from Site A does not abolish CL binding, but instead increases protein affinity for CL. MD calculations suggest this unexpected behavior results from the mutation-induced severe weakening of the H-bond connecting the Nε of His26 with the backbone oxygen of Glu44, which lowers the conformational stability compared to the wt species, overcoming the decreased surface electrostatic interaction.

Computational evidence support the hypothesis of neuroglobin also acting as an electron transfer species.

Neuroglobin (Ngb) is a recently identified hexa-coordinated globin, expressed in the nervous system of humans. Its physiological role is still debated: one hypothesis is that Ngb serves as an electron transfer (ET) species, possibly by reducing cytochrome c and preventing it to initiate the apoptotic cascade. Here, we use the perturbed matrix method (PMM), a mixed quantum mechanics/molecular dynamics approach, to investigate the redox thermodynamics of two neuroglobins, namely the human Ngb and GLB-6 from invertebrate Caenorhabditis elegans. In particular, we calculate the reduction potential of the two globins, resulting in an excellent agreement with the experimental values, and we predict the reorganization energies, λ, which have not been determined experimentally yet. The calculated λ values match well those reported for known ET proteins and thereby support a potential involvement in vivo of the two globins in ET processes.

Voltage Dependence of Conformational Dynamics and Subconducting States of VDAC-1.

The voltage-dependent anion channel 1 (VDAC-1) is an important protein of the outer mitochondrial membrane that transports energy metabolites and is involved in apoptosis. The available structures of VDAC proteins show a wide ß-stranded barrel pore, with its N-terminal α-helix (N-α) bound to its interior. Electrophysiology experiments revealed that voltage, its polarity, and membrane composition modulate VDAC currents. Experiments with VDAC-1 mutants identified amino acids that regulate the gating process. However, the mechanisms for how these factors regulate VDAC-1, and which changes they trigger in the channel, are still unknown. In this study, molecular dynamics simulations and single-channel experiments of VDAC-1 show agreement for the current-voltage relationships of an "open" channel and they also show several subconducting transient states that are more cation selective in the simulations. We observed voltage-dependent asymmetric distortions of the VDAC-1 barrel and the displacement of particular charged amino acids. We constructed conformational models of the protein voltage response and the pore changes that consistently explain the protein conformations observed at opposite voltage polarities, either in phosphatidylethanolamine or phosphatidylcholine membranes. The submicrosecond VDAC-1 voltage response shows intrinsic structural changes that explain the role of key gating amino acids and support some of the current gating hypotheses. These voltage-dependent protein changes include asymmetric barrel distortion, its interaction with the membrane, and significant displacement of N-α amino acids.

Unambiguous Assignment of Reduction Potentials in Diheme Cytochromes.

Perturbed matrix method calculations are performed on a diheme cytochrome c (DHC) protein, in order to assign previously experimentally detemined reduction potentials (E0) to their corresponding heme groups. Very good agreement between calculated values to experimental data prove that the present approach can be used as a predictive tool of redox thermodynamic properties of multicenter redox proteins in the absence of experimental data, or in synergy with state-of-the art spectroscopic and electrochemical approaches to obtain a detailed picture of electron transfer processes within these complex systems.

Effect of motional restriction on the unfolding properties of a cytochrome c featuring a His/Met-His/His ligation switch.

The K72A/K73H/K79A variant of cytochrome c undergoes a reversible change from a His/Met to a His/His axial heme ligation upon urea-induced unfolding slightly below neutral pH. The unfolded form displays a dramatically lower reduction potential than the folded species along with a pseudo-peroxidase activity. We have studied electrochemically the effects of urea-induced unfolding on the protein electrostatically immobilized on an electrode surface functionalized by means of a negatively charged molecular spacer. The latter mimics the electrostatic interaction with the inner mitochondrial membrane. This behavior has been compared with the unfolding of the same species in solution. This system constitutes a model to decipher the role of the above electrostatic interaction in the unfolding of cytochrome c at physiological pH upon interaction with the membrane component phospholipid cardiolipin in the early stages of the apoptosis cascade. We found that immobilization obstacles protein unfolding due to structural constraints at the interface imposed by protein-SAM interaction.

How the dynamics of the metal-binding loop region controls the acid transition in cupredoxins.

Many reduced cupredoxins undergo a pH-dependent structural rearrangement, triggered by protonation of the His ligand belonging to the C-terminal hydrophobic loop, usually termed the acid transition. At variance with several members of the cupredoxin family, the acid transition is not observed for azurin (AZ). We have addressed this issue by performing molecular dynamics simulations of AZ and four mutants, in which the C-terminal loop has been replaced with those of other cupredoxins or with polyalanine loops. All of the loop mutants undergo the acid transition in the pH range of 4.4-5.5. The main differences between AZ and its loop mutants are the average value of the active site solvent accessible surface area and the extent of its fluctuations with time, together with an altered structure of the water layer around the copper center. Using functional mode analysis, we found that these variations arise from changes in nonbonding interactions in the second coordination sphere of the copper center, resulting from the loop mutation. Our results strengthen the view that the dynamics at the site relevant for function and its surroundings are crucial for protein activity and for metal-containing electron transferases.

The Active Site Loop Modulates the Reorganization Energy of Blue Copper Proteins by Controlling the Dynamic Interplay with Solvent.

Understanding the factors governing the rate of electron transfer processes in proteins is crucial not only to a deeper understanding of redox processes in living organisms but also for the design of efficient devices featuring biological molecules. Here, molecular dynamics simulations performed on native azurin and four chimeric cupredoxins allow for the calculation of the reorganization energy and of structure-related quantities that were used to clarify the molecular determinants to the dynamics/function relationship in blue copper proteins. We find that the dynamics of the small, metal-binding loop region controls the outer-sphere reorganization energy not only by determining the exposure of the active site to solvent but also through the modulation of the redox-dependent rearrangement of the whole protein scaffold and of the surrounding water molecules.