Structures of protonated arginine dimer and bradykinin investigated by density functional theory: further support for stable gas-phase salt bridges

The gas-phase structures and energetics of both protonated arginine dimer and protonated bradykinin were investigated using a combination of molecular mechanics with conformational searching to identify candidate low-energy structures, and density functional theory for subsequent minimization and energy calculations. For protonated arginine dimer, a good correlation (R = 0.88) was obtained between the molecular mechanics and EDF1 6-31+G* energies, indicating that mechanics with MMFF is suitable for finding low-energy conformers. For this ion, the salt-bridge or ion-zwitterion form was found to be 5.7 and 7.2 kcal/mol more stable than the simple protonated or ion-molecule form at the EDF1 6-31++G** and B3LYP 6-311++G** levels. For bradykinin, the correlation between the molecular mechanics and DFT energies was poor (R = 0.28), indicating that many low-energy structures are likely passed over in the mechanics conformational searching. This result suggests that structures of this larger peptide ion obtained using mechanics calculations alone are not necessarily reliable. The lowest energy structure of the salt-bridge form of bradykinin is 10.6 kcal/mol lower in energy (EDF1) than the lowest energy simple protonated form at the 6-311G* level. Similarly, the average energy of all salt-bridge structures investigated is 13.6 kcal/mol lower than the average of all the protonated forms investigated. To the extent that a sufficient number of structures are investigated, these results provide some additional support for the salt-bridge form of bradykinin in the gas phase.     

Structure and unimolecular chemistry of protonated sulfur betaines, (CH3)2S(+)(CH2)(n)CO2H (n = 1 and 2)

The fixed charge zwitterionic sulfur betaines dimethylsulfonioacetate (DMSA) (CH(3))(2)S(+)CH(2)CO(2)(-) and dimethylsulfoniopropionate (DMSP) (CH(3))(2)S(+)(CH(2))(2)CO(2)(-) have been synthesized and the structures of their protonated salts (CH(3))(2)S(+)CH(2)CO(2)H···Cl(-) [DMSA.HCl] and (CH(3))(2)S(+)(CH(2))(2)CO(2)H···Pcr(-) [DMSP.HPcr] (where Pcr = picrate) have been characterized using X-ray crystallography. The unimolecular chemistry of the [M+H](+) of these betaines was studied using two techniques; collision-induced dissociation (CID) and electron-induced dissociation (EID) in a hybrid linear ion trap Fourier transform ion cyclotron resonance mass spectrometer. Results from the CID study show a richer series of fragmentation reactions for the shorter chain betaine and contrasting main fragmentation pathways. Thus while (CH(3))(2)S(+)(CH(2))(2)CO(2)H fragments via a neighbouring group reaction to generate (CH(3))(2)S(+)H and the neutral lactone as the most abundant fragmentation channel, (CH(3))(2)S(+)CH(2)CO(2)H fragments via a 1,2 elimination reaction to generate CH(3)S(+)=CH(2) as the most abundant fragment ion. To gain insights into these fragmentation reactions, DFT calculations were carried out at the B3LYP/6-311++G(2d,p) level of theory. For (CH(3))(2)S(+)CH(2)CO(2)H, the lowest energy pathway yields CH(3)S(+)=CH(2)via a six-membered transition state. The two fragment ions observed in CID of (CH(3))(2)S(+)(CH(2))(2)CO(2)H are shown to share the same transition state and ion-molecule complex forming either (CH(3))(2)S(+)H or (CH(2))(2)CO(2)H(+). Finally, EID shows a rich and relatively similar fragmentation channels for both protonated betaines, with radical cleavages being observed, including loss of ˙CH(3).     

Experimental and theoretical studies of potassium cation interactions with the acidic amino acids and their amide derivatives

The binding of K(+) to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four potassium cation-bound amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Formed by electrospray ionization, these complexes have energy-dependent CID cross sections that are analyzed to provide 0 K bond energies after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each K(+)(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311+G(2d,2p) basis set. Theoretical bond dissociation energies are in good agreement with the experimental values. This coordinated examination of both experimental work and quantum chemical calculations allows for a comprehensive understanding of the molecular interactions of K(+) with the Asx and Glx amino acids. K(+) binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 9+/-1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the K(+) binding affinity for the longer-chain amino acids (Glx) is enhanced by 5+/-1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced. Further, a detailed comparison between experimental and theoretical results reveals interesting differences in the binding of K(+) and Na(+) to these amino acids.     

Water molecule adsorption on protonated dipeptides

Equilibrium constants for the adsorption of the first water molecule on six protonated dipeptides (Gly-Gly+H(+), Gly-Ala+H(+), Ala-Gly+H(+), Ala-Ala+H(+), Pro-Gly+H(+), and Gly-Trp+H(+)) have been measured as a function of temperature, and DeltaH(o) and DeltaS(o) determined. Density functional theory calculations were performed for both the unsolvated peptides and the peptide water complexes at the B3LYP/6-311++G level. MP2/6-311++G** calculations were also carried out for Gly/Ala peptides. The calculations suggest that adsorption of a water molecule by these simple dipeptides is a complex process, both the unsolvated peptide and the peptide-water complexes have multiple conformations with similar free energies. Average DeltaH(o) and DeltaS(o) values derived from the calculations are in reasonable agreement with the experimental results. According to the calculations, the dominant water adsorption process involves a significant conformational change to accommodate a bridging water molecule. DeltaH(o) is diminished for Pro-Gly+H(+) mainly because the water interacts with a secondary amine, whereas for Gly-Trp+H(+), DeltaH(o) is significantly decreased by the loss of cation-pi interactions upon water adsorption. For unsolvated peptides the proton affinities of the N-terminus and the backbone carbonyl groups are known to be similar. Addition of a single water molecule causes a significant stabilization of the N-terminus protonation site.     

Protonation processes and electronic spectra of histidine and related ions

A full structural assignment of the neutral, protonated, and deprotonated histidine conformers in the gas phase is presented. A total of 3024 unique trial structures were generated by all combinations of internal single-bond rotamers of these species and optimized at the B3LYP/6-311G* level and further optimized at the B3LYP/6-311++G** level. A set of unique conformers is found, and their relative energies, free energies, dipole moments, rotational constants, electron affinities, ionization energies, and harmonic frequencies are determined. The population ratio of histidine and its tautomer is 1:0.16 at 298 K. Massive conformational changes are observed due to protonation and deprotonation, and the intramolecular H-bonds are characterized with the atoms in molecules theory. The calculated proton dissociation energy, gas-phase acidity, proton affinity, and gas-phase basicity are in excellent agreement with the experiments. The deprotonation and protonation of gaseous histidine both occur on the imidazole ring, explaining the versatile biofunctions of histidine in large biomolecules. The UV spectra of neutral and singly and doubly protonated histidine are investigated with the TDDFT/B3LYP/6-311+G(2df,p) calculations. The S0-S1, S0-S2, and S0-S3 excitations of histidine are mixed pipi*/npi* transitions at 5.37, 5.44, and 5.69 eV, respectively. The three excitation energies for histidine tautomer are 4.85, 5.47, and 5.52 eV, respectively. The three excitations for protonated histidine are mainly npi* transitions at 5.45, 5.67, and 5.82 eV, respectively. The S0-S1 excitation of protonated histidine produces ImH-CbetaH2-CalphaH(COOH)-NH2+, while the S0-S2 and S0-S3 transitions produce ImH-CbetaH2-CalphaH(NH2)-(COOH)+. These data may help to understand the mechanisms of the UV fragmentation of biomolecules.     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

Protonated nitro group: structure, energy and conjugation

Structure of protonated nitro compounds was investigated by calculations at the levels MP2(FC)/6-311++G(2d,2p)//MP2(FC)/6-311++G(2d,2p)(nitromethane and reference compounds) or B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p)(nitrobenzene and its 18 meta- and para-substituted derivatives). The group NO2H+ reveals many similarities with the isoelectronic group CO2H as the preferred conformation, conformational equilibrium, and stabilization by interaction (resonance) within the group quantified by means of isodesmic reactions. However, there is a difference in the interaction with donor groups (for instance in 4-nitroaniline) that is much stronger with NO2H+ than with CO2H. This interaction may be called resonance and may be described by standard resonance formulas, but these formulas predict only partially the geometry and cannot explain the great interaction energy.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

A novel salt bridge mechanism highlights the need for nonmobile proton conditions to promote disulfide bond cleavage in protonated peptides under low-energy collisional activation

The gas-phase fragmentation mechanisms of small models for peptides containing intermolecular disulfide links have been studied using a combination of tandem mass spectrometry experiments, isotopic labeling, structural labeling, accurate mass measurements of product ions, and theoretical calculations (at the MP2/6-311 + G(2d,p)//B3LYP/3-21G(d) level of theory). Cystine and its C-terminal derivatives were observed to fragment via a range of pathways, including loss of neutral molecules, amide bond cleavage, and S-S and C-S bond cleavages. Various mechanisms were considered to rationalize S-S and C-S bond cleavage processes, including charge directed neighboring group processes and nonmobile proton salt bridge mechanism. Three low-energy fragmentation pathways were identified from theoretical calculations on cystine N-methyl amide: (1) S-S bond cleavage dominated by a neighboring group process involving the C-terminal amide N to form either a protonated cysteine derivative or protonated sulfenyl amide product ion (44.3 kcal mol(-1)); (2) C-S bond cleavage via a salt bridge mechanism, involving abstraction of the alpha-hydrogen by the N-terminal amino group to form a protonated thiocysteine derivative (35.0 kcal mol(-1)); and (3) C-S bond cleavage via a Grob-like fragmentation process in which the nucleophilic N-terminal amino group forms a protonated dithiazolidine (57.9 kcal mol(-1)). Interestingly, C-S bond cleavage by neighboring group processes have high activation barriers (63.1 kcal mol(-1)) and are thus not expected to be accessible during low-energy CID experiments. In comparison to the energetics of simple amide bond cleavage, these S-S and C-S bond cleavage reactions are higher in energy, which helps rationalize why bond cleavage processes involving the disulfide bond are rarely observed for low-energy CID of peptides with mobile proton(s) containing intermolecular disulfide bonds. On the other hand, the absence of a mobile proton appears to "switch on" disulfide bond cleavage reactions, which can be rationalized by the salt bridge mechanism. This potentially has important ramifications in explaining the prevalence of disulfide bond cleavage in singly protonated peptides under MALDI conditions.     

Noncovalent interactions of Zn+ with N-donor ligands (pyridine, 4,4'-dipyridyl, 2,2'-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

The binding interactions in complexes of Zn(+) with nitrogen donor ligands, (N-L) = pyridine (x = 1-4), 4,4'-dipyridyl (x = 1-3), 2,2'-dipyridyl (x = 1-2), and 1,10-phenanthroline (x = 1-2), are examined in detail. The bond dissociation energies (BDEs) for loss of an intact ligand from the Zn(+)(N-L)(x) complexes are reported. Experimental BDEs are obtained from thermochemical analyses of the threshold regions of the collision-induced dissociation cross sections of Zn(+)(N-L)(x) complexes. Density functional theory calculations at the B3LYP/6-31G* level of theory are performed to determine stable structures of these species and to provide molecular parameters needed for the thermochemical analysis of experimental data. Relative stabilities of the various conformations of these N-donor ligands and their complexes to Zn(+) as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) and M06/6-311+G(2d,2p) levels of theory using the B3LYP/6-31G* optimized geometries. The experimental BDEs for the Zn(+)(N-L)(x) complexes are in reasonably good agreement with values derived from density functional theory calculations. BDEs derived from M06 calculations provide better agreement with the measured values than those based on B3LYP calculations. Trends in the sequential BDEs are explained in terms of sp polarization of Zn(+) and repulsive ligand-ligand interactions. Comparisons are made to the analogous Cu(+)(N-L)(x) and Ni(+)(N-L)(x) complexes previously studied.     

Experimental and theoretical studies of sodium cation complexes of the deamidation and dehydration products of asparagine, glutamine, aspartic acid, and glutamic acid

The deamidation and dehydration products of Na+(L), where L = asparagine (Asn), glutamine (Gln), aspartic acid (Asp), and glutamic acid (Glu), are examined in detail utilizing collision-induced dissociation (CID) with Xe in a guided ion beam tandem mass spectrometer (GIBMS). Results establish that the Na+(L) complexes decompose upon formation in our dc discharge/flow tube ion source to form a bis-ligand complex, Na+(L-HX)(HX), composed of a sodium cation, the (L-HX) decomposition product, and HX, where HX = NH3 for the amides and H2O for the acids. Analysis of the energy-dependent CID cross sections for the Na+(L-HX)(HX) complexes provides unambiguous identification of the (L-HX) fragmentation products as 3-amino succinic anhydride (a-SA) for Asx and oxo-proline (O-Pro) for Glx. Furthermore, these experiments establish the 0 K sodium cation affinities for these five-membered ring decomposition products and the H2O and NH3 binding affinities of the Na+(a-SA) and Na+(O-Pro) complexes after accounting for unimolecular decay rates, the internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations are determined for a number of geometric conformations of all reaction species as well as a number of candidate species for (L-HX) at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both the experimental work and quantum chemical calculations allows for a complete characterization of the products of deamidation and dehydration of Asx and Glx, as well as the details of Na+, H2O, and NH3 binding to the decomposition species.     

Experimental and theoretical studies of sodium cation interactions with D-arabinose, xylose, glucose, and galactose

The binding of Na (+) to arabinose (Ara), xylose (Xyl), glucose (Glc), and galactose (Gal) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated monosaccharide complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross-sections provides 0 K sodium cation affinities for experimental complexes after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-neutral collisions. Quantum chemical calculations for a number of geometric conformations of each Na (+)(L) complex with a comprehensive analysis of the alpha and beta anomeric forms are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows for determination of the bond energy for both the alpha and beta forms of each monosaccharide studied here. An understanding of the energetic contributions of individual structural characteristics as well as the energetic trends in binding among the monosaccharides is developed. Structural characteristics that affect the energetics of binding involve multidentate sodium cation coordination, ring sterics, and hydrogen bonding schemes. The overall trend in sodium binding affinities for the eight ligands follows beta-Ara < alpha-Ara < beta-Xyl < beta-Glc < alpha-Glc < alpha;-Xyl < alpha-Gal < beta-Gal.     

Effects of the residues on the excitation energies of protonated Schiff base of retinal (PSBR) in bR: a TD-DFT study

Effects of the residues on the excitation energies of protonated Schiff base of retinal (PSBR) in bacteriorhodopsin have been investigated by means of time-dependent density functional theory. The residues around PSBR are replaced by the point charges on atoms. The structures of PSBR and residues are referred from X-ray data. The atomic charges on the each residue were calculated the B3LYP/6-311G(d,p) level. The excitation energy of PSBR perturbed by the point charges on atoms of each residue was calculated at the B3LYP/6-31G(d,p) level. A total of 23 residues and five water molecules around PSBR were considered in the calculations. The large spectral shifts were caused by the Asp212 and Asp85. The origin of the spectral shifts was discussed on the basis of theoretical results.     

Fragmentation studies and electrospray ionization mass spectrometry of lapachol: protonated,deprotonated and cationized species

Electrospray ionization mass spectrometric analysis of lapachol (2‐hydroxy‐3‐(3‐methyl‐2‐butenyl)‐1,4‐naphthoquinone) was accomplished in order to elucidate the gas‐phase dissociation reactions of this important biologically active natural product. The occurrence of protonated and cationized species in the positive mode and of deprotonated species in the negative mode was explored by means of collision‐induced dissociation (CID) experiments. For the protonated molecule, the H₂O and C₄H₈ losses occur by two competitive channels. For the deprotonated molecule, the even‐electron rule is not conserved, and the radicalar species are eliminated by formation of distonic anions. The fragmentation mechanism for each ion was suggested on the basis of computational thermochemistry. Atomic charges, relative energies, and frontier orbitals were employed aiming at a better understanding of the gas‐phase reactivity of lapachol. Potential energy surfaces for fragmentation reactions were obtained by the B3LYP/6‐31+G(d,p) model. Copyright © 2010 John Wiley & Sons, Ltd.     

Dimers of boroglycine and methylamine boronic acid: a computational comparison of the relative importance of dative versus hydrogen bonding

Boronic acids are widely used in materials science, pharmacology, and the synthesis of biologically active compounds. In this Article, geometrical structures and relative energies of dimers of boroglycine, H2N-CH2-B(OH)2, and its constitutional isomer H3C-NH-B(OH)2, were computed using second-order M?ller-Plesset perturbation theory and density functional theory; Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the MP2 calculations, and the Pople 6-311++G(d,p) basis set was employed for a majority of the DFT calculations. Effects of an aqueous environment were incorporated into the results using PCM and COSMO-RS methodology. The lowest-energy conformer of the H2N-CH2-B(OH)2 dimer was a six-membered ring structure (chair conformation; Ci symmetry) with two intermolecular B:N dative-bonds; it was 14.0 kcal/mol lower in energy at the MP2/aug-cc-pVDZ computational level than a conformer with the classic eight-centered ring structure (Ci symmetry) in which the boroglycine monomers are linked by a pair of H-O...H bonds. Compared to the results of MP2 calculations with correlation-consistent basis sets, DFT calculations using the PBE1PBE and TPSS functionals with the 6-311++G(d,p) basis set were significantly better at predicting relative conformational energies of the H2N-CH2-B(OH)2 and H3C-NH-B(OH)2 dimers than corresponding calculations using the BLYP, B3LYP, OLYP, and O3LYP functionals, particularly with respect to dative-bonded structures.     

类锗烯H2C=GeLiCl与RH(R=F,OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了不饱和类锗烯H2C=GeLiCl与RH(R=F, OH, NH2)的插入反应. 在B3LYP/6-311+G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C=GeLiCl与HF、H2O 或NH3发生插入反应的机理相同. QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为173.53、194.48和209.05 kJ·mol-1, 反应热分别为60.18、72.93和75.34 kJ·mol-1. 相同条件下发生插入反应时, 反应活性顺序都是H—F>H—OH>H—NH2.