Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Crystalline ethane-1,2-diol does not have intra-molecular hydrogen bonding: Experimental and theoretical charge density studies

The X-ray structure and electron density distribution of ethane-1,2-diol (ethylene glycol), obtained at a resolution extending to 1.00 Å−1 in sin θ/λ (data completion = 100% at 100 K) by in situ cryocrystallization technique is reported. The diol is in the gauche (g′Gt) conformation with the crystal structure stabilised by a network of inter-molecular hydrogen bonds. In addition to the well-recognized O–H···O hydrogen bonds there is topological evidence for C–H···O inter-molecular interactions. There is no experimental electron density based topological evidence for the occurrence of an intra-molecular hydrogen bond. The O···H spacing is not, vert, similar0.45 Å greater than in the gas-phase with an O–H···O angle close to 90°, calling into question the general assumption that the gauche conformation of ethane-1,2-diol is stabilised by the intra-molecular oxygen–hydrogen interaction.

Molecular and crystal structure of 4-ethynylcyanobenzene. An example of Csingle bondHcdots, three dots, centeredN hydrogen bonding

The molecular and crystal structures of 4-ethynylcyanobenzene arereported. The packing of molecules in the crystal is found to be homologous with the crystal structures of HCN, cyanoacetylene and 4-cyano-4'-ethynylbiphenyl. Alternatively, these four crystals could be said to constitute a structural homologous series. The influence of C-H center dot center dot center dot N hydrogen bonding in directing a linear supramolecular arrangement of molecules with ethynyl and cyano groups at opposite ends, is illustrated. (C) 2010 Elsevier B.V. All rights reserved.

Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

Molecular self-assembly of cadmium-triazolate complexes via hydrogen bonding: Synthesis, structures and photoluminescent properties

Two new cadmium coordination polymers namely Cd(HAmTrz-COO)(4)(NH4+)(2)] 1; and Cd(HAmTrz)(2)I-2](n) 2; (HAmTrz-COOH = 3-amino-1,2,4-triazole-5-carboxylic acid), have been prepared based on HAmTrz-COOH as ligand. The crystal structures of 1 and 2 have been determined by single-crystal X-ray diffraction technique. In coordination-complex 1 four triazole ligands coordinate via N1 nitrogen leading to a tetrahedral geometry around cadmium ion, while in 2 the ligand prefers to coordinate to the metal in a bidentate bridging mode. The structures of both the coordination polymers can be envisaged as 3D hydrogen bonded networks. Thermogravimetric analysis shows that 2 is more stable than 1 owing to different coordination numbers of cadmium atoms. Photoluminescence properties of both the compounds have been investigated in the solid state. (C) 2011 Elsevier B.V. All rights reserved.

Intra-molecular hydrogen bonding with organic fluorine in the solution state: Deriving evidence by a two dimensional NMR experiment

The direct evidence for the existence of intra-molecular C-F center dot center dot center dot H-N hydrogen bond in organofluorine molecules, in the liquid state, is derived using NMR spectroscopy by the detection of long range interactions among fluorine, nitrogen and hydrogen atoms. The present study reports the determination of the relative signs and magnitudes of through space and through bond couplings to draw unambiguous evidence on the existence of weak molecular interactions involving organic fluorine. It is a simple, easy to implement, N-15 natural abundant two dimensional heteronuclear N-15-H-1 double quantum-single quantum correlation experiment. The existence of intra-molecular hydrogen bond is conclusively established in the investigated molecules. (C) 2011 Elsevier B.V. All rights reserved.

Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Herringbone to cofacial solid state packing via H-bonding in diketopyrrolopyrrole (DPP) based molecular crystals: influence on charge transport

Themono-alkylation of DPP derivatives leads to cofacial pi-pi stacking via H-bonding unlike their di-alkylated counterparts, which exhibit a classical herringbone packing pattern. Single crystal organic field-effect transistor (OFET) measurements reveal a significant enhancement of charge carrier mobility for mono-hexyl DPP derivatives.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.

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.

Stability of Dimeric Interface in Banana Lectin: Insight from Molecular Dynamics Simulations

Banana lectin (Banlec) is a homodimeric non-glycosylated protein. It exhibits the b-prism I structure. High-temperature molecular dynamics simulations have been utilized to monitor and understand early stages of thermally induced unfolding of Banlec. The present study elucidates the behavior of the dimeric protein at four different temperatures and compares the structural and conformational changes to that of the minimized crystal structure. The process of unfolding was monitored by following the radius of gyration, the rms deviation of each residue, change in relative solvent accessibility and the pattern of inter- and intra-subunit interactions. The overall study demonstrates that the Banlec dimer is a highly stable structure, and the stability is mostly contributed by interfacial interactions. It maintains its overall conformation during high-temperature (400–500 K) simulations, with only the unstructured loop regions acquiring greater momentum under such condition. Nevertheless, at still higher temperatures (600 K) the tertiary structure is gradually lost which later extends to loss of secondary structural elements. The pattern of hydrogen bonding within the subunit and at the interface across different stages has been analyzed and has provided rationale for its intrinsic high stability.

Zeeman effect of chlorine nuclear quadrupole resonance and hydrogen bonding in 2,6-, 2,5- and 3,5-dichlorophenols

The Zeeman effect of chlorine nuclear quadrupole resonance in polycrystalline samples of 2,6-, 2,5 and 3,5-dichlorophenol has been investigated at room temperature in order to study the effect of hydrogen bonding on the electric field gradient asymmetry parameter n. While the two n.q.r. lines in 3,5-dichlorophenol gave an asymmetry parameter of 10%, those in 2,6- and 2,5-dichlorophenol gave different values of n for the two chlorines. The chlorine atom which is ortho to the OH group and involved in hydrogen bonding (i.e., corresponding to the low frequency line) gave an asymmetry parameter of 0.21 in 2,6-dichlorophenol and 0.17 in 2,5-dichlorophenol while the other chlorine (i.e., corresponding to the high frequency line) gave a lower value of 0.12 in 2,6-dichlorophenol and 0.11 in 2,5-dichlorophenol. These values of n are discussed in terms of hydrogen bonding and bond parameters.

Molecular conformation of N,N -diarylthioureas: An assessment by 1H NMR and infrared spectroscopy

Several N,N -dipyridyl- and N-phenyl-N -pyridyl-thioureas were examined in different solvents at various temperatures by 1H NMR in order to study their conformational properties. The influence of concentration and the methyl substituent in the pyridine ring on the chemical shifts of the NH and pyridine groups was investigated. The observed chemical shifts are analysed in terms of the conformational properties of the molecules. Free energy barriers to the internal rotation about the C N bonds have been determined. Infrared spectra have been measured to supplement the NMR studies. Intramolecular hydrogen bonding played a major role in the preferred conformation of pyridylthioureas. The data further revealed an interesting dynamic exchange phenomenon occurring in symmetric N,N -dipyridylthioureas between two intramolecularly hydrogen bonded conformers.

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W I and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the Ob atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.