Effect of solvent upon CH...O hydrogen bonds with implications for protein folding

Computational study of the dimerization of glyphosate: mechanism and effect of solvent

A computational study on the structure and stability of different series of glyphosate (Glyph) dimers comprising nonionized (N) and zwitterionic structures (Z) for neutral monomers, followed by an analysis of energetics of Glyph dimerization process have been performed by means of quantum chemical calculations in different media. Optimized geometries for energy minima, as well as relative potential and free energies of the possible various conformers of each series of Glyph dimer were computed as a function of the medium at B3LYP-D3/6-311++G(2d,2p) level. The solvation model based on density (SMD) is employed for all solution phase computations. Non-ionized dimers (DN), anion-cation (AC) and either zwitterion-zwitterion (DZP and DZC) or non-ionized-phosphonate zwitterion (NZP) ionized neutral forms of Glyph dimer are predicted to exist in the gas phase and in solution in large contrast to Glyph monomers. The DZC dimer form exhibiting a centrosymmetric arrangement of two carboxylate zwitterion units was found to be the most stable dimer structure in all media. In aqueous solution, the DZP and AC dimer type structures are significantly stabilized by hydration. The tautomerisms between DZC, DZP and AC dimer type structures have been investigated in the gas phase and in solution. The DZC type structures are more prone to experience proton transfer in water than in the gas phase and in cyclohexane. The mechanism for the tautomerization process in neutral ionized Glyph dimers proceeds via two direct proton transfer paths: DZP ⇋ AC ⇋ DZC. Results show that solvents play a key role in modulating the energetics of the dimerization process of Glyph. Solvation in cyclohexane, favors the dimerization process however, hydration opposed it. In aqueous solution, the mechanism of the dimerization of Glyph from its phosphonate zwitterionic monomer form (ZP) could be described by a set of equilibria including direct proton transfer paths as follows: 2ZP ⇋ DZP ⇋ AC ⇋ DZC. According to our results, in aqueous solution, DZC Glyph dimers and their corresponding DZP and AC tautomers should be present in higher concentration than phosphonate zwitterionic Glyph monomers for high Glyph concentration, a fact that seems controversial in the literature.

Enhancement of Halogen Bond Strength by Intramolecular H-Bonds

Quantum calculations study the potential of an intramolecular H-bond between the halogen atom (X) of a halobenzene and a substituent placed ortho to it, to amplify the ability of X to engage in a halogen bond (XB) with a Lewis base. H-bonding substituents NH₂, CH₂CH₂OH, CH₂OH, OH, and COOH were added to halobenzenes (X = Cl, Br, I). The amino group had little effect, but those containing OH increased the CX···N XB energy to a NH₃ nucleophile by about 0.5 kcal/mol; the increment associated with COOH is larger, nearly 2 kcal/mol. These energy increments were approximately doubled if two such H-bonding substituents are present. Combining a pair of ortho COOH groups with an electron-withdrawing NO₂ group in the para position has a particularly large effect, raising the XB energy by about 4 kcal/mol, which can amount to as much as a 4-fold magnification.

Transition State Stabilization by SCF2 -H⋅⋅⋅O Bifurcated Hydrogen Bond

The difluoromethylthio group (SCF₂ H), which is generally considered a highly lipophilic weak hydrogen bonding donor, has attracted special interest from the pharmaceutical and agrochemical industry. Remarkably, there have been relatively few literature investigations of SCF₂ H hydrogen bonding interactions. Here, we report the determination of the hydrogen bond acidity parameter A of the SCF₂ H in the most popularly used electrophilic difluoromethylthiolating reagent. We present kinetic and computational evidence of the RSCF₂ -H⋅⋅⋅O2 bifurcated hydrogen bond for stabilizing the SCF₂ H-transferring transition state, which could cause a reversal of apparent electrophilic reactivity of difluoromethylthiolating and trifluoromethylthiolating reagents. Solvent effects on the RSCF₂ -H⋅⋅⋅O2 bifurcated hydrogen bonds will also be discussed.

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁ YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.

The unusual stability of H-bonded complexes in solvent caused by greater solvation energy of complex compared to those of isolated fragments

Here, the effect of solvent on the stability of non-covalent complexes, was studied. These complexes were from previously published S22, S66, and X40 datasets, which include hydrogen-, halogen- and dispersion-bonded complexes. It was shown that the charge transfer in the complex determines whether the complex is stabilized or destabilized in solvent.

Proline C-H Bonds as Loci for Proline Assembly via C-H/O Interactions

Proline residues within proteins lack a traditional hydrogen bond donor. However, the hydrogens of the proline ring are all sterically accessible, with polarized C-H bonds at Hα and Hδ that exhibit greater partial positive character and can be utilized as alternative sites for molecular recognition. C-H/O interactions, between proline C-H bonds and oxygen lone pairs, have been previously identified as modes of recognition within protein structures and for higher-order assembly of protein structures. In order to better understand intermolecular recognition of proline residues, a series of proline derivatives was synthesized, including 4R-hydroxyproline nitrobenzoate methyl ester, acylated on the proline nitrogen with bromoacetyl and glycolyl groups, and Boc-4S-(4-iodophenyl)hydroxyproline methyl amide. All three derivatives exhibited multiple close intermolecular C-H/O interactions in the crystallographic state, with H⋅⋅⋅O distances as close as 2.3 Å. These observed distances are well below the 2.72 Å sum of the van der Waals radii of H and O, and suggest that these interactions are particularly favorable. In order to generalize these results, we further analyzed the role of C-H/O interactions in all previously crystallized derivatives of these amino acids, and found that all 26 structures exhibited close intermolecular C-H/O interactions. Finally, we analyzed all proline residues in the Cambridge Structural Database of small-molecule crystal structures. We found that the majority of these structures exhibited intermolecular C-H/O interactions at proline C-H bonds, suggesting that C-H/O interactions are an inherent and important mode for recognition of and higher-order assembly at proline residues. Due to steric accessibility and multiple polarized C-H bonds, proline residues are uniquely positioned as sites for binding and recognition via C-H/O interactions.

Molecular duplexes featuring NH···N, CH···O and CH···π interactions in solid-state self-assembly of triazine-based compounds

Synthetic supramolecular structures constructed through the cooperative action of numerous non-covalent forces are highly desirable as models to unravel and understand the complexity of systems created in nature via self-assembly. Taking advantage of the low cost of 2,4,6-trichloro-1,3,5-triazine (cyanuric chloride) and the sequential nucleophilic substitution reactions with almost all types of nucleophiles, a series of six structurally related novel s-triazine derivatives 1-6 were synthesized and structurally characterized based on their physical, spectral and crystallographic data. The solid-state structures of all the six compounds showed intriguing and unique molecular duplexes featuring NH···N, CH···O and CH···π interactions. Careful analysis of different geometric parameters of the involved H-bonds indicates that they are linear, significant and are therefore responsible for guiding the three-dimensional structure of these compounds in the solid state. The prevalence of sextuple hydrogen bond array-driven molecular duplexes and the possibility of structural modifications on the s-triazine ring render these novel triazine derivatives 1-6 attractive as a platform to create heteroduplex constructs and their subsequent utility in the field of supramolecular chemistry and crystal engineering.

Solvation stabilizes intercarbonyl n→π* interactions and polyproline II helix

n→π* interactions between consecutive carbonyls stabilize the α-helix and polyproline II helix (PPII) conformations in proteins. n→π* interactions have been suggested to provide significant conformational biases to the disordered states of proteins. To understand the roles of solvation on the strength of n→π* interactions, computational investigations were conducted on a model n→π* interaction, the twisted-parallel-offset formaldehyde dimer, as a function of explicit solvation of the donor and acceptor carbonyls, using water and HF. In addition, the effects of urea, thiourea, guanidinium, and monovalent cations on n→π* interaction strength were examined. Solvation of the acceptor carbonyl significantly strengthens the n→π* interaction, while solvation of the donor carbonyl only modestly weakens the n→π* interaction. The n→π* interaction strength was maximized with two solvent molecules on the acceptor carbonyl. Urea stabilized the n→π* interaction via simultaneous engagement of both oxygen lone pairs on the acceptor carbonyl. Solvent effects were further investigated in the model peptides Ac-Pro-NMe₂, Ac-Ala-NMe₂, and Ac-Pro₂-NMe₂. Solvent effects in peptides were similar to those in the formaldehyde dimer, with solvation of the acceptor carbonyl increasing n→π* interaction strength and resulting in more compact conformations, in both the proline endo and exo ring puckers, as well as a reduction in the energy difference between these ring puckers. Carbonyl solvation leads to an energetic preference for PPII over both the α-helix and β/extended conformations, consistent with experimental data that protic solvents and protein denaturants both promote PPII. Solvation of the acceptor carbonyl weakens the intraresidue C5 hydrogen bond that stabilizes the β conformation.

Biotin's Lessons in Drug Design

At the heart of drug design is the discovery of molecules that bind with high affinity to their drug targets. Biotin forms the strongest known noncovalent ligand-protein interactions with avidin and streptavidin, achieving femtomolar and picomolar affinities, respectively. This is made even more exceptional because biotin achieves this with a meagre molecular weight of 240 Da. Surprisingly, the approaches by which biotin achieves this are not in the standard repertoire of current medicinal chemistry practice. Biotin's biggest lesson is the importance of nonclassical H-bonds in protein-ligand complexes. Most of biotin's affinity stems from its flexible valeric acid side chain that forms CH-π, CH-O, and classical H-bonds with the lipophilic region of the binding pocket. Biotin also utilizes an oxyanion hole, a sulfur-centered H-bond, and water solvation in the bound state to achieve its potency. The facets and advantages of biotin's approach to binding should be more widely adopted in drug design.

Exploring the Role of Consecutive Addition of Nitrogen Atoms on Stability and Reactivity of Hydrogen-Bonded Azine-Water Complexes

The second-order Møller-Plesset perturbation theory (MP2) and density functional theory with dispersion function calculations have been applied to investigate the hydrogen-bonding interaction between azines and water. The study suggests that the ability of nitrogen present in azine to act as a hydrogen-bond acceptor decreases in the order of pyridine (PY) > diazine (DZ) > triazine (TZ) > tetrazine (TTZ) > pentazine (PZ) > hexazine (HZ). Natural bond orbital (NBO) analysis, atoms in molecules, symmetry-adapted perturbation theory (SAPT), and molecular electrostatic potential studies reflect the factors important for hydrogen-bond strength as well as for the structural, electronic, and vibrational changes occurring during complexation. NBO analysis reflects that upon gradual addition of nitrogen atoms, hyperconjugation leads to an increase in the population of antibonding O-H bond, thus causing elongation and weakening of O-H bond in complexes incorporating N···H-O_W interaction, whereas rehybridization leads to an increase in the s character of the carbon hybrid orbital in C-H bond, thus causing contraction and shortening of C-H bond in complexes having C-H···O_W interactions. From the topological analysis, an excellent linear correlation is found to exist between stabilization energy (ΔE _BSSE), electron density (ρ_c), and its Laplacian (∇²ρ_c) at the bond critical points.

Ligands of Therapeutic Utility for the Liver X Receptors

Liver X receptors (LXRs) have been increasingly recognized as a potential therapeutic target to treat pathological conditions ranging from vascular and metabolic diseases, neurological degeneration, to cancers that are driven by lipid metabolism. Amidst intensifying efforts to discover ligands that act through LXRs to achieve the sought-after pharmacological outcomes, several lead compounds are already being tested in clinical trials for a variety of disease interventions. While more potent and selective LXR ligands continue to emerge from screening of small molecule libraries, rational design, and empirical medicinal chemistry approaches, challenges remain in minimizing undesirable effects of LXR activation on lipid metabolism. This review provides a summary of known endogenous, naturally occurring, and synthetic ligands. The review also offers considerations from a molecular modeling perspective with which to design more specific LXRβ ligands based on the interaction energies of ligands and the important amino acid residues in the LXRβ ligand binding domain.

Assessment of the Presence and Strength of H-Bonds by Means of Corrected NMR

The downfield shift of the NMR signal of the bridging proton in a H-bond (HB) is composed of two elements. The formation of the HB causes charge transfer and polarization that lead to a deshielding. A second factor is the mere presence of the proton-accepting group, whose electron density and response to an external magnetic field induce effects at the position of the bridging proton, exclusive of any H-bonding phenomenon. This second positional shielding must be subtracted from the full observed shift in order to assess the deshielding of the proton caused purely by HB formation. This concept is applied to a number of H-bonded systems, both intramolecular and intermolecular. When the positional shielding is removed, the remaining chemical shift is in much better coincidence with other measures of HB strength.

Frequent side chain methyl carbon-oxygen hydrogen bonding in proteins revealed by computational and stereochemical analysis of neutron structures

The propensity of backbone Cα atoms to engage in carbon-oxygen (CH · · · O) hydrogen bonding is well-appreciated in protein structure, but side chain CH · · · O hydrogen bonding remains largely uncharacterized. The extent to which side chain methyl groups in proteins participate in CH · · · O hydrogen bonding is examined through a survey of neutron crystal structures, quantum chemistry calculations, and molecular dynamics simulations. Using these approaches, methyl groups were observed to form stabilizing CH · · · O hydrogen bonds within protein structure that are maintained through protein dynamics and participate in correlated motion. Collectively, these findings illustrate that side chain methyl CH · · · O hydrogen bonding contributes to the energetics of protein structure and folding.

C–H/O interactions of nucleic bases with a water molecule: a crystallographic and quantum chemical study

The C–H/O interactions of nucleic bases are substantially stronger than the C–H/O interactions of benzene and pyridine. These results can be very important for molecular recognition of DNA and RNA.

Can a C-H···O interaction be a determinant of conformation?

Whether nonconventional hydrogen bonds, such as the C-H···O interaction, are a consequence or a determinant of conformation is a long-running and unresolved issue. Here we outline a solid-state and quantum mechanical study designed to investigate whether a C-H···O interaction can override the significant trans-planar conformational preferences of α-fluoroamide substituents. A profound change in dihedral angle from trans-planar((OCCF)) to cis-planar((OCCF)) observed on introducing an acceptor group for a C-H···O hydrogen bond is consistent with this interaction functioning as a determinant of conformation in certain systems. This testifies to the potential influence of the C-H···O hydrogen bond and is consistent with the assignment of this interaction as a contributor to overall conformation in both model and natural systems.

Molecular Self Assembly: Solvent Guests Tune the Conformation of a Series of 2,6-Bis(2-anilinoethynyl)pyridine-Based Ureas

The conformations of 2,6-bis(2-anilinoethynyl)pyridine-based urea receptors were studied by single crystal X-ray diffraction methods and revealed a rich conformational flexibility influenced by solvents. Whereas receptor L(1) in DMSO prefers an "S" conformation, receptor L(1) crystallizes in an "O" conformation from DMSO/CH(3)OH binary solvent system, and a "W" conformation in the ternary solvent mixture DMSO/toluene/1,4-dioxane. In the case of L(2), the molecule adopts an "S" conformation where water molecules are sandwiched between two molecules of L(2) to form a dimer. Similar to L(2), L(3) also forms a dimer where water molecules are sandwiched between L(3) molecules, which are capped with two molecules of DMSO. Such a capping DMSO solvate is lacking in the case of L(2). Taken together, these results demonstrate that the conformation of 2,6-bis(2-anilinoethynyl) pyridine-based urea receptors can be dramatically manipulated and tuned by the choice of crystallization solvents.

Weak H-bonds. Comparisons of CH···O to NH···O in proteins and PH···N to direct P···N interactions

Whereas CH···O H-bonds are usually weaker than interpeptide NH···O H-bonds, this is not necessarily the case within proteins. The nominally weaker CH···O are surprisingly strong, comparable to, and in some cases stronger than, the NH···O H-bonds in the context of the forces that hold together the adjacent strands in protein β-sheets. The peptide NH is greatly weakened as proton donor in certain conformations of the protein backbone, particularly extended structures, and forms correspondingly weaker H-bonds. The PH group is a weak proton donor, but will form PH···N H-bonds. However, there is a stronger interaction in which P can engage, in which the P atom, not the H, directly approaches the N electron donor to establish a direct P···N interaction. This approach is stabilized by the same sort of electron transfer from the N lone pair to the P-H σ* antibond that characterizes the PH···N H-bond.

Identification of spectroscopic patterns of CH...O H-bonds in proteins

Ab initio calculations are used to identify characteristics of vibrational and NMR spectra that signal the involvement of a protein backbone in a CH...O H-bond and that distinguish this sort of interaction from other H-bonds in which a protein might participate. Glycine and alanine dipeptides, in both their C7 and C5 minimum-energy structures, are paired with formamide in a number of different H-bonding arrangements. The CH...O H-bond is characterized by a small contraction of the C-H bond length, along with a blue shift in its stretching frequency, accompanied by an intensification of this vibrational band. In the context of NMR spectra, the bridging CH proton's chemical shift is moved downfield by 1-2 ppm. The aforementioned features are not produced by other H-bonds in which the protein backbone might participate, such as NH proton donation or accepting a proton via the peptide C=O.

Chemical mechanisms of histone lysine and arginine modifications

Histone lysine and arginine residues are subject to a wide array of post-translational modifications including methylation, citrullination, acetylation, ubiquitination, and sumoylation. The combinatorial action of these modifications regulates critical DNA processes including replication, repair, and transcription. In addition, enzymes that modify histone lysine and arginine residues have been correlated with a variety of human diseases including arthritis, cancer, heart disease, diabetes, and neurodegenerative disorders. Thus, it is important to fully understand the detailed kinetic and chemical mechanisms of these enzymes. Here, we review recent progress towards determining the mechanisms of histone lysine and arginine modifying enzymes. In particular, the mechanisms of S-adenosyl-methionine (AdoMet) dependent methyltransferases, FAD-dependent demethylases, iron dependent demethylases, acetyl-CoA dependent acetyltransferases, zinc dependent deacetylases, NAD(+) dependent deacetylases, and protein arginine deiminases are covered. Particular attention is paid to the conserved active-site residues necessary for catalysis and the individual chemical steps along the catalytic pathway. When appropriate, areas requiring further work are discussed.