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.

Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

Solvent effects on the electronic absorption spectrum of camphor using continuum, discrete or explicit approaches

We address the effect of solvation on the lowest electronic excitation energy of camphor. The solvents considered represent a large variation in-solvent polarity. We consider three conceptually different ways of accounting for the solvent using either an implicit, a discrete or an explicit solvation model. The solvatochromic shifts in polar solvents are found to be in good agreement with the experimental data for all three solvent models. However, both the implicit and discrete solvation models are less successful in predicting solvatochromic shifts for solvents of low polarity. The results presented suggest the importance of using explicit solvent molecules in the case of nonpolar solvents. (C) 2009 Elsevier B.V. All rights reserved.

Expression-system-dependent modulation of HIV-1 envelope glycoprotein antigenicity and immunogenicity.

Recombinant expression systems differ in the type of glycosylation they impart on expressed antigens such as the human immunodeficiency virus type 1 (HIV-1) envelope glycoproteins, potentially affecting their biological properties. We performed head-to-head antigenic, immunogenic and molecular profiling of two distantly related Env surface (gp120) antigens produced in different systems: (a) mammalian (293 FreeStyle cells; 293F) cells in the presence of kifunensine, which impart only high-mannose glycans; (b) insect cells (Spodoptera frugiperda, Sf9), which confer mainly paucimannosidic glycans; (c) Sf9 cells recombinant for mammalian glycosylation enzymes (Sf9 Mimic), which impart high-mannose, hybrid and complex glycans without sialic acid; and (d) 293F cells, which impart high-mannose, hybrid and complex glycans with sialic acid. Molecular models revealed a significant difference in gp120 glycan coverage between the Sf9-derived and wild-type mammalian-cell-derived material that is predicted to affect ligand binding sites proximal to glycans. Modeling of solvent-exposed surface electrostatic potentials showed that sialic acid imparts a significant negative surface charge that may influence gp120 antigenicity and immunogenicity. Gp120 expressed in systems that do not incorporate sialic acid displayed increased ligand binding to the CD4 binding and CD4-induced sites compared to those expressed in the system that do, and imparted other more subtle differences in antigenicity in a gp120 subtype-specific manner. Non-sialic-acid-containing gp120 was significantly more immunogenic than the sialylated version when administered in two different adjuvants, and induced higher titers of antibodies competing for CD4 binding site ligand-gp120 interaction. These findings suggest that non-sialic-acid-imparting systems yield gp120 immunogens with modified antigenic and immunogenic properties, considerations that should be considered when selecting expression systems for glycosylated antigens to be used for structure-function studies and for vaccine use.

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.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Exploring Cyclopentadienone Antiaromaticity: Charge Density Studies of Various Tetracyclones

A systematic study of six tetracyclones has been carried out using experimental and theoretical charge density analysis. A three pronged approach based on quantum theory of atoms in molecules (QTAIM), nucleus independent chemical shifts (NICS) criterion, and source function (SF) contributions has been performed to establish the degree of antiaromaticity of the central five-membered ring in all the derivatives. Electrostatic potentials mapped on the isodensity surface show that electron withdrawing substituents turn both C and O atoms of the carbonyl group more electropositive while retaining the direction of polarity.

Studies of Organic Molecule Recognition by Synthetic Cyclophane Receptors in Aqueous Media

The behaviors of six new cyclophane receptors for organic guest molecules in aqueous media are reported. These new hosts are modifications of more basic parent structures, and the main goal of their examination has been to determine how the modifications affect host selectivity for cationic guests. In particular, we have been interested in determining how additional non-covalent binding interactions can complement the cation-π interactions active in the parent systems. Three types of modifications were made to these systems. Firstly, neutral methoxy and bromine substituents were added to produce four of the six new macrocycles. Secondly, two additional aromatic rings (relative to the parent host) capable of making cation-π interactions with charged guest species were appended. Thirdly, a negatively charged carboxyl group was attached to produce a cavity in which electrostatic interactions should enhance cationic guest binding. ^1H-NMR and circular dichroic techniques were employed to determine the binding affinities of a wide variety of organic guests for the parent and modified structures in aqueous media.

Bromination of the parent host greatly enhances its binding in a general fashion, primarily as the result of hydrophobic interactions. The addition of methoxy groups does not enhance binding, apparently as a result of a collapse of the hosts into a conformation that is not suitable for binding. The appendage of extra aromatic rings enhances the binding of positively charged guests, most likely in response to more complete encapsulation of guest species. The addition of a negatively charged carboxylate enhances the binding to only selective groups of cationic guests. AM1 calculations of the electrostatic potentials of several guests molecules suggests that the enhancements seen with the modified receptor compared to the parent are most likely the result of close contact between regions of highest potential on the guest and the appended carboxylate.

OPTICAL PHONON MODES IN QUANTUM WIRES

Based on numerical results, the characteristics of each type of optical phonon mode in one-dimensional (1D) quantum wires are illustrated. The analytical formulae for 1D bulk-like optical displacements and associated electrostatic potentials are presented. Compared with the two-dimensional (2D) case, both the optical displacements and Frohlich potentials for the bulk-like modes vanish at the interfaces, but the finite dispersion of bulk phonons has a more pronounced effect on the 1D phonon modes.

OPTICAL-PHONON MODES AND FROHLICH POTENTIAL IN ONE-DIMENSIONAL QUANTUM-WELL WIRES

By extending our microscopic model on optical-phonon modes in quantum wells to one-dimensional (1D) quantum-well wires (QWW), the optical displacements and associated electrostatic potentials of optical-phonon modes in 1D QWW are calculated. The modes can be clearly divided into confined LO bulklike, TO bulklike modes, and extended interfacelike modes provided the bulk phonon dispersion is ignored. The character of each type of mode is illustrated with special attention to the interfacelike modes, which are hybrids of longitudinal- and transverse-optical waves from the corresponding bulk materials. Based on the numerical results, approximate analytical formulas for bulklike modes are presented. As in 2D wells, both the optical displacements and Frohlich potentials for the bulklike modes vanish at the interfaces. The finite dispersion of bulk phonons has a more pronounced effect on the 1D phonon modes because interfacelike modes show mixed characteristics of 2D interface and bulklike modes.

Chemical structural effects on γ-ray spectra of positron annihilation in fluorobenzenes

Spectra of ?-ray Doppler shifts for positron annihilation in benzene and its fluoro-derivatives are simulated using low energy plane wave positron (LEPWP) approximation. The results are compared with available measurements. It is found that the Doppler shifts in these larger aromatic compounds are dominated by the contributions of the valence electrons and that the LEPWP model overestimates the measurements by approximately 30%, in agreement with previous findings in noble gases and small molecules. It is further revealed that the halogen atoms not only switch the sign of the charges on carbon atoms that they bond to, but that they also polarize other C-H bonds in the molecule leading to a redistribution of the molecular electrostatic potentials. As a result, it is likely that the halogen atoms contribute more significantly to the annihilation process. The present study also suggests that, while the Doppler shifts are sensitive to the number of valence electrons in the molecules, they are less sensitive to the chemical structures of isomers that have the same numbers and type of atoms and, hence, the same numbers of electrons. Further investigation of this effect is warranted. © EDP Sciences, Società Italiana di Fisica, Springer-Verlag 2012.

Quantitative structure-hydrophobicity relationships of molecular fragments and beyond

Quantitative structure-property relationship (QSPR) models were firstly established for the hydrophobic substituent constant (πX) using the theoretical descriptors derived solely from electrostatic potentials (EPSs) at the substituent atoms. The descriptors introduced are found to be related to hydrogen-bond basicity, hydrogen-bond acidity, cavity, or dipolarity/polarizability terms in linear solvation energy relationship, which endows the models good interpretability. The predictive capabilities of the models constructed were also verified by rigorous Monte Carlo cross-validation. Then, eight groups of meta- or para- disubstituted benzenes and one group of substituted pyridines were investigated. QSPR models for individual systems were achieved with the ESP-derived descriptors. Additionally, two QSPR models were also established for Rekker's fragment constants (foct), which is a secondary-treatment quantity and reflects average contribution of the fragment to logP. It has been demonstrated that the descriptors derived from ESPs at the fragments, can be well used to quantitatively express the relationship between fragment structures and their hydrophobic properties, regardless of the attached parent structure or the valence state. Finally, the relations of Hammett σ constant and ESP quantities were explored. It implies that σ and π, which are essential in classic QSAR and represent different type of contributions to biological activities, are also complementary in interaction site.

Cyclocreatine, an anticancer and neuroprotective agent. Spectroscopic, structural and theoretical study

The structural, spectroscopic and theoretical study of cyclocreatine (1-carboxymethyl-2-iminoimidazolidine, CyCre) has been performed prompted by the biological relevance of the molecule and its potential role as a ligand in biometalic compounds. The crystal structure of CyCre has been determined by X-ray diffraction methods. The compound crystallizes as a zwitterion in the monoclinic system, space group P2(1)/c. The crystal is further stabilized by a network of N-H center dot center dot center dot O bonds. Infrared and Raman spectra of the solid, electronic spectra of aqueous solutions at different pH values and (1)H and (13)C NMR spectra have been recorded and analyzed. Band assignments were accomplished with the help of theoretical calculations. Optimized molecular geometries, harmonic vibrational frequencies and molecular electrostatic potentials were calculated using methods based on the density functional theory. (C) 2010 Elsevier B.V. All rights reserved.