Comparison of various types of hydrogen bonds involving aromatic amino acids

Ab initio calculations are used to compare the abilities of the aromatic groups of the Phe, Tyr, Trp, and His amino acids (modeled respectively by benzene, phenol, indole, and imidazole) to form H-bonds of three different types. Strongest of all are the conventional H-bonds (e.g., OH..O and OH..N). His forms the strongest such H-bond, followed by Tyr, and then by Trp. Whereas OH..phi bonds formed by the approach of a proton donor to the pi electron cloud above the aromatic system are somewhat weaker, they nonetheless represent an important class of stabilizing interactions. The strengths of H-bonds in this category follow the trend Trp > His > Tyr approximately Phe. CH.O interactions are weaker still, and only those involving His and Trp are strong enough to make significant contributions to protein structure. A protonated residue such as HisH(+) makes for a very powerful proton donor, such that even its CH..O H-bonds are stronger than the conventional H-bonds formed by neutral groups.     

Synthesis of pyridazines functionalized with amino acid side chains

A simple route for the preparation of 3,4,6-substituted pyridazines is described by using Tebbe olefination of esters then Diels–Alder reaction of the resulting enol ethers with tetrazine.     

Statistically based reduced representation of amino acid side chains

Preferred conformations of amino acid side chains have been well established through statistically obtained rotamer libraries. Typically, these provide bond torsion angles allowing a side chain to be traced atom by atom. In cases where it is desirable to reduce the complexity of a protein representation or prediction, fixing all side-chain atoms may prove unwieldy. Therefore, we introduce a general parametrization to allow positions of representative atoms (in the present study, these are terminal atoms) to be predicted directly given backbone atom coordinates. Using a large, culled data set of amino acid residues from high-resolution protein crystal structures, anywhere from 1 to 7 preferred conformations were observed for each terminal atom of the non-glycine residues. Side-chain length from the backbone C(alpha) is one of the parameters determined for each conformation, which should itself be useful. Prediction of terminal atoms was then carried out for a second, nonredundant set of protein structures to validate the data set. Using four simple probabilistic approaches, the Monte Carlo style prediction of terminal atom locations given only backbone coordinates produced an average root mean-square deviation (RMSD) of approximately 3 A from the experimentally determined terminal atom positions. With prediction using conditional probabilities based on the side-chain chi(1) rotamer, this average RMSD was improved to 1.74 A. The observed terminal atom conformations therefore provide reasonable and potentially highly accurate representations of side-chain conformation, offering a viable alternative to existing all-atom rotamers for any case where reduction in protein model complexity, or in the amount of data to be handled, is desired. One application of this representation with strong potential is the prediction of charge density in proteins. This would likely be especially valuable on protein surfaces, where side chains are much less likely to be fixed in single rotamers. Prediction of ensembles of structures provides a method to determine the probability density of charge and atom location; such a prediction is demonstrated graphically.     

Hydrogen bonding between phosphate and amino acid side chains

Hydrogen bonds between polar groups of amino acid side chains (histidine, lysine, glutamic acid) and phosphate ions have been studied by infrared spectroscopy. Proton transfer from amino acid groups to phosphate occur mainly in case that tribasic and dibasic phosphate ions take part in hydrogen bonds. Conformational changes and continuum are strongly related to the degree of proton transfer and hydration. It is pointed out that the aforementioned properties should be of great significance for nucleation and growth of prostatic and renal stones.     

Potentials of mean force between ionizable amino acid side chains in water

Potentials of mean force (PMF) between all possible ionizable amino acid side chain pairs in various protonation states were calculated using explicit solvent molecular dynamics simulations with umbrella sampling and the weighted histogram analysis method. The side chains were constrained in various orientations inside a spherical cluster of 200 water molecules. Beglov and Roux's Spherical Solvent Boundary Potential was used to account for the solvent outside this sphere. This approach was first validated by calculating PMFs between monatomic ions (K(+), Na(+), Cl(-)) and comparing them to results from the literature and results obtained using Ewald summation. The strongest interaction (-4.5 kcal/mol) was found for the coaxial Arg(+).Glu(-) pair. Many like-charged side chains display a remarkable lack of repulsion, and occasionally a weak attraction. The PMFs are compared to effective energy curves obtained with common implicit solvation models, namely Generalized Born (GB), EEF1, and uniform dielectric of 80. Overall, the EEF1 curves are too attractive, whereas the GB curves in most cases match the minima of the PMF curves quite well. The uniform dielectric model, despite some fortuitous successes, is grossly inadequate.     

9-BBN: An amino acid protecting group for functionalization of amino acid side chains in organic solvents

9-Borabicyclononane (9-BBN) has been utilized to protect functionalized amino acids for potential chemoselective side chain manipulation. The 9-BBN group imparts organic solubility to otherwise hydrophilic molecules and is tolerant of a wide range of reaction conditions. The high degree of solubility of these molecules in THF is particularly noteworthy. It is cleaved with either aqueous HCl or by exchange with ethylenediamine in methanol. [reaction: see text]     

Role of pyridine hydrogen-bonding sites in recognition of basic amino acid side chains

A series of three, new artificial receptors for guanidinium and ammonium guests has been synthesized. All three receptors have highly preorganized clefts bearing two carboxylate groups. They differ in the number of nitrogen atoms contained in their clefts, as follows: four N atoms in receptor 3, three N atoms in 4, and two nitrogens in 5. Crystallographic studies have produced the solid-state structures of the following guanidinium complexes of each receptor: 3.2CH(3)CH(2)NHC(NH(2))(2)(+), 4.2CH(3)NHC(NH(2))(2)(+), and 5.2C(NH(2))(3)(+). The conformations of the receptor molecules in all three complexes are very similar. N-Alkylguanidinium guests are bound in the clefts of 3 and 4 in similar manners, despite the loss of one hydrogen-bond acceptor nitrogen in 4 and the possible hindrance of the cavity by a CH group. In the guanidinium complex of 5, neither guest enters the cavity containing two CH groups. Complexation studies were conducted in methanol by (1)H NMR titration for several guanidinium and ammonium guests, including derivatives of the amino acids arginine and lysine. Receptor 5 binds all such guests weakly (K(s) < 4000), while 3 binds most guests very strongly (K(s) > 100 000). Receptor 3 is selective for arginine versus lysine, while 4 binds lysine better than does 3. The results generally underscore the importance of receptor preorganization and hydrogen-bonding complementarity in the design of receptors that can serve as probes for biomolecules.     

The influence of Na+ on neighbouring hydrogen bonds of aliphatic amino acid

Summary The influence of Na⁺ on hydrogen bonds of the OH O and NH O type between an aliphatic amino acid (glycine zwitterion) and water is investigated byab initio calculations with minimal Gaussian basis sets. Distortion of the hydration shell caused by Na⁺, and interaction energies contributing to the over-all stabilization are discussed.

---

Der Einfluß von Na⁺ auf die benachbarten Wasserstoffbindungen in aliphatischen Aminosäuren

Zusammenfassung Der Einfluß von Na⁺ auf die Wasserstoffbindungen vom OH O- und NH O-Typ in aliphatischen Aminosäuren (Glycin-Zwitterion) und Wasser wurde mittelsab initio Berechnungen mit einem minimalen Gausschen Basissatz untersucht. Die durch Na⁺-Ionen hervorgerufenen Verzerrungen der Hydratationsschale und die zur Gesamtstabilisierung beitragenden Wechselwir-kungsenergien werden diskutiert.

     

Electronic properties of amino acid side chains: quantum mechanics calculation of substituent effects

Background  

Electronic properties of amino acid side chains such as inductive and field effects have not been characterized in any detail. Quantum mechanics (QM) calculations and fundamental equations that account for substituent effects may provide insight into these important properties. PM3 analysis of electron distribution and polarizability was used to derive quantitative scales that describe steric factors, inductive effects, resonance effects, and field effects of amino acid side chains.     

Cooperativity of conventional and unconventional hydrogen bonds involving imidazole

Correlated ab initio calculations are used to investigate the cooperativity of H‐bonds between imidazole and a pair of water molecules. H‐bonds comprise not only the conventional NH … O and OH … N types, but also CH … O and OH … φ (wherein a proton is donated to the delocalized π cloud lying above the aromatic ring). Conventional and OH … φ H‐bonds obey the normal principles of cooperativity, wherein these bonds are strengthened when a central molecule serves simultaneously as both proton donor and acceptor. In contrast, CH … O bonds do not appear to be amenable to such positive cooperativity. When placed in a polarizable medium, all H‐bonds weaken as the dielectric constant of the solvent grows. The qualitative aspects of the cooperativity are not affected by the medium, although some weakening is observed. Calculations also consider the effects of cooperativity on other aspects of the complexes, including the intermolecular distance, the effect on the covalent X‐H bond length, and IR and NMR spectral data. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Calculation of the free energy of solvation for neutral analogs of amino acid side chains

The ability of the GROMOS96 force field to reproduce partition constants between water and two less polar solvents (cyclohexane and chloroform) for analogs of 18 of the 20 naturally occurring amino acids has been investigated. The estimations of the solvation free energies in water, in cyclohexane solution, and chloroform solution are based on thermodynamic integration free energy calculations using molecular dynamics simulations. The calculations show that while the force field reproduces the experimental solvation free energies of nonpolar analogs with reasonable accuracy the solvation free energies of polar analogs in water are systematically overestimated (too positive). The dependence of the calculated free energies on the atomic partial charges was also studied.     

Array of aromatic amino acid side chains located near the chromophore of photoactive yellow protein

The role of the array of aromatic amino acid side chains located close to the chromophore binding loop of photoactive yellow protein (PYP) was studied using the alanine-substitution mutagenesis. Phe92, Tyr94, Phe96 and Tyr98 were replaced with alanine (F92A, Y94A, F96A and Y98A, respectively), then these mutants were characterized by UV-visible absorption spectra, circular dichroism (CD) spectra, thermal stability and photocycle kinetics. Absorption maxima of F92A, Y94A, F96A and Y98A were 444, 442, 439 and 447 nm, respectively, different to wild type (WT) at 446 nm. Far-UV CD spectra of mutants other than F92A were different from WT, indicating that Tyr94, Phe96 and Tyr98 maintain the native secondary structure of PYP. Mid-point temperatures of thermal denaturation of F92A, Y94A and F96A, estimated by the CD signal at 222 nm, were 5-10 degrees C lower than WT. Time constants of the photocycle estimated by flash-induced absorbance change were 0.36 s for WT and 1.4 s for Y98A, however, 100, 30 and 3000 times slower than WT for F92A, Y94A and F96A, respectively. Tyr98 is located in the loop region, whereas Phe92, Tyr94 and Phe96 are incorporated in the beta4 strand, showing that aromatic amino acid residues in the beta-sheet regulate the absorption spectrum, thermal stability and photocycle of PYP. Aromatic rings of Phe92, Tyr94 and Phe96 lie nearly perpendicular to the aromatic ring of Phe75 or chromophore. Possible weak hydrogen bonds between the aromatic ring hydrogen and pi-electrons of these residues are discussed.     

Ab initio study of hydrogen-bond formation between cyclic ethers and selected amino acid side chains

Binding energies for hydrogen-bonded complexes of six cyclic ethers with five hydrogen-bond donor molecules that mimic selected amino acid side chains have been calculated at the MP2/6-31G*, MP2/6-31+G*, MP2/6-311++G**(single point), and MP2/aug-cc-pvtz levels, using geometries obtained with or without counterpoise corrections throughout the geometry optimization. The calculated basis set superposition error (BSSE) amounts to 10-20% and 5-10% of the uncorrected binding energies for the neutral and ionic species, respectively, at the MP2/aug-cc-pvtz level. The authors conclude that the O...H distances in the hydrogen bonds and binding energies for the studied systems may be determined with uncertainties of up to 0.08 A and 1-2 kcal/mol, respectively, in comparison with the MP2/aug-cc-pvtz values at a reasonable computational cost by performing standard geometry optimization at the MP2/6-31+G* level. Hydrogen-bond formation energies are more negative for cyclic ethers compared to their counterparts with a C=C double bond in the ring next to the oxygen atom. The less negative hydrogen-bonding energy and the increased O...H separation have been attributed to the reduced basicity of the ether oxygen when the lone pairs can enter conjugation with the pi-electrons of the Calpha=Cbeta double bond. The present study is the first step toward the development of an affordable computational level for estimating the binding energies of natural product, fused ring ether systems to the human estrogen receptor.     

Biaryl-bridged macrocyclic peptides: conformational constraint via carbogenic fusion of natural amino acid side chains

A general method for constraining peptide conformations via linkage of aromatic sidechains has been developed. Macrocyclization of suitably functionalized tri-, tetra- and pentapeptides via Suzuki-Miyaura cross-coupling has been used to generate side chain to side chain, biaryl-bridged 14- to 21-membered macrocyclic peptides. Biaryl bridges possessing three different configurations, meta-meta, meta-ortho, and ortho-meta, were systematically explored through regiochemical variation of the aryl halide and aryl boronate coupling partners, allowing fine-tuning of the resultant macrocycle conformation. Suzuki-Miyaura macrocyclizations were successfully achieved both in solution and on solid phase for all three sizes of peptide. This approach constitutes a means of constraining peptide conformation via direct carbogenic fusion of side chains of naturally occurring amino acids such as phenylalanine and tyrosine, and so is complementary to strategies involving non-natural, for example, hydrocarbon, bridges.     

Intermolecular stretching vibration of hydrogen bonds involving a sulfur atom

The intermolecular stretching vibration of hydrogen bonded complexes between phenol derivatives and tetramethylthiourea or dimethylthioformamide is observed in the range 113 – 131 cm⁻¹. The corresponding force constants k_σ are computed using the Lippincott-Schroeder potential function and the matrix method of Wilson. The k_σ values are higher than these previously reported for the complexes of the same phenols with aliphatic amides. These results are in qualitative agreement with the theoretical predictions.     

Infrared study of hydrogen bonds involving N-heterocyclic bases and phenols

The hydrogen bond complexes between phenols and N-heteroaromatic bases 2,4,6-tri(2-pyridyl)-1,3,5-triazine, 2,2′,2′-terpyridine, quinoxaline, pyrido[2,3-b]pyrazine, pyzazino [2,3f]quinoxaline and 5-nitrozphenanthroline are investigated by infrared spectroscopy in 1,2-dichloroethane. The stability constants of the complexes involving N-heteroaromatic bases characterized by tow vicinal nitrogen atoms having lone pairs pointing to each other are higher than predicted from their basicity. Possible differences between protonation and hydrogen bond formation are discussed. Nheteroaromatic bases such as tri(2-pyridyl)-1,3,5-triazine or phenanthrolines cannot be considered as proton sponges but their behaviour is intermediate between that of the classical heteroaromatic bases and the proton sponges.     

Syntheses and structural studies of Schiff bases involving hydrogen bonds

New Schiff bases have been prepared by reacting 3-hydroxy-4-pyridine- carboxaldehyde with various amines. NMR spectroscopic methods provided clear evidence that the Schiff bases exist in the solid state and in solution as hydroxyimino tautomers with the E-configuration. A study of the stabilities of the tautomeric forms and the different conformers has been carried out using density functional calculations at the B3LYP/6-31G** level.     

Understanding hydrogen bonding of hydroxamic acids with some amino acid side chain model molecules

The hydrogen-bonding abilities of a few amino acid side chains have been studied through aggregation of methylamine, methanol, and acetic acid (as model molecules) with formo- and thioformo- hydroxamic acids using ab initio calculations. Forty six aggregates representing all possible H-bond interactions between these amino acid side chain groups and two most stable keto and enol tautomeric forms of both hydroxamic acids have been optimized. Although participation of conventional H-bond donors and acceptors leads to significant stabilization energies, yet C–H···O, C–H···N, S–H···O, and S–H···N etc. unconventional H-bonds also contribute to stabilize interactions in many aggregates. Strength of H-bonds of the molecules with formo- and thioformo- hydroxamic acid studied follows the order acetic acid > methylamine > methanol. A comparative study of atomic charges and orbital interactions employing NBO analysis has been carried out to explore the role of bond polarizations, charge transfer, and electron delocalizations as contributors to stabilization energy.