Confinamento das moléculas H2 e O2 como modelo para o estudo da ligação ligante-sítio ativo em macromoléculas biológicas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Quantum confinement in hydrogen bond

In this work, the quantum confinement effect is proposed as the cause of the displacement of the vibrational spectrum of molecular groups that involve hydrogen bonds. In this approach, the hydrogen bond imposes a space barrier to hydrogen and constrains its oscillatory motion. We studied the vibrational transitions through the Morse potential, for the NH and OH molecular groups inside macromolecules in situation of confinement (when hydrogen bonding is formed) and nonconfinement (when there is no hydrogen bonding). The energies were obtained through the variational method with the trial wave functions obtained from supersymmetric quantum mechanics formalism. The results indicate that it is possible to distinguish the emission peaks related to the existence of the hydrogen bonds. These analytical results were satisfactorily compared with experimental results obtained from infrared spectroscopy. (c) 2015 Wiley Periodicals, Inc.

Quantum confinement in hydrogen bond of DNA and RNA

The hydrogen bond is a fundamental ingredient to stabilize the DNA and RNA macromolecules. The main contribution of this work is to describe quantitatively this interaction as a consequence of the quantum confinement of the hydrogen. The results for the free and confined system are compared with experimental data. The formalism to compute the energy gap of the vibration motion used to identify the spectrum lines is the Variational Method allied to Supersymmetric Quantum Mechanics.

Low-energy electron collisions with glycine

We report cross sections for elastic electron scattering by gas phase glycine (neutral form), obtained with the Schwinger multichannel method. The present results are the first obtained with a new implementation that combines parallelization with OpenMP directives and pseudopotentials. The position of the well known pi* shape resonance ranged from 2.3 eV to 2.8 eV depending on the polarization model and conformer. For the most stable isomer, the present result (2.4 eV) is in fair agreement with electron transmission spectroscopy assignments (1.93 +/- 0.05 eV) and available calculations. Our results also point out a shape resonance around 9.5 eV in the A' symmetry that would be weakly coupled to vibrations of the hydroxyl group. Since electron attachment to a broad and lower lying sigma* orbital located on the OH bond has been suggested the underlying mechanism leading to dissociative electron attachment at low energies, we sought for a shape resonance around similar to 4 eV. Though we obtained cross sections with the target molecule at the equilibrium geometry and with stretched OH bond lengths, least-squares fits to the calculated eigenphase sums did not point out signatures of this anion state (though, in principle, it could be hidden in the large background). The low energy (similar to 1 eV) integral cross section strongly scales as the bond length is stretched, and this could indicate a virtual state pole, since dipole supported bound states are not expected at the geometries addressed here. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3687345]

Positron scattering from methane

We report on measurements of total cross sections for positron scattering from the fundamental organic molecule methane (CH(4)). The energy range of these measurements was 0.1-50 eV, whereas the energy resolution was similar to 100 meV when our Ni moderator was used and similar to 260 meV when the W moderator was employed. To assist us in interpreting these data, Schwinger multichannel calculations were performed at both static and static plus polarization levels of approximation for elastic positron scattering from 0.001 to 10 eV. These calculations are found to be in quite good qualitative agreement with our measured data, and they clearly educe the crucial role played by the target polarization in the low energy positron-CH(4) scattering dynamics.

Positron scattering from the cyclic ethers oxirane, 1,4-dioxane, and tetrahydropyran

In this paper we report original measurements of total cross sections (TCSs) for positron scattering from the cyclic ethers oxirane (C2H4O), 1,4-dioxane (C4H8O2), and tetrahydropyran (C5H10O). The present experiments focus on the low energy range from similar to 0.2 to 50 eV, with an energy resolution smaller than 300 meV. This study concludes our systematic investigation into TCSs for a class of organic compounds that can be thought of as sub-units or moieties to the nucleotides in living matter, and which as a consequence have become topical for scientists seeking to simulate particle tracks in matter. Note that as TCSs specify the mean free path between collisions in such simulations, they have enjoyed something of a recent renaissance in interest because of that application. For oxirane, we also report original Schwinger multichannel elastic integral cross section (ICS) calculations at the static and static plus polarisation levels, and with and without Born-closure that attempts to account for the permanent dipole moment of C2H4O. Those elastic ICSs are computed for the energy range 0.5-10 eV. To the best of our knowledge, there are no other experimental results or theoretical calculations against which we can compare the present positron TCSs. However, electron TCSs for oxirane (also known as ethylene oxide) and tetrahydropyran do currently exist in the literature and a comparison to them for each species will be presented. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3696378]

Shape resonance spectra of lignin subunits

We report integral cross sections for elastic electron scattering by the lignin subunits phenol, guaiacol, and p-coumaryl alcohol. Our calculations employed the Schwinger multichannel method with pseudopotentials and indicate three to four pi* shape resonances for each of these systems, suggesting that low-energy electrons could efficiently transfer energy into the lignin matrix. We also discuss dissociation mechanisms based on the calculated cross sections, available experimental data, virtual orbital analysis, and the knowledge on electron interactions with biomolecules. Our results point out a physical-chemical basis for electron-driven biomass delignification. The latter would be an essential step for efficient biofuel production from lignocellulosic materials.

Positron collisions with ethene

We present experimental and theoretical cross sections for positron collisions with ethene molecules. The experimental total cross sections (TCSs) were obtained with a linear transmission technique, for energies from 0.1 eV up to 70 eV. The calculations employed the Schwinger multichannel method and were performed in the static plus polarization approximation for energies up to 10 eV. Our calculated elastic cross sections indicate a Ramsauer-Townsend minimum around 2.8 eV and a virtual state, in agreement with previous calculations by da Silva et al. [Phys. Rev. Lett. 77, 1028 (1996)]. We found reasonable agreement between the calculated elastic integral cross section and the measured total cross section below the positronium formation threshold. The present results are also in quite good agreement with available theoretical and experimental data, although for the experiments this is only true for TCSs above about 7 eV.

Electron collisions with the HCOOH...(H2O)n (n=1, 2 and 3) complexes in liquid phase: the influence on the π* shape resonance of formic acid

In this conference we report cross sections for elastic collisions of low-energy electrons with the HCOOH…(H2O)n complexes, with n = 1, 2 and 3. The scattering cross sections were computed with the Schwinger multichannel method [K. Takatsuka and V. McKoy, Phys. Rev. A 24 , 2473 (1981); Phys. Rev. A 30 , 1734 (1984)] with pseudopotentials [M. H. F. Bettega, L. G. Ferreira, and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993)] in the static-exchange and static-exchange plus polarization approximations, for energies from 0.5 eV to 6 eV. We considered some diÆerent hydrogen-bonded structures for the complexes that were generated with classical Monte Carlo simulations [K. Coutinho and S. Canuto, J. Chem. Phys. 113, 9132, (2000)]. The aim of this work is to investigate the effect of the surrounding water molecules on the π* shape resonance of the solute. Previous theoretical and experimental studies carried out in the gas phase reported a π* state for HCOOH at around 1.9 eV. For the n = 1 case and for all complexes, the stabilization of the resonance was observed (it appears at lower energy compared to the value obtained in the gas phase), as reported previously for the CH2O…H2O complexes [T. C. Freitas, M. A. P. Lima, S. Canuto, and M. H. F. Bettega, Phys. Rev. A 80, 062710 (2009)]. This result indicates that the presence of the solvent may affect the processes related to the π* state, such as the molecular dissociation by electron impact. For the n = 2 case we have observed both stabilization and destabilization of the π* resonance, that is associated with the hydrogen bond donor or acceptor role of the water molecules in the complexes. For the n = 3 case, preliminary static-exchange results show the stabilization of the π* state. We propose an explanation of the stabilization/destabilization of the π* state in terms of the polarization of the solute due to the surrounding water molecules and the net charge in the solute.

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

Shape resonance spectrum of cytosine-guanine pairs.

Single and double strand breaks in DNA can be caused by low-energy electrons, the most abundant secondary products of the interaction of ionizing radiation to the biological matter. Attachment of these electrons to biomolecules lead to the formation of transient negative ions (TNIs) [1], often referred to as resonances, a process that may lead to significant vibrational excitation and dissociation. In the present study, we employ the parallel version [2] of the Schwinger Multichannel Method implemented with pseudopotentials [3] to obtain the shape resonance spectrum of cytosine-guanine (CG) pairs, with special attention to π* transient anion states. Recent experimental studies pointed out a quasi-continuum vibrational excitation spectrum for electron collisions against formic acid dimers [4], suggesting that electron attachment into π* valence orbitals could induce proton transfer in these dimers. In addition, our previous studies on the shape resonance spectra of the hydrogen-bonded complexes comprising formic acid and formamide units indicated interesting electron delocalization (localization) effects arising from the presence (absence) of inversion symmetry centers in the complexes [5]. In the present work, we extend the studies on hydrogen-bonded complexes to the CG pair, where localization of ¼¤ anions would be expected, based on the previous results. References [1]. B. Boudaïffa, P. Cloutier, D. Hunting, M. A. Huels, L. Sanche, Science 287, 1658 (2000). [2]. J. S. dos Santos, R. F. da Costa , M. T. do N. Varella, J. Chem. Phys. 136, 084307 (2012). [3]. M. H. F. Bettega, L. G. Ferreira, M. A. P. Lima, Phys. Rev. A 47, 1111 (1993). [4]. M. Allan, Phys. Rev. Lett. 98, 123201 (2007). [5]. T. C. Freitas, S. dA. Sanchez, M. T. do N. Varella, M. H. F. Bettega, Phys. Rev. A 84, 062714 (2011).

Shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil.

We report on the shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil, as obtained from fixed-nuclei elastic scattering calculations performed with the Schwinger multichannel method with pseudopotentials. Our results are in good agreement with the available electron transmission spectroscopy data, and support the existence of three π* resonances in uracil and 5-fluorouracil. As expected, the anion states are more stable in the substituted molecules than in uracil. Since the stabilization is stronger in 5-chlorouracil, the lowest π* resonance in this system becomes a bound anion state. The present results also support the existence of a low-lying σ ∗ CCl shape resonance in 5- chlorouracil. Exploratory calculations performed at selected C–Cl bond lengths suggest that the σ ∗ CCl resonance could couple to the two lowest π* states, giving rise to a very rich dissociation dynamics. These facts would be compatible with the complex branching of the dissociative electron attachment cross sections, even though we cannot discuss any details of the vibration dynamics based only on the present fixed-nuclei results.

Low energy elastic and electronically inelastic electron scattering from biomolecules.

Reactions initiated by collisions with low-energy secondary electrons has been found to be the prominent mechanism toward the radiation damage on living tissues through DNA strand breaks. Now it is widely accepted that during the interaction with these secondary species the selective breaking of chemical bonds is triggered by dissociative electron attachment (DEA), that is, the capture of the incident electron and the formation of temporary negative ion states [1,2,3]. One of the approaches largely used toward a deeper understanding of the radiation damage to DNA is through modeling of DEA with its basic constituents (nucleotide bases, sugar and other subunits). We have tried to simplify this approach and attempt to make it comprehensible at a more fundamental level by looking at even simple molecules. Studies involving organic systems such as carboxylic acids, alcohols and simple ¯ve-membered heterocyclic compounds are taken as starting points for these understanding. In the present study we investigate the role played by elastic scattering and electronic excitation of molecules on electron-driven chemical processes. Special attention is focused on the analysis of the in°uence of polarization and multichannel coupling e®ects on the magnitude of elastic and electronically inelastic cross-sections. Our aim is also to investigate the existence of resonances in the elastic and electronically inelastic channels as well as to characterize them with respect to its type (shape, core-excited or Feshbach), symmetry and position. The relevance of these issues is evaluated within the context of possible applications for the modeling of discharge environments and implications in the understanding of mutagenic rupture of DNA chains. The scattering calculations were carried out with the Schwinger multichannel method (SMC) [4] and its implementation with pseudopotentials (SMCPP) [5] at di®erent levels of approximation for impact energies ranging from 0.5 eV to 30 eV. References [1] B. Boudai®a, P. Cloutier, D. Hunting, M. A. Huels and L. Sanche, Science 287, 1658 (2000). [2] X. Pan, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 90, 208102 (2003). [3] F. Martin, P. D. Burrow, Z. Cai, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 93, 068101 (2004). [4] K. Takatsuka and V. McKoy, Phys. Rev. A 24, 2437 (1981); ibid. Phys. Rev. A 30, 1734 (1984). [5] M. H. F. Bettega, L. G. Ferreira and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993).

Electron Interactions with Disulfide Bridges

Because of its electronic properties, sulfur plays a major role in a variety of metabolic processes and, more in general, in the chemistry of life. In particular, S-S bridges between cysteines are present in the amino acid backbone of proteins. Protein disulfur radical anions may decay following different paths through competing intra and intermolecular routes, including bond cleavage, disproportionation, protein-protein cross linking, and electron transfer. Indeed, mass spectrometry ECD (electron capture dissociation massspectroscopy) studies have shown that capture of low-energy (<0.2 eV) electrons by multiply protonated proteins is followed by dissociation of S-S bonds holding two peptide chains together. In view of the importance of organic sulfur chemistry, we report on electron interactions with disulphide bridges. To study these interactions we used as prototypes the molecules dimethyl sulfide [(CH3)2S] and dimethyl disulfide [(H3C)S2(CH3)]. We seek to better understand the electron-induced cleavage of the disulfide bond. To explore dissociative processes we performed electron scattering calculations with the Schwinger Multichannel Method with pseudopotentials (SMCPP), recently parallelized with OpenMP directives and optimized with subroutines for linear algebra (BLAS) and LAPACK routines. Elastic cross sections obtained for different S-S bond lengths indicate stabilization of the anion formed by electron attachment to a σ*SS antibonding orbital, such that dissociation would be expected.

Horn-schunck optical flow with a multi-scale strategy

[EN] The seminal work of Horn and Schunck [8] is the first variational method for optical flow estimation. It introduced a novel framework where the optical flow is computed as the solution of a minimization problem. From the assumption that pixel intensities do not change over time, the optical flow constraint equation is derived. This equation relates the optical flow with the derivatives of the image. There are infinitely many vector fields that satisfy the optical flow constraint, thus the problem is ill-posed. To overcome this problem, Horn and Schunck introduced an additional regularity condition that restricts the possible solutions. Their method minimizes both the optical flow constraint and the magnitude of the variations of the flow field, producing smooth vector fields. One of the limitations of this method is that, typically, it can only estimate small motions. In the presence of large displacements, this method fails when the gradient of the image is not smooth enough. In this work, we describe an implementation of the original Horn and Schunck method and also introduce a multi-scale strategy in order to deal with larger displacements. For this multi-scale strategy, we create a pyramidal structure of downsampled images and change the optical flow constraint equation with a nonlinear formulation. In order to tackle this nonlinear formula, we linearize it and solve the method iteratively in each scale. In this sense, there are two common approaches: one that computes the motion increment in the iterations, like in ; or the one we follow, that computes the full flow during the iterations, like in. The solutions are incrementally refined ower the scales. This pyramidal structure is a standard tool in many optical flow methods.