The effect of the protein environment on the structure and charge distribution of the retinal Schiff base in bacteriorhodopsin

The effects of the amino acid side chains of the binding pocket of bacteriorhodopsin (bR) and of a water molecule on the structure of the retinal Schiff base have been studied using Becke3LYP/6-31G* level of density functional theory. A model protonated Schiff base structure including six conjugated double bonds and methyl substituents was optimized in the presence of several amino acid side chains and of a water molecule, separately. The Schiff base structure was also calculated in the form of a neutral species. At each optimized complex geometry the atomic charges of the model Schiff base were calculated using Mulliken population analysis. In agreement with previously proposed counterion(s) of the protonated retinal Schiff base in bR, the results show that Asp₈₅ and Asp₂₁₂, which are present in the form of negatively charged groups, have significantly large effects on the structure and electronic configuration of both unprotonated and protonated model Schiff bases. The presence of a water molecule in the vicinity of the Schiff base demonstrates significant effects which are comparable to those of aspartate groups. Other side chains studied did not show any significant effect in this direction. Apart from the aspartate groups and the water molecule, in none of the other complexes studied are the atomic charges and the bond alternation of the model Schiff base significantly influenced by the presence of the neighboring amino acids. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

Accuracy of free energies of hydration for organic molecules from 6-31g*-derived partial charges

Absolute free energies of hydration have been computed for 13 diverse organic molecules using partial charges derived from ab initio 6-31G* wave functions. Both Mulliken charges and charges fit to the electrostatic potential surface (EPS) were considered in conjunction with OPLS Lennard–Jones parameters for the organic molecules and the TIP4P model of water. Monte Carlo simulations with statistical perturbation theory yielded relative free energies of hydration. These were converted to absolute quantities through perturbations to reference molecules for which absolute free energies of hydration had been obtained previously in TIP4P water. The average errors in the computed absolute free energies of hydration are 1.1 kcal/mol for the 6-31G* EPS charges and 4.0 kcal/mol for the Mulliken charges. For the EPS charges, the largest individual errors are under 2 kcal/mol except for acetamide, in which case the error is 3.7 kcal/mol. The hydrogen bonding between the organic solutes and water has also been characterized. © John Wiley & Sons, Inc.     

Various atomic charge calculation schemes of CoMFA on HIF‐1 inhibitors of moracin analogs

Selection of appropriate partial charges in a molecule is crucial to derive good quantitative structure–activity relationship models. In this work, several partial atomic charges were assigned and tested in a comparative molecular field analysis (CoMFA) models. Many CoMFA models were generated for a series of hypoxia inducible factor 1 (HIF‐1) inhibitors using various partial atomic charges including charge equalization, Mülliken population analysis (MPA), natural population analysis, and electrostatic potential (ESP)‐derived charges. These atomic charges were investigated at various theoretical levels such as empirical, semiempirical, Hartree–Fock (HF), and density functional theory (DFT). Among them, Merz‐Singh‐Kollman (MK) ESP‐derived charges at the level of HF/6‐31G* gave the highest predictive q² with experimental pIC₅₀ values. With this charge scheme, a detailed analysis of CoMFA model was performed to understand the electrostatic interactions between ligand and receptor. More elaborate charge calculation schemes such as HF and DFT correlated more strongly with activity than empirical or semiempirical schemes. The choice of optimization methods was important. As geometries were fully optimized at the given levels of theory, the aligned structures were different. They differed considerably, especially for the flexible parts. This was likely the source of the substantial variation of q² values, even when the same steric factor was considered without electrostatic parameters. ESP‐derived charges were most appropriate to describe CoMFA electrostatic interactions among MPA, NBA, and ESP charges. Overall q² values vary considerably (0.8–0.5) depending on the charge schemes applied. The results demonstrate the need to consider more appropriate atomic charges rather than default CoMFA charges. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Cluster calculations of the H2O/Pt(111) system

The electronic interaction between water and a Pt(111) surface as evaluated for different Pt_x(H₂O)_y clusters is discussed. Hartree–Fock–Slater (HFS ) one-electron ground state energies, ionization potentials, partial densities of states, and Mulliken occupation numbers are related to bonding shifts, as well as initial and final state screening for different orientations of the molecule. The formation of Pt? H₂O bonds are sensitive to the orientation since surface oriented H atoms bridge the spatial separation between O 2p and Pt 5d orbitals and thus increase the intermixing of metal and adsorbate orbitals. The dipole moment and the net charge of the H₂O molecule is also discussed. Finally, approximations of the metal–H₂O potential for use in statistical models of the liquid–metal interface are suggested.     

A POSSIBLE MECHANISM FOR THE ROLE OF N-RETINYLIDENE PHOSPHATIDYL ETHANOLAMINE IN THE REGENERATION OF RHODOPSIN

Abstract— A series of molecular orbital calculations on a model Schiff base comparable to protonated N -retinylidene phosphatidyl ethanolamine isomers has been made. The effect of the charged oxygen atoms of the phosphate moiety on the distribution of positive charge along the polyene chain of these isomers has been calculated. The stabilizing coulombic energy of interaction of these opposite charges and the possibility of free rotation around carbonxarbon double bonds in the electronically excited state has led to the conclusion that an 11- cis Schiff base isomer is the most probable product of the photoisomerization of an all-trans Schiff base.
The formation of a stable unprotonated all-trans Schiff base in aqueous detergent dispersion and its subsequent conversion to the protonated form, both with absorption spectra in conformity with the literature, is demonstrated.     

Interaction of chromophore, 11-<Emphasis Type="Italic">cis</Emphasis>-retinal,with amino acid residues of the visual pigment rhodopsin in the region of protonated Schiff base: A molecular dynamics study

A molecular dynamics study of the dark adapted visual pigment rhodopsin molecule was carried out. The interaction of the chromophore group, 11-cis-retinal, with the nearest amino acid residues in the chromophore center of the molecule, namely, in the region of the protonated Schiff base linkage, was analyzed. Most likely, the interaction of the CH=NH bond with the negatively charged amino acid residue Glu113 cannot be described as a simple electrostatic interaction of two oppositely charged groups. One can propose that not only Glu113 but also Glu181 and Ser186 are involved in stabilization of the protonated Schiff base linkage. Accord-ing to calculations, Glu181 interacts, as the counter-ion, with the Schiff base indirectly via Ser186. The intramolecular mechanisms of protonated Schiff base stabilization in rhodopsin are discussed. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 19–27, January, 2007.     

Semiempirical treatment of electrostatic potentials and partial charges in combined quantum mechanical and molecular mechanical approaches

A semiempirical treatment of electrostatic potentials and partial charges is presented. These are the basic components needed for the evaluation of electrostatic interaction energies in combined quantum mechanical and molecular mechanical approaches. The procedure to compute electrostatic potentials uses AM1 and MNDO wave functions and is based on one previously suggested by Ford and Wang. It retains the NDDO approximation and is thus both easy to implement and computationally efficient. Partial atomic charges are derived from a semiempirical charge equilibration model, which is based on the principle of electronegativity equalization. Large sets of ab initio restricted Hartee-Fock (RHF/6-31G*) reference data have been used to calibrate the semiempirical models. Applying the final parameters (C, H, N, O), the ab initio electrostatic potentials are reproduced with an average accuracy of 20% (AM1) and 25% (MNDO), respectively, and the ab initio potential derived charges normally to within 0.1 e. In most cases our parameterized models are more accurate than the much more expensive quasi ab initio techniques, which employ deorthogonalized semiempirical wave functions and have generally been preferred in previous applications. © 1996 John Wiley & Sons, Inc.     

Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations

Ab initio MP2/6-31G* interaction energies were calculated for more than 80 geometries of stacked cytosine dimer. Diffuse polarization functions were used to properly cover the dispersion energy. The results of ab initio calculations were compared with those obtained from three electrostatic empirical potential models, constructed as the sum of a Lennard-Jones potential (covering dispersion and repulsion contributions) and the electrostatic term. Point charges and point multipoles of the electrostatic term were also obtained at the MP2/6-31G* level of theory. The point charge MEP model (atomic charges derived from molecular electrostatic potential) satisfactorily reproduced the ab initio data. Addition of π-charges localized below and above the cytosine plane did not affect the calculated energies. The model employing the distributed multipole analysis gave worse agreement with the ab initio data than the MEP approach. The MP2 MEP charges were also derived using larger sets of atomic orbitals: cc-pVDZ, 6-311 + G(2d, p), and aug-cc-pVDZ. Differences between interaction energies calculated using these three sets of point charges and the MP2/6-31G* charges were smaller than 0.8 kcal/mol. The correlated ab initio calculations were also compared with the density functional theory (DFT) method. DFT calculations well reproduced the electrostatic part of interaction energy. They also covered some nonelectrostatic short-range effects which were not reproduced by the empirical potentials. The DFT method does not include the dispersion energy. This energy, approximated by an empirical term, was therefore added to the DFT interaction energy. The resulting interaction energy exhibited an artifact secondary minimum for a 3.9-4.0 vertical separation of bases. This defect is inherent in the DFT functionals, because it is not observed for the Hartree-Fock + dispersion interaction energy.© 1996 John Wiley & Sons, Inc.     

Quantum-chemical prediction of the scandium extraction power of phosphinic acids

The extraction power of phosphinic acids (PAs) toward scandium was predicted on the basis of PA characteristics (dissociation enthalpies and energies, L?wdin and Mulliken effective charges on the oxygen and phosphorus atoms of the phosphoryl group) and the energy of scandium complexation with PAs during extraction, which were calculated by quantum-chemical methods with the use of the HyperChem Release 8.0.4 (Hypercube Inc.) and PC GAMESS version 6.4 software. It was demonstrated that the use of the minimal basis set STO-3G provided the highest multiple correlation coefficient of log K _Sc (K _Sc is the concentration constant of scandium extraction) with the L?wdin effective charge of the oxygen atom and the dissociation energy of PAs (R = 0.963).     

Accurate partial atomic charges for high-energy molecules using class IV charge models with the MIDI! basis set

We have recently developed a new class IV charge model for calculating partial atomic charges in molecules. The new model, called charge model 3 (CM3), was parameterized for calculations on molecules containing H, Li, C, N, O, F, Si, S, P, Cl, and Br by Hartree–Fock theory and by hybrid density functional theory (HDFT) based on the modified Perdew–Wang density functional with several basis sets. In the present article, we extend CM3 for calculating partial atomic charges by Hartree–Fock theory with the economical but well balanced MIDI! basis set. Then, using a test set of accurate dipole moments for molecules containing nitramine functional groups (which include many high-energy materials), we demonstrate the utility of several parameters designed to improve the charges in molecules containing both N and O atoms. We also show that one of our most recently developed CM3 models that is designed for use with wave functions calculated at the mPWXPW91/MIDI! level of theory (where X denotes a variable percentage of Hartree–Fock exchange) gives accurate charge distributions in nitramines without additional parameters for N and O. To demonstrate the reliability of partial atomic charges calculated with CM3, we use these atomic charges to calculate polarization free energies for several nitramines, including the commonly used explosives 1,3,5-trinitro-s-triazine (RDX) and 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (HNIW), in nitromethane. These polarization energies are large and negative, indicating that electrostatic interactions between the charge distribution of the molecule and the solvent make a large contribution to the free energy of solvation of nitramines. By extension, the same conclusion should apply to solid-state condensation. Also, in contrast to some other charge models, CM3 yields atomic charges that are relatively insensitive to the presence of buried atoms and small conformational changes in the molecule, as well as to the level of treatment of electron correlation. This type of charge model should be useful in the future development of solvation models and force fields designed to estimate intramolecular interactions of nitramines in the condensed phase.     

Absorption studies of neutral retinal Schiff base chromophores

The neutral retinal Schiff base is connected to opsin in UV sensing pigments and in the blue-shifted meta-II signaling state of the rhodopsin photocycle. We have designed and synthesized two model systems for this neutral chromophore and have measured their gas-phase absorption spectra in the electrostatic storage ring ELISA with a photofragmentation technique. By comparison to the absorption spectrum of the protonated retinal Schiff base in vacuo, we found that the blue shift caused by deprotonation of the Schiff base is more than 200 nm. The absorption properties of the UV absorbing proteins are thus largely determined by the intrinsic properties of the chromophore. The effect of approaching a positive charge to the Schiff base was also studied, as well as the susceptibility of the protonated and unprotonated chromophores to experience spectral shifts in different solvents.     

Theoretical study of 3d‐metal mononitrides using DFT method

3d‐Metal mononitrides are studied using the density functional theory method. The lowest spin state for these dimers is obtained using the B3LYP hybrid functional with the 6‐311+G* basis set. The equilibrium geometries, vibrational frequencies, binding energies, Mulliken, and natural orbital population analysis charges, natural orbital electronic configuration, electron affinity, and ionization potential are obtained. Mulliken as well as natural orbital population analysis charges indicate that for all dimers, in cations most of the positive charge localized on the transition metal atom where in anions most of the negative charge localized on nitrogen atom. The binding energies for 3d‐metal mononitrides are higher than those for monocarbides and monoxides. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Microscopic Protonation/Deprotonation Equilibria of the Anti-Inflammatory Agent Piroxicam

The microscopic ionization behavior of piroxicam was investigated using two different approaches, i.e., direct UV spectroscopy and an indirect analogue approach (deductive method). The best microscopic pK_a values (pK_a12 = 4.60, pK_a21 = 5.40, pK_a22 = 2.72, and pK_a11 = 1.92) were obtained by the deductive method using as pK_a22 the pK_a of the enolic O-methylated piroxicam 2 . The results show remarkable electrostatic effects in the protonation/deprotonation equilibria, a marked increase in the acidity of the enolic function (2.68 pK_a units) being caused by the pyridinium group. The electronic structure of piroxicam was studied based on ¹H-NMR chemical shifts at various ionization states, indicating an extended electron conjugation through the molecule. The partition measurements in octan-1-ol/H₂O of zwitterionic compound 3 (the pyridyl N-methyl derivative of piroxicam ( 1 )) suggest that the two opposite charges in zwitterionic piroxicam are indeed in a close intramolecular proximity.     

Computation of the influence of chemical substitution on the pK a of pyridine using semiempirical and ab initio methods

The physical properties of chemicals are strongly influenced by their protonation state, affecting, for example, solubility or hydrogen-bonding characteristics. The ability to accurately calculate protonation states in the form of pK _as is, therefore, desirable. Calculations of pK _a changes in a series of substituted pyridines are presented. Computations were performed using both ab initio and semiempirical approaches, including free energies of solvation via reaction-field models. The selected methods are readily accessible with respect to both software and computational feasibility. Comparison of calculated and experimental pK _as shows the experimental trends to be reasonably reproduced by the computations with root-mean-square differences ranging from 1.22 to 4.14 pK _a units. Of the theoretical methods applied the best agreement occurred using the second-order M?ller–Plesset/6-31G(d)/isodensity surface polarized continuum solvation model, while the more computationally accessible Austin model 1/Solvent model 2 (SM2) approach yielded results similar to the ab initio methods. Analysis of component contributions to the calculated pK _as indicates the largest source of error to be associated with the free energies of solvation of the protonated species followed by the gas-phase protonation energies; while the latter may be improved via the use of higher levels of theory, enhancements in the former require improvements in the solvation models. The inclusion of alternate minimum in the computation of pK _as is also indicated to contribute to differences between experimental and calculated pK _a values. Received: 27 April 1999 / Accepted: 27 July 1999 / Published online: 2 November 1999     

A comparative study of effective core potential and full-electron calculations in Mo compounds. I. An analysis of topological properties of bond charge distribution

Effective core potential (ECP) and full-electron (FE) calculations for MoS₄^?2, MoO₄^?2, and MoOCl₄ compounds were analyzed. Geometry parameters, binding energies, charge distributions, and topological properties of the electronic density were studied for Mo? L bonds (L = S, O, Cl). Results clearly indicate that those approaches that include valence plus 4s and 4p electrons (ECP2 methods) are able to reproduce the topological properties of Mo? L bonds, charge distributions, and geometries with respect to those obtained by FE methods. ECP methods that consider only the 4d and 5s valence electrons (ECP1) fail in the calculation of molecular properties. The use of 5p functions in ECP1 approaches produces a negative Mulliken charge on Mo. Bader's charges give more consistent results than Mulliken's ones. A new parameter for measuring the degree of ionicity is proposed. © 1994 by John Wiley & Sons, Inc.     

Importance of electrostatic polarizability in calculating cysteine acidity constants and copper(I) binding energy of Bacillus subtilis CopZ

CopZ is a copper chaperone from Bacillus subtilis. It is an important part of Cu(I) trafficking. We have calculated pK_a values for the CXXC motif of this protein, which is responsible for the Cu(I) binding, and the Cu(I) binding constants. Polarizable and fixed‐charges formalisms were used, and solvation parameters for the both models have been refitted. We had to partially redevelop parameters for the protonated and deprotonated cysteine residues. We have discovered that the polarizable force field (PFF) is qualitatively superior and allows a uniformly better level of energetic results. The PFF pK_a values for cysteine are within about 0.8–2.8 pH units of the experimental data, while the fixed‐charges OPLS formalism yields errors of up to tens of units. The PFF magnitude of the copper binding energy is about 10 kcal/mol or 50% higher than the experimental value, while the using the refitted OPLS parameters leads to an overall positive binding energy, thus predicting no thermodynamically stable complex. At the same time, the agreement of the polarizable S···Cu(I) distances with the experimental results is within 0.08 Å range, and the nonpolarizable calculations lead to an error of about 0.4 Å. Moreover, the accuracy of the PFF has been achieved without any explicit fitting to either pK_a or CopZ···Cu(I) binding energies. We believe that this makes our polarizable technique a choice method in reproducing protein—copper binding and further supports the notion that explicit treatment of electrostatic polarization is crucial in many biologically relevant studies, especially ion binding and transport. © 2012 Wiley Periodicals, Inc.     

Spectroscopic study of rigid-rod polymers. II. Protonation effect

The protonation of a heterocyclic rigid-rod polymer poly(p-phenylene benzbisthiazole) and its model compound has been studied by UV-visible and Raman spectroscopy. Because of the two nitrogens on the heterocyclic ring, spectroscopic features of unprotonated, half-protonated, and fully protonated structures have been identified. For the fully protonated molecule, there is also an increase in conjugation between the phenyl ring and the heterocyclic ring.     

On the representation of electrostatic fields around ab initio charge distributions

Summary We compare two methods (Mulliken charges and a distributed multipole analysis, DMA) of representing an ab initio charge distribution for calculating the electrostatic field and potential outside the molecule, using pyrimidine and the RNA base uracil as examples. This is done using a 3-D graphical display of the electrostatic fields, which, when used with real-time rotation, zooming and clipping, has many advantages for qualitatively assessing the electrostatic interactions of a molecule. The errors involved in using Mulliken point charges may be of similar magnitude to the total electrostatic field in regions which are important in recognition processes. The DMA representation automatically includes the anisotropic electrostatic effects of non-spherical features in the charge distribution of each atom, and yet the displayed electrostatic fields around the atoms which have lone-pair density do not show marked anisotropy.