Costantino Zazza

UV-Vis Spectra of the Anticancer Campothecin Family Drugs in Aqueous Solution: Specific Spectroscopic Signatures Unraveled by a Combined Computational and Experimental Study

The ultraviolet-visible absorption spectrum of camptothecin (CPT) has been been recorded in aqueous solution at pH 5.3, where the equilibrium among the different CPT forms is shifted toward the lactonic one. Time-dependent density functional theory (TD-DFT) computations lead to a remarkable reproduction of the experimental spectrum only upon addition of explicit water molecules in interaction with specific moieties of the camptothecin molecule. Molecular dynamics (MD) simulations enforcing boundary periodic conditions for CPT embedded with 865 water molecules, with a force field derived from DFT computations, show that the experimental spectrum is due to the contributions of CPT molecules with different solvation patterns. A similar solvent effect is observed for several CPT derivatives, including the clinically relevant SN-38 and topotecan drugs. The quantitative agreement between TD-DFT/MD computations and experimental data allow us to identify specific spectroscopic signatures diagnostic of the drug environment and to develop procedures that can be used to monitor the drug-DNA/protein interaction.

A single-site mutation (F429H) converts the enzyme CYP 2B4 into a heme oxygenase: a QM/MM study.

The intriguing deactivation of the cytochrome P450 (CYP) 2B4 enzyme induced by mutation of a single residue, Phe429 to His, is explored by quantum mechanical/molecular mechanical calculations of the O-OH bond activation of the (Fe(3+)OOH)(-) intermediate. It is found that the F429H mutant of CYP 2B4 undergoes homolytic instead of heterolytic O-OH bond cleavage. Thus, the mutant acquires the following characteristics of a heme oxygenase enzyme: (a) donation by His429 of an additional NH---S H-bond to the cysteine ligand combined with the presence of the substrate retards the heterolytic cleavage and gives rise to homolytic O-OH cleavage, and (b) the Thr302/water cluster orients nascent OH(•) and ensures efficient meso hydroxylation.

Theoretical modeling of enzyme reactions: the thermodynamics of formation of compound 0 in horseradish peroxidase.

In this paper, by using the perturbed matrix method (PMM) in combination with basic statistical mechanical relations both based on nanosecond time-scale molecular dynamics (MD) simulations, we quantitatively address the thermodynamics of compound 0 (Cpd 0) formation in horseradish peroxidase (HRP) enzyme. Our results, in the same trend of low-temperature experimental data, obtained in cryoenzymology studies indicate that such a reaction can be described essentially as a stepwise spontaneous process: a first step mechanically constrained, strongly exothermic proton transfer from the heme-H2O2 complex to the conserved His42, followed by a solvent-protein relaxation involving a large entropy increase. Critical evaluation of PMM/MD data also reveals the crucial role played by specific residues in the reaction pocket and, more in general, by the conformational fluctuations of the overall environment in physiological conditions.

Theoretical modeling of the valence UV spectra of 1,2,3-triazine and uracil in solution.

Assessment of the perturbed matrix method (PMM) ability in reproducing valence UV absorption spectra is carried out on two model systems: 1,2,3-triazine in methanol solution and uracil in water solution. Results show that even using the simplest definition of the quantum center, i.e. the portion of the system explicitly treated quantum mechanically, PMM provides rather good results. This paper further confirms the possibility of using PMM as a theoretical-computational tool, complementary to other methodologies, for addressing the electronic properties in molecular systems of high complexity.

Oxygen Adsorption on β-Cristobalite Polymorph: Ab Initio Modeling and Semiclassical Time-Dependent Dynamics

The adsorption dynamics of atomic oxygen on a model beta-cristobalite silica surface has been studied by combining ab initio electronic structure calculations with a molecular dynamics semiclassical approach. We have evaluated the interaction potential of atomic and molecular oxygen interacting with an active Si site of a model beta-cristobalite surface by performing DFT electronic structure calculations. As expected, O is strongly chemisorbed, E(b) = 5.57 eV, whereas molecular oxygen can be weakly adsorbed with a high-energy barrier to the adsorption state of approximately 2 eV. The binding energies calculated for silica clusters of different sizes have revealed the local nature of the O,O(2)-silica interaction. Semiclassical collision dynamic calculations show that O is mainly adsorbed in single-bounce collisions, with a smaller probability for adsorption via a multicollision mechanism. The probability for adsorption/desorption (reflected) collisions at the three impact energies is small but not negligible at the higher energy considered in the trajectory calculations, about P(r) = 0.2 at E(kin) = 0.8 eV. The calculations give evidence of a complex multiphonon excitation-deexcitation mechanism underlying the dynamics of stable adsorption and inelastic reflection collisions.

Oxygen Adsorption on β-Quartz Model Surfaces: Some Insights from Density Functional Theory Calculations and Semiclassical Time-Dependent Dynamics

The O/β-quartz interaction is described by combining our time-dependent semiclassical approach to atom-molecule/surface scattering with first-principles electronic structure calculations at the DFT (PBE0) level of accuracy. In particular, the O, O(2) interaction potentials with an on-top Si atom and its nearest O atom both localized over three different silica clusters have been calculated as a function of the oxygen-silica approaching distance. The calculated DFT potential energy surface has been used in semiclassical trajectory calculations to investigate the sticking and inelastic reflection of oxygen atoms from a model β-quartz surface. The collisional mechanism, including the role played by the phonon dynamics, is brought to light and accurate sticking probabilities are calculated at five impact energies in the range [0.05-0.8] eV and T(S) = 1000 K. The different catalytic response of β-quartz and β-cristobobalite to the atomic oxygen flux is also discussed and highlighted.

Boundary condition effects on the dynamic and electric properties of hydration layers.

Water solvation has a central role in several biochemical processes ranging from protein folding to biomolecular recognition and enzyme catalysis. Because of its importance, the structure and dynamics of hydration layers around biological macromolecules have been the targets of a great number of experimental and computational studies. In the present contribution, we have investigated the effects of periodic boundary conditions (PBCs), as used in conjunction with molecular dynamics (MD) simulations, on the dynamic and electric properties of water layers. In particular, we have systematically performed MD simulations of neat water and biomolecules in aqueous solutions by imposing a different external dielectric constant, a generally overlooked parameter in PBC simulations. The effect of the system size has also been addressed. Overall, our results consistently indicate that the dipole moment properties of water layers, and specifically the dipole moment fluctuations and the reorientational correlation functions, can be sensitive to the choice of the external boundary conditions, whereas other molecular properties, such as the self-diffusion coefficient and the reorientational relaxation times, are not affected. We think that our investigation may help to assess appropriate simulation conditions for modeling the aqueous environment of relevant biochemical systems and processes.

Computational study on compound I redox-active species in horseradish peroxydase enzyme: conformational fluctuations and solvation effects.

Molecular dynamics simulations for the compound I species in horseradish peroxidase were carried out over a nanoseconds time-scale. Results indicate that the supramolecular assembly composed of compound I in interaction with highly conserved distal residues (His42 and Arg38) exists in two well-defined conformations basically differing in the local position of the distal histidine (i.e., His42). Furthermore, we observe the presence of a biological channel in the distal side of the heme cavity, between Arg38 and Pro139 residues, that represents a direct link connecting the compound I (CpdI) species to bulk molecules. Our investigation supports the idea that when CpdI is formed, the biological machinery relaxes the local electrostatic forces, opening a structural channel through which an exchange of water molecules with the bulk solvent takes place without any significant kinetic barrier. Interestingly, we also show that the combined effect of enzyme and solvent, modulated by thermal fluctuations, affects the order and the energy difference between the lowest doublet and quartet magnetic states of the CpdI-His42-Arg38 complex. Consequently, when passing from the gas phase to the biological environment, the doublet spin state becomes slightly more stable than the higher spin multiplicity state. The distribution of the perturbed low-lying states energy variation, induced by surrounding fluctuations, yields results rather close to that obtained by other state of the art quantum-mechanics/molecular mechanics calculations, and still in line with previous polarizable continuum calculations.

On the catalytic role of structural fluctuations in enzyme reactions: computational evidence on the formation of compound 0 in horseradish peroxidase

In this study the question as to whether and to what extent horseradish peroxidase (HRP) enzyme structural flexibility affects the free energy barrier for the formation of the key intermediate compound 0 (Cpd0) is addressed by the use of a combined application of molecular dynamics simulations and perturbed matrix method calculations. Results are intriguing and indicate that, within the simulated conditions, free energy profiles are substantially affected by structural fluctuations of the whole surrounding biological environment (i.e. HRP enzyme and solvent). In this respect our results show that the combined effect of enzyme and solvent provides a substantial lowering of the free energy barrier to the formation of Cpd0, with respect to both gas-phase and QM/MM results carried out at a comparable level of theory. A careful inspection of such observations and their general implications in currently employed methodological approaches to the modelling of enzyme reactions, is also discussed.

In Silico Characterization of a Fourfold Magnesium Organometallic Compound in PTCDA Thin Films

In this contribution, using first principles calculations within a density functional theory framework, we report, for the first time, evidence for the formation of a fourfold magnesium organometallic compound upon metal deposition on perylene-3,4,9,10-tetracarboxyl dianhydride (PTCDA) organic semiconductor. Current investigation clearly indicates that in the bulk of the organic crystallographic structure the magnesium atom mainly interacts with three PTCDA molecules. The reactive metal is bound both to carboxylic oxygen atoms of the anhydride-end moieties and to a perylene carbon atom which changes its hybridization state, from sp(2) to sp(3), in the presence of metal impurities. In turn, the analysis of the electronic structure of the reacted system prevalently reveals the formation of four covalent bonds, as a consequence of a weak charge transfer toward the organic material. Such a result confirms the capability of the PTCDA thin films to host metal atoms providing, inside their structural empty channels, a rather accessible and soft chemical environment. Interestingly, in the light of these findings and of previous works, a relationship between first ionization potential of the doping metal and the character of the newly formed chemical bonds is confirmed.