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.

Hydrogen-Bond-Driven Controlled Molecular Marriage in Covalent Cages

A supramolecular approach that uses hydrogen-bonding interaction as a driving force to accomplish exceptional self-sorting in the formation of imine-based covalent organic cages is discussed. Utilizing the dynamic covalent chemistry approach from three geometrically similar dialdehydes (A, B, and D) and the flexible triamine tris(2-aminoethyl)amine (X), three new 3+2] self-assembled nanoscopic organic cages have been synthesized and fully characterized by various techniques. When a complex mixture of the dialdehydes and triamine X was subjected to reaction, it was found that only dialdehyde B (which has OH groups for H-bonding) reacted to form the corresponding cage B3X2 selectively. Surprisingly, the same reaction in the absence of aldehyde B yielded a mixture of products. Theoretical and experimental investigations are in complete agreement that the presence of the hydroxyl moiety adjacent to the aldehyde functionality in B is responsible for the selective formation of cage B3X2 from a complex reaction mixture. This spectacular selection was further analyzed by transforming a nonpreferred (non-hydroxy) cage into a preferred (hydroxy) cage B3X2 by treating the former with aldehyde B. The role of the H-bond in partner selection in a mixture of two dialdehydes and two amines has also been established. Moreover, an example of unconventional imine bond metathesis in organic cage-to-cage transformation is reported.

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.

Dynamics of the chemical bond: inter- and intra-molecular hydrogen bond

In this discussion, we show that a static definition of a `bond' is not viable by looking at a few examples for both inter-and intra-molecular hydrogen bonding. This follows from our earlier work (Goswami and Arunan, Phys. Chem. Chem. Phys. 2009, 11, 8974) which showed a practical way to differentiate `hydrogen bonding' from `van der Waals interaction'. We report results from ab initio and atoms in molecules theoretical calculations for a series of Rg center dot center dot center dot HX complexes (Rg = He/Ne/Ar and X = F/Cl/Br) and ethane-1,2-diol. Results for the Rg center dot center dot center dot HX/DX complexes show that Rg center dot center dot center dot DX could have a `deuterium bond' even when Rg center dot center dot center dot HX is not `hydrogen bonded', according to the practical criterion given by Goswami and Arunan. Results for ethane-1,2-diol show that an `intra-molecular hydrogen bond' can appear during a normal mode vibration which is dominated by the O center dot center dot center dot O stretching, though a `bond' is not found in the equilibrium structure. This dynamical `bond' formation may nevertheless be important in ensuring the continuity of electron density across a molecule. In the former case, a vibration `breaks' an existing bond and in the later case, a vibration leads to `bond' formation. In both cases, the molecule/complex stays bound irrespective of what happens to this `hydrogen bond'. Both these cases push the borders on the recent IUPAC recommendation on hydrogen bonding (Arunan et al. Pure. Appl. Chem. 2011, 83 1637) and justify the inclusive nature of the definition.

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.

Intramolecular hydrogen bond: Can it be part of the basis set of valence internal coordinates in normal mode analysis?

It has been shown earlier1] that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent' basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.

Hydrogen Bonds Computing Server (HBCS): an online web server to compute hydrogen-bond interactions and their precision

Hydrogen bonds in biological macromolecules play significant structural and functional roles. They are the key contributors to most of the interactions without which no living system exists. In view of this, a web-based computing server, the Hydrogen Bonds Computing Server (HBCS), has been developed to compute hydrogen-bond interactions and their standard deviations for any given macromolecular structure. The computing server is connected to a locally maintained Protein Data Bank (PDB) archive. Thus, the user can calculate the above parameters for any deposited structure, and options have also been provided for the user to upload a structure in PDB format from the client machine. In addition, the server has been interfaced with the molecular viewers Jmol and JSmol to visualize the hydrogen-bond interactions. The proposed server is freely available and accessible via the World Wide Web at http://bioserver1.physics.iisc.ernet.in/hbcs/.

Fe as Hydrogen/Halogen Bond Acceptor in Square Pyramidal Fe(CO)(5)

Hydrogen bonded complexes formed between the square pyramidal Fe(CO)(5) with HX (X = F, Cl, Br), showing X-H center dot center dot center dot Fe interactions, have been investigated theoretically using density functional theory (DFT) including dispersion correction. Geometry, interaction energy, and large red shift of about 400 cm(-1) in the FIX stretching frequency confirm X-H center dot center dot center dot Fe hydrogen bond formation. In the (CO)(5)Fe center dot center dot center dot HBr complex, following the significant red shift, the HBr stretching mode is coupled with the carbonyl stretching modes. This clearly affects the correlation between frequency shift and binding energy, which is a hallmark of hydrogen bonds. Atoms in Molecule (AIM) theoretical analyses show the presence of a bond critical point between the iron and the hydrogen of FIX and significant mutual penetration. These X-H center dot center dot center dot Fe hydrogen bonds follow most but not all of the eight criteria proposed by Koch and Popelier (J. Phys. Chem. 1995, 99, 9747) based on their investigations on C-H center dot center dot center dot O hydrogen bonds. Natural bond orbital (NBO) analysis indicates charge transfer from the organometallic system to the hydrogen bond donor. However, there is no correlation between the extent of charge transfer and interaction,energy, contrary to what is proposed in the recent IUPAC recommendation (Pure Appl.. Chem. 2011, 83, 1637). The ``hydrogen bond radius'' for iron has been determined to be 1.60 +/- 0.02 angstrom, and not surprisingly it is between the covalent (127 angstrom) and van der Waals (2.0) radii of Fe. DFT and AIM theoretical studies reveal that Fe in square pyramidal Fe(CO)(5) can also form halogen bond with CIF and ClH as ``halogen bond donor''. Both these complexes show mutual penetration as well, though the Fe center dot center dot center dot Cl distance is closer to the sum of van der Waals radii of Fe and Cl in (CO)5Fe center dot center dot center dot ClH, and it is about 1 angstrom less in (CO)(5)Fe center dot center dot center dot ClF.

Hydrogen bonding in peptide helices - Analysis of two independent helices in the crystal structure of a peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe

The crystal structure determination of the heptapeptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe reveals two peptide helices in the asymmetric unit, Crystal parameters are: space group P2(1), a = 10.356(2) Angstrom, b = 19.488(5) Angstrom, c = 23.756(6) Angstrom, beta = 102.25(2)degrees), V = 4685.4 Angstrom(3), Z = 4 and R = 5.7% for 7615 reflections [I>3 sigma(I)]. Both molecules adopt largely alpha-helical conformations with variations at the C-terminus, Helix type Is determined by analysing both 4-->1 and 5-->1 hydrogen-bond interactions and comparison with the results of analysis of protein structures. The presence of two 4-->1 hydrogen-bond interactions, besides four 5-->1 interact ions in both the conformations provides an opportunity to characterize bifurcated hydrogen bonds at high resolution, Comparison of the two helical conformations with related peptide structures suggests that distortions at the C-terminus are more facile than at the N-terminus.

Stereochemical studies on cyclic peptides. XIII. Energy minimization studies on cyclic hexapeptides having hydrogen bonds

The basic cyclic hexapeptide conformations which accommodate hydrogen bonded β and γ turns in the backbone have been worked out using stereochemical criteria and energy minimization procedures. It was found that cyclic hexapeptides can be made up of all possible combinations of 4 ± 1 hydrogen bonded types I, I', II and II' β turns, giving rise to symmetric conformations having twofold and inversion symmetries as well as nonsymmetric structures. Conformations having exclusive features of 3 ± 1 hydrogen bonded γ turns were found to be possible in threefold and S6 symmetric cyclic hexapeptides. The results show that the cyclic hexapeptides formed by the linking of two β turn tripeptide fragments differ mainly in (a) the hydrogen bonding scheme present in the β turn tripeptides and (b) the conformation at the α-carbon atoms where the two tripeptide fragments link. The different hydrogen bonding schemes found in the component β turns are: 1) a β turn with only a 4 ± 1 hydrogen bond, 2) a type I or I' β turn with 4 ± 1 and 3 ± 1 hydrogen bonds occurring in a bifurcated form and 3) a type II or II' β turn having both the 4 ± 1 and the 3 ± 1 hydrogen bonds with the same acceptor oxygen atom. The conformation at the linking α-carbon atoms was found to lie either in the extended region or in the 3 ± 1 hydrogen bonded γ turn or inverse γ turn regions. Further, the threefold and the S6 symmetric conformations have three γ turns interleaved by three extended regions or three inverse γ turns, respectively. The feasibility of accommodating alanyl residues of both isomeric forms in the CHP minima has been explored. Finally, the available experimental data are reviewed in the light of the present results.

Conformation of the LL and LD hairpin bends with internal hydrogen bonds in proteins and peptides

The conformation of three linked peptide units having an internal 4 → 1 type of hydrogen bond has been studied in detail, and the low energy conformations are listed. These conformations all lead to the reversal of the chain direction, and may therefore be called as “hairpin bends” or “U-bends”. Since this bend can occur at the end of two chains hydrogen-bonded in the antiparallel β-conformation, it is also known as the “β-bend”. Two types of conformation are possible when the residues at the second and third Cα atoms are both of type L (the LL bend), while only one type is possible for the LD and the DL bend. The LL bend can also accommodate the sequences LG, GL, GG (G = glycine), while the LD bend can accommodate the sequences LG, GD and GG. The conformations for the sequences DD and DL are exact inverses (or mirror images) of those for the sequences LL and LD, respectively, and have dihedral angles (phi2, ψ2), (phi3, ψ3) of the same magnitudes, but of opposite signs as those for the former types, which are listed, along with the characteristics (length, angle and energy) of the hydrogen bonds. A comparison of the theoretical predictions with experimental data (from X-ray diffraction and NMR studies) on proteins and peptides, show reasonably good agreement. However, a systematic trend is observable in the experimental data, slightly deviating from theory, which indicates that some deformations occur in the shapes of the peptide units forming the bend, differing from that of the standard planar peptide unit.

The 310 Helical Conformation of the Amino Terminal Decapeptide of Suzukacillin. 270 MHz 'H NMR Evidence for Eight Intramolecular Hydrogen Bonds

The amino terminal suzukacillin decapeptide fragment, Boc-Aib-Pro-Val-Aib-Val-Ala-Aib-Ala-Aib-Aitbh-eO Me, two pentapeptides Boc-AibPrc-Val-AibVal-OMe and Boc-Ala-AibAla-AibAibOMe, and the tripeptide Boc-Ala-AibAibOMe have been studied by 270-MHz 'H NMR spectroscopy. By use of solvent dependence of chemical shifts in a CDC13-(CD3),S0 system and temperature dependence of amide NH chemical shifts in (CD3),S0, the intramolecularly hydrogen bonded NH groups in these peptides have been identified. The tripeptide possesses one hydrogen bond, both pentapeptides show evidence for three intramolecular hydrogen bonds, and the decapeptide has eight NH groups participating in hydrogen bonding. An Ala( 1)-Aib(2) @ turn is proposed for the tripeptide. Both pentapeptides favor 310 helical conformations composed of three consecutive B turns. The decapeptide adopts a 310 helical conformation with some flexibility at the Va1(5)-Ala(6) segment. The proposed conformations are consistent with the known stereochemical preferences of Aib residues.

Studies On Hydrogen Bonds: Part IV. Proposed Working Criteria for Assessing Qualitative Strength of Hydrogen Bonds

The data obtained in the earlier parts of this series for the donor and acceptor end parameters of N-H. O and O-H. O hydrogen bonds have been utilised to obtain a qualitative working criterion to classify the hydrogen bonds into three categories: "very good" (VG), "moderately good" (MG) and weak (W). The general distribution curves for all the four parameters are found to be nearly of the Gaussian type. Assuming that the VG hydrogen bonds lie between 0 and ± la, MG hydrogen bonds between ± 1 and ± 2, W hydrogen bonds beyond ± 2 (where is the standard deviation), suitable cut-off limits for classifying the hydrogen bonds in the three categories have been derived. These limits are used to get VG and MG ranges for the four parameters 1 and θ (at the donor end) and ± and ± (at the acceptor end). The qualitative strength of a hydrogen bond is decided by the cumulative application of the criteria to all the four parameters. The criterion has been further applied to some practical examples in conformational studies such as α-helix and can be used for obtaining suitable location of hydrogen atoms to form good hydrogen bonds. An empirical approach to the energy of hydrogen bonds in the three categories has also been presented.

Conformations of Heterochiral and Homochiral Proline-Pseudoproline Segments in Peptides: Context Dependent Cis-Trans Peptide Bond Isomerization

The pseudoproline residue (Psi Pro, L-2,2-dimethyl-1,3-thiazolidine-4-carboxylic acid) has been introduced into heterochiral diproline segments that have been previously shown to facilitate the formation of beta-hairpins, containing central two and three residue turns. NMR studies of the octapeptide Boc-Leu-Phe-Val-(D)Pro-Psi Pro-Leu-Phe-Val-OMe (1), Boc-Leu-Val-Val-(D)Pro-Psi Pro-Leu-Val-Val-OMe (2), and the nonapeptide sequence Boc-Leu-Phe-Val-(D)Pro-Psi Pro-(D)Ala-Leu-Phe-Val-OMe (3) established well-registered beta-hairpin structures in chloroform solution, with the almost exclusive population of the trans conformation for the peptide bond preceding the Psi Pro residue. The beta-hairpin conformation of 1 is confirmed by single crystal X-ray diffraction. Truncation of the strand length in Boc-Val-(D)Pro-Psi Pro-Leu-OMe (4) results in air increase in the population of the cis conformer, with a cis/trans ratio of 3.65. Replacement of Psi Pro in 4 by (L)Pro in 5, results in almost exclusive population of the trans form, resulting in an incipient beta-hairpin conformation, stabilized by two intramolecular hydrogen bonds. Further truncation of the sequence gives an appreciable rise in the population of cis conformers in the tripeptide piv-(D)Pro-Psi Pro-Leu-OMe (6). In the homochiral segment Piv-Pro Psi Pro-Leu-OMe (7) only the cis form is observed with the NMR evidence strongly supporting a type VIa beta-turn conformation, stabilized by a 4 -> 1 hydrogen bond between the Piv (CO) and Leu (3) NH groups. The crystal structure of the analog peptide 7a (Piv-Pro-Psi(H,CH3)Pro-Leu-NHMe) confirms the cis peptide bond geometry for the Pro-Psi(H,CH3)pro peptide bond, resulting in a type VIa beta-turn conformation.

Structure of the hydrogen-bonded water molecule in crystals

A correlation of the structural data on IS hydrates obtained by x-ray diffraction, neutron diffraction, and proton magnetic resonance reveals that when a water molecule is hydrogen bonded into a crystal structure and the angle subtended at the donor water oxygen by the acceptor atoms deviates from the vapor H-O-H angle, bent hydrogen bonds are formed in preference to distortion of the H-O-H angle. Theoretical justification for this result is obtained from energy considerations by calculating the energy of formation of bent hydrogen bonds on the basis of the Lippincott-Schroeder potential function model for the hydrogen bond and the energy of deformation of the H-O-H angle from spectroscopic force constants.