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.

Adsorption of quadrupolar, diatomic nitrogen onto graphitized thermal carbon black and in slit-shaped carbon pores. Effects of surface mediation

In this paper, we study the effect of solid surface mediation on the intermolecular potential energy of nitrogen, and its impact on the adsorption of nitrogen on a graphitized carbon black surface and in carbon slit-shaped pores. This effect arises from the lower effective interaction potential energy between two particles close to the surface compared to the potential energy of the same two particles when they are far away from the surface. A simple equation is proposed to calculate the reduction factor and this is used in the Grand Canonical Monte Carlo (GCMC) simulation of nitrogen adsorption on graphitized thermal carbon black. With this modification, the GCMC simulation results agree extremely well with the experimental data over a wide range of pressure; the simulation results with the original potential energy (i.e. no surface mediation) give rise to a shoulder in the neighbourhood of monolayer coverage and a significant over-prediction of the second and higher layer coverages. The influence of this surface mediation on the dependence of the pore-filling pressure on the pore width is also studied. It is shown that such surface mediation has a significant effect on the pore-filling pressure. This implies that the use of the local isotherms obtained from the potential model without surface mediation could give rise to a serious error in the determination of the pore-size distribution.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Manipulating double-decker molecules at the liquid-solid interface

We have used a scanning tunneling microscope to manipulate heteroleptic phthalocyaninato, naphthalocyaninato, porphyrinato double-decker molecules at the liquid/solid interface between 1-phenyloctane solvent and graphite. We employed nano-grafting of phthalocyanines with eight octyl chains to place these molecules into a matrix of heteroleptic double-decker molecules; the overlayer structure is epitaxial on graphite. We have also used nano-grafting to place double-decker molecules in matrices of single-layer phthalocyanines with octyl chains. Rectangular scans with a scanning tunneling microscope at low bias voltage resulted in the removal of the adsorbed doubledecker molecular layer and substituted the double-decker molecules with bilayer-stacked phthalocyanines from phenyloctane solution. Single heteroleptic double-decker molecules with lutetium sandwiched between naphthalocyanine and octaethylporphyrin were decomposed with voltage pulses from the probe tip; the top octaethylporphyrin ligand was removed and the bottom naphthalocyanine ligand remained on the surface. A domain of decomposed molecules was formed within the double-decker molecular domain, and the boundary of the decomposed molecular domain self-cured to become rectangular. We demonstrated a molecular “sliding block puzzle” with cascades of double-decker molecules on the graphite surface.

Molecules, morphology, and ecology indicate a recent, amphibious ancestry for echidnas

The semiaquatic platypus and terrestrial echidnas (spiny anteaters) are the only living egg-laying mammals (monotremes). The fossil record has provided few clues as to their origins and the evolution of their ecological specializations; however, recent reassignment of the Early Cretaceous Teinolophos and Steropodon to the platypus lineage implies that platypuses and echidnas diverged >112.5 million years ago, reinforcing the notion of monotremes as living fossils. This placement is based primarily on characters related to a single feature, the enlarged mandibular canal, which supplies blood vessels and dense electrosensory receptors to the platypus bill. Our reevaluation of the morphological data instead groups platypus and echidnas to the exclusion of Teinolophos and Steropodon and suggests that an enlarged mandibular canal is ancestral for monotremes (partly reversed in echidnas, in association with general mandibular reduction). A multigene evaluation of the echidna–platypus divergence using both a relaxed molecular clock and direct fossil calibrations reveals a recent split of 19–48 million years ago. Platypus-like monotremes (Monotrematum) predate this divergence, indicating that echidnas had aquatically foraging ancestors that reinvaded terrestrial ecosystems. This ecological shift and the associated radiation of echidnas represent a recent expansion of niche space despite potential competition from marsupials. Monotremes might have survived the invasion of marsupials into Australasia by exploiting ecological niches in which marsupials are restricted by their reproductive mode. Morphology, ecology, and molecular biology together indicate that Teinolophos and Steropodon are basal monotremes rather than platypus relatives, and that living monotremes are a relatively recent radiation.

Diversity of dicotyledenous-infecting geminiviruses and their associated DNA molecules in Southern Africa, including the South-west Indian Ocean Islands

The family Geminiviridae comprises a group of plant-infecting circular ssDNA viruses that severely constrain agricultural production throughout the temperate regions of the world, and are a particularly serious threat to food security in sub-Saharan Africa. While geminiviruses exhibit considerable diversity in terms of their nucleotide sequences, genome structures, host ranges and insect vectors, the best characterised and economically most important of these viruses are those in the genus Begomovirus. Whereas begomoviruses are generally considered to be either monopartite (one ssDNA component) or bipartite (two circular ssDNA components called DNA-A and DNA-B), many apparently monopartite begomoviruses are associated with additional subviral ssDNA satellite components, called alpha- (DNA-αs) or betasatellites (DNA-βs). Additionally, subgenomic molecules, also known as defective interfering (DIs) DNAs that are usually derived from the parent helper virus through deletions of parts of its genome, are also associated with bipartite and monopartite begomoviruses. The past three decades have witnessed the emergence and diversification of various new begomoviral species and associated DI DNAs, in southern Africa, East Africa, and proximal Indian Ocean islands, which today threaten important vegetable and commercial crops such as, tobacco, cassava, tomato, sweet potato, and beans. This review aims to describe what is known about these viruses and their impacts on sustainable production in this sensitive region of the world. © 2012 by the authors licensee MDPI, Basel, Switzerland.

Chimaeric HIV-1 subtype C Gag molecules with large in-frame C-terminal polypeptide fusions form virus-like particles

HIV-1 Pr55 Gag virus-like particles (VLPs) are strong immunogens with potential as candidate HIV vaccines. VLP immunogenicity can be broadened by making chimaeric Gag molecules: however, VLPs incorporating polypeptides longer than 200 aa fused in frame with Gag have not yet been reported. We constructed a range of gag-derived genes encoding in-frame C-terminal fusions of myristoylation-competent native Pr55Gag and p6-truncated Gag (Pr50Gag) to test the effects of polypeptide length and sequence on VLP formation and morphology, in an insect cell expression system. Fused sequences included a modified reverse transcriptase-Tat-Nef fusion polypeptide (RTTN, 778 aa), and truncated versions of RTTN ranging from 113 aa to 450 aa. Baculovirus-expressed chimaeric proteins were examined by western blot and electron microscopy. All chimaeras formed VLPs which could be purified by sucrose gradient centrifugation. VLP diameter increased with protein MW, from ∼100 nm for Pr55Gag to ∼250 nm for GagRTTN. The presence or absence of the Gag p6 region did not obviously affect VLP formation or appearance. GagRT chimaeric particles were successfully used in mice to boost T-cell responses to Gag and RT that were elicited by a DNA vaccine encoding a GagRTTN polypeptide, indicating the potential of such chimaeras to be used as candidate HIV vaccines. © 2008 Elsevier B.V. All rights reserved.

Automatic particle picking of biological molecules imaged by electron microscopy

One of the next great challenges of cell biology is the determination of the enormous number of protein structures encoded in genomes. In recent years, advances in electron cryo-microscopy and high-resolution single particle analysis have developed to the point where they now provide a methodology for high resolution structure determination. Using this approach, images of randomly oriented single particles are aligned computationally to reconstruct 3-D structures of proteins and even whole viruses. One of the limiting factors in obtaining high-resolution reconstructions is obtaining a large enough representative dataset ($>100,000$ particles). Traditionally particles have been manually picked which is an extremely labour intensive process. The problem is made especially difficult by the low signal-to-noise ratio of the images. This paper describes the development of automatic particle picking software, which has been tested with both negatively stained and cryo-electron micrographs. This algorithm has been shown to be capable of selecting most of the particles, with few false positives. Further work will involve extending the software to detect differently shaped and oriented particles.

Deposition of hydrocarbon molecules on diamond (001) surfaces : atomic scale modeling

The impact induced chemisorption of hydrocarbon molecules (CH3 and CH2) on H-terminated diamond (001)-(2x1) surface was investigated by molecular dynamics simulation using the many-body Brenner potential. The deposition dynamics of the CH3 radical at impact energies of 0.1-50 eV per molecule was studied and the energy threshold for chemisorption was calculated. The impact-induced decomposition of hydrogen atoms and the dimer opening mechanism on the surface was investigated. Furthermore, the probability for dimer opening event induced by chemisorption of CH, was simulated by randomly varying the impact position as well as the orientation of the molecule relative to the surface. Finally, the energetic hydrocarbons were modeled, slowing down one after the other to simulate the initial fabrication of diamond-like carbon (DLC) films. The structure characteristic in synthesized films with different hydrogen flux was studied. Our results indicate that CH3, CH2 and H are highly reactive and important species in diamond growth. Especially, the fraction of C-atoms in the film having sp(3) hybridization will be enhanced in the presence of H atoms, which is in good agreement with experimental observations. (C) 2002 Elsevier Science B.V. All rights reserved.

Insights into the mechanism of the reaction between ttrachloro-p-bnzoquinone and hydrogen peroxide and their implications in the catalytic role of water molecules in producing the hydroxyl radial

Detailed mechanisms for the formation of hydroxyl or alkoxyl radicals in the reactions between tetrachloro-p-benzoquinone (TCBQ) and organic hydroperoxides are crucial for better understanding the potential carcinogenicity of polyhalogenated quinones. Herein, the mechanism of the reaction between TCBQ and H2O2 has been systematically investigated at the B3LYP/6-311++G** level of theory in the presence of different numbers of water molecules. We report that the whole reaction can easily take place with the assistance of explicit water molecules. Namely, an initial intermediate is formed first. After that, a nucleophilic attack of H2O2 onto TCBQ occurs, which results in the formation of a second intermediate that contains an OOH group. Subsequently, this second intermediate decomposes homolytically through cleavage of the O-O bond to produce a hydroxyl radical. Energy analyses suggest that the nucleophilic attack is the rate-determining step in the whole reaction. The participation of explicit water molecules promotes the reaction significantly, which can be used to explain the experimental phenomena. In addition, the effects of F, Br, and CH3 substituents on this reaction have also been studied.

Evidence for the adsorption of molecules at special sites located at copper/zinc oxide interfaces. Part 2. -- A fourier-transform infrared spectroscopy study of methanol adsorption on reduced and oxidised Cu/ZnO/SiO2 catalysts

FTIR spectra are reported of methanol adsorbed at 295 K on ZnO/SiO 2, on reduced Cu/ZnO/SiO2 and on Cu/ZnO/SiO2 which had been preoxidised by exposure to nitrous oxide. Methanol on ZnO/SiO2 gave methoxy species on ZnO and SiO, in addition to both strongly and weakly physisorbed methanol on SiO2. The corresponding adsorption of methanol on reduced Cu/ZnO/SiO2 also gave methoxy species on Cu and a small amount of bridging formate. Reaction of methanol with a reoxidised Cu/ZnO/SiO2 catalyst resulted in an enhanced quantity of methoxy species on Cu. Heating adsorbed species on Cu/ZnO/SiO2 at 393 K led to the loss of methoxy groups on Cu and the concomitant formation of formate species on both ZnO and Cu. The comparable reaction on a reoxidised Cu/ZnO/SiO2 catalyst gave an increased amount of formate species on ZnO and this correlated with an increased quantity of methoxy groups lost from Cu. An explanation is given in terms of adsorption of formate and formaldehyde species at special sites located at the copper/zinc oxide interface.

Evidence for the adsorption of molecules at special sites located at copper/zinc oxide interfaces: Part 1. -- A Fourier-transform infrared study of formic acid and formaldehyde adsorption on reduced and oxidised Cu/ZnO/SiO2 catalysts

Fourier-transform infrared (FTIR) spectra are reported of formic acid and formaldehyde on ZnO/SiO2, reduced Cu/ZnO/SiO2 and reoxidised Cu/ZnO/SiO2 catalyst. Formic acid adsorption on ZnO/SiO2 produced mainly bidentate zinc formate species with a lesser quantity of unidentate zinc formate. Formic acid on reduced Cu/ZnO/SiO2 catalyst resulted not only in the formation of bridging copper formate structures but also in an enhanced amount of formate relative to that for ZnO/SiO2 catalyst. Formic acid on reoxidised Cu/ZnO/SiO2 gave unidentate formate species on copper in addition to zinc formate moieties. The interaction of formaldehyde with ZnO/SiO2 catalyst resulted in the formation of zinc formate species. The same reaction on reduced Cu/ZnO/SiO2 catalyst gave bridging formate on copper and a remarkable increase in the quantity of formate species associated with the zinc oxide. Adsorption of formaldehyde on a reoxidised Cu/ZnO/SiO2 catalyst produced bridging copper formate and again an apparent increase in the concentration of zinc formate species. An explanation in terms of the adsorption of molecules at special sites located at the interface between copper and zinc oxide is given.