Potential functions for hydrogen bond interactions - IV. Minimum Energy Conformation of the α-Helical Structure of PoIy-L-Alanine

Making use of the empirical potential functions for peptide NH .. O bonds, developed in this laboratory, the relative stabilities of the rightand left-handed α-helical structures of poly-L-alanine have been investigated, by calculating their conformational energies (V). The value of Vmin of the right-handed helix (αP) is about - 10.4 kcal/mole, and that of the left-handed helix (αM) is about - 9.6 kcal/mole, showing that the former is lower in energy by 0.8 kcal/mole. The helical parameters of the stable conformation of αP are n ∼ 3.6 and h ∼ 1.5 Å. The hydrogen bond of length 2.85 Å and nonlinearity of about 10° adds about 4.0 kcal/ mole to the stabilising energy of the helix in the minimum enregy region. The energy minimum is not sharply defined, but occurs over a long valley, suggesting that a distribution of conformations (φ{symbol}, ψ) of nearly the same energy may occur for the individual residues in a helix. The experimental data of a-helical fibres of poly-L-alanine are in good agreement with the theoretical results for αP. In the case of proteins, the mean values of (φ{symbol}, ψ) for different helices are distributed, but they invariably occur within the contour for V = Vmin + 2 kcal/mole for αP.

Vibration spectra of ferroelectric alkali trihydrogen selenites and the nature of the hydrogen bond potential in these crystals

Raman spectra of the ferroelectric LiH3 (SeO3)2 and NaH3(SeO3)2 and the anti-ferroelectric KH3 (SeO3)2 have been recorded at room temperature using a He-Ne and also an Ar-ion laser source. The infrared absorption spectra of these crystals and their deuterated analogues have been recorded in the region 400–4000 cm−1 both below and above the Curie temperature. From an analysis of the spectrum in the region 400–900 cm−1 it is concluded that (i) in LiH3 (SeO3)2 the protons are ordered in an asymmetric double minimum potential with a low barrier and the spectrum can be interpreted in terms of HSeO3− and H2SeO3 vibrations, (ii) in NaH3 (SeO3)2 all three protons occupy a single minimum potential at room temperature and below the transition temperature the groups HSeO3− and H2SeO3 are present, (iii) the proton at the inversion centre in KH3(SeO3)2 is in a broad troughed potential well and the low temperature spectrum is more likely to be due to H3SeO3+ and SeO32− species. This deviation of the spectrum from that of the previous two crystals is attributed to the difference in H-bond scheme and hence the absence of any cooperative motion of protons in this crystal.

Raman spectra of single crystals of adipic and sebacic acids and a study of the hydrogen bond vibrations in carboxylic acid dimers

Raman spectra of single crystals of adipic and sebacic acids have been photographed for the first time using λ 2537 excitation. The spectra have been divided into four regions: (a) internal frequencies; (b) summations and overtones; (c) external vibrations; and (d) low-frequency hydrogen bond oscillations. Tentative correlations have been given for all the internal frequencies and summations and overtones. A series of diffuse weak bands observed in the spectra of both these acids in the not, vert, similar2400–2800 cm−1 have been explained as a superposition of O---H frequencies lowered due to hydrogen bond formation over the summations and overtones of fundamentals mainly in the not, vert, similar1000–1500 cm−1 region. Rotatory type of external oscillations of the two formula units of these molecules in their unit cells have been identified at 76, 99, 118 and 165 cm−1 in adipic acid and 66, 95, 117 and 177 cm−1 in the spectrum of sebacic acid. A brief discussion of the low frequency hydrogen bond vibrations in these acids has been made. Making use of the Lippincott—Schroeder potential and assuming a highly anharmonic potential curve for the hydrogen bond, the vibrational frequencies of the bond have been theoretically evaluated. There is very good agreement between these and the experimental values. The results for adipic acid in cm−1 are: 304 (0 → 1), 270 (1 → 2), 241 (2 → 3), 222 (3 → 4) 201 (4 → 5), 183 (5 → 6). In the case of sebacic acid some of the intermediate and higher transitions are absent in the spectrum recorded by the author. From the above data for adipic acid the dissociation energy of the hydrogen bond was evaluated as 5·9 kcal/mole in fair agreement with the values derived from conventional methods.

Potential functions for hydrogen bond interactions - III. Empirical Potential Function for the Peptide N-H...O=C Hydrogen Bond

A careful comparison of the distribution in the (R, θ)-plane of all NH ... O hydrogen bonds with that for bonds between neutral NH and neutral C=O groups indicated that the latter has a larger mean R and a wider range of θ and that the distribution was also broader than for the average case. Therefore, the potential function developed earlier for an average NH ... O hydrogen bond was modified to suit the peptide case. A three-parameter expression of the form {Mathematical expression}, with △ = R - Rmin, was found to be satisfactory. By comparing the theoretically expected distribution in R and θ with observed data (although limited), the best values were found to be p1 = 25, p3 = - 2 and q1 = 1 × 10-3, with Rmin = 2·95 Å and Vmin = - 4·5 kcal/mole. The procedure for obtaining a smooth transition from Vhb to the non-bonded potential Vnb for large R and θ is described, along with a flow chart useful for programming the formulae. Calculated values of ΔH, the enthalpy of formation of the hydrogen bond, using this function are in reasonable agreement with observation. When the atoms involved in the hydrogen bond occur in a five-membered ring as in the sequence[Figure not available: see fulltext.] a different formula for the potential function is needed, which is of the form Vhb = Vmin +p1△2 +q1x2 where x = θ - 50° for θ ≥ 50°, with p1 = 15, q1 = 0·002, Rmin = 2· Å and Vmin = - 2·5 kcal/mole. © 1971 Indian Academy of Sciences.

Charge Density Analysis of Ferulic Acid: Robustness of a Trifurcated C-H center dot center dot center dot O Hydrogen Bond

The toplogical features of a sporadic trifurcated C-H center dot center dot center dot O interaction region, where an oxygen atom acts as an acceptor of three weak hydrogen bonds, has been investigated by experimental and theoretical charge density analysis of ferulic acid. The interaction energy of the asymmetric molecular dimer formed by the trifurcated C-H center dot center dot center dot O motif, based on the multipolar model, is shown to be greater than the corresponding asymmetric O-H center dot center dot center dot O dimer in this crystal structure. Further, the hydrogen bond energies associated with these interaction motifs have been estimated from the local kinetic and potential energy densities at the bond critical points. The trends suggest that the interaction energy of the trifurcated C-H center dot center dot center dot O region is comparable to that of a single O-H center dot center dot center dot O hydrogen bond.

Estimation of the 2.0(5) helix type i -> i hydrogen bond energy at Aib*-Oxa motif: an isodesmic approach

Bending at the valence angle N-C-alpha-C' (tau) is a known control feature for attenuating the stability of the rare intramolecular hydrogen bonded pseudo five-membered ring C-5 structures, the so called 2.0(5) helices, at Aib. The competitive 3(10)-helical structures still predominate over the C5 structures at Aib for most values of tau. However at Aib*, a mimic of Aib where the carbonyl 0 of Aib is replaced with an imidate N (in 5,6-dihydro-4H-1,3-oxazine = Oxa), in the peptidomimic Piv-Pro-Aib*-Oxa (1), the C(5)i structure is persistent in both crystals and in solution. Here we show that the i -> i hydrogen bond energy is a more determinant control for the relative stability of the C5 structure and estimate its value to be 18.5 +/- 0.7 kJ/mol at Aib* in 1, through the computational isodesmic reaction approach, using two independent sets of theoretical isodesmic reactions. (C) 2014 Elsevier Ltd. All rights reserved.

Tautomeric Preference and Conformation Locking in Fenobam, Thiofenobam, and Their Analogues: The Decisive Role of Hydrogen Bond Hierarchy

The crystal and molecular structures of the potential antidepressant drug fenobam and its derivatives are examined in terms of the preferred form among the two possible tautomeric structures. In this study, chemical derivatization has been utilized as a means to ``experimentally simulate'' the tautomeric preference and conformational variability in fenobam. Eight new derivatives of fenobam have been synthesized, and structural features have been characterized by single-crystal X-ray diffraction and NMR spectroscopy. The specific tautomeric preference found in all of these compounds and their known crystal forms have been construed in terms of the stabilizing intramolecular N-H center dot center dot center dot O and N-H center dot center dot center dot S hydrogen bonding. The hierarchy of intramolecular hydrogen bonds evidenced as the preference of the C-H center dot center dot center dot O hydrogen bond over C-H center dot center dot center dot N and that of the C-H center dot center dot center dot N hydrogen bond over C-H center dot center dot center dot S explains the two distinct conformations adopted by fenobam and thiofenobam derivatives. The relative energy values of different molecular conformations have been calculated and compared.

Isotopic fractionation in proteins as a measure of hydrogen bond length

If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent, it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor (I) is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds, we calculate (I) as a function of the proton donor-acceptor distance R. For numerical results, we use a parameterization of the model for symmetric 0-H. ``.0 bonds R. H. McKenzie, Chem. Phys. Lett. 535, 196 (2012)]. We consider the relative contributions of the 0-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunneling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total (I) as a function of R with NMR experimental results for enzymes, and in particular with an earlier model parametrization (D(R), used previously to determine bond lengths. (C) 2015 AIP Publishing LLC.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Hydrogen-bond dynamics near a micellar surface: Origin of the universal slow relaxation at complex aqueous interfaces

The dynamics of hydrogen bonds among water molecules themselves and with the polar head groups (PHG) at a micellar surface have been investigated by long molecular dynamics simulations. The lifetime of the hydrogen bond between a PHG and a water molecule is found to be much longer than that between any two water molecules, and is likely to be a general feature of hydrophilic surfaces of organized assemblies. Analyses of individual water trajectories suggest that water molecules can remain bound to the micellar surface for more than 100 ps. The activation energy for such a transition from the bound to a free state for the water molecules is estimated to be about 3.5 kcal/mol.

Trifurcated (four-center) hydrogen bond in solid state crystal structure of 5'-amino-5'-deoxyadenosine p-toluenesulfonate

The crystal structure of 5'-amino-5'-deoxyadenosine (5'-Am.dA) p-toluenesulfonate has been determined by X-ray crystallographic methods. It belongs to the orthorhombic space group P2(1)2(1)2(1) with a = 7.754(3)Angstrom, b = 8.065(1)Angstrom and c = 32.481(2)Angstrom. This nucleoside side shows a syn conformation about the glycosyl bond and C2'-endo-C3'-exo puckering for the ribose sugar. The orientation of N5' atom is gauche-trans about the exocyclic C4'-C5' bond. The amino nitrogen N5' forms a trifurcated hydrogen bond with N3, O9T and O4' atoms. Adenine bases form A.A.A triplets through hydrogen bonding between N6, N7 and N1 atoms of symmetry related nucleoside molecules.

Cytotoxicity of half sandwich ruthenium(II) complexes with strong hydrogen bond acceptor ligands and their mechanism of action

Neutral and cationic organometallic ruthenium(II) piano stool complexes of the type [(eta(6)-cymene)R-uCl(X)(Y)] (complexes R1-R8) has been synthesized and characterized. In cationic complexes, X, Y is either a eta(2) phosphorus ligand such as 1,1-bis(diphenylphosphino)methane (DPPM) and 1,2-bis(diphenylphosphino)ethane (DPPE) or partially oxidized ligands such as 1,2-bis(diphenylphosphino)methane monooxide (DPPMO) and 1,2-bis(diphenylphosphino)ethane monooxide (DPPEO) which are strong hydrogen bond acceptors. In neutral complexes. X is chloride and Y is a monodentate phosphorous donor. Complexes with DPPM and DPPMO ligands ([(eta(6)-cymene)Ru(eta(2)-DPPM)Cl]PF6 (R2), [(eta(6)-cymene)Ru(eta(2)-DPPMO)Cl]PF6 (R3), [(eta(6)-cymene)Ru(eta(1)-DPPM)Cl-2] (R5) and [(eta(6)-cymene)Ru(eta(1)-DPPMO)Cl-2] (R6) show good cytotoxicity. Growth inhibition study of several human cancer cell lines by these complexes has been carried out. Mechanistic studies for R5 and R6 show that inhibition of cancer cell growth involves both cell cycle arrest and apoptosis induction. Using an apoptosis PCR array, we identified the sets of antiapoptotic genes that were down regulated and pro-apoptotic genes that were up regulated. These complexes were also found to be potent metastasis inhibitors as they prevented cell invasion through matrigel. The complexes were shown to bind DNA in a non intercalative fashion and cause unwinding of plasmid DNA in cell-free medium by competitive ethidium bromide binding, viscosity measurements, thermal denaturation and gel mobility shift assays.

Peptide mimics for structural features in proteins. Crystal structures of three heptapeptide helices with a C-terminal 6-->1 hydrogen bond.

The crystal structure determination of three heptapeptides containing alpha-aminoisobutyryl (Aib) residues as a means of helix stabilization provides a high-resolution characterization of 6-->1 hydrogen-bonded conformations, reminiscent of helix-terminating structural features in proteins. The crystal parameters for the three peptides, Boc-Val-Aib-X-Aib-Ala-Aib-Y-OMe, where X and Y are Phe, Leu (I), Leu, Phe (II) and Leu, Leu (III) are: (I) space group P1, Z = 1, a = 9.903 A, b = 10.709 A, c = 11.969 A, alpha = 102.94 degrees, beta = 103.41 degrees, gamma = 92.72 degrees, R = 4.55%; (II) space group P21, Z = 2, a = 10.052 A, b = 17.653 A, c = 13.510 A, beta = 108.45 degrees, R = 4.49%; (III) space group P1, Z = 2 (two independent molecules IIIa and IIIb in the asymmetric unit), a = 10.833 A, b = 13.850 A, c = 16.928 A, alpha = 99.77 degrees, beta = 105.90 degrees, gamma = 90.64 degrees, R = 8.54%. In all cases the helices form 3(10)/alpha-helical (or 3(10)helical) structures, with helical columns formed by head-to-tail hydrogen bonding. The helices assemble in an all-parallel motif in crystals I and III and in an antiparallel motif in II. In the four crystallographically characterized molecules, I, II, IIIa and IIIb, Aib(6) adopts a left-handed helical (hL) conformation with positive phi, psi values, resulting in 6-->1 hydrogen-bond formation between Aib(2) CO and Leu(7)/Phe(7) NH groups. In addition a 4-->1 hydrogen bond is seen between Aib(3) CO and Aib(6) NH groups. This pattern of hydrogen bonding is often observed at the C-terminus of helices proteins, with the terminal pi-type turn being formed by four residues adopting the hRhRhRhL conformation.

Influence of the hydrogen bond on the N-H stretching frequencies in amino-acids

The spectra of glycine, its addition compounds and other amino-acids exhibit Raman lines in the region from 3250 cm.−1 to 2500 cm.−1 It has been shown that these lines cannot be assigned to N-H...O stretching vibrations, where the N atom has the covalency of three, but to N+-H...O stretching vibration where the N atom has the covalency of four. Using the data obtained with triglycine sulphate which has the largest number of N+-H...O bonds and whose H bond lengths are known, the correlation curve giving the relation between the N+-H...O stretching frequencies and the corresponding H bond lengths has been drawn. Using this correlation curve, the N+-H...O stretching frequencies appearing inα-glycine,γ-glycine, diglycine hydrochloride, diglycine hydrobromide,l-asparagine monohydrate anddl-alanine have been satisfactorily accounted for on the basis of the known hydrogen bond lengths in these substances.