Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Syntheses, x-ray structures, photochemistry, redox properties, and DFT calculations of interconvertible fac- and mer-[Mn(SPS)(CO)3] isomers containing a flexible SPS-based pincer ligand

The lithium salt of the anionic SPS pincer ligand composed of a central hypervalent lambda(4)-phosphinine ring bearing two ortho-positioned diphenylphosphine sulfide side arms reacts with [Mn(CO)(5)Br] to give fac-[Mn(SPS)(CO)(3)], This isomer can be converted photochemicaily to mer-[Mn(SPS)(CO)(3)], with a very high quantum yield (0.80 +/- 0.05). The thermal backreaction is slow (taking ca. 8 h at room temperature), in contrast to rapid electrodecatalyzed mer-to-fac isomerization triggered by electrochemical reduction of mer-[Mn(SPS)(CO)(3)]. Both geometric isomers of [Mn(SPS)(CO)(3)] have been characterized by X-ray crystallography. Both isomers show luminescence from a low-lying (IL)-I-3 (SPS-based) excited state. The light emission of fac-[Mn(SPS)(CO)(3)] is largely quenched by the efficient photoisomerization occurring probably from a low-lying Mn-CO dissociative excited state. Density functional theory (DFT) and time-dependent DFT calculations describe the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of fac- and mer-[Mn(CO)(3)(SPS)] as ligand-centered orbitals, largely localized on the phosphinine ring of the SPS pincer ligand. In line with the ligand nature of its frontier orbitals, fac-[Mn(SPS)(CO)(3)] is electrochemically reversibly oxidized and reduced to the corresponding radical cation and anion, respectively. The spectroscopic (electron paramagnetic resonance, IR, and UV-vis) characterization of the radical species provides other evidence for the localization of the redox steps on the SIPS ligand. The smaller HOMO-LUMO energy difference in the case of mer-[Mn(CO)(3)(SPS)], reflected in the electronic absorption and emission spectra, corresponds with its lower oxidation potential compared to that of the fac isomer. The thermodynamic instability of mer-[Mn(CO)(3)(SPS)], confirmed by the DFT calculations, increases upon one-electron reduction and oxidation of the complex.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Surface and Quantum Confinement Effects in ZnO Nanocrystals

ZnO nanocrystals are studied using theoretical calculations based on the density functional theory. The two main effects related to the reduced size of the nanocrystals are investigated: quantum confinement and a large surface:volume ratio. The effects of quantum confinement are studied by saturating the surface dangling bonds of the nanocrystals with hypothetical H atoms. To understand the effects of the surfaces of the nanocrystals, all saturation is removed and the system is relaxed to its minimum energy position. Several different surface motifs are reported, which should be observed experimentally. Spin-polarized calculations are performed in the nonsaturated nanocrystals, leading to different magnetic moments. We propose that this magnetic moment can be responsible for the intrinsic magnetism observed in ZnO nanostructures.

Realistic calculations of carbon-based disordered systems

Carbon nanotubes rank amongst potential candidates for a new family of nanoscopic devices, in particular for sensing applications. At the same time that defects in carbon nanotubes act as binding sites for foreign species, our current level of control over the fabrication process does not allow one to specifically choose where these binding sites will actually be positioned. In this work we present a theoretical framework for accurately calculating the electronic and transport properties of long disordered carbon nanotubes containing a large number of binding sites randomly distributed along a sample. This method combines the accuracy and functionality of ab initio density functional theory to determine the electronic structure with a recursive Green`s functions method. We apply this methodology on the problem of nitrogen-rich carbon nanotubes, first considering different types of defects and then demonstrating how our simulations can help in the field of sensor design by allowing one to compute the transport properties of realistic nanotube devices containing a large number of randomly distributed binding sites.

Study of the carbon dioxide chemical fixation-activation by guanidines

Fixation of CO(2) is one of the most important priorities of the scientific community dedicated to reduce global warming. In this work, we propose new methods for the fixation of CO2 using the guanidine bases tetramethylguanidine (TMG) and 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a]-pyrimidine (TBD). In order to understand the reactions occurring during the CO(2) fixation and release processes, we employed several experimental methods, including solution and solid-state NMR, FTIR, and coupled TGA-FTIR. Quantum mechanical NMR calculations were also carried out. Based on the results obtained, we concluded that CO(2) fixation with both TMG and TBD guanidines is a kinetically reversible process, and the corresponding fixation products have proved to be useful as transcarboxylating compounds. Afterward, CO(2) thermal releasing from this fixation product with TBD was found to be an interesting process for CO(2) capture and isolation purposes. (C) 2008 Elsevier Ltd. All rights reserved.

Crystal structure and theoretical calculations of Julocrotine, a natural product with antileishmanial activity

Julocrotine, N-(2,6-dioxo-l-phenethyl-piperidin-3-yl)-2-methyl-butyramide, is a potent antiproliferative agent against the promastigote and amastigote forms of Leishmania amazonensis (L.). In this work, the crystal structure of Julocrotine was solved by X-ray diffraction, and its geometrical parameters were compared with theoretical calculations at the B3LYP and HF level of theory. IR and NMR spectra also have been obtained and compared with theoretical calculations. IR absorptions calculated with the B3LYP level of theory employed together with the 6-311G+(d,p) basis set, are close to those observed experimentally. Theoretical NMR calculations show little deviation from experimental results. The results show that the theory is in accordance with the experimental data. (C) 2007 Wiley Periodicals, Inc.

Study of the carbon dioxide chemical fixation-activation by guanidines

Fixation of CO2 is one of the most important priorities of the scientific community dedicated to reduce global warming. In this work, we propose new methods for the fixation of CO2 using the guanidine bases tetramethylguanidine (TMG) and 1,3,4,6,7,8-hexahydro-2H-pyrimido[1,2-a]-pyrimidine (TBD). In order to understand the reactions occurring during the CO2 fixation and release processes, we employed several experimental methods, including solution and solid-state NMR, FTIR, and coupled TGA-FTIR. Quantum mechanical NMR calculations were also carried out. Based on the results obtained, we concluded that CO2 fixation with both TMG and TBD guanidines is a kinetically reversible process, and the corresponding fixation products have proved to be useful as transcarboxylating compounds. Afterward, CO2 thermal releasing from this fixation product with TBD was found to be an interesting process for CO2 capture and isolation purposes. (C) 2008 Elsevier Ltd. All rights reserved.

Optical and structural properties of GaAs/GaInP quantum wells grown by Chemical Beam Epitaxy

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

Band structure calculations of InP wurtzite/zinc-blende quantum wells

Semiconductor nanowhiskers (NWs) made of III-V compounds exhibit great potential for technological applications. Controlling the growth conditions, such as temperature and diameter, it is possible to alternate between zinc-blende (ZB) and wurtzite (WZ) crystalline phases, giving origin to the so called polytypism. This effect has great influence in the electronic and optical properties of the system, generating new forms of confinement to the carriers. A theoretical model capable to accurately describe electronic and optical properties in these polytypical nanostructures can be used to study and develop new kinds of nanodevices. In this study, we present the development of a wurtzite/zinc-blende polytypical model to calculate the electronic band structure of nanowhiskers based on group theory concepts and the k.p method. Although the interest is in polytypical superlattices, the proposed model was applied to a single quantum well of InP to study the physics of the wurtzite/zinc-blende polytypism. By the analysis of our results, some trends can be predicted: spatial carriers' separation, predominance of perpendicular polarization (xy plane) in the luminescence spectra, and interband transition blueshifts with strain. Also, a possible range of values for the wurtzite InP spontaneous polarization is suggested. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4767511]

Electrostatic potentials and polarization effects in proton-molecule interactions by means of multipoles from the quantum theory of atoms in molecules

Some atomic multipoles (charges, dipoles and quadrupoles) from the Quantum Theory of Atoms in Molecules (QTAIM) and CHELPG charges are used to investigate interactions between a proton and a molecule (F2, Cl2, BF, AlF, BeO, MgO, LiH, H2CO, NH3, PH3, BF3, and CO2). Calculations were done at the B3LYP/6-311G(3d,3p) level. The main aspect of this work is the investigation of polarization effects over electrostatic potentials and atomic multipoles along a medium to long range of interaction distances. Large electronic charge fluxes and polarization changes are induced by a proton mainly when this positive particle approaches the least electronegative atom of diatomic heteronuclear molecules. The search for simple equations to describe polarization on electrostatic potentials from QTAIM quantities resulted in linear relations with r-4 (r is the interaction distance) for many cases. Moreover, the contribution from atomic dipoles to these potentials is usually the most affected contribution by polarization what reinforces the need for these dipoles to a minimal description of purely electrostatic interactions. Finally, CHELPG charges provide a description of polarization effects on electrostatic potentials that is in disagreement with physical arguments for certain of these molecules. (c) 2012 Wiley Periodicals, Inc.

Persistent currents in quantum rings

In this thesis, three different types of quantum rings arestudied. These are quantum rings with diamagnetic,paramagnetic or spontaneous persistent currents. It turns out that the main observable to characterizequantum rings is the Drude weight. Playing a key role inthis thesis, it will be used to distinguish betweendiamagnetic (positive Drude weight) and paramagnetic(negative Drude weight) ring currents. In most models, theDrude weight is positive. Especially in the thermodynamiclimit, it is positive semi-definite. In certain modelshowever, intuitivelysurprising, a negative Drude weight is found. This rareeffect occurs, e.g., in one-dimensional models with adegenerate ground state in conjunction with the possibilityof Umklapp scattering. One aim of this thesis is to examineone-dimensional quantum rings for the occurrence of anegative Drude weight. It is found, that the sign of theDrude weight can also be negative, if the band structurelacks particle-hole symmetry. The second aim of this thesis is the modeling of quantumrings intrinsically showing a spontaneous persistentcurrent. The construction of the model starts from theextended Hubbard model on a ring threaded by anAharonov-Bohm flux. A feedback term through which thecurrent in the ring can generate magnetic flux is added.Another extension of the Hamiltonian describes the energystored in the internally generated field. This model isevaluated using exact diagonalization and an iterativescheme to find the minima of the free energy. The quantumrings must satisfy two conditions to exhibit a spontaneousorbital magnetic moment: a negative Drude weight and aninductivity above the critical level. The magneticproperties of cyclic conjugated hydrocarbons likebenzene due to electron delocalization [magnetic anisotropy,magnetic susceptibility exaltation, nucleus-independent chemical shift (NICS)]---that have become important criteriafor aromaticity---can be examined using this model. Corrections to the presented calculations are discussed. Themost substantial simplification made in this thesis is theneglect of the Zeeman interaction of the electron spins withthe magnetic field. If a single flux tube threads a quantumring, the Zeeman interaction is zero, but in mostexperiments, this situation is difficult to realize. In themore realistic situation of a homogeneous field, the Zeemaninteraction has to be included, if the electrons have atotal spin component in the direction of the magnetic field,or if the magnetic field is strong.

Ab initio quantum mechanical investigation of structural and chemical-physical properties of selected minerals for minero-petrological, structural ceramic and biomaterial applications

The purpose of this thesis is the atomic-scale simulation of the crystal-chemical and physical (phonon, energetic) properties of some strategically important minerals for structural ceramics, biomedical and petrological applications. These properties affect the thermodynamic stability and rule the mineral-environment interface phenomena, with important economical, (bio)technological, petrological and environmental implications. The minerals of interest belong to the family of phyllosilicates (talc, pyrophyllite and muscovite) and apatite (OHAp), chosen for their importance in industrial and biomedical applications (structural ceramics) and petrophysics. In this thesis work we have applicated quantum mechanics methods, formulas and knowledge to the resolution of mineralogical problems ("Quantum Mineralogy”). The chosen theoretical approach is the Density Functional Theory (DFT), along with periodic boundary conditions to limit the portion of the mineral in analysis to the crystallographic cell and the hybrid functional B3LYP. The crystalline orbitals were simulated by linear combination of Gaussian functions (GTO). The dispersive forces, which are important for the structural determination of phyllosilicates and not properly con-sidered in pure DFT method, have been included by means of a semi-empirical correction. The phonon and the mechanical properties were also calculated. The equation of state, both in athermal conditions and in a wide temperature range, has been obtained by means of variations in the volume of the cell and quasi-harmonic approximation. Some thermo-chemical properties of the minerals (isochoric and isobaric thermal capacity) were calculated, because of their considerable applicative importance. For the first time three-dimensional charts related to these properties at different pressures and temperatures were provided. The hydroxylapatite has been studied from the standpoint of structural and phonon properties for its biotechnological role. In fact, biological apatite represents the inorganic phase of vertebrate hard tissues. Numerous carbonated (hydroxyl)apatite structures were modelled by QM to cover the broadest spectrum of possible biological structural variations to fulfil bioceramics applications.