Biomedical subjects

Ivano Tavernelli

Publications and source records attributed to Ivano Tavernelli.

10 recordsLinked to original sources

PubMed2006

Microsolvation effects on the excited-state dynamics of protonated tryptophan.

To better understand the complex photophysics of the amino acid tryptophan, which is widely used as a probe of protein structure and dynamics, we have measured electronic spectra of protonated, gas-phase tryptophan solvated with a controlled number of water molecules and cooled to approximately 10 K. We observe that, even at this temperature, the bare molecule exhibits a broad electronic spectrum, implying ultrafast, nonradiative decay of the excited state. Surprisingly, the addition of two water molecules sufficiently lengthens the excited-state lifetime that we obtain a fully vibrationally resolved electronic spectrum. Quantum chemical calculations at the RI-CC2/aug-cc-pVDZ level, together with TDDFT/pw based first-principles MD simulations of the excited-state dynamics, clearly demonstrate how interactions with water destabilize the photodissociative states and increase the excited-state lifetime.

Models, Molecular↗

PubMed2006

Role of protein frame and solvent for the redox properties of azurin from Pseudomonas aeruginosa.

We have coupled hybrid quantum mechanics (density functional theory; Car-Parrinello)/molecular mechanics molecular dynamics simulations to a grand-canonical scheme, to calculate the in situ redox potential of the Cu(2+) + e(-) --> Cu(+) half reaction in azurin from Pseudomonas aeruginosa. An accurate description at atomistic level of the environment surrounding the metal-binding site and finite-temperature fluctuations of the protein structure are both essential for a correct quantitative description of the electronic properties of this system. We report a redox potential shift with respect to copper in water of 0.2 eV (experimental 0.16 eV) and a reorganization free energy lambda = 0.76 eV (experimental 0.6-0.8 eV). The electrostatic field of the protein plays a crucial role in fine tuning the redox potential and determining the structure of the solvent. The inner-sphere contribution to the reorganization energy is negligible. The overall small value is mainly due to solvent rearrangement at the protein surface.

Azurin↗

PubMed2006

Diabatic free energy curves and coordination fluctuations for the aqueous Ag+Ag2+ redox couple: a biased Born-Oppenheimer molecular dynamics investigation.

Biased Born-Oppenheimer molecular dynamics simulations are performed to compute redox potential and free energy curves for the redox half reaction Ag(+)-->Ag(2+)+e(-) in aqueous solution. The potential energy surfaces of reactant and product state are linearly coupled and the system transferred from the reduced state to the oxidized state by variation of the coupling parameter from 0 to 1. The redox potential is obtained by thermodynamic integration of the average ionization energy of Ag(+). Diabatic free energy curves of reduced (R) and oxidized (O) states are obtained to good statistical accuracy by reweighting and combining the set of biased distributions of the ionization energy. The diabatic free energy curves of Ag(+) and Ag(2+) are parabolic over a wide range of the reaction coordinate in agreement with the linear response assumption that underlies Marcus theory. However, we observe deviations from parabolic behavior in the equilibrium region of Ag(+) and find different values for the reorganization free energy of R (1.4 eV) and O (0.9 eV). The computed reorganization free energy of Ag(2+) is in good agreement with the experimental estimate of 0.9-1.2 eV obtained from photoelectron spectroscopy. As suggested by our calculations, the moderate deviation from linear response behavior found for Ag(+) is likely related to the highly fluxional solvation shell of this ion, which exhibits water exchange reactions on the picosecond time scale of the present molecular dynamics simulation.

Journal Article↗

PubMed2005

Density-functional molecular-dynamics study of the redox reactions of two anionic, aqueous transition-metal complexes.

The thermochemistry of the RuO(4)(2-)+MnO(4)(-)-->RuO(4)(-)+MnO(4)(2-) redox reaction in aqueous solution is studied by separate density-functional-based ab initio molecular-dynamics simulations of the component half reactions RuO(4)(2-)-->RuO(4)(-)+e(-) and MnO(4)(2-)-->MnO(4)(-)+e(-). We compare the results of a recently developed grand-canonical method for the computation of oxidation free energies to the predictions by the energy-gap relations of the Marcus theory that can be assumed to apply to these reactions. The calculated redox potentials are in good agreement. The subtraction of the half-reaction free energies gives an estimate of the free energy of the full reaction. The result obtained from the grand-canonical method is -0.4 eV, while the application of the Marcus theory gives -0.3 eV. These should be compared to the experimental value of 0.0 eV. Size effects, in response to increasing the number of water molecules in the periodic model system from 30 to 48, are found to be small ( approximately 0.1 eV). The link to the Marcus theory also has enabled us to compute reorganization free energies for oxidation. For both the MnO(4)(2-) and RuO(4)(2-) redox reactions we find the same reorganization free energy of 0.8 eV (1.0 eV in the larger system). The results for the free energies and further analysis of solvation and electronic structure confirm that these two tetrahedral oxoanions show very similar behavior in solution in spite of the central transition-metal atoms occupying a different row and column in the periodic table.

Journal Article↗

PubMed2005

Variational optimization of effective atom centered potentials for molecular properties.

In plane wave based electronic structure calculations the interaction of core and valence electrons is usually represented by atomic effective core potentials. They are constructed in such a way that the shape of the atomic valence orbitals outside a certain core radius is reproduced correctly with respect to the corresponding all-electron calculations. Here we present a method which, in conjunction with density functional perturbation theory, allows to optimize effective core potentials in order to reproduce ground-state molecular properties from arbitrarily accurate reference calculations within standard density functional calculations. We demonstrate the wide range of possible applications in theoretical chemistry of such optimized effective core potentials (OECPs) by means of two examples. We first use OECPs to tackle the link atom problem in quantum mechanics/molecular mechanics (QM/MM) schemes proposing a fully automatized procedure for the design of link OECPs, which are designed in such a way that they minimally perturb the electronic structure in the QM region. In the second application, we use OECPs in two sample molecules (water and acetic acid) such as to reproduce electronic densities and derived molecular properties of hybrid (B3LYP) quality within general gradient approximated (BLYP) density functional calculations.

Journal Article↗

PubMed2004

Optimization of effective atom centered potentials for london dispersion forces in density functional theory.

We add an effective atom-centered nonlocal term to the exchange-correlation potential in order to cure the lack of London dispersion forces in standard density functional theory. Calibration of this long-range correction is performed using density functional perturbation theory and an arbitrary reference. Without any prior assignment of types and structures of molecular fragments, our corrected generalized gradient approximation density functional theory calculations yield correct equilibrium geometries and dissociation energies of argon-argon, benzene-benzene, graphite-graphite, and argon-benzene complexes.

Biopolymers↗

PubMed2004

Electronic structure and solvation of copper and silver ions: a theoretical picture of a model aqueous redox reaction.

Electronic states and solvation of Cu and Ag aqua ions are investigated by comparing the Cu(2+) + e(-)--> Cu(+) and Ag(2+) + e(-)--> Ag(+) redox reactions using density functional-based computational methods. The coordination number of aqueous Cu(2+) is found to fluctuate between 5 and 6 and reduces to 2 for Cu(+), which forms a tightly bound linear dihydrate. Reduction of Ag(2+) changes the coordination number from 5 to 4. The energetics of the oxidation reactions is analyzed by comparing vertical ionization potentials, relaxation energies, and vertical electron affinities. The model is validated by a computation of the free energy of the full redox reaction Ag(2+) + Cu(+) --> Ag(+) + Cu(2+). Investigation of the one-electron states shows that the redox active frontier orbitals are confined to the energy gap between occupied and empty states of the pure solvent and localized on the metal ion hydration complex. The effect of solvent fluctuations on the electronic states is highlighted in a computation of the UV absorption spectrum of Cu(+) and Ag(+).

Journal Article↗

PubMed2003

Dynamics of RNase-A and S-protein: a molecular dynamics simulation of the transition toward a folding intermediate.

The description at atomic level of protein folding is an ambitious goal in biophysics, particularly because of the difficulty in obtaining structural information on unfolded states. Computer simulations can contribute in achieving this goal. Here we report the results of a 10-ns comparative simulation on bovine ribonuclease A and its S-protein, obtained by removal from the native molecule of the first 20 residues, the so-called S-peptide. The atomic trajectories have been analyzed by standard procedures and by applying concepts previously developed for disordered systems. Furthermore, we used a novel approach, described in the preceding paper, to represent graphically the energy landscape of the simulated systems. Relative to RNase-A, the S-protein, while largely maintaining its structural organization, displays an increased structural flexibility, it gains ergodicity and its core loses order, thus indicating that the removal of the S-peptide from ribonuclease A triggers the transition to a folding intermediate with reduced compactness. This finding also has biochemical relevance since the S-protein is recognized as not properly folded by the machinery responsible for the control of the folding quality in the endoplasmic reticulum.

Computer Simulation↗

PubMed2003

Protein dynamics, thermal stability, and free-energy landscapes: a molecular dynamics investigation.

Proteins have a complex free-energy landscape because of their rich topology and the nature of their nonbonded interaction potential. This has important consequences because the roughness of the landscape affects the ease with which a chain folds and also determines the dynamic behavior of the folded structure, thus influencing its functional and stability properties. A detailed description of the free-energy landscape is therefore of paramount importance for a quantitative understanding of the relationships between structure, dynamics, stability, and functional behavior of proteins. The free-energy landscape of a protein is a high-dimensional hypersurface, difficult to rationalize. Therefore, achieving its detailed graphical representation in a way that goes beyond the familiar funnel-like free-energy model is still a big challenge. We describe here an approach based on global structural parameters that allows a two-dimensional representation of the free-energy landscape from simulated atomic trajectories. As shown in this and in the accompanying article, our representation of the landscape, combined with other conformational analyses, provides valuable information on its roughness and on how atomic trajectories evolve with time.

Algorithms↗

PubMed2002

Ab initio molecular dynamics for molecules with variable numbers of electrons.

The ab initio molecular dynamics method is extended to treat exchange of electrons between molecules and a reservoir at fixed chemical potential. The method is based on a rigorously grand-canonical density functional approach using separate potential energy surfaces for each oxidation state. It is shown that the resulting discontinuous dependency of excess charge on chemical potential is consistent with the statistical thermodynamics of equilibrium gas-phase reactions. The method is illustrated by an application to the adiabatic redox dynamics of an aniline molecule.

Journal Article↗