Laura Zanetti-Polzi

Surface packing determines the redox potential shift of cytochrome c adsorbed on gold

Thermodynamic and dynamic properties of iso-1-cytochrome c covalently bound to a bare gold surface are here investigated by large scale atomistic simulations. The reduction potential of the protein for low and high surface concentrations is calculated showing a good agreement with experimental estimates. The origin of the dependence of the reduction potential on the surface concentration is investigated and is demonstrated to stem from the changing polarizability of the environment surrounding the protein, a mechanism reminiscent of crowding effects. Moreover, structural analyses are performed revealing relevant changes induced by the presence of the electrode on the dynamic properties of cytochrome c. In particular, one of the two cavities previously identified on the protein surface [Bortolotti et al. J. Am. Chem. Soc., 2012, 134, 13670], and that reversibly open in cytochrome c freely diffusing in solution, is found to be deformed when the protein is adsorbed on gold. This modification exemplifies a mechanism that potentially leads to changes in the protein properties by surface-induced modification of its dynamical behavior.

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.

A few key residues determine the high redox potential shift in azurin mutants

The wide range of variability of the reduction potential (E(0)) of blue-copper proteins has been the subject of a large number of studies in the past several years. In particular, a series of azurin mutants have been recently rationally designed tuning E(0) over a very broad range (700 mV) without significantly altering the redox-active site [Marshall et al., Nature, 2009, 462, 113]. This clearly suggests that interactions outside the primary coordination sphere are relevant to determine E(0) in cupredoxins. However, the molecular determinants of the redox potential variability are still undisclosed. Here, by means of atomistic molecular dynamics simulations and hybrid quantum/classical calculations, the mechanisms that determine the E(0) shift of two azurin mutants with high potential shifts are unravelled. The reduction potentials of native azurin and of the mutants are calculated obtaining results in good agreement with the experiments. The analysis of the simulations reveals that only a small number of residues (including non-mutated ones) are relevant in determining the experimentally observed E(0) variation via site-specific, but diverse, mechanisms. These findings open the path to the rational design of new azurin mutants with different E(0).

Modeling triplet flavin-indole electron transfer and interradical dipolar interaction: a perturbative approach

A benchmark biochemical reaction is here theoretically investigated by means of a perturbative approach in order to model each reaction step. The reaction is the flavin-indole electron transfer, involving also a spin-state relaxation of the ionic complex. The whole reaction path is modeled and the kinetics of the process is studied. The dipolar interaction between the two radicals is explicitly considered during the dynamic evolution of the system in order to investigate the proper conditions for the triplet-to-singlet transition to occur.

Structure of the lipodepsipeptide syringomycin e in phospholipids and sodium dodecylsulphate micelle studied by circular dichroism, NMR spectroscopy and molecular dynamics

Syringomycin E (SRE) is a member of a family of lipodepsipeptides that characterize the secondary metabolism of the plant-associated bacteria Pseudomonas syringae pv. syringae. It displays phytotoxic, antifungal and haemolytic activities, due to the membrane interaction and ion channel formation. To gain an insight into the conformation of SRE in the membrane environment, we studied the conformation of SRE bound to SDS micelle, a suitable model for the membrane-bound SRE. In fact, highly similar circular dichroism (CD) spectra were obtained for SRE bound to sodium dodecylsulphate (SDS) and to a phospholipid bilayer, indicating the conformational equivalence of SRE in these two media, at difference with the CD spectrum of SRE in water solution. The structure of SDS-bound SRE was determined by NMR spectroscopy combined with molecular dynamics calculations in octane environment. The results of this study highlight the influence of the interaction with lipids in determining the three-dimensional structure of SRE and provide the basis for further investigations on structural determinants of syringomycin E-membrane interaction.

Alternative Electron Transfer Channels Ensure Ultrafast Deactivation of Light-Induced Excited States in Riboflavin Binding Protein

Flavoproteins, containing flavin chromophores, are enzymes capable of transferring electrons at very high speeds. The ultrafast photoinduced electron-transfer (ET) kinetics of riboflavin binding protein to the excited riboflavin was studied by femtosecond spectroscopy and found to occur within a few hundred femtoseconds [ Zhong and Zewail, Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 11867-11872 ]. This ultrafast kinetics was attributed to the presence of two aromatic rings that could transfer the electron to riboflavin: the side chains of tryptophan 156 and tyrosine 75. However, the underlying ET mechanism remained unclear. Here, using a hybrid quantum mechanical-molecular dynamics approach, we perform ET dynamics simulations taking into account the motion of the protein and the solvent upon ET. This approach reveals that ET occurs via a major reaction channel involving tyrosine 75 (83%) and a minor one involving tryptophan 156 (17%). We also show that the protein environment is designed to ensure the fast quenching of the riboflavin excited state.

Parallel folding pathways of Fip35 WW domain explained by infrared spectra and their computer simulation

We present a calculation of the amide I′ infrared (IR) spectra of the folded, unfolded, and intermediate states of the WW domain Fip35, a model system for β‐sheet folding. Using an all‐atom molecular dynamics simulation in which multiple folding and unfolding events take place we identify six conformational states and then apply perturbed matrix method quantum‐mechanical calculations to determine their amide I′ IR spectra. Our analysis focuses on two states previously identified as Fip35 folding intermediates and suggests that a three‐stranded core similar to the folded state core is the main source of the spectroscopic differences between the two intermediates. In particular, we propose a hypothesis for why folding via one of these intermediates was not experimentally observed by IR T‐jump.

Structural, Thermodynamic, and Kinetic Properties of Gramicidin Analogue GS6 Studied by Molecular Dynamics Simulations and Statistical Mechanics

Gramicidin S (GS) analogues belong to an important class of cyclic peptides, characterized by an antiparallel double‐stranded β‐sheet structure with Type II′ β‐turns. Such compounds can be used as model systems to understand the folding/unfolding process of β‐hairpins and more in general of β‐structures. In the present study, we specifically investigate the folding/unfolding behavior of the hexameric Gramicidin S analogue GS6 by using all‐atoms molecular dynamics (MD) simulations at different temperatures, coupled to a statistical mechanical model based on the Quasi Gaussian Entropy theory. Such an approach permits to describe the structural, thermodynamic, and kinetic properties of the peptide and to quantitatively characterize its folding/unfolding transitions.

On the nature of solvatochromic effect: The riboflavin absorption spectrum as a case study

We present here the calculation of the absorption spectrum of riboflavin in acetonitrile and dimethyl sulfoxide using a hybrid quantum/classical approach, namely the perturbed matrix method, based on quantum mechanical calculations and molecular dynamics simulations. The calculated spectra are compared to the absorption spectrum of riboflavin previously calculated in water and to the experimental spectra obtained in all three solvents. The experimentally observed variations in the absorption spectra upon change of the solvent environment are well reproduced by the calculated spectra. In addition, the nature of the excited states of riboflavin interacting with different solvents is investigated, showing that environment effects determine a recombination of the gas-phase electronic states and that such a recombination is strongly affected by the polarity of the solvent inducing significant changes in the absorption spectra.