Resonance description of S-nitrosothiols: insights into reactivity

A resonance representation of the electronic structure of S-nitrosothiols as a combination of the conventional R-S-N=O structure, a zwitterionic structure R-S+=N-O-, and a RS-/NO+ ion pair is proposed. The resonance forms are employed to predict and rationalize the structural and conformational properties of RSNOs, their interaction with Lewis acids, and their overall chemical reactivity.     

Ionization Behavior and Ionization-Dependent Conformation of Raclopride,a Dopamine D2 Receptor Antagonist

Raclopride, an antipsychotic 6-methoxysalicylamide (= 2-hydroxy-6-methoxybenzamide) derivative, was shown by titrimetry and UV-pholometry to exist in zwitterionic form at physiological pH. Calculations revealed that the neutral and zwitterionic forms differ considerably in their conformational behavior, the latter form being energetically favored by an intramolecular phenolate–ammonium ionic bond. These findings indicate that raciopride and other halogenated 6-melhoxysalicylamides with a highly acidic phenolic group may not resemble other ortho-methoxybenzamides in their stereoelectronic structure and mode of binding to the dopamine D₂ receptor.     

Molecular dynamics simulation study on zwitterionic structure to maintain the natural behavior of polyalanine13 in aqueous environment

Molecular dynamics simulations are applied to the initial stage of polyalanine13 conformational transi- tion from α-helix to random coil in aqueous environment and the interaction of polyalanine13 with zwitterionic and hydrophobic surfaces respectively in the same condition. The analysis of secondary structure, hydrogen bonds, RMSD, dihedral distribution, and the degree of adsorption are performed. The results show that zwitterionic structure maintains the natural behavior of polyalanine13 in water to a better extent, which should be an indirect proof of the hypothesis of "maintain of normal structure."     

A semi-empirical and ab initio combined approach for the full conformational searches of gaseous lysine and lysine–H2O complex

A full structural search of the canonical, zwitterionic, protonated and deprotonated lysine conformers in gas phase is presented. A total of 17,496 canonical, 972 zwitterionic, 11,664 protonated and 1458 trial deprotonated structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were initially optimized at the AM1 level, and the resulting structures were determined at the B3LYP/6-311G^* level. A total of 927 canonical, 730 protonated and 193 deprotonated conformers were found, but there were no stable zwitterionic structures in the gas phase. The most stable conformers of the canonical, protonated and deprotonated lysine were further optimized at the B3LYP/6-311++G^** level. The energies of the most stable structures were determined at the MP2/6-311G(2df,p) level and the vibrational frequencies were calculated at the B3LYP/6-311++G^** level. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, vertical ionization energies, enthalpies, Gibbs free energies and conformational distributions of gaseous lysine were presented. Numerous new structures are found and the lowest-energy lysine conformer is more stable than the existing one by 1.1 kcal/mol. Hydrogen bonds are classified and may cause significant red-shifts to the associated vibrational frequencies. The calculated proton affinity/dissociation energy and gas-phase basicity/acidity are in good agreement with the experiments. Calculations are also presented for the canonical lysine–H₂O and zwitterionic lysine–H₂O clusters. Interaction between lysine and H₂O significantly affects the relative conformational stabilities. Only one water molecule is sufficient to produce the stable zwitterionic structures in gas phase. The lowest-energy structure is found to be zwitterions when applying the conductor-like polarized continuum solvent model (CPCM) to the lysine–H₂O complexes.     

Stepwise solvation of an amino acid: the appearance of zwitterionic structures

How many solvent molecules are required to solvate an amino acid? This apparently simple question, which relates to the number of solvent molecules necessary to change the amino acid from its gas-phase neutral structure to the zwitterionic solvated structure, remains unanswered to date. Here we present experimental and theoretical (density functional theory: B3LYP/6-31+G**) infrared spectra for tryptophan-watern complexes where n = 1-6, which suggest that the zwitterionic structure becomes competitive in energy at the high end of the series. Compelling evidence for a gradual transition to zwitterionic structures comes from tryptophan-methanol complexes up to n = 9. Starting from n = 5, the infrared spectra show increasing intensity in the diagnostic asymmetric COO- stretch and in the weaker NH3+ bending modes as the cluster size increases. Moreover, convergence toward the Fourier transform infrared spectrum of a solution of tryptophan in methanol is clearly observed. For small solvent complexes (n = 1-4), the microsolvation by methanol and water is shown to behave very similarly. A detailed comparison of the experimental and the theoretical spectra allows us to determine both the preferred solvent binding sites on the amino acid and the evolution of conformational structures of tryptophan as the number of attached solvent molecules increases.     

Monte Carlo simulation of water solvent with biomolecules: Serine and the corresponding zwitterion

Monte Carlo simulation results are reported for clusters consisting of 250 water molecules surrounding serine, both in the neutral form and in two zwitterionic conformers (in order to gain some insight into conformational effects). Calculations were carried out at 300 K and using two-body potentials obtained by means of quantum-mechanical calculations. The spatial dependence of the average interaction energies was investigated. The solvation structure was investigated by means of radial distribution functions and probability density maps, which showed a few water molecules directly solvating the hydrophilic groups and, beyond them, a more or less rich and complex hydrogen-bonded network of solvent molecules.     

Ab initio studies of aspartic acid conformers in gas phase and in solution

Systematic and extensive conformational searches of aspartic acid in gas phase and in solution have been performed. For the gaseous aspartic acid, a total of 1296 trial canonical structures and 216 trial zwitterionic structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 139 canonical conformers were found, but no stable zwitterionic structure was found. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, and vertical ionization energies of the canonical conformers were determined. Single-point energies were also calculated at the MP2/6-311++G** and CCSD/6-311++G** levels. The equilibrium distributions of the gaseous conformers at various temperatures were calculated. The proton affinity and gas phase basicity were calculated and the results are in excellent agreement with the experiments. The conformations in the solution were studied with different solvation models. The 216 trial zwitterionic structures were first optimized at the B3LYP/6-311G* level using the Onsager self-consistent reaction field model (SCRF) and then optimized at the B3LYP/6-311++G** level using the conductorlike polarized continuum model (CPCM) SCRF theory. A total of 22 zwitterions conformers were found. The gaseous canonical conformers were combined with the CPCM model and optimized at the B3LYP/6-311++G** level. The solvated zwitterionic and canonical structures were further examined by the discrete/SCRF model with one and two water molecules. The incremental solvation of the canonical and zwitterionic structures with up to six water molecules in gas phase was systematically examined. The studies show that combining aspartic acid with at least six water molecules in the gas phase or two water molecules and a SCRF solution model is required to provide qualitatively correct results in the solution.     

A robust force field based method for calculating conformational energies of charged drug-like molecules

The binding affinity of a drug-like molecule depends among other things on the availability of the bioactive conformation. If the bioactive conformation has a significantly higher energy than the global minimum energy conformation, then the molecule is unlikely to bind to its target. Determination of the global minimum energy conformation and calculation of conformational penalties of binding is a prerequisite for prediction of reliable binding affinities. Here, we present a simple and computationally efficient procedure to estimate the global energy minimum for a wide variety of structurally diverse molecules, including polar and charged compounds. Identifying global energy minimum conformations of such compounds with force field methods is problematic due to the exaggeration of intramolecular electrostatic interactions. We demonstrate that the global energy minimum conformations of zwitterionic compounds generated by conformational analysis with modified electrostatics are good approximations of the conformational distributions predicted by experimental data and with molecular dynamics performed in explicit solvent. Finally the method is used to calculate conformational penalties for zwitterionic GluA2 agonists and to filter false positives from a docking study.     

Stabilization of zwitterionic proline by DMSO

Zwitterionic stabilization and metal‐free organocatalysis are two emerging topics. In this work, the numbers of DMSO molecules required to render zwitterionic proline geometrically stable, energetically preferential, and conformationally predominant have been determined, as one, three, and three, respectively. Conformations are analyzed for proline conformers interacted with one, two, and three DMSO molecules, and three DMSO molecules are enough to fill up the first shell of proline. Relative stabilities of two selected canonical structures are dependent on the DMSO contents, while zwitterionic stabilities improve monotonously with increase of DMSO contents. DMSO causes a conformational diversity and good zwitterionic stabilization effects, which result from the synergetic effects of two types of H‐bonding interactions. With increase of DMSO contents, type‐2 H‐bonding (CH as donors) contributes more to zwitterionic stabilization. At any DMSO content, zwitterionic proline is facile to form because of low activation energies, and this study helps to understand proline‐catalyzed processes. © 2015 Wiley Periodicals, Inc.     

聚氨酯表面构建磺铵两性离子结构及其抗血小板粘附性的研究

通过己二异氰酸酯(HDI)在聚醚氨酯(PU)表面构建磺铵两性离子结构,以改善其不凝血性能,首先用HDI活化PU表面,生成PU-NCO衍生物;然后通过N,N-二甲基乙醇胺(DMEEA)中的-OH和PU表面的-NCO反应生成PU-N(CH₃)₂;最后用丙磺酸内酯(PS)进行开环.生成磺铵两性离子结构,用ATR-IR表征了各步反应,对构建前后材料的抗血小板粘附性能进行了比较,结果表明,磺铵两性离子结构具有优异的抗血小板粘附性.     

Hydrophilic interaction liquid chromatography coupled with MS/MS to detect and quantify dicarboxyethyl glutathione,a metabolic biomarker of the fumarate hydratase deficient cancer cell

The identification of a glutathione (GSH) fumarate conjugate, dicarboxyethyl glutathione, formed during the nonenzymatic succination of GSH by fumarate was confirmed in fumarate hydratase deficient cells using a product ion scan approach followed by hydrophilic interaction liquid chromatography coupled with MS/MS. GSH and its conjugates, including dicarboxyethyl glutathione and glutathione disulfide, were successfully separated on a zwitterionic stationary phase and detected by MS/MS operated under negative ESI mode. The relative quantitation of the analytes in cell extracts was carried out and a correction model was established to determine correction factors under matrix effects and the response mismatch between the analytes. These factors were calculated and iteratively used to measure all analytes in cell extracts, based on calibration curves constructed in neat solution. The model was a closed‐loop calculation, consisting of two sides with each side of the loop presenting a calculation pathway. Deviation of the correction factors obtained from these pathways manifested the model accuracy. The model was evaluated and there was no significant difference between the two pathways.     

Effects of microsolvation and aqueous solvation on the tautomers of histidine: a computational study on energy, structure and IR spectrum

Microsolvation and combined microsolvation-continuum approaches are employed to investigate the structures and energies of canonical and zwitterionic histidine conformers. The effect of hydration on the relative conformational stability is examined. The strategy of exploring singly and doubly hydrated structures and the possible microsolvation patterns are described. We find that bonding water molecule may significantly change the relative conformational stabilities. In gas phase, the singly and doubly hydrated canonical forms are more stable than their zwitterionic counterparts. In solution, the continuum solvent model shows that bare zwitterionic form is more stable than bare canonical form by 1.1 kcal/mol. This energy separation is increased to 2.2 and 3.4 kcal/mol with inclusion of one and two explicit water molecules, respectively. We have also observed that the doubly hydrated structures obtained by combining two water molecules simultaneously to the solute molecule are preferred over the stepwise hydration. Hydrogen bond energies for the most stable hydrated histidine tautomers are determined by the atoms in molecules theory. The infrared (IR) spectra for the most stable singly and doubly hydrated structures of both histidine tautomers in gas phase are characterized. The stretching frequencies for NH of imidazole ring and OH of COOH are red shifted due to the hydrations. The IR spectra for the most stable zwitterionic tautomers in solution are also presented and discussed in connection with the comparison to the experiments. The pKa values obtained for the ring protonated zwitterions with two explicit water molecules appear to be in good agreement with the experiments.     

Photoredox Reaction of Naphthoquinone C-Glycoside Revisited: Insight into Stereochemical Aspect

Insights into the naphthoquinone photoredox reactions have been gained from the reactions of a C-glycosyl naphthoquinone. The retentiveness comes from the balance of the lifetime of intermediary zwitterionic species relative to the conformational change to bring about the stereo-mutating C−C bond rotation.     

Gaseous arginine conformers and their unique intramolecular interactions

Extensive ab initio calculations were employed to characterize stable conformers of gaseous arginine, both the canonical and zwitterionic tautomers. Step-by-step geometry optimizations of possible single-bond rotamers at the B3LYP/6-31G(d), B3LYP/6-31++G(d,p), and MP2/6-31++G(d,p) levels yield numerous structures that are more stable than any known ones. The final electronic energies of the conformers were determined at the CCSD/6-31++G(d,p) level. The lowest energies of the canonical and zwitterionic structures are lower than the existing values by 2.0 and 2.3 kcal/mol, respectively. The relative energies, rotational constants, dipole moments, and harmonic frequencies of the stable conformers remain for future experimental verification. The conformational distributions at various temperatures, estimated according to thermodynamic principles, consist almost exclusively of the newly found structures. One striking feature is the occurrence of blue-shifting hydrogen bonds in all six of the most stable conformers. A unique feature of important conformations is the coexistence of dihydrogen and blue- and red-shifting hydrogen bonds. In addition to the hydrogen bonds, the stereoelectronic effects were also found to be important stabilization factors. The calculated and measured proton affinities agree within the theoretical and experimental uncertainties, affirming the high quality of our conformational search. The theoretical gas-phase basicity of 245.9 kcal/mol is also in good agreement with the experimental value of 240.6 kcal/mol. The extensive searches establish firmly that gaseous arginine exists primarily in the canonical and not the zwitterionic form.     

Intrinsic Neutral and Anionic Structures of Glutathione

Gas‐phase intrinsic structures of intact neutral and anionic glutathione (GSH) have been determined by means of a combination of negative ion photo‐electron spectroscopy and quantum chemistry calculations. The inferred structures of the neutral parents of those peptide anions are canonical (non‐zwitterionic). These intrinsic structures are compared to those already known in aqueous solution or determined by crystallography in binding sites of enzymes.     

Sheet-like assemblies of charged amphiphilic α/β-peptides at the air-water interface

There is growing interest in the design of molecules that undergo predictable self-assembly. Bioinspired oligomers with well-defined conformational propensities are attractive from this perspective, since they can be constructed from diverse building blocks, and self-assembly can be directed by the identities and sequence of the subunits. Here we describe the structure of monolayers formed at the air-water interface by amphiphilic α/β-peptides with 1:1 alternation of α- and β-amino acid residues along the backbone. Two of the α/β-peptides, one a dianion and the other a dication, were used to determine differences between self-assemblies of the net negatively and positively charged oligomers. Two additional α/β-peptides, both zwitterionic, were designed to favor assembly in a 1:1 molar ratio mixture with parallel orientation of neighboring strands. Monolayers formed by these α/β-peptides at the air-water interface were characterized by surface pressure-area isotherms, grazing incidence X-ray diffraction (GIXD), atomic force microscopy and ATR-FTIR. GIXD data indicate that the α/β-peptide assemblies exhibited diffraction features similar to those of β-sheet-forming α-peptides. The diffraction data allowed the construction of a detailed model of an antiparallel α/β-peptide sheet with a unique pleated structure. One of the α/β-peptide assemblies displayed high stability, unparalleled among previously studied assemblies of α-peptides. ATR-FTIR data suggest that the 1:1 mixture of zwitterionic α/β-peptides assembled in a parallel arrangement resembling that of a typical parallel β-sheet secondary structure formed by α-peptides. This study establishes guidelines for design of amphiphilic α/β-peptides that assemble in a predictable manner at an air-water interface, with control of interstrand orientation through manipulation of Coulombic interactions along the backbone.     

A thermoreversible conformational change in two poly(zwitterions) evidenced by thermally stimulated depolarization currents

A thermally induced conformational change of zwitterionic groups of two polymers has been observed by thermally stimulated depolarization currents. Different dipole moments of the two conformations is a prerequisite. The appearance of peaks of opposite sign confirms the process. Giant dielectric constants can be obtained at room temperature using high-temperature poling and rapid cooling.     

Glutathione transferase A1-1: catalytic role of water

For more than almost 30 years now, glutathione transferases (GSTs) have been known as xenobiotic/endobiotic detoxification enzymes. GSTs catalyze the nucleophilic addition of glutathione (GSH) sulphur thiolate to a wide range of electrophilic substrates, building up a less toxic and more soluble compound, which can then be removed from the cell. Recently we proposed a consistent GSH activation mechanism. By performing QM/MM calculations, we demonstrated that a water molecule, following a first conformational rearrangement of GSH, is capable of assisting a proton transfer between the GSH thiol and alpha carboxylic groups. In this study we go further in the analysis of the water role in GSH activation by performing a long Molecular Dynamics (MD) study on glutathione transferase A1-1 Thr68 mutants complexed with GSH and the GSH decarboxylated analogue (dGSH), for which experimental kinetic data are available.     

Interpretation of synchrotron radiation circular dichroism spectra of anionic, cationic, and zwitterionic dialanine forms

Electronic absorption and synchrotron radiation circular dichroism (SRCD) spectra of the anionic, cationic, and zwitterionic forms of L-alanyl-L-alanine (AA) in aqueous solutions were measured and interpreted by molecular dynamics (MD) and ab initio computations. Time-dependent density functional theory (TD DFT) was applied to predict the electronic excited states. The modeling enabled the assessment of the role of molecular conformation, charge, and interaction with the polar environment in the formation of the spectral shapes. Particularly, inclusion of explicit solvent molecules in the computations appeared to be imperative because of the participation of water orbitals in the amide electronic structure. Implicit dielectric continuum solvent models gave inferior results for clusters, especially at low-energy transitions. Because of the dispersion of transition energies, tens of water/AA clusters had to be averaged in order to obtain reasonable spectral shapes with a more realistic inhomogeneous broadening. The modeling explained most of the observed differences, as the anionic and zwitterionic SRCD spectra were similar and significantly different from the cationic spectrum. The greatest deviation between the experimental and theoretical curves observed for the lowest-energy negative anion signal can be explained by the limited precision of the TD DFT method, but also by the complex dynamics of the amine group. The results also indicate that differences in the experimental spectral shapes do not directly correlate with the peptide main-chain conformation. Future peptide and protein conformational studies based on circular dichroic spectroscopy can be reliable only if such effects of molecular dynamics, solvent structure, and polar solvent-solute interactions are taken into account.