Extreme stability of an unsolvated alpha-helix

High-temperature ion mobility measurements have been performed for alpha-helical Ac-A15K+H+ and globular Ac-KA15+H+ peptides. The alpha-helical and globular conformations do not melt into random coils as the temperature is raised. Instead, both conformations survive to the point where the peptide signals vanishes due to fragmentation. This occurs at 600 K for the globular Ac-KA15+H+ peptide and at 725 K for the alpha-helical Ac-A15K+H+. For the helical Ac-A15K+H+ peptide it appears that fragmentation is triggered by disruption of the helical conformation.     

Contribution of nonlocal interactions to DNA elasticity

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.     

Effect of thioxopeptide bonds on alpha-helix structure and stability

Thioxoamide (thioamide) bonds are nearly isosteric substitutions for amides but have altered hydrogen-bonding and photophysical properties. They are thus well-suited backbone modifications for physicochemical studies on peptides and proteins. The effect of thioxoamides on protein structure and stability has not been subject to detailed experimental investigations up to date. We used alanine-based model peptides to test the influence of single thioxoamide bonds on alpha-helix structure and stability. The results from circular dichroism measurements show that thioxoamides are strongly helix-destabilizing. The effect of an oxo-to-thioxoamide backbone substitution is of similar magnitude as an alanine-to-glycine substitution resulting in a helix destabilization of about 7 kJ/mol. NMR characterization of a helical peptide with a thioxopeptide bond near the N-terminus indicates that the thioxopeptide moiety is tolerated in helical structures. The thioxoamide group is engaged in an i, i+4 hydrogen bond, arguing against the formation of a 3(10)-helical structure as suggested for the N-termini of alpha-helices in general and for thioxopeptides in particular.     

Salt-specific stability and denaturation of a short salt-bridge-forming alpha-helix

The structure of a single alanine-based Ace-AEAAAKEAAAKA-Nme peptide in explicit aqueous electrolyte solutions (NaCl, KCl, NaI, and KF) at large salt concentrations (3-4 M) is investigated using approximately 1 mus molecular dynamics (MD) computer simulations. The peptide displays 71% alpha-helical structure without salt and destabilizes with the addition of NaCl in agreement with experiments of a somewhat longer version. It is mainly stabilized by direct and indirect (" i + 4")EK salt bridges between the Lys and Glu side chains and a concomitant backbone shielding mechanism. NaI is found to be a stronger denaturant than NaCl, while the potassium salts hardly show influence. Investigation of the molecular structures reveals that consistent with recent experiments Na (+) has a much stronger affinity to side chain carboxylates and backbone carbonyls than K (+), thereby weakening salt bridges and secondary structure hydrogen bonds. At the same time, the large I (-) has a considerable affinity to the nonpolar alanine in line with recent observations of a large propensity of I (-) to adsorb to simple hydrophobes, and thereby "assists" Na (+) in its destabilizing action. In the denatured states of the peptide, novel long-lived (10-20 ns) "loop" configurations are observed in which single Na (+) ions and water molecules are hydrogen-bonded to multiple backbone carbonyls. In an attempt to analyze the denaturation behavior within the preferential interaction formalism, we find indeed that for the strongest denaturant, NaI, the protein is least hydrated. Additionally, a possible indication for protein denaturation might be a preferential solvation of the peptide backbone by the destabilizing cosolute (sodium). The mechanisms found in this work may be of general importance to understand salt effects on protein secondary structure stability.     

Hydrophobic side-chain interactions in a family of dimeric amide foldamers-potential alpha-helix mimetics

A series of new alpha-helix mimetics based on a benzamide scaffold and potentially able to disrupt protein-protein interactions have been synthesized and characterized by X-ray analysis. Inspection of the solid state structures of aromatic amide dimers confirmed that the molecules adopt a curved conformation with intramolecular H-bonding between the amide NH and the alkoxy oxygen of the neighboring aromatic fragment (d_NH…_O ∼ 2 Å). Adjacent dimer molecules are prone to form supramolecular assemblies due to both hydrophobic alkyl side-chain/side-chain interactions and intermolecular H-bonding.     

Deciphering aromatic oligoamide foldamer-DNA interactions

Finest selection: Side-chain selective, end-group selective, diastereoselective, and RNA- vs. DNA-selective interactions have been revealed between multiturn helical aromatic amide foldamers having cationic side chains and G-quadruplex aptamers.     

Electrostatic control of aromatic stacking interactions

A supramolecular approach has been used to investigate the free energies of intermolecular aromatic stacking interactions. Chemical double mutant cycles have been used to measure the effect of a range of substituents on face-to-face stacking interactions with phenyl and pentafluorophenyl rings. Electrostatic effects dominate the trends in interaction energy.     

Nativelike structure in designed four alpha-helix bundles driven by buried polar interactions

We show that a single internal polar interaction per helix is sufficient to engender structural specificity in that helix in helical bundle proteins. Furthermore, we use histidine-binding cofactors of different shapes which bind directly into the core, demonstrating that this structural specificity is not the result of a prescribed complimentary, "knobs in holes" core packing. We show that we can switch structural specificity of individual helices on and off by ligating cofactors, singly and in pairs, which bind either one or two histidine ligands. To our knowledge, this is the first demonstration of such extensive manipulation of protein structure by ligand binding, an important result of general interest to those working with self-assembled molecular systems. Finally, as these proteins were designed without the use of computational modeling, we not only demonstrate that designing a uniquely structured cofactor binding protein is not as difficult as is generally believed, we have determined why this is so: hydrophobic core complementarity, which is very difficult to design, is not necessary. Instead, a much simpler design process entails the creation of core polar interactions which themselves can drive conformational specificity.     

On the stability of a single-turn alpha-helix: the single versus multiconformation problem

The pentapeptide Ac-HAAAH-NH2, cyclized through its imidazoles by PdII to give [Pd(en)(peptide)]2+, has recently been evaluated by 2-D NMR and simulated annealing as a single alpha-helix conformation in solution. In the present work, we have questioned this assumption by developing Pd2+ parameters for AMBER*, performing an extensive conformational search for the [Pd(en)(peptide)]2+, and deconvoluting the averaged NMR data into eight rapidly equilibrating conformations with populations ranging from 2 to 55%. None of the latter correspond to the alpha-helix, although a 3% form possesses a related structure. As a critical component of interpreting an averaged NMR spectrum in terms of a single conformation, we advise testing this assumption with a method that permits conformational deconvolution.     

Substituent effects on edge-to-face aromatic interactions

Chemical double mutant cycles have been used to measure the magnitude of edge-to-face aromatic interactions in hydrogen-bonded zipper complexes as a function of substituents on both aromatic rings. The interaction energies vary depending on the combination of substituents from +1.0 kJ mol-1 (repulsive), to -4.9 kJ mol-1 (attractive). The results correlate with the Hammett substituent constants which indicates that electrostatic interactions are responsible for the observed differences in interaction energy. The experiments can be rationalised based on local electrostatic interactions between the protons on the edge ring and the pi-electron density on the face ring as well as global electrostatic interactions between the overall dipoles on the two aromatic groups.     

Selective aromatic interactions in beta-hairpin peptides

To probe the selectivity possible in hydrophobic clusters, we have compared the cross-strand interactions of phenylalanine (Phe) and cyclohexylalanine (Cha) in a beta-hairpin peptide. We have found a preference for self-association among the aromatic residues, which provides 0.55 kcal/mol in stability relative to Cha-Cha cross-strand pair. NMR analysis of the Phe-Phe cross-strand pair indicates that it interacts in an edge-face interaction, despite the fact that it is highly solvent-exposed. The interaction geometry as well as the enthalpic and entropic values for the peptide containing the Phe-Phe cross-strand pair suggest that the preference for self-association arises from inherent differences in the nature of aromatic and aliphatic interactions in water.     

Substituent effects on aromatic stacking interactions

Synthetic supramolecular zipper complexes have been used to quantify substituent effects on the free energies of aromatic stacking interactions. The conformational properties of the complexes have been characterised using NMR spectroscopy in CDCl(3), and by comparison with the solid state structures of model compounds. The structural similarity of the complexes makes it possible to apply the double mutant cycle method to evaluate the magnitudes of 24 different aromatic stacking interactions. The major trends in the interaction energy can be rationalised using a simple model based on electrostatic interactions between the pi-faces of the two aromatic rings. However, electrostatic interactions between the substituents of one ring and the pi-face of the other make an additional contribution, due to the slight offset in the stacking geometry. This property makes aromatic stacking interactions particularly sensitive to changes in orientation as well as the nature and location of substituents.     

Intermolecular interactions in nitrogen-containing aromatic systems

π–π and CH···N interactions are vital in biological systems. In this study, stacking and hydrogen-bonded interactions in pyrazine and triazine dimers were investigated by density functional theory combined with symmetry-adapted perturbation theory (DFT-SAPT) and counterpoise (CP)-corrected supermolecular MP2, SCS-MP2, B3LYP-D and CCSD(T) calculations. All interaction energies were computed using the optimized structures at the CP-corrected SCS/aug-cc-pVDZ level, which gave 1–2 kJ/mol lower interaction energies than the ones computed at the MP2 level. For both dimers, doubly hydrogen-bonded and cross-(displaced) stacked orientations were found to be the lowest energy ones. The reference CCSD(T) calculations favored the former structure in both dimer systems, whereas MP2 and SCS-MP2 located the latter as the lowest energy isomer. In particular, the former was found to be lower in energy than the latter by 2.28 and 1.01 kJ/mol at the CCSD(T)/aug-cc-pVDZ level for pyrazine and triazine, respectively. B3LYP-D produced interaction energies in agreement with the CCSD(T) at the equilibrium geometries, but it overestimates them at the short range and underestimates at the long intermonomer separations. Furthermore, it tends to give smaller equilibrium distances compared to the CCSD(T). DFT-SAPT method was in a good agreement with the reference CCSD(T) calculations. This suggests that DFT-SAPT can be employed to compute the full potential energy surface of these dimers. Moreover, DFT-SAPT calculations showed that the electrostatic and dispersion contributions are the most important energy components stabilizing these dimers. The present study aims to show which theoretical method is the most promising one for the investigation of intermolecular interactions dominated by π–π and CH···N. Therefore, the findings obtained in this study can be used to unravel the structures of nucleic acid bases and other systems stabilized by π–π and CH···N interactions.     

Metal ligand aromatic cation-pi interactions in metalloproteins: ligands coordinated to metal interact with aromatic residues

Cation-pi interactions between aromatic residues and cationic amino groups in side chains and have been recognized as noncovalent bonding interactions relevant for molecular recognition and for stabilization and definition of the native structure of proteins. We propose a novel type of cation-pi interaction in metalloproteins; namely interaction between ligands coordinated to a metal cation--which gain positive charge from the metal--and aromatic groups in amino acid side chains. Investigation of crystal structures of metalloproteins in the Protein Data Bank (PDB) has revealed that there exist quite a number of metalloproteins in which aromatic rings of phenylalanine, tyrosine, and tryptophan are situated close to a metal center interacting with coordinated ligands. Among these ligands are amino acids such as asparagine, aspartate, glutamate, histidine, and threonine, but also water and substrates like ethanol. These interactions play a role in the stability and conformation of metalloproteins, and in some cases may also be directly involved in the mechanism of enzymatic reactions, which occur at the metal center. For the enzyme superoxide dismutase, we used quantum chemical computation to calculate that Trp163 has an interaction energy of 10.09 kcal mol(-1) with the ligands coordinated to iron.     

Investigations on the chromatographic behavior of hybrid reversed-phase materials containing electron donor-acceptor systems. II. Contribution of pi-pi aromatic interactions

The chromatographic behavior of three naphthalimide-type stationary phases were elucidated in terms of hydrophobic, silanophilic and pi-pi interaction properties, employing besides common chromatographic column test methods from Engelhardt and Tanaka, also new test mixtures of geometrical and functional aromatic isomers. It was found that the presence of electron donor/acceptor moieties within a reversed phase system did not only increase the overall retention times for aromatic solutes, but also lead to an enhanced shape selectivity of the hybrid stationary phase. In this context, shape discrimination is primarily based on the number of accessible pi-electrons for pi-pi interaction with the embedded electron deficient ligand moieties. The most outstanding results were obtained for the 1,4,5,8-naphthalenediimidic selector with its horizontal arrangement on the silica surface, which enables a direct face-to-face pi-pi interaction with aromatic solutes, with only little hydrophobic contribution.     

An evaluation of force-field treatments of aromatic interactions

Experimental measurements of edge-to-face aromatic interactions have been used to test a series of molecular mechanics force fields. The experimental data were determined for a range of differently substituted aromatic rings using chemical double mutant cycles on hydrogen-bonded zipper complexes. These complexes were truncated for the purposes of the molecular mechanics calculations so that problems of conformational searching and the optimisation of large structures could be avoided. Double-mutant cycles were then carried out in silico using these truncated systems. Comparison of the experimental aromatic interaction energies and the X-ray crystal structures of these truncated complexes with the calculated data show that conventional molecular mechanics force fields (MM2, MM3, AMBER and OPLS) do not perform well. However, the XED force field which explicitly represents electron anisotropy as an expansion of point charges around each atom reproduces the trends in interaction energy and the three-dimensional structures exceedingly well. Collapsing the XED charges onto atom centres or the use of semi-empirical atom-centred charges within the XED force field gives poor results. Thus the success of XED is not related to the methods used to assign the atomic charge distribution but can be directly attributed to the use of off-atom centre charges.     

Interplay between edge-to-face aromatic and hydrogen-bonding interactions

The interplay between two important noncovalent interactions involving aromatic rings is studied by means of MP2/6-31++G** ab initio calculations. They indicate that synergistic effects are present in complexes where edge-to-face aromatic interactions and hydrogen-bonding interactions coexist. These synergistic effects have been studied bu using the atoms in molecules theory and the molecular interaction potential with polarization partition scheme. Experimental evidence for such interactions has been obtained from the Cambridge Structural Database.     

Substituent effects on the edge-to-face aromatic interactions

The edge-to-face interactions for either axially or facially substituted benzenes are investigated by using ab initio calculations. The predicted maximum energy difference between substituted and unsubstituted systems is approximately 0.7 kcal/mol (approximately 1.2 kcal/mol if substituents are on both axially and facially substituted positions). In the case of axially substituted aromatic systems, the electron density at the para position is an important stabilizing factor, and thus the stabilization/destabilization by substitution is highly correlated to the electrostatic energy. This results in its subsequent correlation with the polarization and charge transfer. Thus, the stabilization/destabilization by substitution is represented by the sum of electrostatic energy and induction energy. On the other hand, the facially substituted aromatic system depends on not only the electron-donating ability responsible for the electrostatic energy but also the dispersion interaction and exchange repulsion. Although the dispersion energy is the most dominating interaction in both axial and facial substitutions, it is almost canceled by the exchange repulsion in the axial substitution, whereas in the facial substitution, together with the exchange repulsion it augments the electrostatic energy. The systems with electron-accepting substituents (NO2, CN, Br, Cl, F) favor the axial substituent conformation, while those with electron-donating substituents (NH2, CH3, OH) favor the facial substituent conformation. The interactions for the T-shape complex systems of an aromatic ring with other counterpart such as H2, H2O, HCl, and HF are also studied.