Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.

Predicting spinor condensate dynamics from simple principles

We study the spin dynamics of quasi-one-dimensional F=1 condensates both at zero and finite temperatures for arbitrary initial spin configurations. The rich dynamical evolution exhibited by these nonlinear systems is explained by surprisingly simple principles: minimization of energy at zero temperature and maximization of entropy at high temperature. Our analytical results for the homogeneous case are corroborated by numerical simulations for confined condensates in a wide variety of initial conditions. These predictions compare qualitatively well with recent experimental observations and can, therefore, serve as a guidance for ongoing experiments.

Asymmetric stochastic switching driven by intrinsic molecular noise

Low-copy-number molecules are involved in many functions in cells. The intrinsic fluctuations of these numbers can enable stochastic switching between multiple steady states, inducing phenotypic variability. Herein we present a theoretical and computational study based on Master Equations and Fokker-Planck and Langevin descriptions of stochastic switching for a genetic circuit of autoactivation. We show that in this circuit the intrinsic fluctuations arising from low-copy numbers, which are inherently state-dependent, drive asymmetric switching. These theoretical results are consistent with experimental data that have been reported for the bistable system of the gallactose signaling network in yeast. Our study unravels that intrinsic fluctuations, while not required to describe bistability, are fundamental to understand stochastic switching and the dynamical relative stability of multiple states.

Oxygen adsorption on Cu(510)

Tässä työssä on tutkittu kuparin (510)-askelpinnan reaktiivisuutta käyttäen apuna kvanttimekaanisia ab initio laskentamenetelmiä. Tutkimus on toteutettu laskemalla happiatomin adsorptioenergia ja tilatiheys erilaisissa potentiaalisissa adsorptiopaikoissa pinnalla. Myös happimolekyylin adsorptiota ja hajoamista ontarkasteltu laskemalla pintaa lähestyvälle molekyylille potentiaalienergiapintoja. Energiapintojen tuloksia on myös täydennetty kvanttimekaanisilla molekyylidynamiikkalaskuilla. Metallisia askelpintoja pidetään yleisesti sileitä pintoja reaktiivisempina happea kohtaan, johtuen askeleen reunan pienentävästä vaikutuksesta molekyylin hajoamisen tiellä olevaan energiavaliin. On kuitenkin olemassa myös tuloksia, jotka osoittavat hapen tarttumisprosessin olevan hallitseva juuri terassialueella, askeleen reunan sijasta. Tässä työssä on todettu hapen adsorboituvan Cu(510)-pinnalla tehokkaimmin juuri terassilla olevaan hollow-paikkaan. Myös adsorptioenergiat ovat tällä pinnalla pienempiä kuin sileällä (100)-pinnalla. Potentiaalienergiapintojen perusteella Cu(510)-pinnan todetaan myös olevan vähemmän reaktiivinen kuin askelpintojen yleisesti odotetaan olevan, vaikka askeleen reunan todetaankin pienentävän happiatominhajoamisen esteenä olevaa energiavallia.

Asymmetric stochastic switching driven by intrinsic molecular noise

Low-copy-number molecules are involved in many functions in cells. The intrinsic fluctuations of these numbers can enable stochastic switching between multiple steady states, inducing phenotypic variability. Herein we present a theoretical and computational study based on Master Equations and Fokker-Planck and Langevin descriptions of stochastic switching for a genetic circuit of autoactivation. We show that in this circuit the intrinsic fluctuations arising from low-copy numbers, which are inherently state-dependent, drive asymmetric switching. These theoretical results are consistent with experimental data that have been reported for the bistable system of the gallactose signaling network in yeast. Our study unravels that intrinsic fluctuations, while not required to describe bistability, are fundamental to understand stochastic switching and the dynamical relative stability of multiple states.

Um potencial de interação para o estudo de materiais e simulações por dinâmica molecular

The Vashishta-Rahman effective interaction potential, based on the Pauling's concept of "ionic radii", has been successfully employed to investigate structural and dynamical properties of different classes of material. By celebrating Pauling's birth centenary, we review the building up of the Vashishta-Rahman potential and we present molecular-dynamics simulation results for structure and dynamics of superionic materials, chalcogenide glasses and metallic oxides.

Estudo da estrutura e da solvatação do HNP-3, um antibiótico natural, por dinâmica molecular

The structure and hydration of the HNP-3 have been derived from molecular dynamics data using root mean square deviation, radial and energy distributions. Three antiparallel beta sheets were found to be preserved. 15 intramolecular hydrogen bonds were identified together with 36 hydrogen bonds on the backbone and 35 on the side chain atoms. From the point of view of the hydration dynamics, the analysis shows a high solvent accessibility of the monomer and attractive interactions with water molecules.

Dependência das vias de dissociação unimolecular com a energia interna da molécula de éter dimetílico

Internal energy dependence of the competitive unimolecular dissociation channels of dimethyl ether were studied with the statistical RRKM formalism. The C-O and C-H fission reactions and the 1,2-H and 1,3-H shifts, and 1,1-H2 and 1,3-H2 molecular eliminations are discussed as a function of energy dependence of k a(E*), the microcanonical rate constant for production of transition states. C-O fission is the dominant process while reaction channels involving C-H fission, 1,1-H2 and 1,3-H2 elimination and production of MeOH should be competitive at energies around 400 kJ mol-1. The less favorable process is the channel of CH4 formation.

Análise estrutural de ciclodextrinas: um estudo comparativo entre métodos teóricos clássicos e quânticos

In the present work, we analyzed the accuracy of distinct theoretical methods to reproduce the solid state structures of cyclodextrins. The a, b and g-cyclodextrins (CD) were considered and also their hydrates with included water molecules: a-CD.2H2O, b-CD.10H2O and g-CD.12H2O. The geometries were fully optimized using Molecular Mechanics (MM2), semiempirical (AM1 and PM3) and ab initio (HF/3-21G) methods and quantitatively compared with experimental data from X ray diffraction. The results obtained from the classical MM2 method were in best agreement with the experiment. The semiempirical and ab initio structures were also in satisfactory accordance with the experimental data. In general, the PM3 method was found to be more suitable than the AM1 to describe the CD geometries, mainly when the intramolecular hydrogen bonds are considered.

Estrutura e propriedades de elipticinas

The ellipticines constitute a broad class of molecules with antitumor activity. In the present work we analyzed the structure and properties of a series of ellipticine derivatives in the gas phase and in solution using quantum mechanical and Monte Carlo methods. The results showed a good correlation between the solvation energies in water obtained with the continuum model and the Monte Carlo simulation. Molecular descriptors were considered in the development of QSAR models using the DNA association constant (log Kapp) as biological data. The results showed that the DNA binding is dominated by electronic parameters, with small contributions from the molecular volume and area.

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

Implementação computacional do modelo carga-fluxo de carga-fluxo de dipolo para cálculo e interpretação das intensidades do espectro infravermelho

The first computational implementation that automates the procedures involved in the calculation of infrared intensities using the charge-charge flux-dipole flux model is presented. The atomic charges and dipoles from the Quantum Theory of Atoms in Molecules (QTAIM) model was programmed for Morphy98, Gaussian98 and Gaussian03 programs outputs, but for the ChelpG parameters only the Gaussian programs are supported. Results of illustrative but new calculations for the water, ammonia and methane molecules at the MP2/6-311++G(3d,3p) theoretical level, using the ChelpG and QTAIM/Morphy charges and dipoles are presented. These results showed excellent agreement with analytical results obtained directly at the MP2/6-311++G(3d,3p) level of theory.