Cation-pi Interactions and oxidative effects on Cu+ and Cu2+ binding to Phe, Tyr, Trp, and His amino acids in the gas phase. Insights from first-principles calculations

The coordination properties of the four natural aromatic amino acids (AA(arom) = Phe, Tyr, Trp, and His) to Cu+ and Cu2+ have been exhaustively studied by means of ab initio calculations. For Cu+-Phe, Cu+-Tyr and Cu+-Trp, the two charge solvated tridentate N/O/ring and bidentate N/ring structures, with the metal cation interacting with the pi system of the ring, were found to be the lowest ones, relative DeltaG(298K) energies being less than 0.5 kcal/mol. The Cu+-His ground-state structure has the metal cation interacting with the NH2 group and the imidazole N. For these low-lying structures vibrational features are also discussed. Unlike Cu+ complexes, the ground-state structure of Cu2+-Phe, Cu2+-Tyr, and Cu2+-Trp does not present cation-pi interactions due to the oxidation of the aromatic ring induced by the metal cation. The ground-state structure of Cu2+-His does not present oxidation of the amino acid, the coordination to Cu2+ being tridentate with the oxygen of the carbonyl group, the nitrogen of the amine, and the N of the imidazole. Other less stable isomers, however, show oxidation of His, particularly of the imidazole ring, which can induce spontaneous proton-transfer reactions from the NH of the imidazole to the NH2 of the backbone. Finally, the computed binding energies for Cu+-AA(arom) and Cu2+-AA(arom) systems have been computed, the order found for the single charged systems being Cu+-His > Cu+-Trp > Cu+-Tyr > Cu+-Phe, in very good agreement with the experimental data.     

A tryptophan-containing open-chain framework for tuning a high selectivity for Ca2+ and 13C NMR observation of a Ca2+-indole interaction in aqueous solution

Two molecular architectures featuring the cation-responsive tryptophan indole were designed and investigated for the development of a novel fluorescent chemosensor for Ca2+. We observed that the Trp-based open-framework chemosensor EW2 exhibits remarkable selectivity for Ca2+ over Mg2+, Ba2+, K+, Na+, and Li+ in water between pH 4.6 and 7.0 on the basis of Ca2+-induced high fluorescence enhancement of the Trp residue. A combined 13C NMR and CD spectroscopic study has demonstrated a dynamic reorientation of the indole ring due to the cation-indole interaction accompanying the Ca2+-induced dramatic fluorescence enhancement. The results suggest that the highly sensitive, metal-ion-dependent Trp indolyl C(3) chemical shifts may serve as a promising indicator for monitoring metal ion-indole noncovalent interaction in solution.     

MP2 study of cation-(pi)n-pi interactions (n = 1-4)

Ab initio calculations at the MP2(full)/6-31++G**, RI-MP2(full)/6-31++G**, and RI-MP2(full)/6-311++G(2d,2p) levels of theory demonstrate important synergic effects between two noncovalent interactions that involve aromatic rings, that is, cation-pi and pi-pi interactions. The presence of a cation interacting with the pi cloud of an aromatic ring favors the face-to-face stacking interaction with additional aromatic rings. This effect is extended in the space up to five stacked aromatic rings.     

Biometal binding-site mimicry with modular, hetero-bifunctionally modified architecture encompassing a Trp/His motif: insights into spatiotemporal noncovalent interactions from a comparative spectroscopic study

Metal-site Trp/His interactions are crucial to diverse metalloprotein functions. This paper presents a study using metal-motif mimicry to capture and dissect the static and transient components of physicochemical properties underlying the Trp/His aromatic side-chain noncovalent interactions across the first- and second-coordination spheres of biometal ions. Modular biomimetic constructs, EDTA-(L-Trp, L-His) or EWH and DTPA-(L-Trp, L-His) or DWH, featuring a function-significant Trp/His pair, enabled extracting the putative hydrophobic/hydrophilic aromatic interactions surrounding metal centers. Fluorescence, circular dichroism (CD) spectroscopic titrations and ESI mass spectrometry demonstrated that both the constructs stoichiometrically bind to Ca(2+), Co(2+), Cu(2+), Ni(2+), Mn(2+), Zn(2+), Cd(2+), and Fe(2+), and such binding was strongly coupled to stereospecific side-chain structure reorientations of the Trp indole and His imidazole rings. A mechanistic dichotomy corresponding to the participation of the indole unit in the binding event was revealed by a scaffold-platform correlation of steady-state fluorescence-response landscape, illuminating that secondary-coordination-sphere ligand cation-π interactions were immediately followed by subsequent transient physicochemical processes including through-space energy transfer, charge transfer and/or electron transfer, depending on the type of metals. The fluorescence quenching of Trp side chain by 3d metal ions can be ascribed to through-space d-π interactions. While the fluorescence titration was capable of illuminating a two-component energetic model, clean isosbestic/isodichroic points in the CD titration spectra indicated that the metallo-constructs, such as Cu(2+)-EWH complex, fold thermodynamically by means of a two-state equilibrium. Further, the metal-ion dependence of Trp conformational variation in the modular architecture of metal-bound scaffolds was evidenced unambiguously by the CD spectra and supported by MMFF calculations; both were capable of distinguishing between the coordination geometry and the preference for metal binding mode. The study thus helps understand how aromatic rings around metal-sites have unique capabilities through the control of the spatiotemporal distribution of noncovalent interaction elements to achieve diverse chemical functionality.     

Hydration of ion-biomolecule complexes: ab initio calculations and gas-phase vibrational spectroscopy of K+ (indole)m(H2O)n

In recent years neutral indole(H2O)n clusters have been used to model the hydration of biomolecules containing an indole moiety. Both experimental and theoretical studies of the binary indole...OH2 system show NH...OH sigma hydrogen-bonding. By introducing a cation to the indole...OH2 model, cation...pi and ion...dipole electrostatic interactions are placed in direct competition with conventional indole...OH2 hydrogen-bonding. The effects, arising from a monovalent potassium cation on (indole)m(H2O)n clusters, were investigated using infrared photodissociation spectroscopy in the OH and NH stretching regions. In K+ (indole)1(H2O)(n < or = 4) and K+ (indole)2(H2O)(m < or = 3) clusters, the electrostatic ion...ligand interaction inhibits the formation of an indole NH...OH2 sigma hydrogen-bond. However, indole...H2O pi hydrogen-bonding via the five-membered indole ring is observed with three or more ligands around the ion.     

Relating Trp‐Glu dipeptide fluorescence to molecular conformation: The role of the discrete chi 1 and chi 2 angles

Molecular dynamics (MD), coupled with fluorescence data for charged dipeptides of tryptophanyl glutamic acid (Trp‐Glu), reveal a detailed picture of how specific conformation affects fluorescence. Fluorescence emission spectra and time‐resolved emission measurements have been collected for all four charged species. MD simulations 20 to 30 ns in length have also been carried out for the Trp‐Glu species, as simulation provides aqueous phase conformational data that can be correlated with the fluorescence data. The calculations show that each dipeptide species is characterized by a similar set of six, discrete Chi 1, Chi 2 dihedral angle pairs. The preferred Chi 1 angles—60°, 180°, and 300°—play the significant role in positioning the terminal amine relative to the indole ring. A Chi 1 angle of 60° results in the arching of the backbone over the indole ring and no interaction of the ring with the terminal amine. Chi 1 values of 180° and 300° result in an extension of the backbone away from the indole ring and a NH₃ cation‐π interaction with indole. This interaction is believed responsible for charge transfer quenching. Two fluorescence lifetimes and their corresponding amplitudes correlate with the Chi 1 angle probability distribution for all four charged Trp‐Glu dipeptides. Fluorescence emission band maxima are also consistent with the proposed pattern of terminal amine cation quenching of fluorescence. © 2013 Wiley Periodicals, Inc.     

Structures and energetics of the cation-pi interactions of Li+, Na+, and K+ with cup-shaped molecules: effect of ring addition to benzene and cavity selectivity

The interactions of alkali metal cations (Li (+), Na (+), and K (+)) with the cup-shaped molecules, tris(bicyclo[2.2.1]hepteno)benzene and tris(7-azabicyclo[2.2.1]hepteno)benzene have been investigated using MP2(FULL)/6-311+G(d,p)//MP2/6-31G(d) level of theory. The geometries and interaction energies obtained for the metal ion complexation with the cup-shaped systems trindene and benzotripyrrole are compared with the results for benzene-metal ion complexes to examine the effect of ring addition to the benzene on structural and binding affinities. The cup-shaped molecules exhibit two faces or cavities (top and bottom). Except for one of the conformers of tris(7-azabicyclo[2.2.1]hepteno)benzene), the metal ions prefer to bind with the top face over bottom face of the cup-shaped molecules. The selectivity of the top face is due to strong interaction of the cation with the pi cloud not only from the central six-membered ring but also from the pi electrons of rim C=C bonds. In contrast, the metal ions under study exhibit preference to bind with the bottom face rather than top face of tris(7-azabicyclo[2.2.1]hepteno)benzene) when the lone pair of electrons of three nitrogen atoms participates in binding with metal ions. This bottom face selectivity could be ascribed to the combined effect of the cation-pi and strong cation-lone pair interactions. As evidenced from the values of pyramidalization angles, the host molecule becomes deeper bowl when the lone pair of electrons of nitrogen atoms participates in binding with cation. Molecular electrostatic potential surfaces nicely explain the cavity selectivity in the cup-shaped systems and the variation of interaction energies for different ligands. Vibrational frequency analysis is useful in characterizing different metal ion complexes and to distinguish top and bottom face complexes of metal ions with the cup-shaped molecules.     

Catalytic asymmetric synthesis of protected tryptophan regioisomers

Tryptophan 1 (Trp) is superior to all other naturally occurring peptide residues in its ability to bind cations (the cation-pi interaction). In an effort to expand the toolbox of Trp-like amino acids, in this note we report catalytic asymmetric syntheses of Trp regioisomers 2a-e, where the alanine unit is attached not to C-3 of indole but to C-2, C-4, C-5, C-6, or C-7. Excellent asymmetric induction is obtained in each case (generally >97% ee). Ab initio calculations suggest that the indole nuclei of 2a-e will bind Na(+) with the same affinity as that of Trp.     

A computational study on the interaction of the nitric oxide ions NO+ and NO- with the side groups of the aromatic amino acids

The interaction of the nitric oxide ions NO+ and NO- with benzene (C6H6) and the aromatic R-groups of the amino acids phenylalanine (Phe), tyrosine (Tyr), histidine (His), and tryptophan (Trp) have been examined using the DFT method B3LYP and the conventional electron correlation method MP2. In particular, the structures and complexation energies of the resulting half-sandwich Ar...NO+/- and sandwich [Ar...NO...Ar]+/- complexes have been considered. For the Ar...NO+ complexes, the presence of an electron rich heteroatom within or attached to the ring is found to not preclude the cation...pi bound complex from being the most stable. Furthermore, unlike the anionic complexes, the pi...cation...pi ([Ar...NO...Ar]+) complexes do not correspond to a "doubling" of the parent half-sandwich.     

Cation-pi interactions with a model for the side chain of tryptophan: structures and absolute binding energies of alkali metal cation-indole complexes

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.     

Influence of the water molecule on cation-pi interaction: ab initio second order Møller-Plesset perturbation theory (MP2) calculations

The influence of introducing water molecules into a cation-pi complex on the interaction between the cation and the pi system was investigated using the MP2/6-311++G method to explore how a cation-pi complex changes in terms of both its geometry and its binding strength during the hydration. The calculation on the methylammonium-benzene complex showed that the cation-pi interaction is weakened by introducing H(2)O molecules into the system. For example, the optimized interaction distance between the cation and the benzene becomes longer and longer, the transferred charge between them becomes less and less, and the cation-pi binding strength becomes weaker and weaker as the water molecule is introduced one by one. Furthermore, the introduction of the third water molecule leads to a dramatic change in both the complex geometry and the binding energy, resulting in the destruction of the cation-pi interaction. The decomposition on the binding energy shows that the influence is mostly brought out through the electrostatic and induction interactions. This study also demonstrated that the basis set superposition error, thermal energy, and zero-point vibrational energy are significant and needed to be corrected for accurately predicting the binding strength in a hydrated cation-pi complex at the MP2/6-311++G level. Therefore, the results are helpful to better understand the role of water molecules in some biological processes involving cation-pi interactions.     

Insight into the cation-π interaction at the metal binding site of the copper metallochaperone CusF

The periplasmic Cu(+)/Ag(+) chaperone CusF features a novel cation-π interaction between a Cu(+)/Ag(+) ion and Trp44 at the metal binding site. The nature and strength of the Cu(+)/Ag(+)-Trp44 interactions were investigated using computational methodologies. Quantum-mechanical (QM) calculations showed that the Cu(+) and Ag(+) interactions with Trp44 are of similar strength (~14 kcal/mol) and bond order. Quantum-mechanical/molecular-mechanical (QM/MM) calculations showed that Cu(+) binds in a distorted tetrahedral coordination environment in the Trp44Met mutant, which lacks the cation-π interaction. Molecular dynamics (MD) simulations of CusF in the apo and Cu(+)-bound states emphasized the importance of the Cu(+)-Trp44 interaction in protecting Cu(+) from water oxidation. The protein structure does not change over the time scale of hundreds of nanoseconds in the metal-bound state. The metal recognition site exhibits small motions in the apo state but remains largely preorganized toward metal binding. Trp44 remains oriented to form the cation-π interaction in the apo state and faces an energetic penalty to move away from the metal ion. Cu(+) binding quenches the protein's internal motions in regions linked to binding CusB, suggesting that protein motions play an essential role in Cu(+) transfer to CusB.     

Insight into indole interactions from alkali metal chloride effects on a tryptophan zipper beta-hairpin peptide

Weakly solvated, low charge density, alkali metal cations (K+ and Rb+) destabilize tryptophan zipper (trpzip) peptides with an effectiveness (for Rb+) similar to that of the protein denaturant urea. An analysis of alkali metal cation effects on polypeptides stabilized predominantly either by hydrogen bonds or by the classical hydrophobic effect indicates that the alkali metals attenuate stabilizing interactions involving the tryptophan indole groups. Destabilization does not result from electrolyte screening of the electrostatic component of the indole-indole interaction, but is likely to involve direct interaction of the low charge density cation with the indole group in a cation-pi interaction. The observations highlight a general simplicity in the nature of molecular interactions in solution, in which stabilizing contributions to polypeptide and protein structures are attenuated by solutes of a complementary nature.     

Enantiomeric Excess Determination for Monosaccharides Using Chiral Transmission to Cold Gas-Phase Tryptophan in Ultraviolet Photodissociation

Chiral transmission between monosaccharides and amino acids via photodissociation in the gas phase was examined using a tandem mass spectrometer fitted with an electrospray ionization source and a cold ion trap in order to investigate the origin of the homochirality of biomolecules in molecular clouds. Ultraviolet photodissociation mass spectra of cold gas-phase noncovalent complexes of the monosaccharide enantiomers glucose (Glc) and galactose (Gal) with protonated l-tryptophan H⁺(l-Trp) were obtained by photoexcitation of the indole ring of l-Trp. l-Trp dissociated via C_α–C_β bond cleavage when noncovalently complexed with d-Glc; however, no dissociation of l-Trp occurred in the homochiral H⁺(l-Trp)(l-Glc) noncovalent complex, where the energy absorbed by l-Trp was released through the evaporation of l-Glc. This enantioselective photodissociation of Trp was due to the transmission of chirality from Glc to Trp via photodissociation in the gas-phase noncovalent complexes, and was applied to the quantitative chiral analysis of monosaccharides. The enantiomeric excess of monosaccharides in solution could be determined by measuring the relative abundance of the two product ions in a single photodissociation mass spectrum of the cold gas-phase noncovalent complex with H⁺(l-Trp), and by referring to the linear relationships derived in this work.

Open image in new window

Graphical Abstract ?

     

Nanozipper formation in the solid state from a self-assembling tripeptide with a single tryptophan residue

The self-assembly of a terminally protected tripeptide Boc-γ-Abu(1)-Ala(2)-Trp(3)-OMe (γ-Abu = γ-aminobutyric acid) 1 results in the formation of a nanostructured supramolecular zipper through various non-covalent interactions in the crystal in which the indole side-chain of the Trp(3) residue plays a key role via N-H?π interactions.     

Cation-pi interaction in model alpha-helical peptides

Cation-pi interactions are increasingly recognized as important in chemistry and biology. Here we investigate the cation-pi interaction by determining its effect on the helicity of model peptides using a combination of CD and NMR spectroscopy. The data show that a single Trp/Arg interaction on the surface of a peptide can make a significant net favorable free energy contribution to helix stability if the two residues are positioned with appropriate spacing and orientation. The solvent-exposed Trp-->Arg (i, i + 4) interaction in helices can contribute -0.4 kcal/mol to the helix stability, while no free energy gain is detected if the two residues have the reversed orientation, Arg-->Trp (i, i + 4). The derived free energy is consistent with other experimental results studied in proteins or model peptides on cation-pi interactions. However in the same system the postulated Phe/Arg (i, i + 4) cation-pi interaction provides no net free energy to helix stability. Thus the Trp-->Arg interaction is stronger than Phe-->Arg. The cation-pi interactions are not sensitive to the screening effect by adding neutral salt as indicated by salt titration. Our results are in qualitative agreement with theoretical calculations emphasizing that cation-pi interactions can contribute significantly to protein stability with the order Trp > Phe. However, our and other experimental values are significantly smaller than estimates from theoretical calculations.     

RESOLUTION OF THE INDIVIDUAL TRYPTOPHANYL CONTRIBUTIONS TO THE NEAR-ULTRAVIOLET DICHROIC ACTIVITY OF APOMYOGLOBIN

Abstract: The individual tryptophanyl contributions to the near-ultraviolet dichroic activity of apomyoglobin in its native conformation have been resolved. This was accomplished by comparing the spectra of two classes of apomyoglobin with different aromatic residue contents and observing the effect of a specific modification of indole residues. The circular dichroism (CD) spectra of apomyoglobins containing two tryptophanyl residues, i.e. Trp A-5 and A-12, show the presence of a positive peak centered at 292 nm, attributable to indolic chromophore, which is missing in the CD spectrum of tuna apomyoglobin possessing only Trp A-12. Moreover, the specific modification of Trp A-5 by 2-hydroxy-5-nitrobenzyl bromide is shown by the lack of the 292 nm peak and the appearance of a positive band at longer wavelength. The pH dependence of the position of this band suggests that it arises from the 2-hydroxy-5-nitrobenzyl moiety. The results suggest that Trp A-12 does not substantially contribute to the optical activity in the near ultraviolet.     

RESOLUTION OF THE INDIVIDUAL TRYPTOPHANYL CONTRIBUTIONS TO THE NEAR-ULTRAVIOLET DICHROIC ACTIVITY OF APOMYOGLOBIN

Abstract –The individual tryptophanyl contributions to the near-ultraviolet dichroic activity of apomyoglobin in its native conformation have been resolved. This was accomplished by comparing the spectra of two classes of apomyoglobin with different aromatic residue contents and observing the effect of a specific modification of indole residues. The circular dichroism (CD) spectra of apomyoglobins containing two tryptophanyl residues, i.e. Trp A-5 and A-12, show the presence of a positive peak centered at 292 nm, attributable to indolic chromophore, which is missing in the CD spectrum of tuna apomyoglobin possessing only Trp A-12. Moreover, the specific modification of Trp A-5 by 2-hydroxy-5-nitrobenzyl bromide is shown by the lack of the 292 nm peak and the appearance of a positive band at longer wavelength. The pH dependence of the position of this band suggests that it arises from the 2-hydroxy-5-nitrobenzyl moiety. The results suggest that Trp A-12 does not substantially contribute to the optical activity in the near ultraviolet.     

Effect of Ni(ii), Cu(ii) and Zn(ii) association on the keto-enol tautomerism of thymine in the gas phase

The effect of Ni(II), Cu(II) and Zn(II) association on the diketo/keto-enol tautomerism of thymine has been investigated through the use of B3LYP density functional theory calculations. Final energies were obtained at the B3LYP/6-311+G(2df,2p)//B3LYP/6-311+G(d,p) level of theory. Ni(II) and Cu(II) lead to an oxidation of thymine which for Zn(ii) is only partial and catalyze the tautomerization process, this catalytic effect being much larger for Ni(2+) and Zn(2+) than for Cu(2+). One of the most significant consequences of the oxidation of the base is that the calculated BDE's are primarily dictated by the value of the second ionization potential of the metal, and therefore follow the sequence Cu(2+) > Ni(2+) > Zn(2+). Also importantly, metal dication association leads to a stabilization of the keto-enol tautomer, which becomes the most stable form upon interaction with Ni(2+) and Zn(2+). This stabilization enhancement is the consequence of three concomitant factors, namely, (i) a stronger interaction of the metal cation with the carbonyl oxygen, (ii) the interaction of the metal with the dehydrogenated ring nitrogen, (iii) an aromatization of the six-membered ring.     

Competition between pi and non-pi cation-binding sites in aromatic amino acids: a theoretical study of alkali metal cation (Li+, Na+, K+)-phenylalanine complexes

To understand the cation-pi interaction in aromatic amino acids and peptides, the binding of M(+) (where M(+) = Li(+), Na(+), and K(+)) to phenylalanine (Phe) is studied at the best level of density functional theory reported so far. The different modes of M(+) binding show the same order of binding affinity (Li(+)>Na(+)>K(+)), in the approximate ratio of 2.2:1.5:1.0. The most stable binding mode is one in which the M(+) is stabilized by a tridentate interaction between the cation and the carbonyl oxygen (O[double bond]C), amino nitrogen (--NH(2)), and aromatic pi ring; the absolute Li(+), Na(+), and K(+) affinities are estimated theoretically to be 275, 201, and 141 kJ mol(-1), respectively. Factors affecting the relative stabilities of various M(+)-Phe binding modes and conformers have been identified, with ion-dipole interaction playing an important role. We found that the trend of pi and non-pi cation bonding distances (Na(+)-pi>Na(+)-N>Na(+)-O and K(+)-pi>K(+)-N>K(+)-O) in our theoretical Na(+)/K(+)-Phe structures are in agreement with the reported X-ray crystal structures of model synthetic receptors (sodium and potassium bound lariat ether complexes), even though the average alkali metal cation-pi distance found in the crystal structures is longer. This difference between the solid and the gas-phase structures can be reconciled by taking the higher coordination number of the cations in the lariat ether complexes into account.